結果
問題 | No.1112 冥界の音楽 |
ユーザー |
![]() |
提出日時 | 2020-07-10 21:54:48 |
言語 | PyPy3 (7.3.15) |
結果 |
AC
|
実行時間 | 100 ms / 2,000 ms |
コード長 | 1,170 bytes |
コンパイル時間 | 274 ms |
コンパイル使用メモリ | 82,268 KB |
実行使用メモリ | 74,716 KB |
最終ジャッジ日時 | 2024-10-11 09:00:10 |
合計ジャッジ時間 | 3,347 ms |
ジャッジサーバーID (参考情報) |
judge4 / judge2 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 34 |
ソースコード
mod = 1000000007 eps = 10**-9 def main(): import sys input = sys.stdin.readline def matmul(A, B): C = [[0] * len(B[0]) for _ in range(len(A))] for i in range(len(A)): for k in range(len(B)): for j in range(len(B[0])): C[i][j] = (C[i][j] + A[i][k] * B[k][j]) % mod return C def matpow(A, p): n = len(A) B = [[0] * n for _ in range(n)] for i in range(n): B[i][i] = 1 while p > 0: if p & 1: B = matmul(B, A) A = matmul(A, A) p >>= 1 return B K, M, N = map(int, input().split()) T = [] mat = [[0] * (K*K) for _ in range(K*K)] for _ in range(M): p, q, r = map(int, input().split()) p -= 1 q -= 1 r -= 1 pre = p + q * K post = q + r * K mat[post][pre] = 1 vec = [[0] for _ in range(K*K)] for i in range(K): vec[i*K][0] = 1 vec2 = matmul(matpow(mat, N-2), vec) ans = 0 for i in range(K): ans = (ans + vec2[i][0])%mod print(ans) if __name__ == '__main__': main()