hdu2825 Wireless Password（AC自动机上dp）

题意：

给定n，m，k，
有m个字符串，每个字符串都是一个咒语，
现在要构造一个长度为n的字符串，满足这个字符串至少包含k个咒语（咒语是这个串的子串就说明包含）。
问有多少种构造方法，答案对20190717取模。

数据范围：n<=25，m<=10，每条咒语长度不超过10.

解法：

对m个咒语建立ac自动机,
令d[i][j][k]表示从根开始走了i步,当前在节点j上,包含的咒语集合为k的方案数.
因为m<=10,因此k用长度为10的二进制数表示.

在ac自动机上dp就行了,
建立ac自动机的时候需要预处理每条路径的咒语状态集合.

最坏复杂度估计:
n*ac自动机节点数*状态总数*枚举下一步转移
=25*100*1000*26=6e7

code：

#include<bits/stdc++.h>
using namespace std;
const int maxm=1e3+5;
const int mod=20090717;
char s[maxm];
int n,m,k;
int cal(int x){//计算二进制中1的位数
    int ans=0;
    while(x){
        ans++;
        x&=(x-1);
    }
    return ans;
}
struct AC{
    int a[maxm][26];
    int fail[maxm];
    int val[maxm];
    int tot=0;
    void init(){
        for(int i=0;i<=tot;i++){
            memset(a[i],0,sizeof a[i]);
            fail[i]=val[i]=0;
        }
        tot=0;
    }
    void add(char *s,int k){
        int len=strlen(s);
        int node=0;
        for(int i=0;i<len;i++){
            int v=s[i]-'a';
            if(!a[node][v])a[node][v]=++tot;
            node=a[node][v];
        }
        val[node]|=(1<<k);
    }
    void build(){
        queue<int>q;
        for(int i=0;i<26;i++){
            if(a[0][i]){
                q.push(a[0][i]);
                fail[a[0][i]]=0;
            }
        }
        while(!q.empty()){
            int x=q.front();q.pop();
            for(int i=0;i<26;i++){
                if(a[x][i]){
                    fail[a[x][i]]=a[fail[x]][i];
                    val[a[x][i]]|=val[a[fail[x]][i]];
                    q.push(a[x][i]);
                }else{
                    a[x][i]=a[fail[x]][i];
                }
            }
        }
    }
    //
    int d[30][maxm][(1<<10)];
    int solve(){
        //init
        for(int i=0;i<=n;i++){
            for(int j=0;j<=tot;j++){
                for(int k=0;k<(1<<m);k++){
                    d[i][j][k]=0;
                }
            }
        }
        //dp
        d[0][0][0]=1;
        for(int i=0;i<n;i++){//步数
            for(int j=0;j<=tot;j++){//节点位置
                for(int s=0;s<(1<<m);s++){//状态
                    if(!d[i][j][s])continue;
                    for(int x=0;x<26;x++){//枚举下一步
                        int v=a[j][x];
                        d[i+1][v][s|val[v]]=(d[i+1][v][s|val[v]]+1ll*d[i][j][s])%mod;
                    }
                }
            }
        }
        //cal
        int ans=0;
        for(int j=0;j<=tot;j++){
            for(int s=0;s<(1<<m);s++){
                if(cal(s)<k)continue;
                ans=(ans+1ll*d[n][j][s])%mod;
            }
        }
        return ans;
    }
}ac;
signed main(){
    while(scanf("%d%d%d",&n,&m,&k)!=EOF){
        if(n+m+k==0)break;
        ac.init();
        for(int i=0;i<m;i++){
            scanf("%s",s);
            ac.add(s,i);
        }
        ac.build();
        //
        int ans=ac.solve();
        printf("%d\n",ans);
    }
    return 0;
}

---