結果
問題 | No.1340 おーじ君をさがせ |
ユーザー |
![]() |
提出日時 | 2021-05-28 15:56:51 |
言語 | PyPy3 (7.3.15) |
結果 |
AC
|
実行時間 | 744 ms / 2,000 ms |
コード長 | 851 bytes |
コンパイル時間 | 397 ms |
コンパイル使用メモリ | 81,536 KB |
実行使用メモリ | 76,928 KB |
最終ジャッジ日時 | 2024-11-07 04:42:20 |
合計ジャッジ時間 | 14,713 ms |
ジャッジサーバーID (参考情報) |
judge1 / judge4 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 59 |
ソースコード
import sysinput = sys.stdin.buffer.readlinesys.setrecursionlimit(10 ** 7)def mat_mul(A, B):res = [[0] * len(B[0]) for _ in range(len(A))]for i in range(len(A)):for k in range(len(A[0])):for j in range(len(B[0])):res[i][j] |= (A[i][k] & B[k][j])return resdef mat_pow(A, n):size = len(A)res = [[0] * size for _ in range(size)]for i in range(size):res[i][i] = 1while n:if n & 1:res = mat_mul(res, A)A = mat_mul(A, A)n >>= 1return resN, M, T = map(int, input().split())deg = [0] * NA = [[0]*(N+1) for _ in range(N+1)]for _ in range(M):a, b = map(int, input().split())deg[a] += 1A[a][b] = 1for i in range(N):if not deg[i]:A[i][N] = 1B = mat_pow(A, T)ans = sum(B[0][:N])print(ans)