矩阵快速幂求斐波那契数列

#include<iostream>
#include<stdio.h>
#include<string.h>
#include<stdlib.h>

using namespace std;

typedef long long ll;
const int MOD=10000;
struct mat
{
    ll a[2][2];
};

mat mat_mul(mat x,mat y);
void mat_pow(int n);

int main()
{
    int n;
    while(cin >> n &&n != -1)
    {
        mat_pow(n);
    }
    return 0;
}

mat mat_mul(mat x,mat y){
    mat res;
    memset(res.a,0,sizeof(res.a));

    for(int i = 0;i < 2;i++){
        for(int j = 0;j < 2;j++){
            for(int k = 0;k < 2;k++){
                res.a[i][j] = (res.a[i][j] + x.a[i][k]*y.a[k][j])%MOD;
            }
        }
    }
    return res;
}

void mat_pow(int n){
    mat c,res;
    c.a[0][0] = c.a[0][1] = c.a[1][0] = 1;
    c.a[1][1] = 0;
    memset(res.a,0,sizeof(res.a));
    for(int i=0;i<2;i++) res.a[i][i]=1;
    while(n){
        if(n&1) res=mat_mul(res,c);
        c=mat_mul(c,c);
        n = n >> 1;
    }
    printf("%lld\n",res.a[0][1]);
}