結果

問題 No.1813 Magical Stones
ユーザー titiatitia
提出日時 2022-01-25 03:15:32
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 631 ms / 2,000 ms
コード長 2,549 bytes
コンパイル時間 165 ms
コンパイル使用メモリ 82,560 KB
実行使用メモリ 111,360 KB
最終ジャッジ日時 2024-05-09 00:07:17
合計ジャッジ時間 13,123 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge3
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 40 ms
52,736 KB
testcase_01 AC 41 ms
52,864 KB
testcase_02 AC 40 ms
52,480 KB
testcase_03 AC 128 ms
100,516 KB
testcase_04 AC 40 ms
52,736 KB
testcase_05 AC 40 ms
52,992 KB
testcase_06 AC 121 ms
78,088 KB
testcase_07 AC 42 ms
53,248 KB
testcase_08 AC 56 ms
64,384 KB
testcase_09 AC 66 ms
68,608 KB
testcase_10 AC 49 ms
60,032 KB
testcase_11 AC 205 ms
86,904 KB
testcase_12 AC 42 ms
53,248 KB
testcase_13 AC 41 ms
52,736 KB
testcase_14 AC 121 ms
100,632 KB
testcase_15 AC 593 ms
108,544 KB
testcase_16 AC 613 ms
108,400 KB
testcase_17 AC 503 ms
111,360 KB
testcase_18 AC 609 ms
109,200 KB
testcase_19 AC 619 ms
108,676 KB
testcase_20 AC 601 ms
108,936 KB
testcase_21 AC 631 ms
108,724 KB
testcase_22 AC 615 ms
108,992 KB
testcase_23 AC 598 ms
108,424 KB
testcase_24 AC 107 ms
76,928 KB
testcase_25 AC 170 ms
82,632 KB
testcase_26 AC 145 ms
78,632 KB
testcase_27 AC 436 ms
100,356 KB
testcase_28 AC 500 ms
100,864 KB
testcase_29 AC 519 ms
102,336 KB
testcase_30 AC 476 ms
102,172 KB
testcase_31 AC 582 ms
110,504 KB
testcase_32 AC 41 ms
52,352 KB
testcase_33 AC 41 ms
53,120 KB
testcase_34 AC 139 ms
77,800 KB
testcase_35 AC 128 ms
78,208 KB
testcase_36 AC 117 ms
78,044 KB
testcase_37 AC 202 ms
81,368 KB
testcase_38 AC 178 ms
93,616 KB
testcase_39 AC 422 ms
100,612 KB
testcase_40 AC 131 ms
78,592 KB
testcase_41 AC 93 ms
76,752 KB
testcase_42 AC 129 ms
78,208 KB
testcase_43 AC 171 ms
79,164 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

N,M=map(int,input().split())

# UnionFind

Group = [i for i in range(N+1)] # グループ分け
Nodes = [1]*(N+1) # 各グループのノードの数

def find(x):
    while Group[x] != x:
        x=Group[x]
    return x

def Union(x,y):
    if find(x) != find(y):
        if Nodes[find(x)] < Nodes[find(y)]:
            
            Nodes[find(y)] += Nodes[find(x)]
            Nodes[find(x)] = 0
            Group[find(x)] = find(y)
            
        else:
            Nodes[find(x)] += Nodes[find(y)]
            Nodes[find(y)] = 0
            Group[find(y)] = find(x)

# （一般グラフの）トポロジカルソート、SCC

E=[[] for i in range(N)]
E_INV=[[] for i in range(N)]

for i in range(M):
    x,y=map(int,input().split())
    x-=1
    y-=1
    E[x].append(y)
    E_INV[y].append(x)

# DFSして帰り際にTOPに点を放り込んでいる。
# NOWで現在地点、USEINDで、どこの辺まで既に見たか、を調べている。
def Top_sort(E):
    Parent=[-1]*N
    USEIND=[0]*N
    TOP=[]

    for ROOT in range(N):
        if Parent[ROOT]!=-1:
            continue
        Parent[ROOT]=ROOT

        NOW=ROOT

        while NOW!=ROOT or USEIND[ROOT]!=len(E[ROOT]):

            if USEIND[NOW]==len(E[NOW]):
                TOP.append(NOW)
                NOW=Parent[NOW]
            elif E[NOW][USEIND[NOW]]==Parent[NOW]:
                USEIND[NOW]+=1
            else:
                NEXT=E[NOW][USEIND[NOW]]
                USEIND[NOW]+=1
                if Parent[NEXT]==-1:
                    Parent[NEXT]=NOW
                    NOW=NEXT
        TOP.append(ROOT)
        
    return TOP[::-1]


USE=[0]*N
SCC=[]

# SCCを調べるための逆順DFS。
# やっていることはhttps://manabitimes.jp/math/1250　などと同じ。
def dfs2(x):
    Q=[x]
    USE[x]=1
    ANS=[]

    while Q:
        x=Q.pop()
        ANS.append(x)
        for to in E_INV[x]:
            if USE[to]==0:
                USE[to]=1
                Q.append(to)
    return ANS

TOP_SORT=Top_sort(E)

for x in TOP_SORT:
    if USE[x]==0:
        SCC.append(dfs2(x))

for i in range(len(SCC)):
    for j in range(len(SCC[i])-1):
        Union(SCC[i][j],SCC[i][j+1])

if Nodes[find(0)]==N:
    print(0)
    exit()

IN=[0]*N
OUT=[0]*N

for i in range(N):
    for to in E[i]:
        if find(i)==find(to):
            continue
        OUT[find(i)]+=1
        IN[find(to)]+=1

ANS=0
ANS2=0

for i in range(N):
    if find(i)==i:
        if OUT[i]==0:
            ANS+=1
        if IN[i]==0:
            ANS2+=1

print(max(ANS,ANS2))
0