51nod 1242 斐波那契数列的第N项（O(logn)求递推式）

51nod 1242 斐波那契数列的第N项（矩阵快速幂）

斐波那契数列的第N项，我们已知的用递推式求（O（n）），用递归更加慢。
那么我们使用一种（O（logn））的算法，来学习下呗。

它叫做 矩阵快速幂，矩阵加快速幂更配哦。
给出矩阵计算式子：

http://blog.csdn.net/qq_33199236/article/details/52516118
ac代码

#include <iostream>
#include <stdio.h>
#include <string.h>
#include <algorithm>
#include <math.h>
#include <queue>
#include <stack>
#include <vector>
using namespace std;
#define INF 0x3f3f3f
#define pi acos(-1.0)
#define MAX 1000010
#define N 105

struct Mat
{
    long long p[N][N];
};
long long n;
int size = 2,mod = 1e9+9;
Mat mul(Mat a,Mat b)
{
    Mat c;
    memset(c.p,0,sizeof(c.p));
    for(int i = 0; i < size; i++)
        for(int j = 0; j < size; j++)
            for(int k = 0; k < size; k++)
                c.p[i][j] = (c.p[i][j]+a.p[i][k]*b.p[k][j])%mod;
    return c;
}
Mat quick_mi(Mat a)
{
    Mat b;
    for(int i = 0; i < size; i++)
        b.p[i][i] = 1;
    while(n)
    {
        if(n&1) b = mul(b,a);
        a = mul(a,a);
        n>>=1;
    }
    return b;
}
int main()
{
    scanf("%lld",&n);
    Mat A,B,C;
    A.p[0][0] = 0;
    A.p[0][1] = 1;
    A.p[1][0] = 1;
    A.p[1][1] = 1;
    B.p[0][0] = 0;
    B.p[1][0] = 1;
    B.p[1][1] = 0;
    B.p[0][1] = 0;
    C = quick_mi(A);
    C = mul(C,B);
    printf("%lld\n",C.p[0][0]);
}