結果
問題 | No.1112 冥界の音楽 |
ユーザー |
👑 |
提出日時 | 2022-07-10 16:29:14 |
言語 | PyPy3 (7.3.15) |
結果 |
AC
|
実行時間 | 93 ms / 2,000 ms |
コード長 | 899 bytes |
コンパイル時間 | 651 ms |
コンパイル使用メモリ | 82,048 KB |
実行使用メモリ | 70,656 KB |
最終ジャッジ日時 | 2025-01-03 12:27:13 |
合計ジャッジ時間 | 4,415 ms |
ジャッジサーバーID (参考情報) |
judge4 / judge3 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 34 |
ソースコード
MOD = 10 ** 9 + 7 def matpow(A, B, w): l = len(A) while w: if w & 1: C = [0] * l for i in range(l): for j in range(l): C[i] += A[i][j] * B[j] C[i] %= MOD B = C C = [[0] * l for _ in range(l)] for i in range(l): for j in range(l): for k in range(l): C[i][j] += A[i][k] * A[k][j] C[i][j] %= MOD A = C w >>= 1 return B k, m, n = map(int, input().split()) k2 = k * k A = [[0] * k2 for _ in range(k2)] f = lambda i, j:i * k + j B = [0] * k2 for i in range(k): B[f(0, i)] = 1 for _ in range(m): p, q, r = map(int, input().split()) p -= 1 q -= 1 r -= 1 A[f(q, r)][f(p, q)] += 1 B = matpow(A, B, n - 2) ans = 0 for i in range(k): ans += B[f(i, 0)] print(ans % MOD)