結果
問題 | No.309 シャイな人たち (1) |
ユーザー | paruki |
提出日時 | 2016-06-28 11:12:28 |
言語 | C++14 (gcc 12.3.0 + boost 1.83.0) |
結果 |
AC
|
実行時間 | 854 ms / 4,000 ms |
コード長 | 2,499 bytes |
コンパイル時間 | 1,686 ms |
コンパイル使用メモリ | 176,680 KB |
実行使用メモリ | 20,480 KB |
最終ジャッジ日時 | 2024-10-11 23:28:36 |
合計ジャッジ時間 | 13,299 ms |
ジャッジサーバーID (参考情報) |
judge3 / judge4 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 850 ms
20,332 KB |
testcase_01 | AC | 845 ms
20,352 KB |
testcase_02 | AC | 854 ms
20,096 KB |
testcase_03 | AC | 835 ms
20,224 KB |
testcase_04 | AC | 18 ms
20,288 KB |
testcase_05 | AC | 201 ms
20,352 KB |
testcase_06 | AC | 849 ms
20,344 KB |
testcase_07 | AC | 849 ms
20,224 KB |
testcase_08 | AC | 836 ms
20,224 KB |
testcase_09 | AC | 839 ms
20,352 KB |
testcase_10 | AC | 848 ms
20,352 KB |
testcase_11 | AC | 843 ms
20,312 KB |
testcase_12 | AC | 845 ms
20,480 KB |
testcase_13 | AC | 8 ms
19,840 KB |
testcase_14 | AC | 7 ms
20,000 KB |
testcase_15 | AC | 212 ms
20,096 KB |
ソースコード
#include "bits/stdc++.h" using namespace std; #define FOR(i,j,k) for(int (i)=(j);(i)<(int)(k);++(i)) #define rep(i,j) FOR(i,0,j) #define each(x,y) for(auto &(x):(y)) #define mp make_pair #define all(x) (x).begin(),(x).end() #define debug(x) cout<<#x<<": "<<(x)<<endl #define smax(x,y) (x)=max((x),(y)) #define smin(x,y) (x)=min((x),(y)) #define MEM(x,y) memset((x),(y),sizeof (x)) #define sz(x) (int)(x).size() typedef long long ll; typedef pair<int, int> pii; typedef vector<int> vi; typedef vector<ll> vll; int bitcnt(int n){ int res = 0; while(n){ if(n&1) ++res; n>>=1; } return res; } const double EPS = 1e-14; int R, C; double P[11][11], prob[12][1 << 11]; int S[11][11], nextHands[1 << 22], dontKnow[1 << 11]; void solve(){ rep(mask, 1 << C) rep(i, C) if((mask >> i & 1) == 0) dontKnow[mask] |= (3 << (2 * i)); rep(mask, 1 << (C * 2)){ queue<int> Q; vi state(C); rep(x, C){ state[x] = mask >> (2 * x) & 3; if(state[x] == 0)Q.push(x + 1), Q.push(x - 1); } while(sz(Q)){ int x = Q.front(); Q.pop(); if(x < 0 || x >= C || state[x] <= 0)continue; if(--state[x] == 0)Q.push(x - 1), Q.push(x + 1); } rep(x, C)if(state[x] <= 0)nextHands[mask] |= 1 << x; } double ans = 0.0; prob[0][0] = 1; rep(y, R){ vector<double> rowProb(1 << C, 1.0); rep(mask, 1 << C)rep(x, C){ if(mask >> x & 1)rowProb[mask] *= P[y][x]; else rowProb[mask] *= 1 - P[y][x]; } rep(upMask, 1 << C){ int needMask = 0; rep(x, C){ int need = 4 - (upMask >> x & 1) - (4 - S[y][x]); smin(need, 3); smax(need, 0); needMask |= need << (x * 2); } rep(knowMask, 1 << C){ int nextUpMask = nextHands[needMask | dontKnow[knowMask]]; prob[y + 1][nextUpMask] += prob[y][upMask] * rowProb[knowMask]; } } rep(mask, 1 << C)ans += bitcnt(mask)*prob[y + 1][mask]; } printf("%0.20f\n", ans); } int main(){ while(cin >> R >> C){ rep(i, R)rep(j, C){ cin >> P[i][j]; P[i][j] *= 0.01; } rep(i, R)rep(j, C)cin >> S[i][j]; rep(i, R + 1)rep(j, 1 << C)prob[i][j] = 0; MEM(nextHands, 0); MEM(dontKnow, 0); solve(); } }