結果
問題 | No.1112 冥界の音楽 |
ユーザー | marroncastle |
提出日時 | 2020-07-14 13:24:58 |
言語 | Python3 (3.12.2 + numpy 1.26.4 + scipy 1.12.0) |
結果 |
AC
|
実行時間 | 1,005 ms / 2,000 ms |
コード長 | 1,149 bytes |
コンパイル時間 | 318 ms |
コンパイル使用メモリ | 12,928 KB |
実行使用メモリ | 14,080 KB |
最終ジャッジ日時 | 2024-04-28 15:35:50 |
合計ジャッジ時間 | 5,133 ms |
ジャッジサーバーID (参考情報) |
judge3 / judge1 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 35 ms
10,752 KB |
testcase_01 | AC | 27 ms
10,752 KB |
testcase_02 | AC | 27 ms
10,880 KB |
testcase_03 | AC | 25 ms
10,752 KB |
testcase_04 | AC | 26 ms
10,880 KB |
testcase_05 | AC | 27 ms
10,880 KB |
testcase_06 | AC | 49 ms
10,752 KB |
testcase_07 | AC | 26 ms
10,880 KB |
testcase_08 | AC | 26 ms
10,752 KB |
testcase_09 | AC | 31 ms
10,880 KB |
testcase_10 | AC | 26 ms
10,880 KB |
testcase_11 | AC | 35 ms
11,008 KB |
testcase_12 | AC | 27 ms
10,752 KB |
testcase_13 | AC | 32 ms
11,008 KB |
testcase_14 | AC | 247 ms
11,520 KB |
testcase_15 | AC | 87 ms
11,008 KB |
testcase_16 | AC | 26 ms
10,752 KB |
testcase_17 | AC | 106 ms
11,136 KB |
testcase_18 | AC | 100 ms
11,264 KB |
testcase_19 | AC | 27 ms
10,752 KB |
testcase_20 | AC | 45 ms
11,008 KB |
testcase_21 | AC | 231 ms
11,648 KB |
testcase_22 | AC | 237 ms
11,776 KB |
testcase_23 | AC | 76 ms
11,008 KB |
testcase_24 | AC | 26 ms
10,880 KB |
testcase_25 | AC | 26 ms
10,880 KB |
testcase_26 | AC | 26 ms
10,880 KB |
testcase_27 | AC | 28 ms
10,880 KB |
testcase_28 | AC | 40 ms
11,136 KB |
testcase_29 | AC | 26 ms
10,880 KB |
testcase_30 | AC | 301 ms
12,160 KB |
testcase_31 | AC | 38 ms
10,880 KB |
testcase_32 | AC | 320 ms
12,160 KB |
testcase_33 | AC | 28 ms
10,880 KB |
testcase_34 | AC | 27 ms
10,752 KB |
testcase_35 | AC | 26 ms
10,880 KB |
testcase_36 | AC | 1,005 ms
14,080 KB |
ソースコード
def power(x,n,mod): digit = len(str(bin(n)))-2 a = len(x) dp = [x]*digit for i in range(1,digit): dp[i] = dot(dp[i-1],dp[i-1],mod) ans = [[0]*a for _ in range(a)] for i in range(a): ans[i][i] = 1 for i in range(digit): if 1<<i & n: ans = dot(dp[i],ans,mod) return ans def dot(A,B,mod): p = len(A) q = len(A[0]) r = len(B) s = len(B[0]) if q!=r: print('掛け算できません') return ans = [[0]*s for _ in range(p)] for i in range(p): for j in range(s): for k in range(q): ans[i][j] += A[i][k]*B[k][j] ans[i][j] %= mod return ans def solve(): mod = 10**9+7 K, M, N = map(int, input().split()) init = [[0] for _ in range(K**2)] for i in range(K): init[i] = [1] mat = [[0]*(K**2) for _ in range(K**2)] for i in range(M): x,y,z = map(int, input().split()) mat[(y-1)*K+(z-1)][(x-1)*K+(y-1)] = 1 res = dot(power(mat,N-2,mod),init,mod) ans = 0 for i in range(K): ans += res[K*i][0] return ans%mod print(solve())