結果
| 問題 |
No.718 行列のできるフィボナッチ数列道場 (1)
|
| コンテスト | |
| ユーザー |
|
| 提出日時 | 2022-06-20 23:45:47 |
| 言語 | Python3 (3.13.1 + numpy 2.2.1 + scipy 1.14.1) |
| 結果 |
AC
|
| 実行時間 | 29 ms / 2,000 ms |
| コード長 | 623 bytes |
| コンパイル時間 | 137 ms |
| コンパイル使用メモリ | 12,800 KB |
| 実行使用メモリ | 10,880 KB |
| 最終ジャッジ日時 | 2024-10-13 12:51:21 |
| 合計ジャッジ時間 | 1,729 ms |
|
ジャッジサーバーID (参考情報) |
judge4 / judge5 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| sample | AC * 3 |
| other | AC * 20 |
ソースコード
N = int(input())-1
A=[]
mod=10**9+7
while N>0:
A.append(N%2)
N//=2
AA=[[1,1],[1,0]]
AAA = []
AAA.append(AA)
for i in range(len(A)-1):
a,b = AA
aa = [(a[0]*a[0]+a[1]*b[0])%mod,(a[0]*a[1]+a[1]*b[1])%mod]
bb = [(b[0]*a[0]+b[1]*b[0])%mod, (b[0]*a[1]+b[1]*b[1])%mod]
AAA.append([aa,bb])
AA = [aa,bb]
B = [[1,0],[0,1]]
for i in range(len(A)):
if A[i]==1:
a,b = AAA[i]
c,d = B
B = [[(a[0]*c[0]+a[1]*d[0])%mod,(a[0]*c[1]+a[1]*d[1])%mod],[(b[0]*c[0]+b[1]*d[0])%mod,(b[0]*c[1]+b[1]*d[1])%mod]]
Fnp,Fn = sum(B[0])%mod,sum(B[1])%mod
#print(Fn,Fnp)
ans = Fn*Fnp%mod
print(ans)