結果

問題 No.298 話の伝達
ユーザー convexineq
提出日時 2021-01-02 15:29:15
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 219 ms / 5,000 ms
コード長 876 bytes
コンパイル時間 266 ms
コンパイル使用メモリ 82,304 KB
実行使用メモリ 96,256 KB
最終ジャッジ日時 2024-10-12 04:22:05
合計ジャッジ時間 2,495 ms
ジャッジサーバーID
(参考情報)
judge4 / judge3
このコードへのチャレンジ
(要ログイン)
ファイルパターン 結果
other AC * 21
権限があれば一括ダウンロードができます

ソースコード

diff #

n,m = map(int,input().split())
prob = [[0.0]*n for _ in range(n)]
g = [[] for _ in range(n)]
indegree = [0]*n
for _ in range(m):
    a,b,c = map(int,input().split())
    prob[a][b] = c/100
    g[a].append(b)
    indegree[b] += 1
from collections import deque
q = deque(i for i in range(n) if indegree[i]==0)
order = []
while q:
    i = q.popleft()
    order.append(i)
    for b in g[i]:
        indegree[b] -= 1
        if indegree[b]==0: q.append(b)
while order[0] != 0: order.pop(0)
while order[-1] != n-1: order.pop()
if not order: print(0); exit()
n = len(order)
dp = [0.0,1.0]
for i in range(1,n):
    p = order[i]
    N = 1<<i
    ndp = [0.0]*(2*N)
    for mask in range(N):
        x = 1.0
        for j in range(i):
            x *= 1.0-(mask>>j&1)*prob[order[j]][p]
        ndp[mask] += dp[mask]*x
        ndp[mask+N] += dp[mask]*(1-x)
    dp = ndp
print(sum(dp[N:]))
0