結果

問題 No.309 シャイな人たち (1)
ユーザー gew1fw
提出日時 2025-06-12 20:04:20
言語 PyPy3
(7.3.15)
結果
WA  
実行時間 -
コード長 2,310 bytes
コンパイル時間 173 ms
コンパイル使用メモリ 82,344 KB
実行使用メモリ 77,948 KB
最終ジャッジ日時 2025-06-12 20:09:44
合計ジャッジ時間 2,632 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge1
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 4 WA * 9
権限があれば一括ダウンロードができます

ソースコード

diff # 

import sys

def compute_probability(q_list, t):
    m = len(q_list)
    if t > m:
        return 0.0
    total = 0.0
    for mask in range(0, 1 << m):
        bits = bin(mask).count('1')
        if bits < t:
            continue
        prob = 1.0
        for i in range(m):
            if (mask >> i) & 1:
                prob *= q_list[i]
            else:
                prob *= (1.0 - q_list[i])
        total += prob
    return total

def main():
    input = sys.stdin.read().split()
    ptr = 0
    R = int(input[ptr])
    ptr += 1
    C = int(input[ptr])
    ptr += 1
    
    P = []
    for _ in range(R):
        row = list(map(int, input[ptr:ptr+C]))
        ptr += C
        P.append(row)
        
    S = []
    for _ in range(R):
        row = list(map(int, input[ptr:ptr+C]))
        ptr += C
        S.append(row)
    
    def get_neighbors(i, j):
        neighbors = []
        if i > 0:
            neighbors.append( (i-1, j) )
        if j > 0:
            neighbors.append( (i, j-1) )
        if j < C-1:
            neighbors.append( (i, j+1) )
        return neighbors
    
    q = [ [0.0 for _ in range(C)] for _ in range(R) ]
    
    while True:
        new_q = [ [0.0 for _ in range(C)] for _ in range(R) ]
        max_change = 0.0
        for i in range(R):
            for j in range(C):
                neighbors = get_neighbors(i, j)
                q_kl = [ q[k][l] for (k,l) in neighbors ]
                m = len(neighbors)
                
                p_ij = P[i][j] / 100.0
                
                t1 = S[i][j]
                prob1 = 0.0
                if t1 <= m:
                    prob1 = compute_probability(q_kl, t1)
                
                t2 = 4
                prob2 = 0.0
                if t2 <= m:
                    prob2 = compute_probability(q_kl, t2)
                
                new_q_val = p_ij * prob1 + (1.0 - p_ij) * prob2
                new_q[i][j] = new_q_val
                change = abs(new_q_val - q[i][j])
                if change > max_change:
                    max_change = change
        
        if max_change < 1e-9:
            break
        
        q = new_q
    
    expected = sum( q[i][j] for i in range(R) for j in range(C) )
    print("{0:.10f}".format(expected))

if __name__ == "__main__":
    main()
0