結果
| 問題 |
No.657 テトラナッチ数列 Easy
|
| コンテスト | |
| ユーザー |
|
| 提出日時 | 2022-01-08 20:44:06 |
| 言語 | PyPy3 (7.3.15) |
| 結果 |
AC
|
| 実行時間 | 1,935 ms / 2,000 ms |
| コード長 | 1,214 bytes |
| コンパイル時間 | 497 ms |
| コンパイル使用メモリ | 81,944 KB |
| 実行使用メモリ | 81,436 KB |
| 最終ジャッジ日時 | 2024-11-14 09:46:02 |
| 合計ジャッジ時間 | 12,475 ms |
|
ジャッジサーバーID (参考情報) |
judge3 / judge2 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| sample | AC * 3 |
| other | AC * 13 |
ソースコード
def matrix_mul(A,B,mod = None):
nA = len(A)
mA = len(A[0])
mB = len(B[0])
tmp = [[0]*mB for _ in range(nA)]
if mod is None:
for i in range(nA):
for j in range(mB):
tmp[i][j] = sum(A[i][k]*B[k][j] for k in range(mA))
return tmp
for i in range(nA):
for j in range(mB):
tmp[i][j] = sum(A[i][k]*B[k][j]%mod for k in range(mA))%mod
return tmp
def matrix_pow(A,n,mod = None):
nbit = list(str(bin(n))[2:])
nbit = [int(i) for i in nbit]
N = len(A)
C = [[0]*N for _ in range(N)]
B = A
for i in range(N):
C[i][i] = 1
if mod is None:
for i in range(len(nbit)):
if nbit[-1-i] == 1:
C = matrix_mul(C,B)
B = matrix_mul(B,B)
return C
for i in range(len(nbit)):
if nbit[-1-i] == 1:
C = matrix_mul(C,B,mod)
B = matrix_mul(B,B,mod)
return C
A = [[0,1,0,0],[0,0,1,0],[0,0,0,1],[1,1,1,1]]
INIT = [[0],[0],[0],[1]]
Q = int(input())
for _ in range(Q):
n = int(input())
print(*matrix_mul(matrix_pow(A,n-1,mod=17),INIT,mod=17)[0])