結果

問題 No.5013 セクスタプル (open)
コンテスト
ユーザー tails
提出日時 2022-12-29 20:02:04
言語 cLay
(20241019-1)
結果
AC  
実行時間 1,982 ms / 2,000 ms
コード長 1,420 bytes
コンパイル時間 2,408 ms
実行使用メモリ 6,952 KB
スコア 18,547
最終ジャッジ日時 2022-12-29 20:05:34
合計ジャッジ時間 208,571 ms
ジャッジサーバーID
（参考情報） 		judge11 / judge13
純コード判定しない問題か言語
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 100
権限があれば一括ダウンロードができます

ソースコード

diff # 

static int f_row(int y,int p[6],int d[36][6]){
	int z=0;
	int v[6]{};
	int n[6]{};
	rep(x,6){
		int i=p[y*6+x];
		rep(k,6){
			int c=d[i][k];
			v[c]|=1<<x;
			n[c]+=1;
		}
	}
	rep(c,6){
		if(v[c]==63){
			z+=n[c]-3;
		}
	}
	return z;
}

static int f_col(int x,int p[6],int d[36][6]){
	int z=0;
	int v[6]{};
	int n[6]{};
	rep(y,6){
		int i=p[y*6+x];
		rep(k,6){
			int c=d[i][k];
			v[c]|=1<<y;
			n[c]+=1;
		}
	}
	rep(c,6){
		if(v[c]==63){
			z+=n[c]-3;
		}
	}
	return z;
}

{
	Timer timer;
	timer.set();
	
	Rand r;
	
	int@--d[36][6];
	int p[36];
	int s=0;
	int s_row[6],s_col[6];
	rep(i,36)p[i]=i;
	rep(i,6)s+=s_row[i]=f_row(i,p,d);
	rep(i,6)s+=s_col[i]=f_col(i,p,d);
	int best_p[36];
	int best_s=-1;
	double a;
	while(a=timer.get()/1.98,a<1.0){
		int i=r.get(36);
		int j=r.get(35);
		j+=j>=i;
		if(1){
			swap(p[i],p[j]);
			int sir=f_row(i/6,p,d);
			int sic=f_col(i%6,p,d);
			int sjr=0;
			int sjc=0;
			int e=sir+sic-s_row[i/6]-s_col[i%6];
			if(i/6!=j/6){
				sjr=f_row(j/6,p,d);
				e+=sjr-s_row[j/6];
			}
			if(i%6!=j%6){
				sjc=f_col(j%6,p,d);
				e+=sjc-s_col[j%6];
			}
			if(r.get(0.,1.)<exp(e*1.3)){
				s_row[j/6]=sjr;
				s_row[i/6]=sir;
				s_col[j%6]=sjc;
				s_col[i%6]=sic;
				s+=e;
				if(s>best_s){
					best_s=s;
					rep(i,36)best_p[i]=p[i];
				}
			}else{
				swap(p[i],p[j]);
			}
		}
	}
	{
		int o[36];
		rep(i,36){
			o[best_p[i]]=i;
		}
		rep(i,36){
			wt(o[i]/6+1,o[i]%6+1);
		}
	}
}
0