結果

問題 No.309 シャイな人たち (1)
コンテスト
ユーザー tjake
提出日時 2015-12-02 17:10:40
言語 Python2
(2.7.18)
結果
WA  
実行時間 -
コード長 1,136 bytes
コンパイル時間 274 ms
コンパイル使用メモリ 6,912 KB
実行使用メモリ 7,552 KB
最終ジャッジ日時 2024-09-14 07:59:19
合計ジャッジ時間 3,567 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge4
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 4 WA * 9
権限があれば一括ダウンロードができます

ソースコード

diff #
raw source code 

inputs = lambda:map(int, raw_input().split())
r, c = inputs()
raw_input()
p = [inputs() + [0] for i in xrange(r)] + [0]*(c+1)
raw_input()
s = [inputs() + [4] for i in xrange(r)] + [4]*(c+1)

dp = [[0]*(2**c) for i in xrange(r+1)]
dpb = [[0]*c for i in xrange(r+1)]

dp[0][0] = 1.0

bit = [[(i>>k)&1 for k in xrange(c+1)] for i in xrange(2**c)]
ans = 0
for i in xrange(r):
    si = s[i]
    pi = p[i]
    dpp = dp[i]
    dpi = dp[i+1]
    dpbp = dpb[i]
    dpbi = dpb[i+1]
    for S in xrange(2**c):
        bit_s = bit[S]
        res = 1.0
        for j in xrange(c):
            x = bit_s[j-1]; y = bit_s[j+1]
            if x+y >= si[j]:
                res *= pi[j]/100.0 if bit_s[j] else (100-pi[j])/100.0
            elif x+y+1 >= si[j]:
                rp = dpbp[j] * (pi[j]/100.0 if bit_s[j] else (100-pi[j])/100.0)
                if not bit_s[j]:
                    rp += (1-dpbp[j])
                res *= rp
            elif bit_s[j]:
                res = 0.0
        dpi[S] += res
        for j in xrange(c):
            if bit_s[j]:
                dpbi[j] += res
    for j in xrange(c):
        ans += dpbi[j]
print ans
0