結果
問題 | No.658 テトラナッチ数列 Hard |
ユーザー | realDivineJK |
提出日時 | 2020-05-27 20:39:09 |
言語 | PyPy3 (7.3.15) |
結果 |
AC
|
実行時間 | 1,580 ms / 2,000 ms |
コード長 | 896 bytes |
コンパイル時間 | 367 ms |
コンパイル使用メモリ | 81,520 KB |
実行使用メモリ | 79,044 KB |
最終ジャッジ日時 | 2024-04-21 05:17:44 |
合計ジャッジ時間 | 8,888 ms |
ジャッジサーバーID (参考情報) |
judge3 / judge5 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 37 ms
52,552 KB |
testcase_01 | AC | 36 ms
53,344 KB |
testcase_02 | AC | 67 ms
68,928 KB |
testcase_03 | AC | 94 ms
76,288 KB |
testcase_04 | AC | 650 ms
78,896 KB |
testcase_05 | AC | 720 ms
77,764 KB |
testcase_06 | AC | 878 ms
78,656 KB |
testcase_07 | AC | 930 ms
78,356 KB |
testcase_08 | AC | 1,124 ms
79,044 KB |
testcase_09 | AC | 1,571 ms
78,528 KB |
testcase_10 | AC | 1,580 ms
79,040 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)