矩阵快速幂

 

 

 

void mult(int a[2][2],int b[2][2])
{
	int i,j,k,c[2][2]={0,0,0,0};
	for(i=0;i<2;i++)
		for(j=0;j<2;j++)
			for(k=0;k<2;k++)
				c[i][j]=(c[i][j]+(a[i][k]*b[k][j])%mod)%mod;
	for(i=0;i<2;i++)for(j=0;j<2;j++) a[i][j]=c[i][j];
}
int cal(int n)
{
	if(n<1)return 0;
	if(n<20) return F[n]; n--;
	int sum[2][2]={1,0,1,0};
	int m[2][2]={1,1,1,0};
	while(n)
	{
		if(n&1)
			mult(sum,m);
		mult(m,m);
		n>>=1;
	}
	return (sum[0][0]+sum[0][1])%mod;
}