poj3070（矩阵快速幂，矩阵乘法）

求一个巨大的裴波那契数列，

非常不可思议地竟然可以把斐波那契数列优化到O（log n*2^3）的复杂度

同时，用结构体记录矩阵，通过operator对矩阵乘法*，进行重定义

返回结构体，以下是矩阵结构体模版

同时利用快速幂的做法，来对矩阵的乘法来进行复杂度降维！

#include<cstdio>
#include<algorithm>
#include<cstring>
#include<cmath>
#include<cstdlib>

using namespace std;

struct mat
{
int m[2][2];
int x,y;
}ans,d;
mat operator * (mat a,mat b)
{
mat c;
memset(c.m,0,sizeof(c.m));
c.x=a.x,c.y=b.y;
for (int i=0;i<c.x;i++)
for (int j=0;j<c.y;j++)
for (int k=0;k<c.x;k++)
c.m[i][j]=(c.m[i][j]+a.m[i][k]*b.m[k][j])% 10000;
return c;
}///operator 的用法！！！！！！结构体的矩阵乘法定义
int n;

int matmod(int n)
{
d.m[0][1]=d.m[1][1]=d.m[1][0]=1;
d.m[0][0]=0;
ans.m[0][0]=ans.m[1][1]=1;
ans.m[0][1]=ans.m[1][0]=0;
d.x=d.y=ans.x=ans.y=2;
n--;
while (n>0)//这里一定要注意！！！！必须大于0时才能进行，当然对于这题来说这里是可以改的！，不过本题的关键还是矩阵快速幂和矩阵乘法以及定义的学习为关键
{
if (n&1) ans=ans*d;
n>>=1;
d=d*d;
}
if (n!=-1) return ans.m[1][1];else return 0;
}
int main()
{
while (scanf("%d",&n)&&n!=-1)
printf("%d\n",matmod(n));
return 0;
}