結果
問題 | No.658 テトラナッチ数列 Hard |
ユーザー | realDivineJK |
提出日時 | 2020-05-27 19:49:51 |
言語 | Python3 (3.12.2 + numpy 1.26.4 + scipy 1.12.0) |
結果 |
TLE
|
実行時間 | - |
コード長 | 928 bytes |
コンパイル時間 | 121 ms |
コンパイル使用メモリ | 12,544 KB |
実行使用メモリ | 17,700 KB |
最終ジャッジ日時 | 2024-10-13 03:45:31 |
合計ジャッジ時間 | 4,113 ms |
ジャッジサーバーID (参考情報) |
judge2 / judge4 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 28 ms
17,572 KB |
testcase_01 | AC | 27 ms
10,880 KB |
testcase_02 | AC | 35 ms
10,752 KB |
testcase_03 | AC | 70 ms
10,624 KB |
testcase_04 | TLE | - |
testcase_05 | -- | - |
testcase_06 | -- | - |
testcase_07 | -- | - |
testcase_08 | -- | - |
testcase_09 | -- | - |
testcase_10 | -- | - |
ソースコード
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: t = matprod(y, bas, n, n, n) y = t[:][:] t = matprod(bas, bas, n, n, n) bas = t[:][:] 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)