結果
問題 | No.309 シャイな人たち (1) |
ユーザー | chocorusk |
提出日時 | 2019-12-29 04:18:16 |
言語 | C++14 (gcc 12.3.0 + boost 1.83.0) |
結果 |
AC
|
実行時間 | 521 ms / 4,000 ms |
コード長 | 2,266 bytes |
コンパイル時間 | 926 ms |
コンパイル使用メモリ | 112,184 KB |
実行使用メモリ | 36,884 KB |
最終ジャッジ日時 | 2024-10-13 07:20:44 |
合計ジャッジ時間 | 8,036 ms |
ジャッジサーバーID (参考情報) |
judge4 / judge3 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 506 ms
36,736 KB |
testcase_01 | AC | 521 ms
36,788 KB |
testcase_02 | AC | 518 ms
36,864 KB |
testcase_03 | AC | 521 ms
36,736 KB |
testcase_04 | AC | 14 ms
20,608 KB |
testcase_05 | AC | 126 ms
24,576 KB |
testcase_06 | AC | 498 ms
36,736 KB |
testcase_07 | AC | 507 ms
36,884 KB |
testcase_08 | AC | 512 ms
36,864 KB |
testcase_09 | AC | 504 ms
36,864 KB |
testcase_10 | AC | 498 ms
36,864 KB |
testcase_11 | AC | 499 ms
36,864 KB |
testcase_12 | AC | 496 ms
36,744 KB |
testcase_13 | AC | 6 ms
20,224 KB |
testcase_14 | AC | 6 ms
20,352 KB |
testcase_15 | AC | 128 ms
24,576 KB |
ソースコード
#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; }