結果
問題 | No.658 テトラナッチ数列 Hard |
ユーザー |
![]() |
提出日時 | 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 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 8 |
ソースコード
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)