pek （北大oj）3070

思路：矩阵快速幂， 二分加速

#include<cstdio>
#include<cstring>
#define ll long long
#define mod 10000

struct jz
{
    ll num[][];
    jz(){ memset(num, , sizeof(num)); }
    jz operator*(const jz&p)const
    {
        jz ans;
        for (int k = ; k < ;++k)
        for (int i = ; i < ;++i)
        for (int j = ; j < ; ++j)
            ans.num[i][j] = (ans.num[i][j] + num[i][k] * p.num[k][j] % mod) % mod;
        return ans;
    }
};
jz POW(jz x, ll n)
{
    jz ans;
    for (int i = ; i < ; ++i)ans.num[i][i] = ;
    for (; n; n>>=, x=x*x)
    if (n & ){ ans = ans*x; }
    return ans;
}
int main()
{
    ll n;
    while (scanf("%lld", &n) && (n != -))
    {
        if (n <= )printf("%lld\n", n);
        else
        {
            jz pp;
            pp.num[][] = ; pp.num[][] = ;
            pp.num[][] = ;
            pp = POW(pp, n - );
            printf("%lld\n", (pp.num[][]+pp.num[][])%mod);
        }
    }
}