結果
問題 | No.658 テトラナッチ数列 Hard |
ユーザー | realDivineJK |
提出日時 | 2020-05-27 20:39:09 |
言語 | PyPy3 (7.3.15) |
結果 |
AC
|
実行時間 | 1,578 ms / 2,000 ms |
コード長 | 896 bytes |
コンパイル時間 | 178 ms |
コンパイル使用メモリ | 82,440 KB |
実行使用メモリ | 79,216 KB |
最終ジャッジ日時 | 2024-10-13 03:47:34 |
合計ジャッジ時間 | 8,452 ms |
ジャッジサーバーID (参考情報) |
judge4 / judge3 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 39 ms
52,224 KB |
testcase_01 | AC | 36 ms
52,352 KB |
testcase_02 | AC | 54 ms
68,224 KB |
testcase_03 | AC | 96 ms
76,392 KB |
testcase_04 | AC | 650 ms
78,568 KB |
testcase_05 | AC | 723 ms
77,688 KB |
testcase_06 | AC | 884 ms
78,792 KB |
testcase_07 | AC | 909 ms
78,256 KB |
testcase_08 | AC | 1,138 ms
78,848 KB |
testcase_09 | AC | 1,578 ms
78,544 KB |
testcase_10 | AC | 1,560 ms
79,216 KB |
ソースコード
Q = int(input()) mod = 17 def matprod(m1, m2, n, l, m): y = [[0 for _ in range(m)] for __ in range(n)] for i in range(n): for j in range(m): for k in range(l): y[i][j] += m1[i][k] * m2[k][j] % mod y[i][j] %= mod return y def matpow(m1, n, m): y = [[i==j for j in range(n)] for i in range(n)] bas = [[m1[i][j] for j in range(n)] for i in range(n)] tmp = m while tmp != 0: if tmp % 2 == 1: y = matprod(y, bas, n, n, n)[:][:] bas = matprod(bas, bas, n, n, n)[:][:] tmp //= 2 return y vec = [0, 0, 0, 1] for _ in range(Q): n = int(input()) M = [ [0, 1, 0, 0], [0, 0, 1, 0], [0, 0, 0, 1], [1, 1, 1, 1] ] m = matpow(M, 4, n-1) S = 0 for i in range(4): S += m[0][i] * vec[i] % mod S %= mod print(S)