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

提出ソース

結果

問題	No.1365 [Cherry 1st Tune] Whose Fault?
ユーザー	iaNTU
提出日時	2021-01-20 03:32:00
言語	C++17 (gcc 12.3.0 + boost 1.83.0)
結果	AC
実行時間	316 ms / 2,000 ms
コード長	5,154 bytes
コンパイル時間	3,419 ms
コンパイル使用メモリ	212,576 KB
実行使用メモリ	11,348 KB
最終ジャッジ日時	2023-08-28 11:21:27
合計ジャッジ時間	24,793 ms
ジャッジサーバーID （参考情報）	judge11 / judge13

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

テストケース

テストケース表示

入力	結果	実行時間実行使用メモリ
testcase_00	AC	2 ms 9,668 KB
testcase_01	AC	3 ms 9,632 KB
testcase_02	AC	3 ms 9,600 KB
testcase_03	AC	3 ms 9,776 KB
testcase_04	AC	3 ms 9,596 KB
testcase_05	AC	3 ms 9,660 KB
testcase_06	AC	3 ms 9,824 KB
testcase_07	AC	3 ms 9,720 KB
testcase_08	AC	3 ms 9,616 KB
testcase_09	AC	3 ms 9,844 KB
testcase_10	AC	3 ms 9,780 KB
testcase_11	AC	3 ms 9,784 KB
testcase_12	AC	3 ms 9,604 KB
testcase_13	AC	30 ms 9,784 KB
testcase_14	AC	20 ms 9,812 KB
testcase_15	AC	6 ms 9,712 KB
testcase_16	AC	4 ms 9,800 KB
testcase_17	AC	18 ms 9,928 KB
testcase_18	AC	10 ms 9,900 KB
testcase_19	AC	29 ms 9,844 KB
testcase_20	AC	16 ms 9,848 KB
testcase_21	AC	26 ms 9,936 KB
testcase_22	AC	22 ms 9,740 KB
testcase_23	AC	55 ms 10,632 KB
testcase_24	AC	123 ms 10,628 KB
testcase_25	AC	254 ms 10,996 KB
testcase_26	AC	172 ms 10,860 KB
testcase_27	AC	216 ms 10,312 KB
testcase_28	AC	88 ms 9,864 KB
testcase_29	AC	72 ms 10,108 KB
testcase_30	AC	316 ms 11,348 KB
testcase_31	AC	277 ms 11,348 KB
testcase_32	AC	150 ms 10,864 KB
testcase_33	AC	253 ms 10,912 KB
testcase_34	AC	197 ms 10,504 KB
testcase_35	AC	110 ms 10,736 KB
testcase_36	AC	228 ms 11,052 KB
testcase_37	AC	49 ms 10,124 KB
testcase_38	AC	55 ms 10,792 KB
testcase_39	AC	294 ms 11,096 KB
testcase_40	AC	295 ms 11,192 KB
testcase_41	AC	220 ms 10,928 KB
testcase_42	AC	206 ms 10,396 KB
testcase_43	AC	3 ms 9,780 KB
testcase_44	AC	50 ms 10,932 KB
testcase_45	AC	87 ms 10,320 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]);

	// Limit check
	// assert(1 <= n && n <= 100000);
	// assert(1 <= q && q <= 200000);
	// for (int i = 1; i <= n; i++) assert(0 <= sa[i] && sa[i] <= 1000);
	// for (int i = 1; i <= n; i++) assert(0 <= sb[i] && sb[i] <= 1000);
	// for (int i = 1; i <= n; i++) assert(0 <= p[i] && p[i] <= 100);
	// for (int i = 1; i <= q; i++) assert(x[i] != y[i]);
	// for (int i = 1; i <= q; i++) assert(1 <= x[i] && x[i] <= n);
	// for (int i = 1; i <= q; i++) assert(1 <= y[i] && y[i] <= n);

	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]];
#pragma unroll
	for (int i = 1; i <= n; i++) srp[i] = 1.0 / p[i];
#pragma unroll
	for (int i = 1; i <= n; i++) f[i] = i;
	
#pragma unroll
	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

結果

テストケース

ソースコード