矩阵快速幂是一种很有用的方法,具体作用就是让矩阵A^n的运算变成log(n)的时间复杂度;
因为矩阵里面A^2=A*A,A^4=A^2*A^2;所以说可以将一个n拆分成这样的运算,大大减少了时间复杂度,和快速幂差不多;
给一个题目用来理解
题目:http://120.78.128.11/Problem.jsp?pid=3089
#include<math.h>
#include<stdio.h>
#define ll long long
int M[40][2][2]={1, 1, 1, 0}, sav[2][2], sav1[2][2];
int main( )
{
int a, b, x, y;
int i, j, k, l, m, arr[100], top;
scanf("%d%d%d", &x, &y, &m);
for(i=1; i<33; i++)
for(j=0; j<2; j++)
for(k=0; k<2; k++)
for(l=0; l<2; l++)
M[i][j][k]=(M[i][j][k]+(ll)M[i-1][j][l]*M[i-1][l][k])%m;///预先处理好每一个矩阵,A的平方,A平方的平方…………
top=0;
y--;
sav[0][0]=sav[1][1]=1;
for(i=30; i>=0&&y; i--)///运算到结束为止
{
if(y>=1<<i)
y-=1<<i;
else
continue;
for(j=0; j<2; j++)
for(k=0; k<2; k++)
for(l=0; l<2; l++)
sav1[j][k]=(sav1[j][k]+(ll)sav[j][l]*M[i][l][k])%m;
for(j=0; j<2; j++)
for(k=0; k<2; k++)
sav[j][k]=sav1[j][k], sav1[j][k]=0;
}
a=((ll)(x*(1+sqrt(5))/2+x-1))%m;
b=((2*a-x+1)%m+m)%m;
printf("%d\n", (b*sav[1][0]+a*sav[1][1])%m);
}