快速矩阵幂HDU2842

解题思想部分参考：http://blog.csdn.net/q3498233/article/details/5786303

题意：

题目大意：有n个环在一条线上，叫你把他们通通取下来最少要做几步操作，
结果mod 200907.要拿下环来是有条件的，i-th环的前面的环如果通通拿下来了，那么(i+1)-th这个环，你可以放进去，也可以取下来.注意，第一个环，可以随意拿进取下，
其他环就是要按照规则才能进行取下和拿进操作。
用矩阵求斐波那契数列


推出公式：
f(n)=2*f(n-2)+f(n-1)+1
f(1)=1 f(2)=2


设f(n)表示取下n个环的步数，那么f(5)就表示取下5个环要走几步啦。
首先我们想，要想取下最后一个环，我们就必须把前面3个环取下，
第四个环留着，所以就需要f(3)步，取下最后一个环要一步，所以总共是
f（3）+1步之后我们把前面3个环放回去，就要f(3)步，之后我们只要再
取走剩下的4个就好了，那不就是f(4)步吗。
所以  f(5) = 2*f(3)+f(4)+1,通过扩展就可得到上面的f(n)的公式了。




构造初始矩阵A和乘法矩阵E
2 1 1 1 1 0
A= 0 0 0 E= 2 0 0

0 0 0 1 0 1

代码：

/*//
HDU 2842 快速矩阵幂

//*/
#include <cstdio>
#include <cstring>
#include <cmath>
int n;
typedef long long LL;
LL mod=200907;
struct Matrix
{
    LL m[3][3];
	void clear()
	{
		memset(m,0,sizeof(m));
	}
}E, Z;

Matrix Mut(Matrix A, Matrix B)
{
    Matrix ans;
    for (int i = 0; i<3; i++)
        for (int j = 0; j<3; j++)
        {
            ans.m[i][j] = 0;
            for (int k = 0; k<3; k++)
            {
                ans.m[i][j] += ((A.m[i][k])*(B.m[k][j]));
                ans.m[i][j]%=mod;
            }
			while(ans.m[i][j]<0)
				ans.m[i][j]+=mod;
        }
    return ans;
}
Matrix Pow(Matrix A, LL b)
{
    Matrix t = A, ans = E;
    while (b)
    {
        if (b % 2) 
			ans = Mut(ans, t);
        b /= 2;
        t = Mut(t, t);
    }
    return ans;
}

int main()
{
    // freopen("in.txt", "r", stdin);
	Matrix A;
	
	while(scanf("%d",&n)!=EOF&&n!=0)
	{
		if(n==1||n==2)
		{
			if(n==1)
				printf("1\n");
			else 
				printf("2\n");
			continue;
		}
		//初始化A
		A.clear();
		A.m[0][0]=2;A.m[0][1]=1;A.m[0][2]=1;

		//初始化E
		E.clear();
		E.m[0][0]=1;E.m[0][1]=1;
		E.m[1][0]=2;
		E.m[2][0]=1;E.m[2][2]=1;

		Matrix ans=Pow(E,n-3);
		ans=Mut(A,ans);
		printf("%I64d\n",ans.m[0][0]);
	}
    return 0;
}