poj 3070 Fibonacci (快速矩陣乘法)

#include<iostream>
#define P 10000
using namespace std;
void Mul(int a[],int b[]){
	int c[4];
	c[0]=(a[0]*b[0]+a[1]*b[2])%P;
	c[1]=(a[0]*b[1]+a[1]*b[3])%P;
	c[2]=(a[2]*b[0]+a[3]*b[2])%P;
	c[3]=(a[2]*b[1]+a[3]*b[3])%P;
	for(int i=0;i<4;i++)
		a[i]=c[i];
}
void FastMul(int n){
	int a[4]={1,1,1,0};
	int ans[4]={1,0,0,1};
	for(;n;n>>=1,Mul(a,a)){
		if(n&1)
			Mul(ans,a);
	}
	cout<<ans[1]<<endl;
}
int main(){
	int n;
	while(cin>>n && n>=0){
		FastMul(n);
	}
	return 0;
}