poj3070

矩阵快速幂：

#include<iostream>
#include<algorithm>
#include<cstring>
using namespace std;
typedef long long ll;
const ll mod=10000;
const ll maxx=1e9+7;
ll n;
struct node{
	ll p[2][2];
}ans,res;
node pow_node(node ax,node bx){
	node temp;
	for(int i=0;i<2;i++){
		for(int j=0;j<2;j++){
			temp.p[i][j]=0;
		}
	}
	for(int i=0;i<2;i++){
		for(int j=0;j<2;j++){
			for(int k=0;k<2;k++){
				temp.p[i][j]+=ax.p[i][k]*bx.p[k][j]%mod;
				temp.p[i][j]%=mod;
			}
		}
	}
	return temp;
}
void pow_mod(ll b){
	for(int i=0;i<2;i++){
		for(int j=0;j<2;j++){
			if(i==j){
				ans.p[i][j]=1;
			}else{
				ans.p[i][j]=0;
			}
		}
	}
	while(b){
		if(b&1){
			ans=pow_node(ans,res);
		}
		res=pow_node(res,res);
		b>>=1;
	}
}
int main(){
	while(cin>>n){
		if(n==-1)break;
		res.p[0][0]=res.p[0][1]=res.p[1][0]=1;
		res.p[1][1]=0;
		pow_mod(n);
		ll sum=0;
		for(int i=0;i<2;i++){
			for(int j=0;j<2;j++){
				sum+=ans.p[i][j];
			}
		}
		cout<<(ans.p[0][1])%mod<<endl;
	}
	return 0;
}