bzoj1009矩阵快速面+kmp

其实kmp真的很次要，求长度为20的kmp感觉真的有点杀鸡用牛刀

这题思路相当明确：一看题就是数位dp，一看n的大小就是矩阵

矩阵的构造用m*m比较方便，本来想写1*m的矩阵乘m*m的，但是感觉想起来太麻烦就偷懒，没想到1A了

log的速度的确可以，87ms贼快，好久没见这么短的运行时间了

#include <cstdio>
int n,m,mod,k;
char ch;
int a[100],ne[100];
int ans[100][100],t[100][100],zy[100][100];
void mul(bool b)
{
    for(int i=0;i<m;i++)
        for(int j=0;j<m;j++)
            for(k=0,t[i][j]=0;k<m;k++)
                t[i][j]=(t[i][j]+(b?ans[i][k]:zy[i][k])*zy[k][j])%mod;
    for(int i=0;i<m;i++)
        for(int j=0;j<m;j++)
            if(b) ans[i][j]=t[i][j];
                else zy[i][j]=t[i][j];
    if(b) n--;else n/=2;
}
int main()
{
    scanf("%d%d%d",&n,&m,&mod);
    for(ch=getchar();ch<'0' || ch>'9';ch=getchar());
    for(int i=1;i<=m;i++,ch=getchar())
        a[i]=ch-'0';
    for(int i=2;i<=m;i++)
    {
        int k=ne[i-1];
        while(k && (a[k+1]!=a[i])) k=ne[k];
        ne[i]=k+(a[k+1]==a[i]);
    }
    for(int i=0;i<m;i++)
    for(int j=0;j<=9;j++)
    {
        int k=i;
        while(k && a[k+1]!=j) k=ne[k];
        k+=(a[k+1]==j);
        if(k!=m) zy[k][i]++;
    }
    for(int i=0;i<m;i++)
        ans[i][i]=1;
    while(n)
        mul(n&1);
    int sum=0;
    for(int i=0;i<m;i++)
        sum=(sum+ans[i][0])%mod;
    printf("%d",sum);
    return 0;
}