poj 3070 Fibonacci（矩阵乘法快速幂）

来源：互联网发布：南方cass软件编辑：程序博客网时间：2024/06/01 23:17

题目分析：

#include<iostream>#include<cstdio>#include<algorithm>#include<cstring>using namespace std;struct node{int matrix[10][10];};node a;node operator * (node x,node y){node temp;for(int i=0;i<2;i++)for(int j=0;j<2;j++){int t=0;for(int k=0;k<2;k++)t+=x.matrix[i][k]*y.matrix[k][j]%10000;temp.matrix[i][j]=t%10000;}    return temp;}node pow(node x,int k){node temp;//单位矩阵temp.matrix[0][0]=temp.matrix[1][1]=1;temp.matrix[0][1]=temp.matrix[1][0]=0;if(k==0)    return temp;else if(k==1)return x;else {node t=pow(x,k/2);//开始没注意着，TLE；放到前面有stackflow！！！囧if(k%2==0)   return t*t;    else   return t*t*x;}}int main(){int n;for(int i=0;i<2;i++)for(int j=0;j<2;j++){a.matrix[i][j]=1;if(i==1&&j==1)    a.matrix[i][j]=0;}while(scanf("%d",&n)!=EOF  &&  n!=-1){         node temp=pow(a,n); int ans=temp.matrix[0][1]; printf("%d\n",ans);}system("pause");return 0;}