poj 3070 快速幂 矩阵相乘

矩阵相乘模板

#include <cstdio>
#include <cstring>
# define ld __int64
struct Matrix{
	ld a[2][2];
	void init()
	{
        a[0][0]=a[1][0]=a[0][1]=1;
		a[1][1]=0;
	}
};
Matrix matrix_mul(Matrix a,Matrix b)
{
	int i,j,k;
	Matrix ans;
	for(i=0;i<2;i++)
	{
		for(j=0;j<2;j++)
		{
			ans.a[i][j]=0;
			for(k=0;k<2;k++)
				ans.a[i][j]+=a.a[i][k]*b.a[k][j];
			ans.a[i][j]%=10000;
		}
	}
	return ans;
}
Matrix mult(Matrix a,ld b)
{
	Matrix ans;
	ans.init();
	while(b)
	{
		if(b&1)
		    ans=matrix_mul(ans,a);
		b>>=1;
		a=matrix_mul(a,a);
	}
	return ans;
}
int main()
{
	int n ;
    while(scanf("%I64d",&n),n!=-1)
	{
		Matrix ans,tmp;
        tmp.init();
		ans=mult(tmp,n);
		printf("%I64d\n",ans.a[1][1]%10000);
	}
	return 0;
}