結果

問題 No.5007 Steiner Space Travel
ユーザー ebicochinealebicochineal
提出日時 2022-07-30 15:35:00
言語 C++14
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 2 ms / 1,000 ms
コード長 8,581 bytes
コンパイル時間 1,335 ms
実行使用メモリ 3,660 KB
スコア 4,232,757
最終ジャッジ日時 2022-07-30 15:35:03
合計ジャッジ時間 2,702 ms
ジャッジサーバーID
（参考情報） 		judge11 / judge12
純コード判定しない問題か言語
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 2 ms
3,576 KB
testcase_01 AC 1 ms
3,488 KB
testcase_02 AC 2 ms
3,576 KB
testcase_03 AC 2 ms
3,524 KB
testcase_04 AC 1 ms
3,524 KB
testcase_05 AC 1 ms
3,576 KB
testcase_06 AC 1 ms
3,532 KB
testcase_07 AC 2 ms
3,484 KB
testcase_08 AC 1 ms
3,488 KB
testcase_09 AC 1 ms
3,524 KB
testcase_10 AC 2 ms
3,424 KB
testcase_11 AC 2 ms
3,540 KB
testcase_12 AC 2 ms
3,660 KB
testcase_13 AC 2 ms
3,544 KB
testcase_14 AC 2 ms
3,416 KB
testcase_15 AC 2 ms
3,528 KB
testcase_16 AC 1 ms
3,572 KB
testcase_17 AC 2 ms
3,580 KB
testcase_18 AC 2 ms
3,472 KB
testcase_19 AC 1 ms
3,472 KB
testcase_20 AC 2 ms
3,524 KB
testcase_21 AC 2 ms
3,424 KB
testcase_22 AC 2 ms
3,576 KB
testcase_23 AC 1 ms
3,568 KB
testcase_24 AC 2 ms
3,484 KB
testcase_25 AC 2 ms
3,516 KB
testcase_26 AC 1 ms
3,484 KB
testcase_27 AC 1 ms
3,420 KB
testcase_28 AC 1 ms
3,520 KB
testcase_29 AC 2 ms
3,420 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#pragma GCC optimize("O3")
#pragma GCC optimize("unroll-loops")
#pragma GCC optimize("inline")

#include <cmath>
#include <string>
#include <array>
#include <vector>
#include <cstdlib>
#include <sstream>
#include <iostream>
#include <algorithm>
#include <limits>
#include <queue>
#include <set>
#include <unordered_set>
#include <map>
#include <unordered_map>
#include <bitset>
#include <chrono>
using namespace std;

template<class F, class S> string in_v_to_str (const pair<F, S> v);
template<class F, class S> string v_to_str (const pair<F, S> v);
string in_v_to_str (const char v) { return "'" + string{v} + "'"; }
string in_v_to_str (const char* v) { return "\"" + string(v) + "\""; }
string in_v_to_str (const string v) { return "\"" + v + "\""; }
template<class T> string in_v_to_str (const T v) { stringstream ss; ss << v; return ss.str(); }
template<class T> string v_to_str (const T v) { stringstream ss; ss << v; return ss.str(); }
template<class T, size_t N> string v_to_str (const array<T, N>& v) { stringstream ss; if (v.size() > 0) { ss << "["; for (size_t i = 0; i < v.size() - 1; ++i) { ss << in_v_to_str(v[i]) << ", "; } ss << in_v_to_str(v[v.size() - 1]) << "]"; } else { ss << "[]"; } return ss.str(); }
template<class T, size_t N> string v_to_str (const array< array<T, N>, N >& v) { stringstream ss; if (v.size() > 0) { ss << "["; for (size_t i = 0; i < v.size() - 1; ++i) { ss << v_to_str(v[i]) << ", "; } ss << v_to_str(v[v.size() - 1]) << "]"; } else { ss << "[-]"; } return ss.str(); }
template<class T> string v_to_str (const vector<T>& v) { stringstream ss; if (v.size() > 0) { ss << "["; for (size_t i = 0; i < v.size() - 1; ++i) { ss << in_v_to_str(v[i]) << ", "; } ss << in_v_to_str(v[v.size() - 1]) << "]"; } else { ss << "[]"; } return ss.str(); }
template<class T> string v_to_str (const vector< vector<T> >& v) { stringstream ss; if (v.size() > 0) { ss << "["; for (size_t i = 0; i < v.size() - 1; ++i) { ss << v_to_str(v[i]) << ", "; } ss << v_to_str(v[v.size() - 1]) << "]"; } else { ss << "[-]"; } return ss.str(); }
template<class T> string v_to_str (const set<T>& v) { stringstream ss; int len = v.size(); ss << (v.size() > 0 ? "{" : "{}"); for (auto& i : v) { ss << in_v_to_str(i) << (len-- > 1 ? ", " : "}"); } return ss.str(); }
template<class K, class V> string v_to_str (const map<K, V>& v) { stringstream ss; int len = v.size(); ss << (v.size() > 0 ? "{" : "{}"); for (auto& i : v) { ss << in_v_to_str(i.first) << " : " << in_v_to_str(i.second) << (len-- > 1 ? ", " : "}"); } return ss.str(); }
template<class T> string v_to_str (const unordered_set<T>& v) { stringstream ss; int len = v.size(); ss << (v.size() > 0 ? "{" : "{}"); for (auto& i : v) { ss << in_v_to_str(i) << (len-- > 1 ? ", " : "}"); } return ss.str(); }
template<class K, class V> string v_to_str (const unordered_map<K, V>& v) { stringstream ss; int len = v.size(); ss << (v.size() > 0 ? "{" : "{}"); for (auto& i : v) { ss << in_v_to_str(i.first) << " : " << in_v_to_str(i.second) << (len-- > 1 ? ", " : "}"); } return ss.str(); }
template<class F, class S> string in_v_to_str (const pair<F, S> v) { stringstream ss; ss << "<" << v_to_str(v.first) << ", " << v_to_str(v.second) << ">"; return ss.str(); }
template<class F, class S> string v_to_str (const pair<F, S> v) { stringstream ss; ss << "<" << v_to_str(v.first) << ", " << v_to_str(v.second) << ">"; return ss.str(); }
string print () { return ""; }
template<typename F, typename... R> string print (const F& f, const R& ...r) { stringstream ss; ss << v_to_str(f); if (sizeof...(r) > 0) { ss << " " << print(r...); } return ss.str(); }

template<typename F, typename... R> void pdebug (const F& f, const R& ...r) {
    stringstream ss;
    ss << v_to_str(f); if (sizeof...(r) > 0) { ss << " " << print(r...); }
    cerr << "<cerr>" << ss.str() << "</cerr>" << endl;
}

template<typename F, typename... R> void fdebug (const F& f, const R& ...r) {
    stringstream ss;
    ss << v_to_str(f); if (sizeof...(r) > 0) { ss << " " << print(r...); }
    cerr << "<cerrfile>" << ss.str() << "</cerrfile>" << endl;
}

template<typename F, typename... R> void tdebug (const F& f, const R& ...r) {
    stringstream ss;
    ss << v_to_str(f); if (sizeof...(r) > 0) { ss << " " << print(r...); }
    cerr << "<cerr>[time]" << ss.str() << "</cerr>" << endl;
}



struct e512pos {
public:
    int x;
    int y;
    e512pos () { this->x = 0; this->y = 0; }
    e512pos (int x, int y) {
        this->x = x;
        this->y = y;
    }
    e512pos operator + (const e512pos& t) { return e512pos(this->x + t.x, this->y + t.y); }
    e512pos operator - (const e512pos& t) { return e512pos(this->x - t.x, this->y - t.y); }
    bool operator == (const e512pos& t) const { return this->x == t.x && this->y == t.y; }
};
namespace std {
    template <> class hash<e512pos> {
    public:
        size_t operator()(const e512pos& t) const{ return t.x<<16 | t.y; }
    };
}
ostream& operator << (ostream& os, const e512pos& p) {
    os << "(" << p.x << ", " << p.y << ")";
    return os;
};


class StopWatch {
public:
    std::chrono::system_clock::time_point start, tstart, end;
    StopWatch () {
        this->start = std::chrono::system_clock::now();
    }
    inline void stop () {
        this->tstart = std::chrono::system_clock::now();
    }
    inline void resume () {
        this->start += std::chrono::system_clock::now() - this->tstart;
    }
    inline uint64_t get_milli_time () {
        this->end = std::chrono::system_clock::now();
        return std::chrono::duration_cast<std::chrono::milliseconds>(end-start).count();
    }
};



inline uint32_t xrnd() {
    static uint32_t y = 2463534242;
    y = y ^ (y << 13);
    y = y ^ (y >> 17);
    return y = y ^ (y << 5);
}


inline float distance (const e512pos& a, const e512pos& b) {
    return sqrt((a.x - b.x) * (a.x - b.x) + (a.y - b.y) * (a.y - b.y));
}

class MM {
public:
    StopWatch sw;
    int N, M;
    vector<e512pos> planet;
    vector< pair<int, int> > best;
    MM () { this->sw = StopWatch(); }
    void input () {
        cin >> this->N >> this->M;
        
        for (int i = 0; i < 256; ++i) {
            this->planet.emplace_back(0, 0);
        }
        
        for (int i = 0; i < this->N; ++i) {
            cin >> this->planet[i+1].x >> this->planet[i+1].y;
        }
        
        
    }
    
    
    void output () {
        for (int i = 0; i < this->M; ++i) {
            cout << print(0, i+1) << endl;
        }
        
        
        cout << this->best.size() + 2 << endl;
        cout << print(1, 1) << endl;
        
        for (int i = 0; i < this->best.size(); ++i) {
            cout << print(this->best[i].first, this->best[i].second) << endl;
        }
        
        cout << print(1, 1) << endl;
        
        pdebug((int)(1000000000 / (1000+sqrt(this->allscore(this->best)))));
    }
    void solve () {
        vector< pair<int, int> > v;
        
        
        e512pos p(0, 0);
        bitset<256> use;
        
        for (int z = 0; z < this->N; ++z) {
            int best = -1;
            float bestdist = 1000000000;
            for (int i = 1; i < this->N+1; ++i) {
                if (use[i]) { continue; }
                int x = (p.x - this->planet[i].x);
                int y = (p.y - this->planet[i].y);
                float dist = sqrt(x*x+y*y);
                
                if (dist < bestdist) {
                    bestdist = dist;
                    best = i;
                }
                
                
            }
            use[best] = true;
            v.emplace_back(1, best);
            p = this->planet[best];
            
        }
        
        
        
        
        this->best = v;
    }
    
    float allscore (vector< pair<int, int> >& v) {
        static const int mat[3][3] = {
            {0, 0, 0},
            {0, 25, 5},
            {0, 5, 0},
        };
        static const int tt[3] = {0, 0, 105};
        float r = 0;
        int p = 0;
        int t = 1;
        for (int i = 0; i < this->N; ++i) {
            
            float d = distance(this->planet[p+tt[t]], this->planet[v[i].second + tt[v[i].first]]);
            
            r += d*d*mat[t][v[i].first];
            t = v[i].first;
            p = v[i].second;
        }
        float d = distance(this->planet[p+tt[t]], this->planet[0]);
        r += d*d*mat[t][0];
        
        return r;
    }
    
    
};

int main () {
    cin.tie(0);
    ios::sync_with_stdio(false);
    
    MM mm;
    mm.input();
    mm.solve();
    mm.output();
    cout.flush();
    return 0;
}
0