結果

問題 No.5008 [Cherry Alpha] Discrete Pendulum with Air Resistance
ユーザー butsurizukibutsurizuki
提出日時 2022-10-03 15:47:33
言語 C++23
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 1,957 ms / 2,000 ms
コード長 2,995 bytes
コンパイル時間 4,124 ms
実行使用メモリ 6,956 KB
スコア 1,463,111,598,896,013
最終ジャッジ日時 2022-10-14 21:38:36
合計ジャッジ時間 104,154 ms
ジャッジサーバーID
（参考情報） 		judge12 / judge15
純コード判定しない問題か言語
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 1,745 ms
4,904 KB
testcase_01 AC 1,690 ms
4,904 KB
testcase_02 AC 1,865 ms
6,952 KB
testcase_03 AC 1,924 ms
4,904 KB
testcase_04 AC 1,858 ms
4,900 KB
testcase_05 AC 1,926 ms
4,904 KB
testcase_06 AC 1,891 ms
4,900 KB
testcase_07 AC 1,859 ms
4,900 KB
testcase_08 AC 1,794 ms
4,900 KB
testcase_09 AC 1,876 ms
4,900 KB
testcase_10 AC 1,861 ms
4,904 KB
testcase_11 AC 1,723 ms
4,900 KB
testcase_12 AC 1,744 ms
4,904 KB
testcase_13 AC 1,670 ms
4,900 KB
testcase_14 AC 1,779 ms
4,900 KB
testcase_15 AC 1,844 ms
4,900 KB
testcase_16 AC 1,829 ms
6,948 KB
testcase_17 AC 1,892 ms
4,900 KB
testcase_18 AC 1,797 ms
6,948 KB
testcase_19 AC 1,803 ms
4,904 KB
testcase_20 AC 1,915 ms
4,900 KB
testcase_21 AC 1,750 ms
6,952 KB
testcase_22 AC 1,818 ms
4,904 KB
testcase_23 AC 1,817 ms
4,904 KB
testcase_24 AC 1,686 ms
6,948 KB
testcase_25 AC 1,828 ms
4,900 KB
testcase_26 AC 1,895 ms
4,908 KB
testcase_27 AC 1,950 ms
6,952 KB
testcase_28 AC 1,731 ms
4,900 KB
testcase_29 AC 1,796 ms
6,956 KB
testcase_30 AC 1,505 ms
4,904 KB
testcase_31 AC 1,774 ms
4,900 KB
testcase_32 AC 1,780 ms
4,900 KB
testcase_33 AC 1,702 ms
4,904 KB
testcase_34 AC 1,841 ms
4,908 KB
testcase_35 AC 1,788 ms
4,904 KB
testcase_36 AC 1,736 ms
4,904 KB
testcase_37 AC 1,907 ms
4,904 KB
testcase_38 AC 1,886 ms
4,900 KB
testcase_39 AC 1,536 ms
4,900 KB
testcase_40 AC 1,872 ms
6,952 KB
testcase_41 AC 1,957 ms
4,900 KB
testcase_42 AC 1,565 ms
4,900 KB
testcase_43 AC 1,840 ms
6,948 KB
testcase_44 AC 1,870 ms
4,904 KB
testcase_45 AC 1,858 ms
4,900 KB
testcase_46 AC 1,835 ms
4,904 KB
testcase_47 AC 1,705 ms
4,900 KB
testcase_48 AC 1,726 ms
4,904 KB
testcase_49 AC 1,675 ms
6,952 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#pragma GCC target("avx2")
#pragma GCC optimize("Ofast")
#pragma GCC optimize("unroll-loops")
#include<bits/stdc++.h>

using namespace std;

int getpos(int tim,int b){
  int p=tim/b;
  int q=tim%b;

  int res;
  if(p&1){res=b-q;}
  else{res=q;}

  if(p&2){res*=-1;}

  return res;
}

int n,k;
vector<int> t,u;
int sz;

long long eval(vector<int> &cnt,bool debug){
  // tairitsu
  double t_all=0.0;
  for(auto &nx : t){
    double t_cur=0.0;
    vector<int> xlis(sz+1);
    for(int i=1;i<=sz;i++){
      xlis[i]=getpos(nx,i);
    }
    for(int i=1;i<=sz;i++){
      for(int j=i;j<=sz;j++){
        double del;
        if(i==j){del=(cnt[i]*(cnt[i]-1))/2;}
        else{del=cnt[i]*cnt[j];}
        if(del==0.0){continue;}
        del*=abs(xlis[i]-xlis[j]);
        del/=((double)(i+j));
        t_cur+=del;
      }
    }
    t_cur*=2.0e7;
    t_cur/=((double)(n*(n-1)));
    t_all+=round(t_cur);
    if(debug){cerr << ((int)round(t_cur)) << "\n";}
  }
  t_all/=((double)k);
  if(debug){cerr << "\n";}

  // kyoucyou
  double k_all=0.0;
  for(auto &nx : u){
    double k_cur=0.0;
    vector<int> xlis;
    for(int i=1;i<=sz;i++){
      if(cnt[i]>0){xlis.push_back(getpos(nx,i));}
    }
    sort(xlis.begin(),xlis.end());
    k_cur=xlis.back()-xlis[0];
    k_cur/=20.0;
    k_cur+=1.0;
    k_cur=1.0e7/sqrt(k_cur);
    k_all+=round(k_cur);
    if(debug){cerr << ((int)round(k_cur)) << "\n";}
  }
  k_all/=((double)k);
  if(debug){cerr << "\n";}

  long long rt=(t_all+0.5);
  long long rk=(k_all+0.5);
  if(debug){
    cerr << rt << "\n";
    cerr << rk << "\n";
    cerr << rt*rk << "\n";
  }
  return rt*rk;
}

long long get_rand(long long lim,mt19937_64 &eg){
  return (long long)(eg()%lim);
}

int main(){
  std::random_device seed_gen;
  std::mt19937_64 engine(seed_gen());
  ios::sync_with_stdio(false);
  cin.tie(nullptr);

  cin >> n >> k;
  t.resize(k);
  u.resize(k);
  for(auto &nx : t){cin >> nx;}
  for(auto &nx : u){cin >> nx;}

  sz=15;

  vector<int> res(sz+1,0);
  for(int i=0;i<n;i++){
    res[1+(i%sz)]++;
  }
  long long maxsc=eval(res,false);

  for(int tr=0;tr<150;tr++){
    vector<int> cv(sz+1,0);
    int bas=1+tr%15;
    for(int i=0;i<n;i++){
      cv[1+get_rand(bas,engine)]++;
    }
    long long cursc=eval(cv,false);
    for(int i=0;i<350;i++){
      int p=1+get_rand(bas,engine);
      int q=1+get_rand(bas,engine);
      if(p==q){continue;}
      bool act=false;
      while(cv[p]>0){
        cv[p]--;cv[q]++;
        long long nexsc=eval(cv,false);
        if(cursc<=nexsc){act=true;cursc=nexsc;}
        else{cv[p]++;cv[q]--;break;}
      }
      if(act){continue;}
      while(cv[q]>0){
        cv[p]++;cv[q]--;
        long long nexsc=eval(cv,false);
        if(cursc<=nexsc){act=true;cursc=nexsc;}
        else{cv[p]--;cv[q]++;break;}
      }
    }
    if(maxsc<=cursc){
      maxsc=cursc;
      res=cv;
    }
  }

  for(int i=1;i<=sz;i++){
    for(int j=0;j<res[i];j++){
      cout << i << " " << i << " " << 1 << "\n";
    }
  }
  eval(res,true);
  return 0;
}
0