結果

問題 No.519 アイドルユニット
コンテスト
ユーザー kurenai3110
提出日時 2017-05-28 23:07:17
言語 C++11（廃止可能性あり）
(gcc 13.3.0)
結果
AC  
実行時間 902 ms / 1,000 ms
コード長 3,095 bytes
コンパイル時間 692 ms
コンパイル使用メモリ 82,012 KB
実行使用メモリ 6,820 KB
最終ジャッジ日時 2024-10-03 03:52:50
合計ジャッジ時間 33,214 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge1
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 34
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <iostream>
#include <vector>
#include <algorithm>
#include <string>
#include <map>
#include <chrono>
#include <cmath>
using namespace std;

class RandomXor {
	unsigned static int x, y, z, w;
public:
	RandomXor() {
		x = 31103110, y = 123456789, z = 521288629, w = 88675123;
	}

	unsigned int rand()
	{
		unsigned int t;
		t = (x ^ (x << 11)); x = y; y = z; z = w;
		return(w = (w ^ (w >> 19)) ^ (t ^ (t >> 8)));
	}
};
unsigned int RandomXor::x, RandomXor::y, RandomXor::z, RandomXor::w;


int main()
{
	auto start = chrono::system_clock::now();
	RandomXor rx;
	int N; cin >> N;
	vector<vector<int>>LOVE(N);
	for (int i = 0; i < N; i++) {
		LOVE[i].resize(N);
		for (int j = 0; j < N; j++) {
			cin >> LOVE[i][j];
		}
	}

	int P = N / 2;
	vector<pair<int, int>>Pi(P);
	int maxLOVE = 0,tmpLOVE=0;
	for (int i = 0; i < P; i++) {
		Pi[i] = make_pair(2 * i, 2 * i + 1);
		tmpLOVE += LOVE[2 * i][2 * i + 1];
	}
	maxLOVE = tmpLOVE;

	double startTemp = 100;
	double endTemp = 10;
	double T = 0.9;
	while (1) {
		auto end = chrono::system_clock::now();
		double t = chrono::duration<double>(end - start).count();
		if (t > 0.9)break;

		int r1, r2;
		r1 = rx.rand() % P;
		r2 = rx.rand() % P;
		if (r1 == r2)continue;

		double temp = startTemp + (endTemp - startTemp) * t / T;
		double probability = exp(-50. / temp);
		//cerr << probability << endl;

		pair<int, int>p1 = Pi[r1], p2 = Pi[r2];
		if (tmpLOVE < tmpLOVE - LOVE[p1.first][p1.second] - LOVE[p2.first][p2.second] + LOVE[p1.first][p2.first] + LOVE[p1.second][p2.second]) {
			tmpLOVE += -LOVE[p1.first][p1.second] - LOVE[p2.first][p2.second] + LOVE[p1.first][p2.first] + LOVE[p1.second][p2.second];
			maxLOVE = max(maxLOVE, tmpLOVE);
			swap(Pi[r1].second, Pi[r2].first);
		}
		else if (tmpLOVE < tmpLOVE - LOVE[p1.first][p1.second] - LOVE[p2.first][p2.second] + LOVE[p1.first][p2.second] + LOVE[p2.first][p1.second]) {
			tmpLOVE += -LOVE[p1.first][p1.second] - LOVE[p2.first][p2.second] + LOVE[p1.first][p2.second] + LOVE[p2.first][p1.second];
			maxLOVE = max(maxLOVE, tmpLOVE);
			swap(Pi[r1].second, Pi[r2].second);
		}
		else if (tmpLOVE < tmpLOVE - LOVE[p1.first][p1.second] - LOVE[p2.first][p2.second] + LOVE[p2.first][p1.second] + LOVE[p1.first][p2.second]) {
			tmpLOVE += -LOVE[p1.first][p1.second] - LOVE[p2.first][p2.second] + LOVE[p2.second][p1.second] + LOVE[p2.first][p1.first];
			maxLOVE = max(maxLOVE, tmpLOVE);
			swap(Pi[r1].first, Pi[r2].second);
		}
		else if (tmpLOVE < tmpLOVE - LOVE[p1.first][p1.second] - LOVE[p2.first][p2.second] + LOVE[p2.first][p1.second] + LOVE[p1.first][p2.second]) {
			tmpLOVE += -LOVE[p1.first][p1.second] - LOVE[p2.first][p2.second] + LOVE[p2.first][p1.second] + LOVE[p1.first][p2.second];
			maxLOVE = max(maxLOVE, tmpLOVE);
			swap(Pi[r1].first, Pi[r2].first);
		}
		else if ((double)rx.rand() / UINT32_MAX < probability) {
			tmpLOVE += -LOVE[p1.first][p1.second] - LOVE[p2.first][p2.second] + LOVE[p1.first][p2.first] + LOVE[p1.second][p2.second];
			maxLOVE = max(maxLOVE, tmpLOVE);
			swap(Pi[r1].second, Pi[r2].first);
		}
	}

	cout << maxLOVE << endl;


    return 0;
}
0