#733999 (PyPy3) No.1813 Magical Stones

• 問題: No.1813 Magical Stones
• ユーザー: titia
• 提出日時: 2022-01-25 03:15:32
• 言語: PyPy3 (7.3.15)
• 結果: AC
• 実行時間: 593 ms / 2,000 ms
• コード長: 2,549 bytes
• コンパイル時間: 590 ms
• コンパイル使用メモリ: 82,260 KB
• 実行使用メモリ: 111,616 KB
• 最終ジャッジ日時: 2024-12-15 08:05:32
• 合計ジャッジ時間: 12,662 ms
• ジャッジサーバーID （参考情報）: judge3 / judge1

ファイルパターン	結果
sample	AC * 4
other	AC * 40

ソースコード

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
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))
 
 
הההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההה
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX