#955383 (C++23) No.5020 Averaging

提出ソース

結果

問題	No.5020 Averaging
ユーザー	Pechi
提出日時	2024-02-25 14:52:59
言語	C++23 (gcc 12.3.0 + boost 1.83.0)
結果	AC
実行時間	953 ms / 1,000 ms
コード長	3,625 bytes
コンパイル時間	4,361 ms
コンパイル使用メモリ	272,416 KB
実行使用メモリ	6,676 KB
スコア	31,024,788
最終ジャッジ日時	2024-02-25 14:54:11
合計ジャッジ時間	54,813 ms
ジャッジサーバーID （参考情報）	judge10 / judge15
	純コード判定しない問題か言語

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

テストケース

テストケース表示

入力	結果	実行時間実行使用メモリ
testcase_00	AC	952 ms 6,676 KB
testcase_01	AC	953 ms 6,676 KB
testcase_02	AC	953 ms 6,676 KB
testcase_03	AC	953 ms 6,676 KB
testcase_04	AC	953 ms 6,676 KB
testcase_05	AC	953 ms 6,676 KB
testcase_06	AC	953 ms 6,676 KB
testcase_07	AC	952 ms 6,676 KB
testcase_08	AC	953 ms 6,676 KB
testcase_09	AC	953 ms 6,676 KB
testcase_10	AC	953 ms 6,676 KB
testcase_11	AC	953 ms 6,676 KB
testcase_12	AC	953 ms 6,676 KB
testcase_13	AC	953 ms 6,676 KB
testcase_14	AC	952 ms 6,676 KB
testcase_15	AC	953 ms 6,676 KB
testcase_16	AC	953 ms 6,676 KB
testcase_17	AC	953 ms 6,676 KB
testcase_18	AC	953 ms 6,676 KB
testcase_19	AC	953 ms 6,676 KB
testcase_20	AC	953 ms 6,676 KB
testcase_21	AC	953 ms 6,676 KB
testcase_22	AC	953 ms 6,676 KB
testcase_23	AC	953 ms 6,676 KB
testcase_24	AC	953 ms 6,676 KB
testcase_25	AC	953 ms 6,676 KB
testcase_26	AC	953 ms 6,676 KB
testcase_27	AC	953 ms 6,676 KB
testcase_28	AC	952 ms 6,676 KB
testcase_29	AC	952 ms 6,676 KB
testcase_30	AC	953 ms 6,676 KB
testcase_31	AC	953 ms 6,676 KB
testcase_32	AC	953 ms 6,676 KB
testcase_33	AC	953 ms 6,676 KB
testcase_34	AC	953 ms 6,676 KB
testcase_35	AC	953 ms 6,676 KB
testcase_36	AC	953 ms 6,676 KB
testcase_37	AC	953 ms 6,676 KB
testcase_38	AC	953 ms 6,676 KB
testcase_39	AC	952 ms 6,676 KB
testcase_40	AC	953 ms 6,676 KB
testcase_41	AC	953 ms 6,676 KB
testcase_42	AC	953 ms 6,676 KB
testcase_43	AC	953 ms 6,676 KB
testcase_44	AC	953 ms 6,676 KB
testcase_45	AC	953 ms 6,676 KB
testcase_46	AC	953 ms 6,676 KB
testcase_47	AC	953 ms 6,676 KB
testcase_48	AC	953 ms 6,676 KB
testcase_49	AC	953 ms 6,676 KB

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

ソースコード

raw source code

#define _USE_MATH_DEFINES
#pragma GCC target("avx2")
#pragma GCC optimize("O3")
#include<bits/stdc++.h>
#include<unordered_map>
#include<random>
#include <array>
#include <complex>
#include<chrono>
using namespace std;
#define LP(I,S,G) for (long long int I = (S); I < (G); ++I)
#define IN(X) 	for (int in = 0; in < X.size(); in++)cin >> X[in]
#define OUT(X) 	for (int in = 0; in < X.size(); in++)cout << X[in]<<" "
#define SORT(X) sort((X).begin(), (X).end())
#define CSORT(X,Y) sort(X.begin(), X.end(),Y)
#define COPY(X,Y) copy(X.begin(), X.end(), Y.begin())
#define ALL(X,Y) for (auto (X) :(Y))
#define FULL(a)  (a).begin(),(a).end()
#define BFS(Q,S) for(Q.push(S);Q.size()!=0;Q.pop())
typedef long long int ll;
typedef unsigned long long int ull;
long long int M = 998244353;

chrono::system_clock::time_point starttime;
using namespace std::chrono;
#ifndef ONLINE_JUDGE
#define DEBUG
#endif

inline float getTime() {
#ifdef DEBUG
	return duration_cast<milliseconds>(system_clock::now() - starttime).count() /4.0;
#else
	return duration_cast<milliseconds>(system_clock::now() - starttime).count();
#endif
}


int dx[] = { -1,0,1,0 }, dy[] = { 0,1,0,-1 };

ll MAX(ll A, ll B) { return ((A) > (B) ? (A) : (B)); }
ll MIN(ll A, ll B) { return ((A) < (B) ? (A) : (B)); }
inline long long int xor128() {
	static long long int x = 123456789, y = 362436069, z = 521288629, w = 88675123;
	long long int t = (x ^ (x << 11));
	x = y; y = z; z = w;
	return (w = (w ^ (w >> 19)) ^ (t ^ (t >> 8)));
}

long long int pow_n(long long int x, long long int n) {
	if (n == 0)
		return 1;
	if (n % 2 == 0)
		return pow_n(x * x, n / 2);
	else
		return x * pow_n(x, n - 1);
}

ll border = 5 * pow_n(10, 17);

double calc_score(vector<ll> a, vector<ll> b, vector<pair<ll, ll>>& ans) {
	for (int i = 0; i < ans.size(); ++i) {
		int u = ans[i].first, v = ans[i].second;
		a[u] = a[v] = (a[u] + a[v]) / 2;
		b[u] = b[v] = (b[u] + b[v]) / 2;
	}
	return log10(abs(border - a[0])+ abs(border - b[0]) + 1);
}

int main(int argc, char* argv[]) {
	starttime = chrono::system_clock::now();
	ios::sync_with_stdio(false);
	std::cin.tie(nullptr);
	mt19937 rnd(0);
	ll n, t = 50, now = 0;
	cin >> n;
	vector<ll> a(n), b(n);
	LP(i, 0, n)cin >> a[i] >> b[i];
	vector<pair<ll, ll>> ans(t);
	LP(i, 0, t) {
		ans[i].first = xor128() % n;
		ans[i].second = xor128() % (n - 1);
		if (ans[i].second >= ans[i].first)++ans[i].second;
	}
	double score = calc_score(a, b, ans);
	vector<int> table(t);
	for (int i = 0; i < t; ++i)table[i] = i;
	double time;
	shuffle(table.begin(), table.end(), rnd);
	int cnt = 0;
	while ((time=getTime()) <= 950) {
		double temp = 0.75*(950 - time) / 950.0;
		int i = table[now];
		int bx = ans[i].first, by = ans[i].second;
		ans[i].first = xor128() % n;
		ans[i].second = xor128() % (n - 1);
		if (ans[i].first != 0 && ans[i].second != 0 && xor128() % 10 == 0)ans[i].first = 0;
		if (ans[i].second >= ans[i].first)++ans[i].second;
		double new_score = calc_score(a, b, ans);
		if (xor128() % 10000 < exp((score-new_score) / temp) * 10000) {
			score = new_score;
		}
		else {
			ans[i].first = bx;
			ans[i].second = by;
		}
		++cnt;
		if (++now == t) {
			now = 0;
			shuffle(table.begin(), table.end(), rnd);
		}
	}
	cerr << cnt;
	for (int i = 0; i < ans.size(); ++i) {
		int u = ans[i].first, v = ans[i].second;
		a[u] = a[v] = (a[u] + a[v]) / 2;
		b[u] = b[v] = (b[u] + b[v]) / 2;
	}
	cerr << "score:" << ll(2000000 - 100000 * log10(max(abs(border - a[0]), abs(border - b[0])) + 1)) << "\n";
	cerr << a[0] << " " << b[0] << "\n";
	cout << ans.size() << "\n";
	ALL(x, ans)cout << x.first+1 << " " << x.second+1 << "\n";
	exit(0);
}

yukicoder

結果

テストケース

ソースコード