結果

問題 No.1364 [Renaming] Road to Cherry from Zelkova
ユーザー 👑 KazunKazun
提出日時 2021-01-20 04:49:55
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 881 ms / 2,500 ms
コード長 2,283 bytes
コンパイル時間 282 ms
コンパイル使用メモリ 87,312 KB
実行使用メモリ 193,128 KB
最終ジャッジ日時 2023-08-28 11:00:17
合計ジャッジ時間 22,338 ms
ジャッジサーバーID
（参考情報） 		judge11 / judge12
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 94 ms
71,836 KB
testcase_01 AC 92 ms
71,888 KB
testcase_02 AC 93 ms
71,780 KB
testcase_03 AC 95 ms
71,660 KB
testcase_04 AC 96 ms
72,488 KB
testcase_05 AC 94 ms
71,704 KB
testcase_06 AC 96 ms
71,828 KB
testcase_07 AC 96 ms
71,776 KB
testcase_08 AC 193 ms
87,524 KB
testcase_09 AC 149 ms
79,780 KB
testcase_10 AC 185 ms
82,012 KB
testcase_11 AC 193 ms
85,404 KB
testcase_12 AC 188 ms
81,836 KB
testcase_13 AC 607 ms
159,036 KB
testcase_14 AC 691 ms
141,008 KB
testcase_15 AC 698 ms
155,708 KB
testcase_16 AC 522 ms
122,716 KB
testcase_17 AC 386 ms
141,672 KB
testcase_18 AC 864 ms
189,636 KB
testcase_19 AC 832 ms
189,292 KB
testcase_20 AC 852 ms
190,024 KB
testcase_21 AC 877 ms
189,740 KB
testcase_22 AC 881 ms
189,864 KB
testcase_23 AC 348 ms
115,556 KB
testcase_24 AC 166 ms
86,672 KB
testcase_25 AC 509 ms
124,592 KB
testcase_26 AC 681 ms
143,056 KB
testcase_27 AC 467 ms
120,744 KB
testcase_28 AC 429 ms
116,388 KB
testcase_29 AC 417 ms
119,464 KB
testcase_30 AC 409 ms
115,488 KB
testcase_31 AC 368 ms
114,820 KB
testcase_32 AC 331 ms
112,000 KB
testcase_33 AC 665 ms
136,472 KB
testcase_34 AC 757 ms
145,836 KB
testcase_35 AC 878 ms
168,468 KB
testcase_36 AC 878 ms
162,016 KB
testcase_37 AC 314 ms
108,496 KB
testcase_38 AC 275 ms
90,648 KB
testcase_39 AC 276 ms
89,708 KB
testcase_40 AC 272 ms
90,092 KB
testcase_41 AC 270 ms
90,788 KB
testcase_42 AC 272 ms
90,780 KB
testcase_43 AC 348 ms
179,572 KB
testcase_44 AC 328 ms
162,436 KB
testcase_45 AC 347 ms
193,128 KB
testcase_46 AC 209 ms
125,652 KB
testcase_47 AC 95 ms
71,768 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#SCC参考元:https://judge.yosupo.jp/submission/22487
def scc(E):
    n = len(E)
    iE = [[] for _ in range(n)]
    for i, e in enumerate(E):
        for v in e:
            iE[v].append(i)
    T = []
    done = [0] * n # 0 -> 1 -> 2
    ct = 0
    for i0 in range(n):
        if done[i0]: continue
        Q = [~i0, i0]
        while Q:
            i = Q.pop()
            if i < 0:
                if done[~i] == 2: continue
                done[~i] = 2
                T.append(~i)
                ct += 1
                continue
            if i >= 0:
                if done[i]: continue
                done[i] = 1
            for j in E[i]:
                if done[j]: continue
                Q.append(~j)
                Q.append(j)

    done = [0] * n
    SCC = []
    for i0 in T[::-1]:
        if done[i0]: continue
        L = []
        Q = [~i0, i0]
        while Q:
            i = Q.pop()
            if i < 0:
                if done[~i] == 2: continue
                done[~i] = 2
                L.append(~i)
                continue
            if i >= 0:
                if done[i]: continue
                done[i] = 1
            for j in iE[i]:
                if done[j]: continue
                Q.append(~j)
                Q.append(j)
        SCC.append(L)
    return SCC
#================================================
import sys
from collections import defaultdict

input=sys.stdin.readline

N,M=map(int,input().split())
Mod=10**9+7

X=[[] for _ in range(N+1)]
E=[defaultdict(lambda :[0,0]) for _ in range(N+1)]

for _ in range(M):
    u,v,w,a=map(int,input().split())

    X[u].append(v)

    b,l=E[u][v]
    E[u][v]=[b+a,(l+a*w)%Mod]

S=scc(X)

inf=float("inf")
Flag=[0]*(N+1)
Flag[0]=1

DP=[[0,0] for _ in range(N+1)]
DP[0]=[1,0]

for U in S:
    if len(U)>=2:
        F=0
        for v in U:
            F|=Flag[v]

        if F:
            for v in U:
                DP[v]=[inf,inf]

    for u in U:
        x,y=DP[u]
        for v in E[u]:
            Flag[v]|=Flag[u]

            if x==inf:
                DP[v]=[inf,inf]
            else:
                a,w=E[u][v]
                DP[v][0]+=x*a
                DP[v][1]+=y*a+x*w

                DP[v][0]%=Mod
                DP[v][1]%=Mod

print(DP[N][1] if DP[N][1]<inf else "INF")
0