AC自动机

AC自动机(Aho-Corasick automaton)是一种多模式匹配算法。
字典树的典型应用是给定n个单词和一个包含m个单词的文章,统计总共有多少单词出现在文章中。
正常利用KMP算法,该问题的复杂度也达到了O(m^2),但AC自动机利用fail指针避免了多模式匹配下的回退问题,时间复杂度仅为O(n)。
AC自动机的构建
1.构建Trie(字典树)
和正常的Trie结构相比,AC自动机的字典树结构中多了一个fail指针。
C语言模板(仅包含26个小写字母的字典树);

struct Node{
    Node *son[26];
    Node *fail;
    int cnt;
    Node(){
        for(int i=0;i<26;i++){
            son[i]=NULL;
        }
        fail=NULL;
        cnt=0;
    }
    BuildTree(const char* str){
        Node *now=this;
        int len=strlen(str);
        for(int i=0;i<len;i++){
            int num=str[i]-'a';
            if(now->son[num]==NULL){
                now->son[num]=new Node();
            }
            now=now->son[num];
        }
        now->cnt=1;
    }
} *p[MAX_L];

2.fail数组的构建

void BuildAC(){
    int head,tail;
    head=0;tail=0;
    Node *now;
    Node *f;
    p[head++]=root;
    while(head!=tail){
        now=p[tail++];
        if(now==root){
            for(int i=0;i<26;i++){
                if(now->son[i]!=NULL)
                {
                    now->son[i]->fail=root;
                    p[head++]=now->son[i];
                }
            }
        }else{
            for(int i=0;i<26;i++){
                f=now->fail;
                if(now->son[i]!=NULL){
                    while(f!=NULL && f->son[i]==NULL){
                        f=f->fail;
                    }
                    if(f==NULL)
                        now->son[i]->fail=root;
                    else
                    {
                        now->son[i]->fail=f->son[i];
                    }
                    p[head++]=now->son[i];
                }
            }
        }

    }

}

模板

#include<iostream>
#include<cstdio>
#include<cstring>
using namespace std;
const int MAX_L=500500;
const int MAX_N=10050;
const int MAX_M=1000050;
struct Node{
    Node *son[26];
    Node *fail;
    int cnt;
    Node(){
        for(int i=0;i<26;i++){
            son[i]=NULL;
        }
        fail=NULL;
        cnt=0;
    }
    BuildTree(const char* str){
        Node *now=this;
        int len=strlen(str);
        for(int i=0;i<len;i++){
            int num=str[i]-'a';
            if(now->son[num]==NULL){
                now->son[num]=new Node();
            }
            now=now->son[num];
        }
        now->cnt++;
    }
    Show(){
        for(int i=0;i<26;i++){
            if(this->son[i]!=NULL)
            {
                printf("%c",'a'+i);
                this->son[i]->Show();
            }
        }
    }
} *p[MAX_L];
char s[MAX_M];
char d[50];
Node * root;
void BuildAC(){
    int head,tail;
    head=0;tail=0;
    Node *now;
    Node *f;
    p[head++]=root;
    while(head!=tail){
        now=p[tail++];
        if(now==root){
            for(int i=0;i<26;i++){
                if(now->son[i]!=NULL)
                {
                    now->son[i]->fail=root;
                    p[head++]=now->son[i];
                }
            }
        }else{
            for(int i=0;i<26;i++){
                f=now->fail;
                if(now->son[i]!=NULL){
                    while(f!=NULL && f->son[i]==NULL){
                        f=f->fail;
                    }
                    if(f==NULL)
                        now->son[i]->fail=root;
                    else
                    {
                        now->son[i]->fail=f->son[i];
                    }
                    p[head++]=now->son[i];
                }
            }
        }

    }

}
int SearchS(const char* str){
    int ans=0;
    int len=strlen(str);
    Node *now=root;
    Node *tmp;
    for(int i=0;i<len;i++)
    {
        while(now!=NULL&&now->son[str[i]-'a']==NULL){
            now=now->fail;
        }
        if(now==NULL)
            now=root;
        if(now->son[str[i]-'a']!=NULL){
            now=now->son[str[i]-'a'];
            tmp=now;
            while(tmp!=root&&tmp->cnt!=-1){
                ans+=tmp->cnt;
                tmp->cnt=-1;
                tmp=tmp->fail;
            }
        }
    }
    return ans;
}
int main(){
    int m,n,ans;
    scanf("%d",&m);
    while(m--){
        root=new Node();
        scanf("%d",&n);
        while(n--){
            scanf("%s",d);
            root->BuildTree(d);
        }
         scanf("%s",s);
         BuildAC();
         ans=SearchS(s);
         printf("%d\n",ans);
    }

    return 0;
}

最后编辑于
©著作权归作者所有,转载或内容合作请联系作者
  • 序言:七十年代末,一起剥皮案震惊了整个滨河市,随后出现的几起案子,更是在滨河造成了极大的恐慌,老刑警刘岩,带你破解...
    沈念sama阅读 196,099评论 5 462
  • 序言:滨河连续发生了三起死亡事件,死亡现场离奇诡异,居然都是意外死亡,警方通过查阅死者的电脑和手机,发现死者居然都...
    沈念sama阅读 82,473评论 2 373
  • 文/潘晓璐 我一进店门,熙熙楼的掌柜王于贵愁眉苦脸地迎上来,“玉大人,你说我怎么就摊上这事。” “怎么了?”我有些...
    开封第一讲书人阅读 143,229评论 0 325
  • 文/不坏的土叔 我叫张陵,是天一观的道长。 经常有香客问我,道长,这世上最难降的妖魔是什么? 我笑而不...
    开封第一讲书人阅读 52,570评论 1 267
  • 正文 为了忘掉前任,我火速办了婚礼,结果婚礼上,老公的妹妹穿的比我还像新娘。我一直安慰自己,他们只是感情好,可当我...
    茶点故事阅读 61,427评论 5 358
  • 文/花漫 我一把揭开白布。 她就那样静静地躺着,像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上,一...
    开封第一讲书人阅读 46,335评论 1 273
  • 那天,我揣着相机与录音,去河边找鬼。 笑死,一个胖子当着我的面吹牛,可吹牛的内容都是我干的。 我是一名探鬼主播,决...
    沈念sama阅读 36,737评论 3 386
  • 文/苍兰香墨 我猛地睁开眼,长吁一口气:“原来是场噩梦啊……” “哼!你这毒妇竟也来了?” 一声冷哼从身侧响起,我...
    开封第一讲书人阅读 35,392评论 0 254
  • 序言:老挝万荣一对情侣失踪,失踪者是张志新(化名)和其女友刘颖,没想到半个月后,有当地人在树林里发现了一具尸体,经...
    沈念sama阅读 39,693评论 1 294
  • 正文 独居荒郊野岭守林人离奇死亡,尸身上长有42处带血的脓包…… 初始之章·张勋 以下内容为张勋视角 年9月15日...
    茶点故事阅读 34,730评论 2 312
  • 正文 我和宋清朗相恋三年,在试婚纱的时候发现自己被绿了。 大学时的朋友给我发了我未婚夫和他白月光在一起吃饭的照片。...
    茶点故事阅读 36,512评论 1 326
  • 序言:一个原本活蹦乱跳的男人离奇死亡,死状恐怖,灵堂内的尸体忽然破棺而出,到底是诈尸还是另有隐情,我是刑警宁泽,带...
    沈念sama阅读 32,349评论 3 314
  • 正文 年R本政府宣布,位于F岛的核电站,受9级特大地震影响,放射性物质发生泄漏。R本人自食恶果不足惜,却给世界环境...
    茶点故事阅读 37,750评论 3 299
  • 文/蒙蒙 一、第九天 我趴在偏房一处隐蔽的房顶上张望。 院中可真热闹,春花似锦、人声如沸。这庄子的主人今日做“春日...
    开封第一讲书人阅读 29,017评论 0 19
  • 文/苍兰香墨 我抬头看了看天上的太阳。三九已至,却和暖如春,着一层夹袄步出监牢的瞬间,已是汗流浃背。 一阵脚步声响...
    开封第一讲书人阅读 30,290评论 1 251
  • 我被黑心中介骗来泰国打工, 没想到刚下飞机就差点儿被人妖公主榨干…… 1. 我叫王不留,地道东北人。 一个月前我还...
    沈念sama阅读 41,706评论 2 342
  • 正文 我出身青楼,却偏偏与公主长得像,于是被迫代替她去往敌国和亲。 传闻我的和亲对象是个残疾皇子,可洞房花烛夜当晚...
    茶点故事阅读 40,904评论 2 335