結果

問題 No.5007 Steiner Space Travel
コンテスト
ユーザー tanaka
提出日時 2022-07-30 16:54:05
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 993 ms / 1,000 ms
コード長 4,757 bytes
コンパイル時間 1,556 ms
実行使用メモリ 5,164 KB
スコア 5,683,013
最終ジャッジ日時 2022-07-30 16:54:47
合計ジャッジ時間 34,054 ms
ジャッジサーバーID
（参考情報） 		judge13 / judge11
純コード判定しない問題か言語
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 30
権限があれば一括ダウンロードができます
コンパイルメッセージ
main.cpp: 関数 ‘int main()’ 内:
main.cpp:181:24: 警告: ‘dy’ はこの関数内初期化されずに使用されるかもしれません [-Wmaybe-uninitialized]
  181 |                 d[idx1]-=dy;
      |                 ~~~~~~~^~~~
main.cpp:180:24: 警告: ‘dx’ はこの関数内初期化されずに使用されるかもしれません [-Wmaybe-uninitialized]
  180 |                 c[idx1]-=dx;
      |                 ~~~~~~~^~~~

ソースコード

diff # 

#include <bits/stdc++.h>

#define rep(i, n) for (int i = 0; i < (int)(n); i++)
#define rep3(i, n, m) for (int i = m; i < (int)(n); i++)
#define repr(i, n) for (int i = n-1; i >= 0; i--)
#define all(v) v.begin(), v.end()
#define SP << " " 
#define ll long long

using namespace std;

//定数周り
int imax=2147483647;
ll lmax=9223372036854775807;

//焼きなましの定数
double TIME_LIMIT=990;
double start_temp=2000000.0;
double end_temp=10000.0;

//入力など
int n, m;
int a[100], b[100];

// 変数
int alp=5, sz=0;
int c[8], d[8];
vector<int> r(99);
vector<vector<int>> s(100);
ll b_sco, n_sco;

//乱数の準備
//int seed=1;
auto seed=(unsigned)time(NULL);
mt19937_64 mt(seed);

void inpt(){
    cin>> n >> m;
    rep(i, 100){
        cin>> a[i] >> b[i];
        // a[i]=mt()%1001;
        // b[i]=mt()%1001;
    }
}

void init(){
    rep(i, r.size()){
        r[i]=i+2;
    }
    shuffle(all(r), mt);
    r.push_back(1);
    rep(i, 8){
        c[i]=mt()%1001;
        d[i]=mt()%1001;
    }
    //c={}
}

ll score(){
    ll sco=0;
    int px=a[0], py=b[0];
    int nx, ny;
    rep(i, 100){
        int tmp=alp*alp;
        rep(j, s[i].size()){
            nx=c[s[i][j]-1];
            ny=d[s[i][j]-1];
            tmp=1;
            if(j==0 || j==s[i].size()-1) tmp=alp;
            sco+=tmp*((px-nx)*(px-nx)+(py-ny)*(py-ny));
            px=nx;
            py=ny;
        }
        nx=a[r[i]-1];
        ny=b[r[i]-1];
        sco+=tmp*((px-nx)*(px-nx)+(py-ny)*(py-ny));
        px=nx;
        py=ny;
    }
    return sco;
}

void print_ans(){
    rep(i, 8) cout<< c[i] SP << d[i] <<endl;
    cout<< r.size()+sz+1 <<endl;
    cout<< 1 SP << 1 <<endl;
    rep(i, r.size()){
        rep(j, s[i].size()) cout<< 2 SP << s[i][j] <<endl;
        cout<< 1 SP << r[i] <<endl;
    }
}

int main(){

    //開始時間の計測
    std::chrono::system_clock::time_point start, current;
    start = chrono::system_clock::now();

    //入力
    inpt();
    //初期状態の計算
    init();
    // print_ans();
    b_sco=score();

    //焼きなまし
    int lp=0;
    while (true) { // 時間の許す限り回す
        lp++;
        //cout<< lp <<endl;
        current = chrono::system_clock::now(); // 現在時刻
        if (chrono::duration_cast<chrono::milliseconds>(current - start).count() > TIME_LIMIT) break;
        //if(lp==20000) break;

        int type=mt()%4+1;
        int idx1, idx2, dx, dy;
        vector<int> pre_s;
        //cout<< "type:" << type <<endl;
        if(type==1){
            // 惑星訪問の順番入れ替え
            idx1=mt()%99;
            while(1){
                idx2=(idx1+mt()%99)%100;
                if(idx2!=99) break;
            }
            swap(r[idx1], r[idx2]);
        }else if(type==2){
            // ステーション移動
            idx1=mt()%m;
            while(1){
                dx=mt()%21-10;
                dy=mt()%21-10;
                if(0<=c[idx1]+dx && c[idx1]+dx<=1000 && 0<=d[idx1]+dy && d[idx1]+dy<=1000) break;
            }
        }else if(type==3){
            // ステーション挿入
            idx1=mt()%m+1;
            idx2=mt()%99;
            pre_s=s[idx2];
            int flag=1;
            rep(i, s[idx2].size()){
                if(s[idx2][i]==idx1){
                    flag=0;
                    break;
                }
            }
            if(flag){
                int ins=0;
                if(s[idx2].size()>0) ins=mt()%s[idx2].size();
                s[idx2].insert(s[idx2].begin()+ins, idx1);
                sz++;
            }
        }else{
            // ステーション削除
            idx2=mt()%99;
            pre_s=s[idx2];
            if(s[idx2].size()>0){
                idx1=mt()%s[idx2].size();
                s[idx2].erase(s[idx2].begin()+idx1);
                sz--;
            }
        }
        n_sco=score();
        //cout<< n_sco <<endl;

        // 温度関数
        double temp = start_temp + (end_temp - start_temp) * chrono::duration_cast<chrono::milliseconds>(current - start).count() / TIME_LIMIT;
        // 遷移確率関数(最大化の場合)
        double prob = exp((b_sco-n_sco)/temp);

        if (prob > (mt()%imax)/(double)imax) { // 確率probで遷移する
        //if (nsco<psco) {
            b_sco=n_sco;
            // cout<< "lp:" << lp <<endl;
            // cout<< n_sco <<endl;
        }else{
            if(type==1){
                swap(r[idx1], r[idx2]);
            }else if(type==2){
                c[idx1]-=dx;
                d[idx1]-=dy;
            }else if(type==3){
                s[idx2]=pre_s;
                sz--;
            }else{
                s[idx2]=pre_s;
                sz++;
            }
        }
    }

    print_ans();

	return 0;
}
0