結果

問題 No.720 行列のできるフィボナッチ数列道場 (2)
ユーザー ixmelixmel
提出日時 2017-06-24 21:41:43
言語 C++11
(gcc 13.3.0)
結果
AC  
実行時間 27 ms / 2,000 ms
コード長 1,688 bytes
コンパイル時間 977 ms
コンパイル使用メモリ 94,704 KB
実行使用メモリ 6,824 KB
最終ジャッジ日時 2024-10-04 03:05:15
合計ジャッジ時間 2,373 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge2
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 20
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include<iostream>
#include<vector>
#include<string>
#include<algorithm>	
#include<map>
#include<set>
#include<utility>
#include<cmath>
#include<cstring>
#include<queue>
#include<stack>
#include<cstdio>
#include<sstream>
#include<iomanip>
#include<assert.h>
#define loop(i,a,b) for(int i=a;i<b;i++) 
#define rep(i,a) loop(i,0,a)
#define pb push_back
#define all(in) in.begin(),in.end()
#define shosu(x) fixed<<setprecision(x)
using namespace std;
//kaewasuretyuui
typedef long long ll;
typedef ll Def;
typedef pair<Def,Def> pii;
typedef vector<Def> vi;
typedef vector<vi> vvi;
typedef vector<pii> vp;
typedef vector<vp> vvp;
typedef vector<string> vs;
typedef vector<double> vd;
typedef vector<vd> vvd;
typedef pair<Def,pii> pip;
typedef vector<pip>vip;
//#define mt make_tuple
//typedef tuple<int,string> tp;
//typedef vector<tp> vt;
const double PI=acos(-1);
const double EPS=1e-7;
const int inf=1e9;
const ll INF=2e18;
int dx[]={0,1,0,-1};
int dy[]={1,0,-1,0};
#define MOD 1000000007
ll m;
vvi mul(vvi A,vvi B){
	vvi C(A.size(),vi(B[0].size()));
	rep(i,2)rep(j,m+1){
		C[i][j]=(A[i][0]*B[0][j]%MOD+A[i][1]*B[1][j]%MOD)%MOD;
	}
	(C[0][m]+=A[0][m])%=MOD;
	(C[1][m]+=A[1][m])%=MOD;
	return C;
}
vvi pow(vvi A,ll n){
	vvi B;
	bool h=false;
	while(n){
		if(n&1){
			if(h)B=mul(B,A);
			else B=A;
			h=true;
		}
		A=mul(A,A);
		n>>=1;
	}
	return B;
}
int main(){
	ll n;
	cin>>n>>m;
	assert(1<=n&&n<=10000000000ll&&2<=m&&m<=30000);
	vi A(m+1,1);
	vvi B(2,vi(m+1,1));
	A[m]=0;B[0][m-2]=2;B[1][m]=0;
	rep(i,m-2){
		A[m-i-3]=(A[m-i-2]+A[m-i-1])%MOD;
		rep(j,2)B[j][m-i-3]=(B[j][m-i-2]+B[j][m-i-1])%MOD;
	}
	vvi out=pow(B,n);
	cout<<(A[0]*out[0][m]%MOD+A[1]*out[1][m]%MOD)%MOD<<endl;
}
0