poj 3070

题目：http://poj.org/problem?id=3070

算是矩阵快速幂的模板吧

#include <stdio.h>
#include <string.h>
#include <iostream>
#include <algorithm>
#include <vector>
#include <math.h>
#define N 160
#define mod 10000
#define _clr(a,val) (memset(a,val,sizeof(a)))

using namespace std;

typedef long long ll;

void cal(ll t[][2],ll b[][2])
{
    int c[2][2];
    _clr(c,0);
    int i,j,k;
    for(k = 0; k < 2; k++)
    {
        for(i = 0; i < 2; i++)
        {
            if(!t[i][k]) continue;
            for(j = 0; j < 2; j++)
            {
                if(!b[k][j]) continue;
                c[i][j] += ((t[i][k] * b[k][j]) % mod);
            }
        }
    }
    for(i = 0; i < 2; i++)
    {
        for(j = 0; j < 2; j++)
        {
            c[i][j] %= mod;
            t[i][j] = c[i][j];
        }
    }
}
void ca(ll a[][2],ll n,ll b[][2])
{
    while(n)
    {
        if(n & 1) cal(b,a);
        cal(a,a);
        n /= 2;
    }
}
int main()
{
    ll n;
    //freopen("data.txt","r",stdin);
    while(cin>>n)
    {
        if(n < 0) break;
        if(n == 0 || n == 1)
        {
            printf("%d\n",n);
            continue;
        }
        ll a[2][2] = {{0,1},{1,1}};
        ll b[2][2] = {{1,0},{0,1}};
        ca(a,n - 1,b);
        int ans = b[1][1] % mod;

        printf("%d\n",ans);
    }
    return 0;
}