結果

問題 No.1813 Magical Stones
ユーザー ygd.
提出日時 2022-01-15 00:02:41
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 869 ms / 2,000 ms
コード長 4,220 bytes
コンパイル時間 210 ms
コンパイル使用メモリ 82,444 KB
実行使用メモリ 120,704 KB
最終ジャッジ日時 2024-11-20 16:16:07
合計ジャッジ時間 17,316 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge1
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 4
other AC * 40
権限があれば一括ダウンロードができます

ソースコード

diff #
プレゼンテーションモードにする

import sys
#input = sys.stdin.readline #文字列につけてはダメ
input = sys.stdin.buffer.readline #文字列につけてはダメ
#sys.setrecursionlimit(1000000)
#import bisect
#import itertools
#import random
from heapq import heapify, heappop, heappush
#from collections import defaultdict 
#from collections import deque
#import copy
#import math
#from functools import lru_cache
#MOD = pow(10,9) + 7
#MOD = 998244353
class UnionFind(object):
    def __init__(self, n=1):
        self.par = [i for i in range(n)]
        self.rank = [0 for _ in range(n)]
        self.size = [1 for _ in range(n)]
    def find(self, x):
        """
        x が属するグループを探索して親を出す。
        """
        if self.par[x] == x:
            return x
        else:
            self.par[x] = self.find(self.par[x])
            return self.par[x]
    def union(self, x, y):
        """
        x と y のグループを結合
        """
        x = self.find(x)
        y = self.find(y)
        if x != y:
            if self.rank[x] < self.rank[y]:
                x, y = y, x
            if self.rank[x] == self.rank[y]:
                self.rank[x] += 1
            self.par[y] = x
            self.size[x] += self.size[y]
    def is_same(self, x, y):
        """
        x と y が同じグループか否か
        """
        return self.find(x) == self.find(y)
    def get_size(self, x):
        """
        x が属するグループの要素数
        """
        x = self.find(x)
        return self.size[x]
def topological(graph, deg):
    start = []
    for i in range(len(deg)):
        if deg[i] == 0:
            start.append(i)
    topo = []
    while start:
        v = start.pop()
        topo.append(v)
        for u in graph[v]:
            deg[u] -= 1
            if deg[u] == 0:
                start.append(u)
    return topo
#https://atcoder.jp/contests/arc030/submissions/20496099より拝借
def scc_dfs1(s, links, status, postorder):
    stack = [s]
    status[s] = 0
    while stack:
        v = stack[-1]
        limit = len(links[v])
        while status[v] < limit:
            u = links[v][status[v]]
            status[v] += 1
            if status[u] != -1:
                continue
            stack.append(u)
            status[u] = 0
            break
        else:
            stack.pop()
            postorder.append(v)
    return
def scc_dfs2(s, rev_links, status):
    stack = [s]
    status[s] = -1
    scc = [s]
    while stack:
        v = stack.pop()
        for u in rev_links[v]:
            if status[u] != -1:
                stack.append(u)
                status[u] = -1
                scc.append(u)
    scc.reverse() #ループの順番が正の順番になるようにする。#要注意
    return scc
def strongly_connected_components(n, links, rev_links):
    status = [-1] * n
    postorder = []
    for v in range(n):
        if status[v] != -1:
            continue
        scc_dfs1(v, links, status, postorder)
    postorder.reverse()
    sccs = []
    for v in postorder:
        if status[v] == -1:
            continue
        sccs.append(scc_dfs2(v, rev_links, status))
    # label。不要なら消しても良い。
    label = [-1] * n
    for i in range(len(sccs)):
        for x in sccs[i]:
            label[x] = i
    return sccs, label #DAG順になったグループごとのリスト, 頂点が何番目のグループにいるか
def main():
    N,M = map(int,input().split())
    G = [[] for _ in range(N)]
    RG = [[] for _ in range(N)]
    uf = UnionFind(N)
    for i in range(M):
        a,b = map(int,input().split())
        a -= 1; b -= 1
        G[a].append(b)
        RG[b].append(a)
        uf.union(a,b)
    
    L, label = strongly_connected_components(N,G,RG)
    if len(L) == 1:
        print(0);exit()
    #print(label)
    IN = [0]*len(L)
    OUT = [0]*len(L)
    for v in range(N):
        for u in G[v]:
            if label[v] == label[u]: continue
            IN[label[v]] += 1
        for u in RG[v]:
            if label[v] == label[u]: continue
            OUT[label[v]] += 1
    #print(IN,OUT)
    ans = max(IN.count(0), OUT.count(0))
    print(ans)
    
if __name__ == '__main__':
    main()
הההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההה
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX
0