1. 程式人生 > >[洛谷P1962]斐波那契數列

[洛谷P1962]斐波那契數列

數列 ont brush quad str times ron light -1

題目大意:求斐波那契數列第n項。第一項和第二項為1。

解題思路:矩陣快速冪模板題。

公式:

$$\begin{bmatrix} F[n] \\F[n-1] \end{bmatrix} \quad =\begin{bmatrix} 1 & 1 \\ 1 & 0 \end{bmatrix} ^n \quad × \begin{bmatrix} F[2] \\ F[1] \end{bmatrix} \quad$$

C++ Code:

#include<cstdio>
#include<cstring>
using namespace std;
struct mat{
	long long a[2][2];
};
mat mul(mat x,mat y){
	mat ans;
	memset(&ans,0,sizeof ans);
	for(int i=0;i<2;++i)
	for(int j=0;j<2;++j)
	for(int k=0;k<2;++k)
	ans.a[i][j]=(ans.a[i][j]+x.a[i][k]*y.a[k][j])%1000000007;
	return ans;
}
void pow(long long n){
	mat p;
	p.a[0][0]=p.a[0][1]=p.a[1][0]=1;
	p.a[1][1]=0;
	mat ans=p;
	while(n){
		if(n&1){
			ans=mul(ans,p);
		}
		p=mul(p,p);
		n>>=1;
	}
	printf("%d\n",(int)ans.a[0][1]);
}
int main(){
	long long n;
	scanf("%lld",&n);
	pow(n-1);
	return 0;
}

[洛谷P1962]斐波那契數列