結果

問題 No.1813 Magical Stones
ユーザー 👑 tamatotamato
提出日時 2022-01-14 23:07:55
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 569 ms / 2,000 ms
コード長 3,227 bytes
コンパイル時間 576 ms
コンパイル使用メモリ 87,468 KB
実行使用メモリ 113,168 KB
最終ジャッジ日時 2023-08-12 21:57:21
合計ジャッジ時間 14,587 ms
ジャッジサーバーID
（参考情報） 		judge12 / judge7
このコードへのチャレンジ（β）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 94 ms
72,092 KB
testcase_01 AC 89 ms
72,092 KB
testcase_02 AC 90 ms
72,112 KB
testcase_03 AC 152 ms
98,220 KB
testcase_04 AC 89 ms
71,744 KB
testcase_05 AC 88 ms
72,148 KB
testcase_06 AC 180 ms
81,024 KB
testcase_07 AC 93 ms
72,388 KB
testcase_08 AC 112 ms
78,072 KB
testcase_09 AC 120 ms
78,216 KB
testcase_10 AC 103 ms
77,252 KB
testcase_11 AC 256 ms
91,240 KB
testcase_12 AC 90 ms
72,200 KB
testcase_13 AC 90 ms
71,740 KB
testcase_14 AC 150 ms
98,120 KB
testcase_15 AC 569 ms
112,888 KB
testcase_16 AC 553 ms
112,016 KB
testcase_17 AC 455 ms
109,240 KB
testcase_18 AC 549 ms
112,120 KB
testcase_19 AC 551 ms
112,196 KB
testcase_20 AC 561 ms
112,816 KB
testcase_21 AC 551 ms
112,252 KB
testcase_22 AC 557 ms
113,168 KB
testcase_23 AC 546 ms
112,132 KB
testcase_24 AC 119 ms
78,304 KB
testcase_25 AC 217 ms
86,532 KB
testcase_26 AC 186 ms
81,852 KB
testcase_27 AC 442 ms
105,480 KB
testcase_28 AC 450 ms
99,644 KB
testcase_29 AC 457 ms
104,716 KB
testcase_30 AC 448 ms
102,996 KB
testcase_31 AC 521 ms
112,660 KB
testcase_32 AC 91 ms
72,084 KB
testcase_33 AC 91 ms
71,940 KB
testcase_34 AC 145 ms
79,260 KB
testcase_35 AC 147 ms
79,416 KB
testcase_36 AC 167 ms
80,664 KB
testcase_37 AC 178 ms
83,168 KB
testcase_38 AC 230 ms
96,356 KB
testcase_39 AC 435 ms
108,592 KB
testcase_40 AC 175 ms
81,260 KB
testcase_41 AC 119 ms
78,092 KB
testcase_42 AC 140 ms
79,492 KB
testcase_43 AC 174 ms
81,520 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