結果

問題 No.5020 Averaging
ユーザー e869120e869120
提出日時 2024-02-20 13:56:55
言語 C++14
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 859 ms / 1,000 ms
コード長 3,298 bytes
コンパイル時間 1,146 ms
コンパイル使用メモリ 95,948 KB
実行使用メモリ 6,676 KB
スコア 79,383,898
最終ジャッジ日時 2024-02-25 12:42:35
合計ジャッジ時間 46,673 ms
ジャッジサーバーID
（参考情報） 		judge14 / judge15
純コード判定しない問題か言語
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 856 ms
6,676 KB
testcase_01 AC 857 ms
6,676 KB
testcase_02 AC 856 ms
6,676 KB
testcase_03 AC 858 ms
6,676 KB
testcase_04 AC 858 ms
6,676 KB
testcase_05 AC 856 ms
6,676 KB
testcase_06 AC 857 ms
6,676 KB
testcase_07 AC 856 ms
6,676 KB
testcase_08 AC 855 ms
6,676 KB
testcase_09 AC 857 ms
6,676 KB
testcase_10 AC 857 ms
6,676 KB
testcase_11 AC 857 ms
6,676 KB
testcase_12 AC 856 ms
6,676 KB
testcase_13 AC 857 ms
6,676 KB
testcase_14 AC 857 ms
6,676 KB
testcase_15 AC 857 ms
6,676 KB
testcase_16 AC 857 ms
6,676 KB
testcase_17 AC 857 ms
6,676 KB
testcase_18 AC 857 ms
6,676 KB
testcase_19 AC 857 ms
6,676 KB
testcase_20 AC 858 ms
6,676 KB
testcase_21 AC 858 ms
6,676 KB
testcase_22 AC 857 ms
6,676 KB
testcase_23 AC 855 ms
6,676 KB
testcase_24 AC 857 ms
6,676 KB
testcase_25 AC 858 ms
6,676 KB
testcase_26 AC 856 ms
6,676 KB
testcase_27 AC 856 ms
6,676 KB
testcase_28 AC 857 ms
6,676 KB
testcase_29 AC 857 ms
6,676 KB
testcase_30 AC 857 ms
6,676 KB
testcase_31 AC 859 ms
6,676 KB
testcase_32 AC 856 ms
6,676 KB
testcase_33 AC 857 ms
6,676 KB
testcase_34 AC 856 ms
6,676 KB
testcase_35 AC 855 ms
6,676 KB
testcase_36 AC 857 ms
6,676 KB
testcase_37 AC 856 ms
6,676 KB
testcase_38 AC 856 ms
6,676 KB
testcase_39 AC 856 ms
6,676 KB
testcase_40 AC 858 ms
6,676 KB
testcase_41 AC 856 ms
6,676 KB
testcase_42 AC 856 ms
6,676 KB
testcase_43 AC 858 ms
6,676 KB
testcase_44 AC 856 ms
6,676 KB
testcase_45 AC 857 ms
6,676 KB
testcase_46 AC 858 ms
6,676 KB
testcase_47 AC 856 ms
6,676 KB
testcase_48 AC 857 ms
6,676 KB
testcase_49 AC 856 ms
6,676 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

// [想定解答 6]
//    - 想定解法 5 は運要素が強く，一発で合わせるのは難しい．
//    - そこで，貪欲法に以下のようなランダム性を加え，プログラムを 0.85 秒動かして一番良いものを解答として出力する．
//        * 想定解法 5 では，max(A[x] - RemainA / 2, B[x] - RemainB / 2) が最小のカードを選択していた．
//        * しかし，今回は上位 min(3, 候補数 / 3) 個の中からランダムに選ぶ．
//    - 79% 程度の得点が得られる．

#include <iostream>
#include <cmath>
#include <vector>
#include <algorithm>
using namespace std;
const long long Target = 500000000000000000LL;

void Operation(vector<long long>& A, vector<long long>& B, int idx1, int idx2) {
	long long avg1 = (A[idx1] + A[idx2]) / 2LL;
	long long avg2 = (B[idx1] + B[idx2]) / 2LL;
	A[idx1] = avg1; B[idx1] = avg2;
	A[idx2] = avg1; B[idx2] = avg2;
}

int main() {
	// Step 1. 入力
	vector<long long> A(45, 0);
	vector<long long> B(45, 0);
	long long N; cin >> N;
	for (int i = 0; i < N; i++) cin >> A[i] >> B[i];

	// Step 2. 最初のステップ
	long long MinV = (1LL << 60);
	long long MinID = -1;
	for (int i = 1; i < N; i++) {
		long long v1 = (A[0] + A[i]) / 2LL;
		long long v2 = (B[0] + B[i]) / 2LL;
		long long err = max(abs(v1 - Target), abs(v2 - Target));
		if (MinV > err) {
			MinV = err;
			MinID = i;
		}
	}
	Operation(A, B, 0, MinID);

	// Step 3. ランダム貪欲
	vector<pair<int, int>> Answer;
	int StartTime = clock();
	long long RecordScore = Target, Loops = 0;
	while (true) {
		vector<long long> C = A;
		vector<long long> D = B;
		long long RemainC = Target * 2 - C[0];
		long long RemainD = Target * 2 - D[0];
		vector<bool> Used(N, false);
		vector<int> Order; Order.push_back(0);
		Loops += 1;

		// 貪欲法スタート
		while (true) {
			vector<pair<long long, int>> Vec;
			for (int i = 1; i < N; i++) {
				if (Used[i] == true) continue;
				long long v1 = abs(RemainC - C[i]); if (v1 >= RemainC) continue;
				long long v2 = abs(RemainD - D[i]); if (v2 >= RemainD) continue;
				Vec.push_back(make_pair(max(v1, v2), i));
			}
			sort(Vec.begin(), Vec.end()); if (Vec.size() == 0) break;
			int idx = Vec[rand() % max(1, min(3, (int)Vec.size() / 3))].second;
			RemainC = min(6LL * Target, 2LL * RemainC - C[idx]);
			RemainD = min(6LL * Target, 2LL * RemainD - D[idx]);
			Order.push_back(idx);
			Used[idx] = true;
		}
		for (int i = 1; i < N; i++) {
			if (Used[i] == false) Order.push_back(i);
		}

		// スコアを出す
		vector<pair<int, int>> CurrentAns;
		CurrentAns.push_back(make_pair(0, MinID));
		for (int i = N - 2; i >= 0; i--) {
			CurrentAns.push_back(make_pair(Order[i], Order[i + 1]));
			Operation(C, D, Order[i], Order[i + 1]);
		}

		// スコアを比較
		long long CandScore = max(abs(Target - C[0]), abs(Target - D[0]));
		if (RecordScore > CandScore) {
			RecordScore = CandScore;
			Answer = CurrentAns;
		}
		if (Loops % 100 == 0 && clock() - StartTime > 85 * CLOCKS_PER_SEC / 100) break; // clock の呼び出しは時間がかかるので 100 回に 1 回呼び出す
	}

	// Step 5. 出力
	cout << Answer.size() << endl;
	for (int i = 0; i < Answer.size(); i++) cout << Answer[i].first + 1 << " " << Answer[i].second + 1 << endl;
	return 0;
}
0