51Nod - 1242 斐波那契（快速幂）

斐波那契数列的定义如下：

 

F(0) = 0

F(1) = 1

F(n) = F(n - 1) + F(n - 2) (n >= 2)

 

(1, 1, 2, 3, 5, 8, 13, 21, 34, 55, 89, 144, 233, 377, ...)

给出n，求F(n)，由于结果很大，输出F(n) % 1000000009的结果即可。

 

Input输入1个数n(1 <= n <= 10^18)。Output输出F(n) % 1000000009的结果。Sample Input

11

Sample Output

89

#include <iostream>
#include<string.h>
#include<stdio.h>
#define ll long long
using namespace std;
const ll mod = ;
struct mat
{
    ll m[][];
    mat()
    {
        memset(m, , sizeof(m));
    }
};
mat mul(mat &A, mat &B)
{
    mat C;
    for (int i = ; i < ; i++)
    {
        for (int j = ; j < ; j++)
        {
            for (int k = ; k < ; k++)
            {
                C.m[i][j] = (C.m[i][j] + A.m[i][k] * B.m[k][j]) % mod;
            }
        }
    }
    return C;
}
mat pow(mat A, ll n)
{
    mat B;
    B.m[][] = B.m[][] = ;
    while (n)
    {
        if (n & )
            B = mul(A, B);
        A = mul(A, A);
        n >>= ;
    }
    return B;
}
int main()
{
    ll n;
    while (cin >> n)
    {
        mat A;
        A.m[][] = A.m[][] = A.m[][] = ;
        mat B = pow(A, n);
        printf("%lld\n", B.m[][]);
    }
    return ;
}