結果
問題 |
No.3228 Very Large Fibonacci Sum
|
ユーザー |
![]() |
提出日時 | 2025-08-10 01:24:07 |
言語 | PyPy3 (7.3.15) |
結果 |
AC
|
実行時間 | 51 ms / 2,000 ms |
コード長 | 952 bytes |
コンパイル時間 | 415 ms |
コンパイル使用メモリ | 82,596 KB |
実行使用メモリ | 63,596 KB |
最終ジャッジ日時 | 2025-08-10 01:24:10 |
合計ジャッジ時間 | 3,005 ms |
ジャッジサーバーID (参考情報) |
judge3 / judge4 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 23 |
ソースコード
import sys input = sys.stdin.readline a,b,c,d,e,N=map(int,input().split()) mod=10**9+7 # 行列の計算(numpyを使えないとき,modを使用) def prod(A,B,k,l,m):# A:k*l,B:l*m C=[[None for i in range(m)] for j in range(k)] for i in range(k): for j in range(m): ANS=0 for pl in range(l): ANS=(ANS+A[i][pl]*B[pl][j])%mod C[i][j]=ANS return C def plus(A,B,k,l):# a,B:k*l C=[[None for i in range(l)] for j in range(k)] for i in range(k): for j in range(l): C[i][j]=(A[i][j]+B[i][j])%mod return C A=[[c,1,0,1],[d,0,0,0],[e,0,1,0],[0,0,0,1]] POWA=[A] # 漸化式を行列累乗で求める(ダブリング) for i in range(60): POWA.append(prod(POWA[-1],POWA[-1],4,4,4)) # ベキを求めて X=[[b,a,1,a]] while N: X=prod(X,POWA[N.bit_length()-1],1,4,4) # n乗の場合 N-=1<<(N.bit_length()-1) print(X[0][3]%mod)