結果
問題 | No.5007 Steiner Space Travel |
ユーザー | Spady_JP |
提出日時 | 2022-07-30 15:16:41 |
言語 | C++17 (gcc 12.3.0 + boost 1.83.0) |
結果 |
TLE
|
実行時間 | - |
コード長 | 3,984 bytes |
コンパイル時間 | 3,976 ms |
スコア | 0 |
最終ジャッジ日時 | 2022-07-30 15:17:20 |
合計ジャッジ時間 | 38,154 ms |
ジャッジサーバーID (参考情報) |
judge10 / judge13 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | TLE | - |
testcase_01 | TLE | - |
testcase_02 | TLE | - |
testcase_03 | TLE | - |
testcase_04 | TLE | - |
testcase_05 | TLE | - |
testcase_06 | TLE | - |
testcase_07 | TLE | - |
testcase_08 | TLE | - |
testcase_09 | TLE | - |
testcase_10 | TLE | - |
testcase_11 | TLE | - |
testcase_12 | TLE | - |
testcase_13 | TLE | - |
testcase_14 | TLE | - |
testcase_15 | TLE | - |
testcase_16 | TLE | - |
testcase_17 | TLE | - |
testcase_18 | TLE | - |
testcase_19 | TLE | - |
testcase_20 | TLE | - |
testcase_21 | TLE | - |
testcase_22 | TLE | - |
testcase_23 | TLE | - |
testcase_24 | TLE | - |
testcase_25 | TLE | - |
testcase_26 | TLE | - |
testcase_27 | TLE | - |
testcase_28 | TLE | - |
testcase_29 | TLE | - |
ソースコード
#ifdef ONLINE_JUDGE #define NDEBUG 1 #endif #include <bits/stdc++.h> //ifブロックはインデントがいる #if __has_include(<atcoder/all>) #include <atcoder/all> 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<int,int>; using PLI = pair<long long,long long>; template<class T> inline bool chmin(T& a, T b) { if (a > b) { a = b; return true; } return false; } template<class T> 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<std::chrono::milliseconds>(dur) .count(); } private: std::chrono::_V2::steady_clock::time_point m_start; }; long long calc(vector<P> &points, vector<P> &station,vector<P> &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<P> points(n); rep(i, n) { cin >> points[i].first >> points[i].second; } vector<P> station(m); rep(i, m) { station[i] = {dist(0, 1000), dist(0, 1000)};} vector<P> 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; }