#607512 (C++17) No.1365 [Cherry 1st Tune] Whose Fault?

提出ソース

結果

問題	No.1365 [Cherry 1st Tune] Whose Fault?
ユーザー	iaNTU
提出日時	2021-01-20 03:15:49
言語	C++17 (gcc 12.3.0 + boost 1.83.0)
結果	AC
実行時間	283 ms / 2,000 ms
コード長	4,573 bytes
コンパイル時間	2,890 ms
コンパイル使用メモリ	217,028 KB
実行使用メモリ	12,400 KB
最終ジャッジ日時	2024-06-09 06:48:27
合計ジャッジ時間	20,720 ms
ジャッジサーバーID （参考情報）	judge3 / judge5

このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示

入力	結果	実行時間実行使用メモリ
testcase_00	AC	3 ms 11,988 KB
testcase_01	AC	3 ms 11,980 KB
testcase_02	AC	3 ms 11,852 KB
testcase_03	AC	2 ms 11,856 KB
testcase_04	AC	3 ms 11,856 KB
testcase_05	AC	3 ms 11,980 KB
testcase_06	AC	2 ms 11,852 KB
testcase_07	AC	3 ms 11,988 KB
testcase_08	AC	2 ms 11,980 KB
testcase_09	AC	2 ms 11,984 KB
testcase_10	AC	3 ms 11,852 KB
testcase_11	AC	2 ms 11,852 KB
testcase_12	AC	2 ms 11,860 KB
testcase_13	AC	25 ms 11,948 KB
testcase_14	AC	19 ms 11,916 KB
testcase_15	AC	6 ms 12,048 KB
testcase_16	AC	4 ms 11,952 KB
testcase_17	AC	16 ms 12,096 KB
testcase_18	AC	9 ms 12,044 KB
testcase_19	AC	25 ms 12,044 KB
testcase_20	AC	15 ms 11,812 KB
testcase_21	AC	24 ms 12,044 KB
testcase_22	AC	20 ms 12,044 KB
testcase_23	AC	49 ms 12,108 KB
testcase_24	AC	109 ms 12,044 KB
testcase_25	AC	225 ms 11,920 KB
testcase_26	AC	155 ms 12,044 KB
testcase_27	AC	200 ms 12,164 KB
testcase_28	AC	79 ms 11,976 KB
testcase_29	AC	64 ms 12,160 KB
testcase_30	AC	283 ms 12,044 KB
testcase_31	AC	247 ms 12,172 KB
testcase_32	AC	135 ms 12,040 KB
testcase_33	AC	231 ms 12,104 KB
testcase_34	AC	174 ms 12,172 KB
testcase_35	AC	94 ms 12,300 KB
testcase_36	AC	199 ms 12,044 KB
testcase_37	AC	44 ms 12,216 KB
testcase_38	AC	48 ms 12,400 KB
testcase_39	AC	262 ms 12,048 KB
testcase_40	AC	261 ms 12,288 KB
testcase_41	AC	195 ms 12,172 KB
testcase_42	AC	179 ms 12,024 KB
testcase_43	AC	3 ms 11,988 KB
testcase_44	AC	43 ms 12,296 KB
testcase_45	AC	76 ms 12,048 KB

権限があれば一括ダウンロードができます

ソースコード

raw source code

// Exported by Exporter.exe

// Included from B.cpp
// Compile flags -Wall -Wextra -Wshadow -D_GLIBCXX_ASSERTIONS -DDEBUG -ggdb3 -fmax-errors=2
#pragma GCC optimize("O3")
#include <bits/stdc++.h>
using namespace std;
#define PB push_back
#define F first
#define S second
#define MP make_pair
#define MTP make_tuple
typedef long long int ll;

constexpr int kN = int(2E5 + 10);
// constexpr int kMod = 998244353;
// constexpr int kInf = 0x3f3f3f3f
// constexpr ll kInf = 0x3f3f3f3f3f3f3f3f;

template <typename T> T ABS(T n) {return n >= 0 ? n : -n;}


// Included from C:\Users\ianli\Desktop\CP\template\Various\Fast_IO.cpp
static inline char Get_Raw_Char() {
	static char buf[1 << 16], *p = buf, *end = buf;
	if (p == end) {
		if ((end = buf + fread(buf, 1, 1 << 16, stdin)) == buf) return '\0';
		p = buf;
	}
	return *p++;
}

static inline int Get_Digit() {
	char c = Get_Raw_Char();
	while (!isdigit(c)) c = Get_Raw_Char();
	return int(c - '0');
}

static inline int Get_PInt() {
	char c = Get_Raw_Char();
	int ret;
	while (!isdigit(c)) c = Get_Raw_Char();
	ret = int(c - '0');
	while (isdigit(c = Get_Raw_Char())) ret = ret * 10 + int(c - '0');
	return ret;
}

static inline int Get_Int() {
	char c = Get_Raw_Char();
	bool neg = false;
	int ret;
	while (!isdigit(c)) {
		if (c == '-') neg = true;
		c = Get_Raw_Char();
	}
	ret = int(c - '0');
	while (isdigit(c = Get_Raw_Char())) ret = ret * 10 + int(c - '0');
	if (neg) return -ret;
	return ret;
}

static inline long long int Get_ll() {
	char c = Get_Raw_Char();
	bool neg = false;
	long long int ret;
	while (!isdigit(c)) {
		if (c == '-') neg = true;
		c = Get_Raw_Char();
	}
	ret = int(c - '0');
	while (isdigit(c = Get_Raw_Char())) ret = ret * 10 + int(c - '0');
	if (neg) return -ret;
	return ret;
}

static inline long long int Get_Pll() {
	char c = Get_Raw_Char();
	long long int ret;
	while (!isdigit(c)) c = Get_Raw_Char();
	ret = int(c - '0');
	while (isdigit(c = Get_Raw_Char())) ret = ret * 10 + int(c - '0');
	return ret;
}

template <typename T> static inline void Read_P(T &n) {
	static_assert(is_integral<T>::value);
	char c = Get_Raw_Char();
	while (!isdigit(c)) c = Get_Raw_Char();
	n = int(c - '0');
	while (isdigit(c = Get_Raw_Char())) n = n * 10 + int(c - '0');
	return ;
}

template <typename T> static inline void Read(T &n) {
	static_assert(is_integral<T>::value);
	char c = Get_Raw_Char();
	bool neg = false;
	while (!isdigit(c)) {
		if (c == '-') neg = true;
		c = Get_Raw_Char();
	}
	n = int(c - '0');
	while (isdigit(c = Get_Raw_Char())) n = n * 10 + int(c - '0');
	if (neg) n = -n;
	return ;
}

template <typename T> static inline void Read_Digit(T &n) {
	static_assert(is_integral<T>::value);
	char c = Get_Raw_Char();
	while (!isdigit(c)) c = Get_Raw_Char();
	n = int(c - '0');
	return ;
}

template <typename T, typename... Targs> static inline void Read(T &n, Targs&... Fargs) {
	Read(n);
	return Read(Fargs...);
}

template <typename T, typename... Targs> static inline void Read_Digit(T &n, Targs&... Fargs) {
	Read_Digit(n);
	return Read_Digit(Fargs...);
}

template <typename T, typename... Targs> static inline void Read_P(T &n, Targs&... Fargs) {
	Read_P(n);
	return Read_P(Fargs...);
}
// End of C:\Users\ianli\Desktop\CP\template\Various\Fast_IO.cpp

int f[kN], p[kN], x[kN], y[kN], sa[kN], sb[kN], idx[kN];
double ans, srp[kN], val[kN], tsp[kN];

double Cal(int n) {
	ll sd = ABS(sa[n] - sb[n]);
	return val[n] = sd * sd / srp[n];
}

int Find(int n) {return n == f[n] ? n : f[n] = Find(f[n]);}
void Merge(int l, int r) {
	l = Find(l), r = Find(r);
	if (l != r) {
		if (idx[l] < idx[r]) swap(l, r);
		ans -= val[l];
		ans -= val[r];
		f[r] = l;
		sa[l] += sa[r];
		sb[l] += sb[r];
		srp[l] += srp[r];
		ans += Cal(l);
	}
	return ;
}

int main() {
	//ios::sync_with_stdio(false);
	//cin.tie(0);
	//freopen("file_name", "r", stdin);
	//freopen("file_name", "w", stdout);
	//fstream input, output;
	//input.open("file_name", ios::in);
	//output.open("file_name", ios::out);
	int n, q;
	mt19937 rng;
	Read_P(n);
	for (int i = 1; i <= n; i++) Read_P(sa[i]);
	for (int i = 1; i <= n; i++) Read_P(sb[i]);
	for (int i = 1; i <= n; i++) Read_P(p[i]);
	Read_P(q);
	for (int i = 1; i <= q; i++) Read_P(x[i], y[i]);
	for (int i = 1; i <= n; i++) idx[i] = rng();

	for (int i = 0; i <= 100; i++) tsp[i] = 1.0 / i;

	for (int i = 1; i <= n; i++) srp[i] = tsp[p[i]];
	for (int i = 1; i <= n; i++) f[i] = i;
	for (int i = 1; i <= n; i++) ans += Cal(i);

	for (int i = 1; i <= q; i++) {
		Merge(x[i], y[i]);
		printf("%.20lf\n", ans);
	}
	//input.close();
	//output.close();
}
// End of B.cpp

yukicoder

結果

テストケース

ソースコード