poj3070 Fibonacci 矩阵快速幂
学了线代之后 终于明白了矩阵的乘法。。
于是 第一道矩阵快速幂。。
实在是太水了。。。
这差不多是个模板了
#include <cstdlib>#include <cstring>#include <cstdio>#include <iostream>using namespace std;int N;struct matrix{ int a[3][3];}origin,res;matrix multiply(matrix x,matrix y){ matrix temp; memset(temp.a,0,sizeof(temp.a)); for(int i=0;i<2;i++) { for(int j=0;j<2;j++) { for(int k=0;k<2;k++) { temp.a[i][j]+=(x.a[i][k]*y.a[k][j]%10000); temp.a[i][j]%=10000; } } } return temp;}void init(){ memset(res.a,0,sizeof(res.a)); res.a[0][0]=res.a[1][1]=1; res.a[1][0]=res.a[0][1]=0; origin.a[0][0]=origin.a[0][1]=origin.a[1][0]=1; origin.a[1][1]=0;}void calc(int n){ while(n) { if(n&1) res=multiply(res,origin); n>>=1; origin=multiply(origin,origin); } printf("%d\n",res.a[0][1]%10000);}int main(){ while(scanf("%d",&N)&&N!=-1) { init(); calc(N); } return 0;}