結果

問題 No.5020 Averaging
コンテスト
ユーザー Pechi
提出日時 2024-02-25 14:36:32
言語 C++17(gcc12)
(gcc 12.3.0 + boost 1.87.0)
結果
AC  
実行時間 954 ms / 1,000 ms
コード長 3,361 bytes
コンパイル時間 3,283 ms
コンパイル使用メモリ 226,852 KB
実行使用メモリ 6,676 KB
スコア 31,289,794
最終ジャッジ日時 2024-02-25 14:37:37
合計ジャッジ時間 53,614 ms
ジャッジサーバーID
（参考情報） 		judge14 / judge12
純コード判定しない問題か言語
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 50
権限があれば一括ダウンロードができます

ソースコード

diff # 

#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() /2.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(max(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);
	while ((time=getTime()) <= 950) {
		double temp = (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;
		}
		if (++now == t) {
			now = 0;
			shuffle(table.begin(), table.end(), rnd);
		}
	}
	cerr << "score:" << ll(2000000 - 100000 * calc_score(a, b, ans)) << "\n";
	cout << ans.size() << "\n";
	ALL(x, ans)cout << x.first+1 << " " << x.second+1 << "\n";
	exit(0);
}
0