結果

問題 No.1340 おーじ君をさがせ
ユーザー ayaoni
提出日時 2021-01-16 02:39:53
言語 PyPy3
(7.3.15)
結果
WA  
実行時間 -
コード長 1,243 bytes
コンパイル時間 312 ms
コンパイル使用メモリ 82,720 KB
実行使用メモリ 80,768 KB
最終ジャッジ日時 2024-11-26 20:58:32
合計ジャッジ時間 15,325 ms
ジャッジサーバーID
(参考情報)
judge4 / judge2
このコードへのチャレンジ
(要ログイン)
ファイルパターン 結果
sample AC * 3
other AC * 39 WA * 20
権限があれば一括ダウンロードができます

ソースコード

diff #

import sys

sys.setrecursionlimit(10**7)
def I(): return int(sys.stdin.readline().rstrip())
def MI(): return map(int,sys.stdin.readline().rstrip().split())
def LI(): return list(map(int,sys.stdin.readline().rstrip().split()))
def LI2(): return list(map(int,sys.stdin.readline().rstrip()))
def S(): return sys.stdin.readline().rstrip()
def LS(): return list(sys.stdin.readline().rstrip().split())
def LS2(): return list(sys.stdin.readline().rstrip())


def matrix_multiplication(A,B):
    l = len(A)
    m = len(B)
    n = len(B[0])
    C = [[0]*n for _ in range(l)]
    for i in range(l):
        for j in range(n):
            for k in range(m):
                C[i][j] |= A[i][k] & B[k][j]
    return C


def matrix_power(A,n):
    l = len(A)
    C = [[0]*l for _ in range(l)]
    for i in range(l):
        C[i][i] = 1
    while n > 0:
        if n % 2 == 1:
            C = matrix_multiplication(C,A)
        A = matrix_multiplication(A,A)
        n >>= 1
    return C


N,M,T = MI()
Graph = [[0]*N for _ in range(N)]
for _ in range(M):
    a,b = MI()
    Graph[a][b] = 1

X = matrix_power(Graph,T)
initial_state = [[1]]+[[0]]*(N-1)

ANS = matrix_multiplication(X,initial_state)
ans = 0
for i in range(N):
    ans += ANS[i][0]

print(ans)
0