【矩阵加速】Fibonacci

POJ 3070
题目：略
题解：推矩阵公式，加速。。。

#include<stdio.h>
#include<iostream>
#include<cmath>
#include<string.h>
using namespace std;
const int N=3,mod=10000;
int n=2;
struct Mt{
	int a[N][N];
	void be(){memset(a,0,sizeof(a));}
	Mt operator * (const Mt x)const
	{
		Mt ans;ans.be();
		for(int i=1;i<=n;i++)
			for(int j=1;j<=n;j++)
				for(int k=1;k<=n;k++)
					ans.a[i][j]=(ans.a[i][j]+(a[i][k]*x.a[k][j])%mod)%mod;
		return ans;
	}
};
int main()
{
	int k;
	while(scanf("%d",&k)&&k!=-1)
	{
		Mt st,ans;st.be(),ans.be();
		st.a[1][1]=st.a[1][2]=st.a[2][1]=1;
		ans.a[1][1]=ans.a[2][2]=1;
		while(k)
		{
			if(k&1)
			{
				ans=ans*st;
			}
			st=st*st;
			k>>=1;
		}
		printf("%d",ans.a[1][2]);
	}
}

//---------------------------
你以为我会这么水过这篇吗真的会
$2\times 2->5\times 5$级矩阵
（之后dp优化可以用）
1.$f[i]=f[i-1]+f[i-2]$
2.$f[i]=f[i-1]+f[i-2]+1$
3.$f[i]=f[i-1]+f[i-2]+n+1$
4. $s[i]=s[i-1]+f[i-1]+f[i-2]+n+1$