hdu2825（状压dp+AC自动机）

题目网址https://vjudge.net/problem/HDU-2825

 

用ac自动机来存储状态 ，用dp来转移

 

dp[i][j][k] 表示当长度为i时，串末尾为 自动机中 j 节点多对应的值，此时串中选择状态为k 的数目

转移时，下一位状态，字符长度加一，i+1，即是可以由，j +一个新的字符，获得如果这个字符正好是串的最后一个字符，那么k的状态变化为 加上该字符的状态；

我发现，百度第一个博客的代码是错的，我wa了两次，原因是，求num数组时j的最大枚举值写错了。

AC代码如下：

# include<iostream>
# include<cstdio>
# include<queue>
# include<cstring>
# include<algorithm>
#include<cstdio>
#include<string>
#include<cstring>
#include<map>
#include<vector>
using namespace std;

const int N=110;
const int mod = 20090717;
int cnt;
int ch[N][30];
int f[N],val[N];
int n,m,k;
int num[1<<11];
int dp[30][N][1<<10];
int ans;
void init()
{
    cnt=1;
    memset(ch,0,sizeof(ch));
    memset(val,0,sizeof(val));
    // memset(last,0,sizeof(last));
}
int idx(char x)
{
    if(x>='a'&&x<='z')
      return x-'a';
    else
        return 26;
}

void insertt(string p,int id)
{
    int n=p.size();
    int root=0;
    for(int i=0; i<n; i++)
    {
        int c=idx(p[i]);
        if(!ch[root][c])
        {
            ch[root][c]=cnt;
            cnt++;
        }
        root=ch[root][c];
    }
    val[root]|=(1<<id);
}

void getFail()
{
    queue<int>q;
    f[0]=0;
    for(int i=0; i<26; ++i)
    {
        if(ch[0][i])
        {
            f[ch[0][i]]=0;
           // last[ch[0][i]]=0;
            q.push(ch[0][i]);
        }
    }

    while(!q.empty())
    {
        int r=q.front();
        q.pop();
        val[r] |= val[f[r]];
        for(int i=0; i<26; ++i)
        {
            int u = ch[r][i];
            if(!u)
            {
                ch[r][i]=ch[f[r]][i];
                continue;
            }
            q.push(u);
            int v = f[r];
            while(v && !ch[v][i]) v = f[v];
            f[u] = ch[f[r]][i];
            //val[u] |= val[f[u]];
           // last[u] = val[f[u]]? f[u] :last[f[u]];
        }
    }
}

void getnum()
{
    for(int i=0; i<(1<<10); i++)
    {
        num[i] =0 ;
        for(int j = 0; j<10; j++)
        {
            if(i&(1<<j))
                num[i]++;
        }
    }
}
int solve()
{
    memset(dp,0,sizeof(dp));
    dp[0][0][0] = 1;
    for(int i=0; i<n ;i++)
    {
        for(int j=0; j<cnt; j++)
        {
            for(int k=0; k<(1<<m);k++)
            {
                if(dp[i][j][k])
                {

                    for(int c=0; c<26;c++)
                    {
                        int nowi = i+1;
                        int nowj = ch[j][c];
                        int nowk = k|val[nowj];
                        dp[nowi][nowj][nowk]+=dp[i][j][k];
                        dp[nowi][nowj][nowk]%=mod;
                    }
                }
            }
        }
    }
    int ans = 0;
    for(int i=0; i<(1<<m); i++)
    {
        if(num[i]<k)continue;
        for(int j=0; j<cnt; j++)
        {
            ans = (ans + dp[n][j][i])%mod;
        }
    }
    return ans%mod;
}
int main()
{
    getnum();

    int t;
    while(~scanf("%d%d%d",&n,&m,&k))
   {
       if(n==0&&m==0&&k==0)
        break;
    init();
    ans = 0;
    string p;
    for(int i=0; i<m; ++i)
    {
        cin>>p;
        insertt(p,i);
    }
    getFail();

    int ans = solve();
    cout<<ans<<endl;
   }
    return 0;
}