結果
問題 | No.5007 Steiner Space Travel |
ユーザー | Spady_JP |
提出日時 | 2022-07-30 14:58:11 |
言語 | C++17 (gcc 12.3.0 + boost 1.83.0) |
結果 |
AC
|
実行時間 | 902 ms / 1,000 ms |
コード長 | 3,756 bytes |
コンパイル時間 | 3,904 ms |
実行使用メモリ | 3,820 KB |
スコア | 1,242,744 |
最終ジャッジ日時 | 2022-07-30 14:58:55 |
合計ジャッジ時間 | 33,150 ms |
ジャッジサーバーID (参考情報) |
judge13 / judge12 |
純コード判定しない問題か言語 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 902 ms
3,680 KB |
testcase_01 | AC | 902 ms
3,640 KB |
testcase_02 | AC | 901 ms
3,796 KB |
testcase_03 | AC | 901 ms
3,752 KB |
testcase_04 | AC | 901 ms
3,740 KB |
testcase_05 | AC | 902 ms
3,764 KB |
testcase_06 | AC | 901 ms
3,728 KB |
testcase_07 | AC | 902 ms
3,692 KB |
testcase_08 | AC | 902 ms
3,748 KB |
testcase_09 | AC | 902 ms
3,764 KB |
testcase_10 | AC | 902 ms
3,820 KB |
testcase_11 | AC | 902 ms
3,740 KB |
testcase_12 | AC | 901 ms
3,688 KB |
testcase_13 | AC | 902 ms
3,644 KB |
testcase_14 | AC | 902 ms
3,724 KB |
testcase_15 | AC | 901 ms
3,688 KB |
testcase_16 | AC | 902 ms
3,636 KB |
testcase_17 | AC | 902 ms
3,640 KB |
testcase_18 | AC | 901 ms
3,744 KB |
testcase_19 | AC | 902 ms
3,640 KB |
testcase_20 | AC | 902 ms
3,636 KB |
testcase_21 | AC | 902 ms
3,740 KB |
testcase_22 | AC | 901 ms
3,644 KB |
testcase_23 | AC | 902 ms
3,772 KB |
testcase_24 | AC | 901 ms
3,640 KB |
testcase_25 | AC | 901 ms
3,748 KB |
testcase_26 | AC | 901 ms
3,820 KB |
testcase_27 | AC | 901 ms
3,636 KB |
testcase_28 | AC | 902 ms
3,636 KB |
testcase_29 | AC | 902 ms
3,820 KB |
ソースコード
#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); while (tm.elapsed() < 900) { int ds = dist(0, out.size()); int dt = dist(0, out.size()); swap(out[ds], out[dt]); auto new_score = calc(points, station, out); if (new_score >= pre_score) { //悪化 swap(out[ds], out[dt]); continue; } pre_score = new_score; 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; }