結果
問題 | No.519 アイドルユニット |
ユーザー |
![]() |
提出日時 | 2017-05-28 23:02:08 |
言語 | C++11 (gcc 13.3.0) |
結果 |
WA
|
実行時間 | - |
コード長 | 3,063 bytes |
コンパイル時間 | 714 ms |
コンパイル使用メモリ | 82,040 KB |
実行使用メモリ | 6,944 KB |
最終ジャッジ日時 | 2024-09-21 15:47:50 |
合計ジャッジ時間 | 33,259 ms |
ジャッジサーバーID (参考情報) |
judge4 / judge2 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
other | AC * 23 WA * 11 |
ソースコード
#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(100. / temp); 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; }