矩阵快速幂板子

const int maxn=2   //2行2列矩阵 
const int mod=1e9+7;
using namespace std;

struct Matrix{
	int m[maxn][maxn];
}ans,x;   //x为要求幂的矩阵 

Matrix Mul(Matrix A,Matrix B)    //矩阵乘法 
{
	Matrix tmp;  
	for(int i=0;i<maxn;i++)
		for(int j=0;j<maxn;j++)
			tmp.m[i][j]=0;
	for(int i=0;i<maxn;i++)
	{
		for(int j=0;j<maxn;j++)
		{
			for(int k=0;k<maxn;k++)
			{
				tmp.m[i][j]+=(A.m[i][k]*B.m[k][j]%mod);
			}
			tmp.m[i][j]=(tmp.m[i][j]+mod)%mod;
		}
	}
	return tmp;
}

int QuickPower(int n) //斐波那契x={1,1,1,0},x的n次方后ans.m[0][1]即为fn 
{
	for(int i=0;i<maxn;i++)
		for(int j=0;j<maxn;j++)
		{
			if(i==j)
			ans.m[i][j]=1;
			else ans.m[i][j]=0;
		}
		while(n)
		{
			if(n&1)
				ans=Mul(ans,x);
			x=Mul(x,x);
			n=n>>1;
		}
		return (ans.m[0][1]%mod+mod)%mod;
}