码迷,mamicode.com
首页 > 其他好文 > 详细

2018/2/22 每日一学 AC自动机

时间:2018-02-22 13:38:20      阅读:135      评论:0      收藏:0      [点我收藏+]

标签:nbsp   amp   ace   return   main   cas   names   print   word   

AC自动机=trie树+KMP。

代码如下:(自行理解)

#include<cstdio>
#include<string.h>
#include<math.h>
#include<queue>
#include<algorithm>
#define N 500006
using namespace std;

char st[1000005];
char keyword[55];
int n,m;
int next[N][26],cnt[N],fail[N],pos;

int newnode()
{
    for(int i=0; i<26; i++) next[pos][i] = 0;
    fail[pos] = cnt[pos] = 0;
    return pos++;
}

void insert(char *s) //构造trie
{
    int i,p = 0;
    for(i=0; s[i]; i++)
    {
        int k = s[i] - a, &x = next[p][k];
        p = x?x : x = newnode();
    }
    cnt[p]++; // 位运算要用
}

void makenext()  // 构造失败指针
{
    int i;
    queue<int>q;
    q.push(0);
    while(!q.empty())  // 这个代码为什么会这么强?
    {
        int u = q.front();
        cnt[u] += cnt[fail[u]];

        q.pop();
        for(i=0; i<26; i++)
        {
            int v = next[u][i];
            if(v == 0)
              next[u][i] = next[fail[u]][i]; // why? 缺少v, 求next
            else
              q.push(v);

            if(u && v)
            {
                fail[v]=next[fail[u]][i]; // 为什么只弄了一下? 如果没有会怎么样? 求fail
            }
        }
    }
}

int query(char *s)
{
    int ret = 0, idx, d = 0;
    for(int i=0; s[i]; i++)
    {
        idx = s[i] - a;
        d = next[d][idx];
        ret += cnt[d];
        cnt[d] = 0;
    }
    return ret;
}

int main()
{
    int cas;
    scanf("%d",&cas);

    while(cas--)
    {
        scanf("%d",&n);
        pos = 0, newnode();
        for(int i=0; i<n; i++)
        {
            scanf("%s",keyword);
            insert(keyword);
        }

        makenext();
        scanf(" %s",st);
        printf("%d\n",query(st));
    }
}

 

2018/2/22 每日一学 AC自动机

标签:nbsp   amp   ace   return   main   cas   names   print   word   

原文地址:https://www.cnblogs.com/Alex-leaves/p/8458179.html

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!