poj_3070 Fibonacci

题目链接：http://poj.org/problem?id=3070

 

//C++代码
#include<iostream>
using namespace std;

const int Mod=10000;

struct Matrix{
	int a[2][2];
}origin,res;

Matrix multiply(Matrix x,Matrix y){
	Matrix temp;
	temp.a[0][0]=temp.a[0][1]=temp.a[1][0]=temp.a[1][1]=0;
	for(int i=0;i<2;i++){
		for(int j=0;j<2;j++){
			for(int k=0;k<2;k++) temp.a[i][j]+=x.a[i][k]*y.a[k][j];
			temp.a[i][j]%=Mod;
		}
	}
	return temp;
}

void calc(int n){
	while(n){
		if(n&1) res=multiply(res,origin);
		n>>=1;
		origin=multiply(origin,origin);
	}
}

int main(){
	int n;
	while(cin>>n,n!=-1){
		if(n==0) cout<<0<<endl;
		else if(n==1) cout<<1<<endl;
		else{
			origin.a[0][0]=origin.a[0][1]=origin.a[1][0]=1;
			origin.a[1][1]=0;
			res.a[0][0]=res.a[1][1]=1;
			res.a[0][1]=res.a[1][0]=0;
			calc(n-1);
			cout<<res.a[0][0]<<endl;
		}
	}
	return 0;
}