#include #include #include #include #include #include #include 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>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>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(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; pairp1 = 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; }