快速矩阵幂 HDU1757

矩阵的快速幂是用来高效地计算矩阵的高次方的。将朴素的o（n）的时间复杂度，降到log（n）。

这篇文章很好得分析了原理：http://www.cnblogs.com/yan-boy/archive/2012/11/29/2795294.html

题意：
if(x<10)
    f(x) = x.
if (x>=10)
f(x) = a0 * f(x-1) + a1 * f(x-2) + a2 * f(x-3) + …… + a9 * f(x-10);
And ai(0<=i<=9) can only be 0 or 1 .


给出 k m  a0-a9
求f(k)%m的值


分析：
由题目可知，我们需要构造一个原始矩阵A和一个
乘法矩阵B，A*B相当于向后计算一位。就是通过
已经计算出f(x-1),..f(x-10)计算出f(x).
现在有a0-a9十个数，也就是10*10的矩阵，直接想很难
想出矩阵如何构造。
那假设只有a0,a1,a2。且a0=3,a1=5,a2=4.
f(0)=2,f(1)=9,f(2)=7.
 n>=3时，f(n)=a0*f(n-1)+a1*f(n-2)+a2*f(n-3).
 则f3=3*f(2)+5*f(1)+4*f(0)=3*7+5*9+4*2=74;
   f4=3*f(3)+5*f(2)+4*f(1)=3*74+5*7+4*9=293


且容易可以得到原始矩阵和乘法矩阵（3*3）分别是
7 9 2       3 1 0 
0 0 0 5 0 1
0 0 0 4 0 0
A的第一行和B的第一列相当于计算新的f(n)，A的
第一行和B的第二列相当于保留7( f(n-1) ),A的第
一行和B的第三列相当于保留9( f(n-2) )

同样的道理：构造出原题的初始矩阵和乘法矩阵

9 8 7 6 5 4 3 2 1 0 a0 1 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 a1 0 1 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 a2 0 0 1 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 a3 0 0 0 1 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 a4 0 0 0 0 1 0 0 0 0
0 0 0 0 0 0 0 0 0 0 a5 0 0 0 0 0 1 0 0 0
0 0 0 0 0 0 0 0 0 0 a6 0 0 0 0 0 0 1 0 0
0 0 0 0 0 0 0 0 0 0 a7 0 0 0 0 0 0 0 1 0
0 0 0 0 0 0 0 0 0 0 a8 0 0 0 0 0 0 0 0 1
0 0 0 0 0 0 0 0 0 0 a9 0 0 0 0 0 0 0 0 0

/*//
HDU 1757 快速矩阵幂
//*/
#include <cstdio>
#include <cstring>

int n,mod;

struct Matrix
{
    int m[12][12];
	void clear()
	{
		memset(m,0,sizeof(m));
	}
}E, Z;

Matrix Mut(Matrix A, Matrix B)
{
    Matrix ans;
    for (int i = 0; i<10; i++)
        for (int j = 0; j<10; j++)
        {
            ans.m[i][j] = 0;
            for (int k = 0; k<10; k++)
            {
                ans.m[i][j] += ((A.m[i][k])*(B.m[k][j]));
                ans.m[i][j] %= mod;
            }
        }
    return ans;
}
Matrix Pow(Matrix A, int 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);
	//A是初始矩阵 E是乘法矩阵
	Matrix A;
	while(scanf("%d%d", &n, &mod)!=EOF)
	{ 
		A.clear();E.clear();
	   for(int i=0;i<10;i++)  //初始化A
	   {
		 A.m[0][i]=9-i;
	   }

	   for(int i=0;i<10;i++)  //初始化E
	   {
		   scanf("%d",&E.m[i][0]);
		   E.m[i][0]%=mod;
	   }
	   if(n<=9)
	   {
		   printf("%d\n",n%mod);
		   continue;
	   }
	   for(int i=1;i<=9;i++)  //填充1
	   {
		   E.m[i-1][i]=1;
	   }

		Matrix ans = Pow(E, n-10);
		ans=Mut(A,ans); //A*(E^(n-10))
		printf("%d\n",ans.m[0][0]);
        
    }
    return 0;
}