結果
問題 |
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)