結果

問題 No.309 シャイな人たち (1)
ユーザー mayoko_mayoko_
提出日時 2015-12-03 10:26:09
言語 C++11
(gcc 11.4.0)
結果
WA  
実行時間 -
コード長 3,145 bytes
コンパイル時間 647 ms
コンパイル使用メモリ 83,496 KB
実行使用メモリ 4,352 KB
最終ジャッジ日時 2023-10-12 09:27:57
合計ジャッジ時間 13,770 ms
ジャッジサーバーID
（参考情報） 		judge11 / judge12
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 WA -
testcase_01 TLE -
testcase_02 -- -
testcase_03 -- -
testcase_04 -- -
testcase_05 -- -
testcase_06 -- -
testcase_07 -- -
testcase_08 -- -
testcase_09 -- -
testcase_10 -- -
testcase_11 -- -
testcase_12 -- -
testcase_13 -- -
testcase_14 -- -
testcase_15 -- -
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<cmath>
//#include<cctype>
#include<climits>
#include<iostream>
#include<string>
#include<vector>
#include<map>
//#include<list>
#include<queue>
#include<deque>
#include<algorithm>
//#include<numeric>
#include<utility>
#include<complex>
//#include<memory>
#include<functional>
#include<cassert>
#include<set>
#include<stack>

const int dx[] = {1, 0, -1, 0};
const int dy[] = {0, 1, 0, -1};
using namespace std;
typedef long long ll;
typedef vector<int> vi;
typedef vector<ll> vll;
typedef pair<int, int> pii;

const int MAX = 11;
int P[MAX][MAX], S[MAX][MAX];
double dp[MAX][MAX];
int cnt[MAX];

int main() {
    int R, C;
    cin >> R >> C;
    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];
    // 最初の行
    for (int s = 0; s < (1<<C); s++) {
        double prob = 1;
        for (int i = 0; i < C; i++) {
            if ((s>>i)&1) prob *= P[0][i]/100.;
            else prob *= (100-P[0][i])/100.;
        }
        if (prob == 0) continue;
        queue<int> que;
        memset(cnt, 0, sizeof(cnt));
        for (int i = 0; i < C; i++) {
            if ((s>>i)&1) cnt[i] = 4-S[0][i];
            if (cnt[i] >= 4) que.push(i);
        }
        while (!que.empty()) {
            int now = que.front(); que.pop();
            if (now-1 >= 0) {
                cnt[now-1]++;
                if (cnt[now-1] == 4) que.push(now-1);
            }
            if (now+1 < C) {
                cnt[now+1]++;
                if (cnt[now+1] == 4) que.push(now+1);
            }
        }
        for (int i = 0; i < C; i++) {
            if (cnt[i] >= 4) dp[0][i] += prob;
        }
    }
    // 以降の行
    for (int i = 1; i < R; i++) {
        for (int s = 0; s < (1<<(2*C)); s++) {
            double prob = 1;
            for (int j = 0; j < C; j++) {
                if ((s>>j)&1) prob *= dp[i-1][j];
                else prob *= 1-dp[i-1][j];
            }
            for (int j = 0; j < C; j++) {
                if ((s>>(j+C))&1) prob *= P[i][j]/100.;
                else prob *= (100-P[i][j])/100.;
            }
            if (prob == 0) continue;
            queue<int> que;
            memset(cnt, 0, sizeof(cnt));
            for (int j = 0; j < C; j++) {
                if ((s>>(j+C))&1) cnt[j] = 4-S[i][j];
                if ((s>>j)&1) cnt[j]++;
                if (cnt[j] >= 4) que.push(i);
            }
            while (!que.empty()) {
                int now = que.front(); que.pop();
                if (now-1 >= 0) {
                    cnt[now-1]++;
                    if (cnt[now-1] == 4) que.push(now-1);
                }
                if (now+1 < 0) {
                    cnt[now+1]++;
                    if (cnt[now+1] == 4) que.push(now+1);
                }
            }
            for (int j= 0; j < C; j++) {
                if (cnt[j] >= 4) dp[i][j] += prob;
            }
        }
    }
    double ans = 0;
    for (int i = 0; i < R; i++) for (int j = 0; j < C; j++) ans += dp[i][j];
    printf("%.10lf\n", ans);
    return 0;
}
0