結果
問題 | No.1112 冥界の音楽 |
ユーザー |
![]() |
提出日時 | 2020-07-10 21:53:48 |
言語 | PyPy3 (7.3.15) |
結果 |
WA
|
実行時間 | - |
コード長 | 1,166 bytes |
コンパイル時間 | 243 ms |
コンパイル使用メモリ | 82,352 KB |
実行使用メモリ | 74,872 KB |
最終ジャッジ日時 | 2024-10-11 08:57:48 |
合計ジャッジ時間 | 3,041 ms |
ジャッジサーバーID (参考情報) |
judge5 / judge3 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 23 WA * 11 |
ソースコード
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]) print(ans) if __name__ == '__main__': main()