分治 C题

#include <iostream>
using namespace std;
void multi(int a[2][2],int c[2][2])
{
    int b[2][2];
    b[0][0] = (a[0][0] *c[0][0] + a[0][1]*c[1][0])%10000;
    b[0][1] = (a[0][0] *c[0][1] + a[0][1]*c[1][1])%10000;
    b[1][0] = (a[1][0] *c[0][0] + a[1][1]*c[1][0])%10000;
    b[1][1] = (a[1][0] *c[0][1] + a[1][1]*c[1][1])%10000;
    a[0][0] = b[0][0];
    a[0][1] = b[0][1];
    a[1][0] = b[1][0];
    a[1][1] = b[1][1];
}
int main()
{
    int t;
    while(cin>>t, t!=-1)
    {
        int ma[2][2] = {{1,1},
            {1,0}
        };
        int an[2][2] = {{1,1},
            {1,0}
        };
        if(t<=1)
        {
            cout<<t<<endl;
            continue;
        }
        t--;
        while(t)
        {
            if(t & 1)
                multi(an,ma);
            t>>=1;
            multi(ma,ma);
        }
        cout<<an[0][1]<<endl;
    }
    return 0;
}