HDU 1757 矩阵快速幂加速递推

题意：

已知：

当x<10时：f(x)=x

否则：f(x) = a0 * f(x-1) + a1 * f(x-2) + a2 * f(x-3) + ……+ a9 * f(x-10);

求：f(x)%m的值。

思路：

矩阵快速幂加速递推。 嗯嗯

// by SiriusRen
#include <cstdio>
#include <cstring>
using namespace std;
int cases,k,ans,a[10][10],mod;
struct matrix{int a[10][10];void init(){memset(a,0,sizeof(a));}}first,cpy,td;
matrix mul(matrix &a,matrix &b,int x){
    matrix temp;temp.init();
    for(int i=0;i<10;i++)
        for(int j=0;j<x;j++)
            for(int k=0;k<10;k++)
            temp.a[i][j]=(a.a[i][k]*b.a[k][j]+temp.a[i][j])%mod;
    return temp;
}
int main(){
    while(~scanf("%d%d",&k,&mod)){
        cpy.init();first.init();
        ans=0;
        for(int i=0;i<10;i++)scanf("%d",&first.a[9][9-i]),cpy.a[9][9-i]=first.a[9][9-i];
        for(int i=0;i<9;i++)first.a[i][i+1]=cpy.a[i][i+1]=1;
        for(int i=0;i<10;i++)td.a[i][0]=i;
        if(k<10)printf("%d\n",k%mod);
        else{k-=10;
            while(k){
                if(k&1)first=mul(cpy,first,10);
                cpy=mul(cpy,cpy,10);k>>=1;
            }
            first=mul(first,td,1);
            printf("%d\n",first.a[9][0]);
        }
    }
}