bzoj1212 [HNOI2004]L语言（DP+AC自动机）

题目

求模式串对主串的最大匹配长度。

题解

DP+AC自动机
设f[i]=true表示前i位可以匹配出来，那么转移方案就是背包的，（要求：主串的后缀与模式串j完全一致），其中len[j]是一个模式串的长度。
如果要是大（bao）力DP的话，显然会很慢，因为我们要配对每一个模式串。
显然这种后缀配前缀的问题应当交由AC自动机来处理。
把匹配串逐一insert到字典树中，同时标记一下结尾，求一个fail指针，然后我们就可以来一波DP了。DP时，对于长度为i的前缀，我们通过跳fail枚举所有 模式串的前缀 与当前主串后缀相同的树节点，一不小心碰到标记过的，说明满足这个主串以i结尾的后缀与这个模式串匹配，此时符合DP方程的转移条件，考虑更新f[i]。

代码

#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
const int maxn=30;
const int maxp=10*maxn;

char s[(const int)1e6+10];int len;
struct Tr{int x,fail,dep,son[25];}tr[maxp];int cnt=1,root=1;

bool flag[maxp];
void insert()
{
    int x=root;
    for(int i=1;i<=len;i++)
    {
        int k=s[i]-'a';
        if(tr[x].son[k]==0)
        {
            tr[x].son[k]=++cnt;
            tr[cnt].dep=tr[x].dep+1;
        }
        x=tr[x].son[k];
    }
    flag[x]=true;
}

int head,tail,q[maxp];
void getfail()
{
    head=0,tail=1;q[0]=root;
    while(head<tail)
    {
        int x=q[head++];
        for(int k=0;k<26;k++)//debug 一共26个字母！ 
        {
            int y=tr[x].son[k];
            if(y==0) continue;
            else if(x==root) tr[y].fail=root;
            else
            {
                int p=tr[x].fail;//debug int p=x;
                while(p!=root && !tr[p].son[k]) p=tr[p].fail;
                //debug tr[y].fail=p;
                if(tr[p].son[k]) tr[y].fail=tr[p].son[k];
                else tr[y].fail=root;
            }
            q[tail++]=y;
        }
    }
}

int f[(const int)1e6+10];//f[i]表示长度为i的前缀能否被翻译 
void solve(int id)
{
    int x=root,ans=0;
    f[0]=id;
    for(int i=1;i<=len;i++)
    {
        int k=s[i]-'a';
        while(x!=root && !tr[x].son[k]) x=tr[x].fail;
        x= !tr[x].son[k]?root:tr[x].son[k] ;//debug x=tr[x].son[k];
        
        for(int p=x;p!=root;p=tr[p].fail)
        {
            if(flag[p] && i-tr[p].dep>=0 && f[i-tr[p].dep]==id){f[i]=id;break;}//DP转移
        }
        if(f[i]==id) ans=i;
    }
    printf("%d\n",ans);
}

int main()
{
    int n,m;
    scanf("%d%d",&n,&m);
    tr[root].dep=0;
    for(int i=1;i<=n;i++)
    {
        scanf("%s",s+1);
        len=strlen(s+1);
        insert();
    }
    getfail();
//    for(int i=1;i<=cnt;i++)
//    {
//        printf("fail %d  : %d\n",i,tr[i].fail);
//    }
    for(int i=1;i<=m;i++)
    {
        scanf("%s",s+1);
        len=strlen(s+1);
        solve(i);
    }
    return 0;
}