AC自动机-洛谷3121 [USACO15FEB]审查（黄金）Censoring (Gold)

https://www.luogu.org/problem/show?pid=3121#sub
首先题目看清楚

FJ注意到列表中的单词不会出现一个单词是另一个单词子串的情况，这意味着每个列表中的单词在S中出现的开始位置是互不相同的

这意味着你在找答案时，当一个节点正常访问的时候，你不用去寻找fail；
你只要在匹配失败的时候取fail就好了；
关于fail的意义，可以看我AC自动机的博客；
然后我们在造trie时，顺便把一个字符串的长度记录下来；
end=一个字符串的长度；
然后我们开2个栈；
a[]记录当前搜到那个点；
b[]记录当前搜那个char；
然后如果发现当前节点是end；
直接弹出栈顶end个元素；
更新now指针，继续搜；
好了；
有个事情要注意；
我原始的AC自动机；

void makezyy(){
    for(int i=0;i<26;i++)if(T[0].nxt[i])q[++r]=T[0].nxt[i];
    while(r>l){
        int x=q[++l];
        for(int i=0;i<26;i++)if(T[x].nxt[i]){
            q[++r]=T[x].nxt[i];
            int y=T[x].fail;
            while(y&&!T[y].nxt[i])y=T[y].fail;
            T[q[r]].fail=T[y].nxt[i];
        }
    }
}
void find(int m){
    int o=0;
    for(int i=0;i<m;i++){
        int x=s[i]-'a';
        int y=o;
        while(y&&!T[y].nxt[x])y=T[y].fail;
        o=T[y].nxt[x];
        a[++tot]=o;
        b[tot]=i;
        if(T[o].E){
            tot-=T[o].E;
            o=a[tot];
        }    
    }
}

大家看，这个好理解，但是用了while循环，会被卡一个点；
那我们看看更优的写法；

void makezyy(){
    for(int i=0;i<26;i++)if(T[0].nxt[i])q[++r]=T[0].nxt[i];
    while(r>l){
        int x=q[++l];
        for(int i=0;i<26;i++)
            if(!T[x].nxt[i])T[x].nxt[i]=T[T[x].fail].nxt[i];
            else{
                q[++r]=T[x].nxt[i];
                T[T[x].nxt[i]].fail=T[T[x].fail].nxt[i];
            }
    }
}
void find(int m){
    int o=0;
    for(int i=0;i<m;i++){
        o=T[o].nxt[s[i]-'a'];
        a[++tot]=o;
        b[tot]=i;
        if(T[o].E){
            tot-=T[o].E;
            o=a[tot];
        }   
    }
}

种方法巧妙地运用地推；
两种方法求出的fail值是相同的；
但是这种方法会无法辨别出某个节点在一开始有没有某个儿子；
因为makezyy之后空的儿子自动变成了fail；
好像还可以优化的吧；

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
using namespace std;
struct trie{
    int nxt[26],fail,E;
}T[100005];
int q[100005],l,r,a[100005],b[100005],tot;
int n,m,ll;
char s[100005],c[100005];
void init(int m){
    int o=0;
    for(int i=0;i<m;i++){
        int x=c[i]-'a';
        if(!T[o].nxt[x])T[o].nxt[x]=++ll;
        o=T[o].nxt[x];
    }T[o].E=m;
}
void makezyy(){
    for(int i=0;i<26;i++)if(T[0].nxt[i])q[++r]=T[0].nxt[i];
    while(r>l){
        int x=q[++l];
        for(int i=0;i<26;i++)
            if(!T[x].nxt[i])T[x].nxt[i]=T[T[x].fail].nxt[i];
            else{
                q[++r]=T[x].nxt[i];
                T[T[x].nxt[i]].fail=T[T[x].fail].nxt[i];
            }
    }
}
void find(int m){
    int o=0;
    for(int i=0;i<m;i++){
        o=T[o].nxt[s[i]-'a'];
        a[++tot]=o;
        b[tot]=i;
        if(T[o].E){
            tot-=T[o].E;
            o=a[tot];
        }   
    }
}
int main()
{
    scanf("%s",s);
    scanf("%d",&n);
    while(n--)scanf("%s",c),init(strlen(c));
    makezyy();
    find(strlen(s));
    for(int i=1;i<=tot;i++)cout<<s[b[i]];
}