結果

問題 No.309 シャイな人たち (1)
ユーザー chocoruskchocorusk
提出日時 2019-12-29 04:18:16
言語 C++14
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 533 ms / 4,000 ms
コード長 2,266 bytes
コンパイル時間 955 ms
コンパイル使用メモリ 111,752 KB
実行使用メモリ 37,048 KB
最終ジャッジ日時 2024-04-21 08:29:54
合計ジャッジ時間 8,097 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge4
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 508 ms
36,812 KB
testcase_01 AC 523 ms
36,864 KB
testcase_02 AC 503 ms
36,980 KB
testcase_03 AC 533 ms
36,864 KB
testcase_04 AC 14 ms
20,608 KB
testcase_05 AC 126 ms
24,704 KB
testcase_06 AC 512 ms
37,028 KB
testcase_07 AC 504 ms
36,840 KB
testcase_08 AC 519 ms
37,048 KB
testcase_09 AC 510 ms
36,892 KB
testcase_10 AC 504 ms
36,992 KB
testcase_11 AC 510 ms
36,892 KB
testcase_12 AC 503 ms
36,716 KB
testcase_13 AC 6 ms
20,352 KB
testcase_14 AC 7 ms
20,608 KB
testcase_15 AC 125 ms
24,704 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <cstdio>
#include <cstring>
#include <iostream>
#include <string>
#include <cmath>
#include <bitset>
#include <vector>
#include <map>
#include <set>
#include <queue>
#include <deque>
#include <algorithm>
#include <complex>
#include <unordered_map>
#include <unordered_set>
#include <random>
#include <cassert>
#include <fstream>
#include <utility>
#include <functional>
#include <time.h>
#include <stack>
#include <array>
#define popcount __builtin_popcount
using namespace std;
typedef long long int ll;
typedef pair<int, int> P;

int main()
{
	int r, c; cin>>r>>c;
	double p[11][11]; int s[11][11]={};
	for(int i=0; i<r; i++) for(int j=0; j<c; j++) cin>>p[i][j];
	for(int i=0; i<r; i++) for(int j=0; j<c; j++) cin>>s[i][j];
	double dp[12][1<<11]={};
	dp[0][0]=1;
	double ans=0;
	double q[11][1<<11];
	for(int i=0; i<r; i++){
		for(int k=0; k<(1<<c); k++){
			q[i][k]=1;
			for(int l=0; l<c; l++){
				if(k&(1<<l)){
					q[i][k]*=p[i][l]/100;
				}else{
					q[i][k]*=(100-p[i][l])/100;
				}
			}
		}
	}
	int nuo[1<<22]={};
	int nue[1<<22];
	for(int i=0; i<(1<<(2*c)); i++){
		int cnt[11];
		bool used[11]={};
		int que[30];
		int t=0, u=0;
		for(int j=0; j<c; j++){
			cnt[j]=((i>>(2*j))&3);
			if(cnt[j]==0){
				que[u]=j; u++;
				used[j]=1;
				nuo[i]^=(1<<j);
			}else if(cnt[j]==3) used[j]=1;
		}
		while(t<u){
			int x=que[t]; t++;
			if(x>0 && !used[x-1]){
				cnt[x-1]--;
				if(cnt[x-1]==0){
					que[u]=x-1; u++;
					used[x-1]=1;
					nuo[i]^=(1<<(x-1));
				}
			}
			if(x<c-1 && !used[x+1]){
				cnt[x+1]--;
				if(cnt[x+1]==0){
					que[u]=x+1; u++;
					used[x+1]=1;
					nuo[i]^=(1<<(x+1));
				}
			}
		}
		nue[i]=popcount(nuo[i]);
	}
	int myo[1<<11]={};
	for(int i=0; i<(1<<c); i++){
		for(int j=0; j<c; j++) if(i&(1<<j)) myo[i]^=(1<<(2*j));
	}
	for(int i=0; i<r; i++){
		for(int k=0; k<(1<<c); k++){
			int x=0, y=k;
			for(int l=0; l<c; l++){
				if(k&(1<<l)){
					if(s[i][l]==0){
						y^=(1<<l);
					}else if(s[i][l]==4){
						y^=(1<<l);
						x^=(3<<(2*l));
					}else{
						x^=(s[i][l]<<(2*l));
					}
				}else{
					x^=(3<<(2*l));
				}
			}
			for(int j=0; j<(1<<c); j++){
				double myon=q[i][k]*dp[i][j];
				dp[i+1][nuo[x-myo[j&y]]]+=myon;
				ans+=myon*nue[x-myo[j&y]];
			}
		}
	}
	printf("%.9lf\n", ans);
	return 0;
}
0