#175790 (C++11) No.519 アイドルユニット

提出ソース

結果

問題	No.519 アイドルユニット
ユーザー	kurenai3110
提出日時	2017-05-28 23:07:17
言語	C++11 (gcc 11.4.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

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

テストケース

テストケース表示

入力	結果	実行時間実行使用メモリ
testcase_00	AC	902 ms 6,816 KB
testcase_01	AC	902 ms 6,816 KB
testcase_02	AC	901 ms 6,820 KB
testcase_03	AC	902 ms 6,816 KB
testcase_04	AC	902 ms 6,820 KB
testcase_05	AC	902 ms 6,820 KB
testcase_06	AC	901 ms 6,816 KB
testcase_07	AC	902 ms 6,816 KB
testcase_08	AC	902 ms 6,820 KB
testcase_09	AC	902 ms 6,816 KB
testcase_10	AC	902 ms 6,820 KB
testcase_11	AC	902 ms 6,816 KB
testcase_12	AC	902 ms 6,816 KB
testcase_13	AC	902 ms 6,820 KB
testcase_14	AC	902 ms 6,820 KB
testcase_15	AC	901 ms 6,820 KB
testcase_16	AC	902 ms 6,820 KB
testcase_17	AC	901 ms 6,816 KB
testcase_18	AC	902 ms 6,816 KB
testcase_19	AC	902 ms 6,816 KB
testcase_20	AC	902 ms 6,816 KB
testcase_21	AC	902 ms 6,820 KB
testcase_22	AC	902 ms 6,816 KB
testcase_23	AC	902 ms 6,820 KB
testcase_24	AC	902 ms 6,816 KB
testcase_25	AC	902 ms 6,816 KB
testcase_26	AC	902 ms 6,820 KB
testcase_27	AC	902 ms 6,816 KB
testcase_28	AC	901 ms 6,816 KB
testcase_29	AC	901 ms 6,816 KB
testcase_30	AC	902 ms 6,816 KB
testcase_31	AC	902 ms 6,820 KB
testcase_32	AC	902 ms 6,816 KB
testcase_33	AC	902 ms 6,816 KB

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

ソースコード

raw source code

#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;
}

yukicoder

結果

テストケース

ソースコード