矩阵快速幂模板（POJ-3070）

模板：

const int MOD = 10000;
struct Node{
	int a[2][2];
};
Node mul(Node x,Node y){ //矩阵乘法 
	Node res;
	for(int i=0;i<2;i++){
		for(int j=0;j<2;j++){
			res.a[i][j] = 0;
			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;
}
int POW(int n){ //矩阵快速幂 
	Node c,res;
	for(int i=0;i<2;i++){
		for(int j=0;j<2;j++){
			c.a[i][j] = 1;
			res.a[i][j] = 0;
		}
		res.a[i][i] = 1;
	}
	c.a[1][1] = 0;
	while(n){
		if(n&1){
			res = mul(res,c);
		}
		c = mul(c,c);
		n>>=1;
	}
	/*for(int i=0;i<2;i++){
		for(int j=0;j<2;j++){
			printf("%d ",res.a[i][j]);
		}
		printf("\n");
	}*/
	return res.a[0][1];
}