2018.09.25 poj3070 Fibonacci（矩阵快速幂）

传送门

矩阵快速幂板题，写一道来练练手。

这一次在poj做题总算没忘了改万能库。

代码：

#include<iostream>
#include<cstdio>
#define mod 10000
#define A a[0][0]
#define B a[0][1]
#define C a[1][0]
#define D a[1][1]
using namespace std;
int n;
struct Matrix{int a[2][2];Matrix(){A=0,B=C=D=1;}};
inline Matrix operator*(Matrix a,Matrix b){
	Matrix ret;
	ret.A=(a.A*b.A+a.B*b.C)%mod;
	ret.B=(a.A*b.B+a.B*b.D)%mod;
	ret.C=(a.C*b.A+a.D*b.C)%mod;
	ret.D=(a.C*b.B+a.D*b.D)%mod;
	return ret;
}
inline int ksm(){
	if(n==0)return 0;
	if(n==1)return 1;
	n-=2;
	Matrix x,ret;
	while(n){
		if(n&1)ret=ret*x;
		x=x*x,n>>=1;
	}
	return ret.D;
}
int main(){
	while(scanf("%d",&n)&&~n)printf("%d\n",ksm());
	return 0;
}