#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
実行時間	315 ms / 2,000 ms
コード長	4,573 bytes
コンパイル時間	2,646 ms
コンパイル使用メモリ	213,152 KB
実行使用メモリ	11,176 KB
最終ジャッジ日時	2023-08-28 11:19:54
合計ジャッジ時間	25,085 ms
ジャッジサーバーID （参考情報）	judge14 / judge15

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

テストケース

テストケース表示

入力	結果	実行時間実行使用メモリ
testcase_00	AC	3 ms 9,592 KB
testcase_01	AC	3 ms 9,544 KB
testcase_02	AC	3 ms 9,724 KB
testcase_03	AC	3 ms 9,544 KB
testcase_04	AC	3 ms 9,796 KB
testcase_05	AC	3 ms 9,756 KB
testcase_06	AC	3 ms 9,672 KB
testcase_07	AC	3 ms 9,648 KB
testcase_08	AC	3 ms 9,584 KB
testcase_09	AC	3 ms 9,544 KB
testcase_10	AC	3 ms 9,592 KB
testcase_11	AC	3 ms 9,740 KB
testcase_12	AC	3 ms 9,664 KB
testcase_13	AC	29 ms 9,744 KB
testcase_14	AC	20 ms 9,884 KB
testcase_15	AC	6 ms 9,708 KB
testcase_16	AC	4 ms 9,848 KB
testcase_17	AC	18 ms 9,912 KB
testcase_18	AC	10 ms 9,940 KB
testcase_19	AC	28 ms 9,720 KB
testcase_20	AC	16 ms 9,656 KB
testcase_21	AC	26 ms 9,728 KB
testcase_22	AC	22 ms 9,860 KB
testcase_23	AC	55 ms 10,500 KB
testcase_24	AC	124 ms 10,608 KB
testcase_25	AC	255 ms 11,056 KB
testcase_26	AC	174 ms 10,556 KB
testcase_27	AC	217 ms 10,380 KB
testcase_28	AC	88 ms 9,792 KB
testcase_29	AC	72 ms 10,020 KB
testcase_30	AC	315 ms 11,048 KB
testcase_31	AC	280 ms 11,176 KB
testcase_32	AC	154 ms 10,748 KB
testcase_33	AC	254 ms 10,788 KB
testcase_34	AC	199 ms 10,580 KB
testcase_35	AC	111 ms 10,840 KB
testcase_36	AC	231 ms 10,936 KB
testcase_37	AC	49 ms 10,172 KB
testcase_38	AC	57 ms 10,720 KB
testcase_39	AC	299 ms 11,140 KB
testcase_40	AC	296 ms 10,968 KB
testcase_41	AC	221 ms 10,796 KB
testcase_42	AC	205 ms 10,448 KB
testcase_43	AC	3 ms 9,540 KB
testcase_44	AC	50 ms 11,012 KB
testcase_45	AC	87 ms 10,148 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

結果

テストケース

ソースコード