poj3070 单位矩阵（转移矩阵构造）+矩阵快速幂

太妙了。。通过矩阵乘法来加速递推

#include<iostream>
#include<cstring>
#include<cstdio>
using namespace std;
#define mod 10000
int n;
void mul(int f[],int a[][]){//一维数组和矩阵相乘
    int c[]={};
    for(int i=;i<;i++)
        for(int j=;j<;j++)
            c[j]=(c[j]+(long long)f[i]*a[i][j])%mod;
    memcpy(f,c,sizeof c);
}
void mulself(int a[][]){//矩阵自乘
    int c[][]={};
    for(int i=;i<;i++)
        for(int j=;j<;j++)
            for(int k=;k<;k++)
                c[i][j]=(c[i][j]+(long long)a[i][k]*a[k][j])%mod;
    memcpy(a,c,sizeof c);
}

int main(){
    while(cin>>n && n>=){
        int f[]={,},a[][]={{,},{,}};
        while(n){
            if(n%)mul(f,a);
            mulself(a);
            n>>=;
        }
        cout<<f[]<<endl;
    }
}
/*
1 1 2 3 5 8 13 21 34
*/