#756739 (C++17) No.5005 3-SAT

提出ソース

結果

問題	No.5005 3-SAT
ユーザー	umezo
提出日時	2022-04-29 17:58:42
言語	C++17 (gcc 12.3.0 + boost 1.83.0)
結果	AC
実行時間	1,928 ms / 2,000 ms
コード長	3,377 bytes
コンパイル時間	2,979 ms
実行使用メモリ	333,648 KB
スコア	87,361
最終ジャッジ日時	2022-04-29 18:02:06
合計ジャッジ時間	201,499 ms
ジャッジサーバーID （参考情報）	judge15 / judge13
	純コード判定しない問題か言語

このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示

入力	結果	実行時間実行使用メモリ
testcase_00	AC	1,919 ms 276,824 KB
testcase_01	AC	1,918 ms 276,184 KB
testcase_02	AC	1,923 ms 278,284 KB
testcase_03	AC	1,921 ms 280,500 KB
testcase_04	AC	1,918 ms 280,688 KB
testcase_05	AC	1,922 ms 283,308 KB
testcase_06	AC	1,922 ms 279,796 KB
testcase_07	AC	1,922 ms 283,732 KB
testcase_08	AC	1,923 ms 277,592 KB
testcase_09	AC	1,926 ms 285,260 KB
testcase_10	AC	1,922 ms 280,156 KB
testcase_11	AC	1,921 ms 278,944 KB
testcase_12	AC	1,922 ms 273,712 KB
testcase_13	AC	1,918 ms 282,620 KB
testcase_14	AC	1,918 ms 279,976 KB
testcase_15	AC	1,920 ms 285,300 KB
testcase_16	AC	1,922 ms 283,816 KB
testcase_17	AC	1,922 ms 282,544 KB
testcase_18	AC	1,915 ms 276,900 KB
testcase_19	AC	1,916 ms 277,776 KB
testcase_20	AC	1,919 ms 281,664 KB
testcase_21	AC	1,927 ms 286,988 KB
testcase_22	AC	1,923 ms 282,384 KB
testcase_23	AC	1,922 ms 282,108 KB
testcase_24	AC	1,917 ms 283,936 KB
testcase_25	AC	1,917 ms 282,888 KB
testcase_26	AC	1,926 ms 284,196 KB
testcase_27	AC	1,921 ms 282,864 KB
testcase_28	AC	1,924 ms 276,832 KB
testcase_29	AC	1,916 ms 275,776 KB
testcase_30	AC	1,923 ms 282,924 KB
testcase_31	AC	1,919 ms 282,312 KB
testcase_32	AC	1,928 ms 284,800 KB
testcase_33	AC	1,921 ms 279,964 KB
testcase_34	AC	1,917 ms 274,312 KB
testcase_35	AC	1,919 ms 283,388 KB
testcase_36	AC	1,917 ms 289,448 KB
testcase_37	AC	1,916 ms 280,984 KB
testcase_38	AC	1,919 ms 286,740 KB
testcase_39	AC	1,915 ms 276,732 KB
testcase_40	AC	1,917 ms 286,060 KB
testcase_41	AC	1,919 ms 283,624 KB
testcase_42	AC	1,920 ms 283,440 KB
testcase_43	AC	1,922 ms 276,384 KB
testcase_44	AC	1,919 ms 288,152 KB
testcase_45	AC	1,918 ms 278,644 KB
testcase_46	AC	1,918 ms 283,544 KB
testcase_47	AC	1,919 ms 284,288 KB
testcase_48	AC	1,916 ms 278,988 KB
testcase_49	AC	1,916 ms 276,588 KB
testcase_50	AC	1,919 ms 281,592 KB
testcase_51	AC	1,919 ms 285,900 KB
testcase_52	AC	1,922 ms 292,264 KB
testcase_53	AC	1,924 ms 279,712 KB
testcase_54	AC	1,924 ms 333,648 KB
testcase_55	AC	1,921 ms 277,400 KB
testcase_56	AC	1,921 ms 280,204 KB
testcase_57	AC	1,920 ms 282,600 KB
testcase_58	AC	1,921 ms 281,392 KB
testcase_59	AC	1,924 ms 283,584 KB
testcase_60	AC	1,920 ms 291,396 KB
testcase_61	AC	1,921 ms 284,296 KB
testcase_62	AC	1,917 ms 280,292 KB
testcase_63	AC	1,917 ms 280,204 KB
testcase_64	AC	1,917 ms 280,500 KB
testcase_65	AC	1,919 ms 288,328 KB
testcase_66	AC	1,918 ms 281,120 KB
testcase_67	AC	1,914 ms 277,452 KB
testcase_68	AC	1,918 ms 286,436 KB
testcase_69	AC	1,915 ms 276,124 KB
testcase_70	AC	1,917 ms 280,620 KB
testcase_71	AC	1,917 ms 286,684 KB
testcase_72	AC	1,916 ms 279,256 KB
testcase_73	AC	1,918 ms 280,512 KB
testcase_74	AC	1,916 ms 281,344 KB
testcase_75	AC	1,917 ms 280,336 KB
testcase_76	AC	1,920 ms 291,620 KB
testcase_77	AC	1,917 ms 287,808 KB
testcase_78	AC	1,920 ms 287,752 KB
testcase_79	AC	1,918 ms 279,272 KB
testcase_80	AC	1,918 ms 281,080 KB
testcase_81	AC	1,918 ms 282,728 KB
testcase_82	AC	1,918 ms 289,240 KB
testcase_83	AC	1,919 ms 284,920 KB
testcase_84	AC	1,921 ms 277,260 KB
testcase_85	AC	1,925 ms 282,316 KB
testcase_86	AC	1,922 ms 287,352 KB
testcase_87	AC	1,924 ms 288,448 KB
testcase_88	AC	1,922 ms 280,232 KB
testcase_89	AC	1,922 ms 288,304 KB
testcase_90	AC	1,923 ms 283,924 KB
testcase_91	AC	1,923 ms 285,048 KB
testcase_92	AC	1,921 ms 284,712 KB
testcase_93	AC	1,922 ms 285,620 KB
testcase_94	AC	1,919 ms 286,876 KB
testcase_95	AC	1,922 ms 280,024 KB
testcase_96	AC	1,922 ms 291,396 KB
testcase_97	AC	1,921 ms 270,380 KB
testcase_98	AC	1,921 ms 288,408 KB
testcase_99	AC	1,918 ms 287,600 KB

権限があれば一括ダウンロードができます

ソースコード

raw source code

#define rep(i, n) for (int i = 0; i < (int)(n); i++)
#define ALL(v) v.begin(), v.end()
typedef long long ll;

#include <bits/stdc++.h>
using namespace std;
using namespace chrono;

const int MAX=2048;

static uint32_t randXor(){
  static uint32_t x=123456789;
  static uint32_t y=362436069;
  static uint32_t z=521288629;
  static uint32_t w=88675123;
  uint32_t t;
 
  t=x^(x<<11);
  x=y;y=z;z=w;
  return w=(w^(w>>19))^(t^(t>>8));
}

static double rand01(){
    return (randXor()+0.5)*(1.0/UINT_MAX);
}

int calc(vector<int>&A,vector<int>&B,vector<int>&C,vector<int>&P,vector<int>&Q,vector<int>&R,vector<int>&S){
  int num=0;
  rep(i,MAX){
    bool b=false;
    if(S[A[i]]==P[i]) b=true;
    if(S[B[i]]==Q[i]) b=true;
    if(S[C[i]]==R[i]) b=true;
    if(b) num++;
    else break;
  }
  return num;
}

map<pair<vector<int>,int>,int> m;

int main(){
  ios::sync_with_stdio(false);
  std::cin.tie(nullptr);
  
  vector<int> A(MAX),B(MAX),C(MAX),P(MAX),Q(MAX),R(MAX);
  vector<int> A0,A1,B0,B1,C0,C1;
  rep(i,MAX){
    cin>>A[i]>>B[i]>>C[i]>>P[i]>>Q[i]>>R[i];
    if(P[i]==0) A0.push_back(i);
    else A1.push_back(i);
    if(Q[i]==0) B0.push_back(i);
    else B1.push_back(i);
    if(R[i]==0) C0.push_back(i);
    else C1.push_back(i);
  }
  A0.push_back(MAX),A1.push_back(MAX),B0.push_back(MAX);
  B1.push_back(MAX),C0.push_back(MAX),C1.push_back(MAX);
  
  auto startClock=system_clock::now();
  int bestscore=0;
  vector<int> S(256),bestS(256);
  
  const static double START_TEMP=1500;
  const static double END_TEMP  = 100;
  const static double END_TIME=1.85;
 
  double time=0.0;
  int cnt=0;
  do{
    cnt++;
    const double progressRatio=time/END_TIME;
    const double temp=START_TEMP+(END_TEMP-START_TEMP)*progressRatio;
    
    int score=-1;
    
    rep(i,2048){
      if(S[A[i]]==P[i] || S[B[i]]==Q[i] || S[C[i]]==R[i]) continue;
      else{
        score=i;
        break;
      }
    }
    
    if(score>bestscore){
      bestscore=score;
      bestS=S;
    }
    
    ll t=randXor();
    if(m[{S,0}]==m[{S,1}] && m[{S,0}]==m[{S,2}]){
      if(t%3==0){
        m[{S,0}]++;
        S[A[score]]^=1;
      }
      else if(t%3==1){
        m[{S,1}]++;
        S[B[score]]^=1;
      }
      else{
        m[{S,2}]++;
        S[C[score]]^=1;
      }
    }
    else if(m[{S,0}]==m[{S,1}] && m[{S,0}]<m[{S,2}]){
      if(t%2==0){
        m[{S,0}]++;
        S[A[score]]^=1;
      }
      else{
        m[{S,1}]++;
        S[B[score]]^=1;
      }
    }
    else if(m[{S,1}]==m[{S,2}] && m[{S,1}]<m[{S,0}]){
      if(t%2==0){
        m[{S,1}]++;
        S[B[score]]^=1;
      }
      else{
        m[{S,2}]++;
        S[C[score]]^=1;
      }
    }  
    else if(m[{S,2}]==m[{S,0}] && m[{S,2}]<m[{S,1}]){
      if(t%2==0){
        m[{S,2}]++;
        S[C[score]]^=1;
      }
      else{
        m[{S,0}]++;
        S[A[score]]^=1;
      }
    }
    else if(m[{S,0}]<m[{S,1}] && m[{S,0}]<m[{S,2}]){
      m[{S,0}]++;
      S[A[score]]^=1;
    }
    else if(m[{S,1}]<m[{S,2}] && m[{S,1}]<m[{S,0}]){
      m[{S,1}]++;
      S[B[score]]^=1;
    }
    else if(m[{S,2}]<m[{S,0}] && m[{S,2}]<m[{S,1}]){
      m[{S,2}]++;
      S[C[score]]^=1;
    }
    
    time=(duration_cast<microseconds>(system_clock::now()-startClock).count()*1e-6);
  }while(time<END_TIME);

  for(int i=255;i>=0;i--) cout<<bestS[i];
  cout<<endl;
  
  cout<<calc(A,B,C,P,Q,R,bestS)<<endl;

  return 0;
}

yukicoder

結果

テストケース

ソースコード