結果
問題 |
No.1112 冥界の音楽
|
ユーザー |
![]() |
提出日時 | 2020-07-10 23:04:09 |
言語 | PyPy3 (7.3.15) |
結果 |
AC
|
実行時間 | 224 ms / 2,000 ms |
コード長 | 901 bytes |
コンパイル時間 | 378 ms |
コンパイル使用メモリ | 82,176 KB |
実行使用メモリ | 76,464 KB |
最終ジャッジ日時 | 2024-10-11 17:33:43 |
合計ジャッジ時間 | 4,228 ms |
ジャッジサーバーID (参考情報) |
judge1 / judge3 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 34 |
ソースコード
import sys def input(): return sys.stdin.readline()[:-1] MOD = 10**9+7 def mat_dot(one, two, mod): return [[sum([(one[i][k]*two[k][j])%mod for k in range(max(len(two), len(one[0])))])%mod for j in range(len(two[0]))] for i in range(len(one))] def mat_pow(mat, exp, mod): size = len(mat) res = [[0 for _ in range(size)] for _ in range(size)] for i in range(size): res[i][i] = 1 cnt = 0 while (1<<cnt) <= exp: if (exp>>cnt)&1: res = mat_dot(res, mat, mod) mat = mat_dot(mat, mat, mod) cnt += 1 return res k, m, n = map(int, input().split()) pre = [0 for _ in range(k*k)] for i in range(k): pre[i] = 1 mat = [[0 for _ in range(k*k)] for _ in range(k*k)] for _ in range(m): p, q, r = map(int, input().split()) mat[(q-1)*k + r-1][(p-1)*k + q-1] += 1 res = mat_pow(mat, n-2, MOD) ans = 0 for i in range(k): for j in range(k*k): ans += pre[j] * res[i*k][j] ans %= MOD print(ans)