【HDU2825】Wireless Password,AC自动机+状态压缩DP

传送门
写在前面:依旧有毒的HDU,读题半天看不懂啊,感觉自动机的精通遥遥无期啊!
思路:
题意——给定m个给定的串(可能相同),要求组建一个给定长度为n的新串,要求新串中至少包含k个给定的串(序号重复算一个),求组成新串的方案数,所有出现的字符都为’a’-‘z’
(看到数据范围这么小,先想到高维DP、状压DP和网络流)建trie树,建fail指针,然后在自动机上进行DP,f[i][j][p]代表DP到了长度为i的新串,此时在自动机上的节点为j,p是用01二进制来表示m个给定串的选取情况,选了就是1,没选就是0(所以建trie树给m个串的标号要注意是1<<(id-1),不是id)。f[i][j][p]能扩展的状态就是 f[i+1][next][p|num[next]] next是自动机上节点j能转移的所有编号(如果没有就回到fail节点找,建立fail的时候把这件事就办了),最后我们要的答案就是f[n][i][j],i是自动机上的任一节点,j在二进制下要求>=k个1,然而j的范围为0-1023,所以预处理出来0-1023所有数的二进制1的数量就可以了
注意:
初始化依旧要注意,同时仔细想想哪些需要清零,哪些能不用memset……
代码:

#include<cstdio>
#include<iostream>
#include<cstring>
#include<algorithm>
#include<queue>
#define mod 20090717
using namespace std;
int n,m,k,root=1,tot; 
char s[15];
int trie[110][26],num[4000],f[30][110][1030],fail[4000],sum[1025];
queue<int>q;
void insert(char s[],int id)
{
    int now=root;
    for (int i=0;i<strlen(s);i++)
    {
        if (!trie[now][s[i]-'a']) trie[now][s[i]-'a']=++tot;
        now=trie[now][s[i]-'a'];
    }
    num[now]|=(1<<id-1);
}
void build()
{
    int now,tmp;
    q.push(root);
    while (!q.empty())
    {
        now=q.front();
        q.pop();
        num[now]|=num[fail[now]];
        for(int i=0;i<26;i++)
        if (trie[now][i])
        {
            tmp=fail[now];
            while (tmp&&!trie[tmp][i]) tmp=fail[tmp];
            if (now!=root&&tmp) fail[trie[now][i]]=trie[tmp][i];
            else fail[trie[now][i]]=root;
            q.push(trie[now][i]);
        }
        else
        {
            if (now==root) trie[now][i]=root;
            else trie[now][i]=trie[fail[now]][i];
        }
    }
}
int main()
{
    for (int i=0;i<1024;i++)
    {
        int t=i;
        while (t)
        {
            if (t&1) sum[i]++;
            t>>=1;
        }
    }
    int x,y,z;
    scanf("%d%d%d",&n,&m,&k);
    while (n||m||k)
    {
        int ans=0;tot=1;
        memset(trie,0,sizeof(trie));
        memset(num,0,sizeof(num));
        for (int i=1;i<=m;i++)
            scanf("%s",s),
            insert(s,i);
        for (int i=0;i<=tot;i++) fail[i]=0;
        build();
        for (int i=0;i<=n;i++)
            for (int j=0;j<=tot;j++)
                for (int p=0;p<(1<<m);p++) f[i][j][p]=0;
        f[0][1][0]=1;
        bool now=0;
        for (int i=0;i<n;i++)
            for (int j=1;j<=tot;j++)
                for (int p=0;p<(1<<m);p++)
                if (f[i][j][p])
                    for (int l=0;l<26;l++)
                        x=i+1,y=trie[j][l],z=num[y]|p,
                        f[x][y][z]=(f[x][y][z]+f[i][j][p])%mod;
        for (int j=0;j<(1<<m);j++)
            if (sum[j]>=k)
                for (int i=1;i<=tot;i++)
                    ans=(ans+f[n][i][j])%mod;
        printf("%d\n",ans);
        scanf("%d%d%d",&n,&m,&k);
    }
}

你可能感兴趣的:(【HDU2825】Wireless Password,AC自动机+状态压缩DP)