题意：

已知：

当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]);

        }

    }

}