GT考试

比较神仙的$dp+KMP+Matrix$综合题目，比较值得一写

$0x00$:首先我打了一个爆搜

不过对正解并无任何启发。。。（逗比发言请忽略）

$0x01$:基础$dp$

状态还是比较好设的，

考虑设$f_{i,j}$表示$dp$到长串的第$i$位，匹配了短串的前缀长度为$j$时的总方案数

直接转移不大好搞，那么我们还需要一个转移数组

设$g_{k,j}$表示在 和短串匹配长度为$k$的串 后面添加一个字符，使新构成的串 和短串匹配长度为$j$，这样的字符的添加方案数

这个东西可以用$KMP$求，具体而言就是枚举每一种长度以及每一种可能的添加字符，

如果这个字符与短串不匹配，那么就跳$nxt$数组找到合法的最长匹配前缀。

那么转移为$f_{i,j}=\sum_{k=0}^{m-1} f_{i-1,k}*g_{k,j}$那么最终的答案就是$\sum_{i=0}^{m-1} f_{n,i}$，然而复杂度$O(nm^2)$直接爆炸

$0x02$:矩阵优化

我们看到这个方程式非常的像矩阵乘法，考虑矩阵加速递推，经过观察，我们发现这个$f$数组貌似可以省去，直接递推$g^n$就可以得到最后的$f_n$

这样，最后的答案就是$\sum_{i=0}^{m-1}g^n_{0,i}$，时间复杂度$O(m^3log_n)$完全可过

#include<bits/stdc++.h>
using namespace std;
namespace AE86{
    inline int read(){
        int x=0,f=1;char ch=getchar();
        while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
        while(ch>='0'&&ch<='9'){x=(x<<1)+(x<<3)+(ch^48);ch=getchar();}return x*f;
    }inline void write(int x,char opt='\n'){
        char ch[20];int len=0;if(x<0)x=~x+1,putchar('-');
        do{ch[len++]=x%10+(1<<5)+(1<<4);x/=10;}while(x);
        for(int i=len-1;i>=0;--i)putchar(ch[i]);putchar(opt);}
}using namespace AE86;
const int NN=1e5+5;
int n,m,mod,nxt[NN],ans;
char s[25];
namespace Matrix{
    struct Ma{
        int p[25][25];
        inline void cle(){memset(p,0,sizeof(p));}
        inline void pre(){for(int i=0;i<m;i++)p[i][i]=1;}
        inline void print(){for(int i=0;i<m;i++){for(int j=0;j<m;j++){cout<<p[i][j]<<" ";} cout<<endl;}}
    };
    inline Ma mul(Ma a,Ma b){
        Ma c;c.cle();
        for(int i=0;i<m;i++)
            for(int j=0;j<m;j++)
                for(int k=0;k<m;k++)
                    (c.p[i][j]+=1ll*a.p[i][k]*b.p[k][j]%mod)%=mod;
        return c;
    }
    inline Ma qmo(Ma a,int b){
        Ma c; c.cle(); c.pre();
        while(b){
            if(b&1) c=mul(c,a);
            b>>=1; a=mul(a,a);
        } return c;
    }
    inline Ma kmp(){
        nxt[1]=0; Ma a; a.cle();
        for(int i=2,j=0;i<=m;i++){
            while(j && s[j+1]!=s[i]) j=nxt[j];
            if(s[j+1]==s[i]) ++j;
            nxt[i]=j;
        }
        for(int i=0;i<m;i++){
            for(int ch='0';ch<='9';ch++){
                int j=i;
                while(j && s[j+1]!=ch) j=nxt[j];
                if(s[j+1]==ch) ++j;
                ++a.p[i][j];
            }
        }
        return a;
    }
}using namespace Matrix;
namespace WSN{
    inline short main(){
        n=read(); m=read(); mod=read(); scanf("%s",s+1);
        Ma a=kmp(); a=qmo(a,n);
        for(int i=0;i<m;i++) (ans+=a.p[0][i])%=mod;
        write(ans);
        return 0;
    }
}
signed main(){return WSN::main();}