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

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