#ifdef ONLINE_JUDGE #define NDEBUG 1 #endif #include //ifブロックはインデントがいる #if __has_include() #include using namespace atcoder; #endif #define repp(i,l,r)for(long long i=(l);i<(r);i++) #define rep(i,n) for (long long i = 0; i < (n); ++i) #define per(i,n) for (long long i = (n); i >= 0; --i) #define all(v) v.begin(), v.end() const int INF = 1<<30; const long long LINF = 1LL << 60; const long long int MOD = 1000000007; using namespace std; using ll = long long; using P = pair; using PLI = pair; template inline bool chmin(T& a, T b) { if (a > b) { a = b; return true; } return false; } template inline bool chmax(T& a, T b) { if (a < b) { a = b; return true; } return false; } unsigned int randxor() { static unsigned int x = 123456789, y = 362436069, z = 521288629, w = 88675123; unsigned int t; t = (x ^ (x << 11)); x = y; y = z; z = w; return (w = (w ^ (w >> 19)) ^ (t ^ (t >> 8))); } int dist(int from, int to) { if (to == 0) return 0; if (to == from) return from; else return randxor() % (to - from) + from; } struct Timer { public: Timer() { restart(); } void restart() { m_start = std::chrono::steady_clock::now(); } auto elapsed() { std::chrono::steady_clock::time_point en = std::chrono::steady_clock::now(); auto dur = en - m_start; return std::chrono::duration_cast(dur) .count(); } private: std::chrono::_V2::steady_clock::time_point m_start; }; long long calc(vector

&points, vector

&station,vector

&out){//type, to int alpha = 5; long double d = 0; if(!(out.front() == out.back() and out.front() == pair(1,1))){ return -LINF; } auto [type, to] = out.front(); for (int i = 1; i < out.size(); i++) { auto x = (type == 1) ? points[out[i - 1].second] : station[out[i - 1].second]; auto y = (out[i].first == 1) ? points[out[i].second] : station[out[i].second]; long double tmp = (x.first - y.first) * (x.first - y.first) + (x.second - y.second) * (x.second - y.second); if (out[i].first == 1 and type == 1) {//ともに惑星 d += alpha * alpha * tmp; }else if(min(out[i].first,type) == 1){ d += alpha * tmp; }else{ d += tmp; } type = out[i].first; to = out[i].second; } return -d; } //ミョ(-ω- ?) int main() { cin.tie(nullptr) ; ios::sync_with_stdio(false) ; Timer tm; int n, m; cin >> n >> m; vector

points(n); rep(i, n) { cin >> points[i].first >> points[i].second; } vector

station(m); rep(i, m) { station[i] = {dist(0, 1000), dist(0, 1000)};} vector

out; rep(i, n) { out.push_back({1, i + 1}); } rep(i, m) { out.push_back({2, i + 1}); }; out.push_back({1, 1}); auto pre_score = calc(points, station, out); double start_temp = 50, end_temp = 10; while (tm.elapsed() < 1000) { int ds = dist(1, out.size() - 1); int dt = dist(1, out.size() - 1); swap(out[ds], out[dt]); auto new_score = calc(points, station, out); double temp = start_temp + (end_temp - start_temp) * (tm.elapsed()) / 900; double prob = exp((new_score-pre_score)/temp); if (prob > (randxor()%INF)/(double)INF) { // 確率probで遷移する pre_score = new_score; }else{ swap(out[ds], out[dt]); cerr << "new_score : " << new_score << "\n"; } } rep(i, m) { cout << station[i].first << " " << station[i].second << "\n"; } cout << out.size() << "\n"; rep(i, out.size()) { cout << out[i].first << " " << out[i].second << "\n"; } cerr << "Judge Score : " << round(1e9 / (1000 + sqrtl(-pre_score))) << "\n"; return 0; }