矩阵连乘求斐波那契数列模板

#include<stdio.h>
#include<string.h>
#include<algorithm>
#include<math.h>
#include<queue>
#define ll long long
#define M 1000000007
#define N 2
using namespace std;
struct mat
{
	ll m[N][N];
};
mat A=
{
	1,1,
	1,0
};
mat I=
{
	1,0,
	0,1
};
mat multi(mat a,mat b)
{
	mat c;
	for(int i=0;i<N;i++)
	{
		for(int j=0;j<N;j++)
		{
			c.m[i][j]=0;
			for(int k=0;k<N;k++)
				c.m[i][j]+=a.m[i][k]*b.m[k][j]%M;
			c.m[i][j]%=M;
		}
	}
	return c;
}
mat power(mat A,int k)//矩阵快速幂 
{
	mat ans=I,p=A;
	while(k)
	{
		if(k&1)
		{
			ans=multi(ans,p);
			k--;
		}
		k>>=1;
		p=multi(p,p);
	}
	return ans;
}
int main()
{
	int n;
	while(scanf("%d",&n)!=EOF)
	{
		mat ans=power(A,n-1);
		printf("%lld\n",ans.m[0][0]);
	}
	return 0;
}