hdu 2825 AC自动机＋状压dp

假设一个字符串长为n，现在在j这个位置，此时已经包含了一个模式串集合（由于m只有10，所有可以用状压来表示模式串），到下一个位置时，一共有26种情况，现在你想知道的是多一个字母后会多几个模式串，由于是多匹配问题便想到了AC自动机，而此时只要纪录走道j这个位置时trie上走到k这个位置。所以状态便是dp[i][j][k]便是长度为i多串走道j这个trie节点同时含有k这个集合的模式串的个数。

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<iostream>
#include<queue>
#include<vector>
#define mod 20090717
using namespace std;
char c[15][15];
int ch[105][26],tot;//从0到tot
int val[105];
int to[105][26];
int fail[105];
vector<int> g[105];
int vv[105];
int dp[26][105][1025];
queue<int >q;
void init(void){
    memset(ch,0,sizeof(ch));
    memset(val,0,sizeof(val));
    memset(fail,0,sizeof(fail));
    memset(vv,0,sizeof(vv));
    memset(to,0,sizeof(to));
}
void get_trie(int n){
    int i,j,f,p;
    tot=0;
    for(i=0;i<n;i++){
        p=(int)strlen(c[i]);
        f=0;
        for(j=0;j<p;j++){
            if(ch[f][c[i][j]-'a']) f=ch[f][c[i][j]-'a'];
            else{
                ch[f][c[i][j]-'a']=++tot;
                f=tot;
            }
        }
        val[f]=i+1;
    }
}
void get_fail(void){
    int i,fa,p;
    for(i=0;i<26;i++){
        if(ch[0][i]) q.push(ch[0][i]);
    }
    while(!q.empty()){
        p=q.front();
        q.pop();
        fa=fail[p];
        for(i=0;i<26;i++){
            if(ch[p][i]){
                q.push(ch[p][i]);
                fa=fail[p];
                while(fa&&!ch[fa][i]) fa=fail[fa];
                fa=ch[fa][i];
                fail[ch[p][i]]=fa;
            }
        }
    }
    for(i=0;i<=tot;i++){
        p=i;
        while(p){
            if(val[p]) vv[i]|=(1<<(val[p]-1));
            p=fail[p];
        }
    }
    for(i=0;i<=tot;i++){
        for(int j=0;j<26;j++){
            p=i;
            while(p&&!ch[p][j]) p=fail[p];
            p=ch[p][j];
            to[i][j]=p;
        }
    }
}
int get(int x){
    int s=0;
    while(x>0){
        if(x&1) s++;
        x=x>>1;
    }
    return s;
}
int solve(int n,int m,int k){
    int i,j,p,f,x;
    p=1<<n;
    for(i=0;i<=m;i++){
        for(j=0;j<=tot;j++){
            for(f=0;f<p;f++) dp[i][j][f]=0;
        }
    }
    dp[0][0][0]=1;
    for(i=0;i<m;i++){
        for(j=0;j<=tot;j++){
            for(f=0;f<p;f++){
                if(dp[i][j][f]==0) continue;
                for(int v=0;v<26;v++){
                    dp[i+1][to[j][v]][vv[to[j][v]]|f]+=dp[i][j][f];
                    dp[i+1][to[j][v]][vv[to[j][v]]|f]%=mod;
                }
            }
        }
    }
    int s=0;
    for(i=0;i<=tot;i++){
        for(j=0;j<p;j++){
            x=get(j);
            if(x>=k) s+=dp[m][i][j];
            s%=mod;
        }
    }
    return s;
}int main()
{
    int i,n,m,k;
    while(1){
        scanf("%d%d%d",&m,&n,&k);
        if(m==0&&n==0&&k==0) break;
        init();
        for(i=0;i<n;i++) scanf("%s",c[i]);
        get_trie(n);
        get_fail();
        printf("%d\n",solve(n,m,k));
    }
}