AC自动机

AC自动机

求解一个目标串匹配多个模式串的问题
下面模板求解的是有多少种模式串在目标串中出现

// Created by CAD
#include <bits/stdc++.h>
#define mst(name, value) memset(name,value,sizeof(name))
using namespace std;

const int maxn=1e6+5;
namespace ac{
    int next[maxn][30],fail[maxn],end[maxn];
    int root,L;
    //新建节点
    int newnode(){
        for(int i=0;i<26;++i)
            next[L][i]=-1;
        end[L++]=0;
        return L-1;
    }
    //初始化
    void init(){
        L=0;
        root=newnode();
    }
    //插入字符串
    void insert(char buf[]){
        int len=strlen(buf);
        int now=root;
        for(int i=0;i<len;i++){
            if(next[now][buf[i]-'a']==-1)
                next[now][buf[i]-'a']=newnode();
            now=next[now][buf[i]-'a'];
        }
        end[now]++;
    }
    //构建AC自动机
    void build(){
        queue<int> Q;
        fail[root]=root;
        for(int i=0;i<26;++i)
            if(next[root][i]==-1)
                next[root][i]=root;
            else{
                fail[next[root][i]]=root;
                Q.push(next[root][i]);
            }
        //求 fail 数组
        while(!Q.empty()){
            int now=Q.front();  Q.pop();
            for(int i=0;i<26;++i)
                if(next[now][i]==-1)
                    next[now][i]=next[fail[now]][i];
                else{
                    fail[next[now][i]]=next[fail[now]][i];
                    Q.push(next[now][i]);
                }
        }
    }
    bool vis[maxn];
    //求目标串中出现了几个模式串（用vis数组标记保证复杂度为 O(n)）
    int query(char buf[]){
        mst(vis,0);
        int len=strlen(buf);
        int now=root,ans=0;
        for(int i=0;i<len;++i){
            now=next[now][buf[i]-'a'];
            int temp=now;
            //可以只对当前 now 进行标记，然后用 dfs 迭代，保证复杂度为 O(n)
            while(temp!=root&&!vis[temp]){
                ans+=end[temp];
//                end[temp]=0;
                vis[temp]=1;
                temp=fail[temp];
            }
        }
        return ans;
    }
}
char buf[maxn];
int main() {
    int n;scanf("%d",&n);
    ac::init();
    for(int i=1;i<=n;++i){
        scanf("%s",buf);
        ac::insert(buf);
    }
    ac::build();
    scanf("%s",buf);
    printf("%d
",ac::query(buf));
    return 0;
}