POJ 3070

链接：http://vjudge.net/contest/126908#problem/A

这是一道很典型的矩阵快速幂求fib数列的题目，算是矩阵快速幂的入门题目了吧。矩阵快速幂和二分求一个数的n次方算法的原理是一样的。原理网上有好多，大家查一查= =。

#include <iostream>
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <cmath>
using namespace std;
typedef long long ll;
struct matrix
{
    int arr[2][2]; //打包矩阵 
};
matrix func(matrix a , matrix b);
int qk(int n);
int main()
{
    int n;
    while(scanf("%d",&n) && n!=-1)
    {
        printf("%d\n",qk(n));
    }
    return 0;
}
int qk(int n)
{
    int i,j,k;
    matrix ans,num;
    num.arr[0][0] = 1; num.arr[0][1] = 1; num.arr[1][0] = 1; num.arr[1][1] = 0;
    for(i=0; i<2; i++)
        for(j=0; j<2; j++)
            ans.arr[i][j] = (i==j); //定义单位矩阵
    while(n > 0)
    {
        if(n&1) //n为奇数的话
            ans = func(ans,num);
        n /= 2;
        num = func(num,num);
    }
    return ans.arr[0][1]%10000;
}
matrix func(matrix a , matrix b)
{
    int i,j,k;
    matrix bb;
    memset(bb.arr,0,sizeof(bb.arr));
    for(i=0; i<2; i++)
    {
        for(j=0; j<2; j++)
            for(k=0; k<2; k++)
                bb.arr[i][j] += (a.arr[i][k]*b.arr[k][j])%10000;
    }
    return bb;
}