結果
問題 | No.1595 The Final Digit |
ユーザー |
|
提出日時 | 2021-07-09 21:44:50 |
言語 | PyPy3 (7.3.8) |
結果 |
AC
|
実行時間 | 84 ms / 2,000 ms |
コード長 | 1,160 bytes |
コンパイル時間 | 256 ms |
使用メモリ | 80,940 KB |
最終ジャッジ日時 | 2023-02-01 22:48:10 |
合計ジャッジ時間 | 3,056 ms |
ジャッジサーバーID (参考情報) |
judge12 / judge13 |
テストケース
テストケース表示入力 | 結果 | 実行時間 使用メモリ |
---|---|---|
testcase_00 | AC | 77 ms
75,848 KB |
testcase_01 | AC | 79 ms
75,600 KB |
testcase_02 | AC | 76 ms
75,620 KB |
testcase_03 | AC | 78 ms
75,792 KB |
testcase_04 | AC | 82 ms
80,508 KB |
testcase_05 | AC | 82 ms
80,476 KB |
testcase_06 | AC | 84 ms
80,688 KB |
testcase_07 | AC | 82 ms
80,940 KB |
testcase_08 | AC | 74 ms
75,576 KB |
testcase_09 | AC | 73 ms
75,984 KB |
testcase_10 | AC | 73 ms
75,768 KB |
testcase_11 | AC | 74 ms
75,524 KB |
testcase_12 | AC | 73 ms
75,648 KB |
testcase_13 | AC | 76 ms
75,764 KB |
testcase_14 | AC | 83 ms
80,648 KB |
testcase_15 | AC | 82 ms
80,676 KB |
testcase_16 | AC | 83 ms
80,428 KB |
testcase_17 | AC | 75 ms
75,740 KB |
testcase_18 | AC | 79 ms
75,740 KB |
testcase_19 | AC | 84 ms
80,896 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 mv(A,x,mod): p = len(A) q = len(A[0]) r = len(x) if q!=r: print('掛け算できません') return ans = [0]*p for i in range(p): for k in range(q): ans[i] += A[i][k]*x[k] ans[i] %= mod return ans def solve(): p,q,r,K = map(int, input().split()) init = [r,q,p] mat = [[1,1,1],[1,0,0],[0,1,0]] ans = mv(power(mat,K-3,10),init,10) return ans[0] print(solve())