結果

問題 No.5020 Averaging
ユーザー FplusFplusFFplusFplusF
提出日時 2024-02-25 15:24:39
言語 C++23
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 952 ms / 1,000 ms
コード長 4,504 bytes
コンパイル時間 3,543 ms
コンパイル使用メモリ 268,216 KB
実行使用メモリ 6,548 KB
スコア 51,553,507
最終ジャッジ日時 2024-02-25 15:25:52
合計ジャッジ時間 53,596 ms
ジャッジサーバーID
（参考情報） 		judge15 / judge14
純コード判定しない問題か言語
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 951 ms
6,548 KB
testcase_01 AC 952 ms
6,548 KB
testcase_02 AC 952 ms
6,548 KB
testcase_03 AC 952 ms
6,548 KB
testcase_04 AC 952 ms
6,548 KB
testcase_05 AC 952 ms
6,548 KB
testcase_06 AC 952 ms
6,548 KB
testcase_07 AC 952 ms
6,548 KB
testcase_08 AC 952 ms
6,548 KB
testcase_09 AC 951 ms
6,548 KB
testcase_10 AC 952 ms
6,548 KB
testcase_11 AC 952 ms
6,548 KB
testcase_12 AC 952 ms
6,548 KB
testcase_13 AC 952 ms
6,548 KB
testcase_14 AC 952 ms
6,548 KB
testcase_15 AC 952 ms
6,548 KB
testcase_16 AC 952 ms
6,548 KB
testcase_17 AC 952 ms
6,548 KB
testcase_18 AC 952 ms
6,548 KB
testcase_19 AC 951 ms
6,548 KB
testcase_20 AC 952 ms
6,548 KB
testcase_21 AC 952 ms
6,548 KB
testcase_22 AC 951 ms
6,548 KB
testcase_23 AC 952 ms
6,548 KB
testcase_24 AC 952 ms
6,548 KB
testcase_25 AC 952 ms
6,548 KB
testcase_26 AC 952 ms
6,548 KB
testcase_27 AC 952 ms
6,548 KB
testcase_28 AC 952 ms
6,548 KB
testcase_29 AC 951 ms
6,548 KB
testcase_30 AC 952 ms
6,548 KB
testcase_31 AC 952 ms
6,548 KB
testcase_32 AC 951 ms
6,548 KB
testcase_33 AC 952 ms
6,548 KB
testcase_34 AC 952 ms
6,548 KB
testcase_35 AC 952 ms
6,548 KB
testcase_36 AC 951 ms
6,548 KB
testcase_37 AC 952 ms
6,548 KB
testcase_38 AC 952 ms
6,548 KB
testcase_39 AC 951 ms
6,548 KB
testcase_40 AC 952 ms
6,548 KB
testcase_41 AC 952 ms
6,548 KB
testcase_42 AC 952 ms
6,548 KB
testcase_43 AC 951 ms
6,548 KB
testcase_44 AC 952 ms
6,548 KB
testcase_45 AC 952 ms
6,548 KB
testcase_46 AC 952 ms
6,548 KB
testcase_47 AC 952 ms
6,548 KB
testcase_48 AC 952 ms
6,548 KB
testcase_49 AC 951 ms
6,548 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>
using namespace std;
using pii=pair<int,int>;
using tii=tuple<int,int,int>;
using qii=tuple<int,int,int,int>;
using ll=long long;
using ld=long double;
const ll INF=1ll<<60;
#define rep(i,n) for (int i=0;i<(int)(n);i++)
#define all(v) v.begin(),v.end()
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;
}
auto program_start=chrono::system_clock::now();
mt19937 mt;
int rand_int(int r){    //[0,r)
    return mt()%r;
}
int rand_int(int l,int r){  //[l,r)
    return l+mt()%(r-l);
}
double rand_double(){   //[0.0,1.0]
    return (double)mt()/(double)mt19937::max();
}
const ll C=5e17;
int N;
vector<ll> A,B;
namespace Solver{
    ll fake_score(vector<int> &v){
        __int128 a_sum=0,b_sum=0;
        rep(i,N){
            a_sum+=(__int128)A[i]*(1ll<<v[i]);
            b_sum+=(__int128)B[i]*(1ll<<v[i]);
        }
        a_sum/=(1ll<<(N-1));
        b_sum/=(1ll<<(N-1));
        return floorl(2000000.0-100000.0*log10l(max(abs(C-(ll)a_sum),abs(C-(ll)b_sum))+1));
    }
    int score(vector<pii> &ans){
        vector<ll> a=A,b=B;
        for(auto [u,v]:ans){
            a[v]=a[u]=(a[u]+a[v])/2;
            b[v]=b[u]=(b[u]+b[v])/2;
        }
        return floorl(2000000.0-100000.0*log10l(max(abs(C-a[0]),abs(C-b[0]))+1));
    }
    void solve(){
        vector<int> k;
        rep(i,N-1){
            k.push_back(i);
            if(i==0) k.push_back(i);
        }
        int pre_fake_score=fake_score(k);
        while(true){
            auto now=chrono::system_clock::now();
            int ms=chrono::duration_cast<chrono::milliseconds>(now-program_start).count();
            if(50<=ms) break;
            int x=rand_int(N),y=rand_int(N);
            if(x==y) continue;
            swap(k[x],k[y]);
            int new_fake_score=fake_score(k);
            if(new_fake_score<pre_fake_score) swap(k[x],k[y]);
            else pre_fake_score=new_fake_score;
        }
        vector<int> ord(N);
        rep(i,N) ord[i]=i;
        sort(all(ord),[&](int i,int j){return k[i]<k[j];});
        vector<pii> start_ans;
        bool zero=false;
        rep(i,N){
            if(zero&&ord[i]==0) continue;
            if(zero||i==N-1) start_ans.emplace_back(0,ord[i]);
            else{
                start_ans.emplace_back(ord[i],ord[i+1]);
                if(ord[i+1]==0) zero=true;
            }
        }
        vector<pii> ans=start_ans;
        int pre_score=score(ans);
        int cnt=0;
        int max_score=0;
        vector<pii> max_score_ans;
        while(true){
            auto now=chrono::system_clock::now();
            int ms=chrono::duration_cast<chrono::milliseconds>(now-program_start).count();
            if(950<=ms) break;
            vector<pii> new_ans=ans;
            int r=mt()%4;
            int x=rand_int(N),y=rand_int(N);
            int sz=new_ans.size();
            if(x==y) continue;
            if(y<x) swap(x,y);
            if(sz==0||(sz<50&&r==0)){
                int d=rand_int(sz+1);
                new_ans.insert(new_ans.begin()+d,{x,y});
            }else if(r==1){
                int d=rand_int(sz);
                new_ans.erase(new_ans.begin()+d);
            }else if(r==2){
                int p=rand_int(sz),q=rand_int(sz);
                if(p==q) continue;
                if(q<p) swap(p,q);
                reverse(new_ans.begin()+p,new_ans.begin()+q+1);
            }else{
                int d=rand_int(sz);
                new_ans[d]={x,y};
            }
            cnt++;
            int new_score=score(new_ans);
            if(pre_score<=new_score){
                if(pre_score<new_score){
                    cnt=0;
                    if(chmax(max_score,new_score)){
                        max_score_ans=new_ans;
                    }
                }
                pre_score=new_score;
                ans=new_ans;
            }
            if(cnt==100000){
                ans=start_ans;
                pre_score=score(ans);
                cnt=0;
            }
        }
        cerr << max_score << endl;
        cout << (int)max_score_ans.size() << '\n';
        for(auto [u,v]:max_score_ans) cout << u+1 << ' ' << v+1 << '\n';
    }
}

int main(){
    ios::sync_with_stdio(false);
    cin.tie(nullptr);
    cin >> N;
    A.resize(N);
    B.resize(N);
    rep(i,N) cin >> A[i] >> B[i];
    Solver::solve();
}
0