結果

問題 No.1813 Magical Stones
ユーザー tamatotamato
提出日時 2022-01-14 23:07:55
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 425 ms / 2,000 ms
コード長 3,227 bytes
コンパイル時間 289 ms
コンパイル使用メモリ 82,640 KB
実行使用メモリ 112,256 KB
最終ジャッジ日時 2024-04-30 16:59:35
合計ジャッジ時間 9,896 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge5
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 39 ms
54,448 KB
testcase_01 AC 37 ms
55,304 KB
testcase_02 AC 36 ms
56,120 KB
testcase_03 AC 93 ms
97,392 KB
testcase_04 AC 37 ms
54,964 KB
testcase_05 AC 37 ms
54,628 KB
testcase_06 AC 118 ms
79,348 KB
testcase_07 AC 39 ms
56,164 KB
testcase_08 AC 58 ms
70,148 KB
testcase_09 AC 64 ms
74,928 KB
testcase_10 AC 45 ms
64,996 KB
testcase_11 AC 180 ms
89,620 KB
testcase_12 AC 39 ms
55,912 KB
testcase_13 AC 38 ms
55,324 KB
testcase_14 AC 90 ms
97,440 KB
testcase_15 AC 406 ms
111,172 KB
testcase_16 AC 425 ms
112,256 KB
testcase_17 AC 327 ms
105,256 KB
testcase_18 AC 414 ms
111,624 KB
testcase_19 AC 419 ms
111,560 KB
testcase_20 AC 417 ms
111,552 KB
testcase_21 AC 408 ms
111,100 KB
testcase_22 AC 421 ms
112,036 KB
testcase_23 AC 414 ms
111,872 KB
testcase_24 AC 65 ms
73,612 KB
testcase_25 AC 156 ms
84,924 KB
testcase_26 AC 123 ms
80,048 KB
testcase_27 AC 326 ms
104,620 KB
testcase_28 AC 334 ms
97,548 KB
testcase_29 AC 341 ms
103,156 KB
testcase_30 AC 329 ms
100,676 KB
testcase_31 AC 386 ms
109,964 KB
testcase_32 AC 38 ms
56,488 KB
testcase_33 AC 37 ms
55,724 KB
testcase_34 AC 85 ms
77,788 KB
testcase_35 AC 88 ms
77,756 KB
testcase_36 AC 101 ms
78,612 KB
testcase_37 AC 106 ms
80,492 KB
testcase_38 AC 158 ms
94,112 KB
testcase_39 AC 319 ms
106,112 KB
testcase_40 AC 117 ms
79,388 KB
testcase_41 AC 64 ms
74,804 KB
testcase_42 AC 85 ms
77,840 KB
testcase_43 AC 115 ms
79,732 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

mod = 1000000007
eps = 10**-9


def main():
    import sys
    input = sys.stdin.readline
    from collections import deque

    # compressed[v]: 縮約後のグラフで縮約前のvが属する頂点
    # num: 縮約後のグラフの頂点数
    # 縮約後のグラフの頂点番号はトポロジカル順
    def SCC(adj, adj_rev):
        N = len(adj) - 1
        seen = [0] * (N + 1)
        compressed = [0] * (N + 1)
        order = []

        for v0 in range(1, N + 1):
            if seen[v0]:
                continue
            st = deque()
            st.append(v0)
            while st:
                v = st.pop()
                if v < 0:
                    order.append(-v)
                else:
                    if seen[v]:
                        continue
                    seen[v] = 1
                    st.append(-v)
                    for u in adj[v]:
                        st.append(u)

        seen = [0] * (N + 1)
        num = 0
        for v0 in reversed(order):
            if seen[v0]:
                continue
            num += 1
            st = deque()
            st.append(v0)
            seen[v0] = 1
            compressed[v0] = num
            while st:
                v = st.pop()
                for u in adj_rev[v]:
                    if seen[u]:
                        continue
                    seen[u] = 1
                    compressed[u] = num
                    st.append(u)

        return num, compressed

    # 縮約後のグラフを構築
    # 先にSCC()を実行してnum, compressedを作っておく
    def construct(adj, num, compressed):
        N = len(adj) - 1
        adj_compressed = [set() for _ in range(num + 1)]
        for v in range(1, N + 1):
            v_cmp = compressed[v]
            for u in adj[v]:
                u_cmp = compressed[u]
                if v_cmp != u_cmp:
                    adj_compressed[v_cmp].add(u_cmp)
        return adj_compressed

    N, M = map(int, input().split())
    adj = [[] for _ in range(N+1)]
    adj_rev = [[] for _ in range(N + 1)]
    for _ in range(M):
        a, b = map(int, input().split())
        adj[a].append(b)
        adj_rev[b].append(a)

    NN, compressed = SCC(adj, adj_rev)
    adj_new = construct(adj, NN, compressed)

    in_num = [0] * (NN + 1)
    out_num = [0] * (NN + 1)
    for v in range(1, NN + 1):
        for u in adj_new[v]:
            in_num[u] += 1
            out_num[v] += 1

    if NN == 1:
        print(0)
        exit()

    seen = [0] * (NN+1)
    v0_list = []
    plus_sum = minus_sum = 0
    for v0 in range(1, NN+1):
        if seen[v0]:
            continue
        v0_list.append(v0)
        seen[v0] = 1
        que = deque()
        que.append(v0)
        seq = []
        while que:
            v = que.popleft()
            seq.append(v)
            for u in adj_new[v]:
                if seen[u] == 0:
                    seen[u] = 1
                    que.append(u)
        p = s = 0
        for v in seq:
            if in_num[v] == 0:
                p += 1
            if out_num[v] == 0:
                s += 1
        plus_sum += p
        minus_sum += s
    print(max(plus_sum, minus_sum))


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