AC自动机模板

最新推荐文章于 2024-04-28 22:37:46 发布

Till_Gloam

最新推荐文章于 2024-04-28 22:37:46 发布

阅读量267

点赞数 1

CC 4.0 BY-SA版权

分类专栏： AC自动机文章标签： AC自动机

本文链接：https://blog.youkuaiyun.com/Monkey_king2017cn/article/details/73693376

AC自动机专栏收录该内容

2 篇文章

订阅专栏

本文介绍了一种基于前缀树（Trie）和Aho-Corasick自动机的模板匹配算法，该算法能够高效地在文本中查找并统计多个模式串出现的频次。通过构建前缀树和设置失败指针，算法可以实现对大量模板串的有效检索。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

大意：
给你一些模板串，然后给你一些文本串，对于每个文本串，输出它中间出现次数最多的模板串和次数（可能有多个）；

这当然是模板。。。。。。。。。。。。。。。。。。

代码：

#include<queue>
#include<cstdio>
#include<cstring>

#define maxn 15000
#define maxc 1000005
#define For(a, b, c) for(a = b; a <= (int)(c); ++a)
using namespace std;
int n;
char T[maxc], P[151][71];
inline int idx(char c){ return c - 'a'; }

struct Trie{
    int sz, ch[maxn][26], val[maxn];
    inline void init(){
        sz = 0;
        memset(ch, 0, sizeof(ch));
        memset(val, 0, sizeof(val));
    }
    inline void insert(char *s, int v){
        int i, u = 0, l = strlen(s);
        For(i, 0, l - 1){
            int c = idx(s[i]);
            if(!ch[u][c]) ch[u][c] = ++sz;
            u = ch[u][c];
        }
        val[u] = v;
    }
}tr;

struct AC_machine{
    int fail[maxn], last[maxn], cnt[maxn];
    inline void init(){
        memset(fail, 0, sizeof(fail));
        memset(last, 0, sizeof(last));
        memset(cnt, 0, sizeof(cnt));
    }
    inline void getfail(){
        int c, u, k, r;
        queue<int> q;
        fail[0] = 0;
        For(c, 0, 25){
            if(tr.ch[0][c]) q.push(tr.ch[0][c]);
        }
        while(!q.empty()){
            r = q.front(), q.pop();
            For(c, 0, 25){
                u = tr.ch[r][c];
                if(!u){ tr.ch[r][c] = tr.ch[fail[r]][c]; continue; }
                q.push(u);
                k = fail[r];
                while(k && !tr.ch[k][c]) k = fail[k];
                fail[u] = tr.ch[k][c];
                last[u] = tr.val[fail[u]] ? fail[u] : last[fail[u]];
            }
        }
    }
    inline void calc(int x){
        while(x){
            ++cnt[tr.val[x]];
            x = last[x];
        }
    }
    inline void find(char *s){
        int i, j, l = strlen(s);
        For(i, 0, l - 1){
            j = tr.ch[j][idx(s[i])];
            if(tr.val[j]) calc(j);
            else if(last[j]) calc(last[j]);
        }
    }
}ac;

int main(){
    int i;
    while(scanf("%d", &n), n){
        tr.init();
        ac.init();
        For(i, 1, n) scanf("%s", P[i]), tr.insert(P[i], i);
        scanf("%s", T);
        ac.getfail();
        ac.find(T);
        int best = -1;
        For(i, 1, n) if(ac.cnt[i] > best) best = ac.cnt[i];
        printf("%d\n", best);
        For(i, 1, n) if(ac.cnt[i] == best){
            printf("%s\n", P[i]);
        }
    }
    return 0;
}