病毒侵袭持续中 HDU - 3065 （AC自动机静态数组版本）

题目链接

病毒侵袭持续中 HDU - 3065

题目

小t非常感谢大家帮忙解决了他的上一个问题。然而病毒侵袭持续中。在小t的不懈努力下，他发现了网路中的“万恶之源”。这是一个庞大的病毒网站，他有着好多好多的病毒，但是这个网站包含的病毒很奇怪，这些病毒的特征码很短，而且只包含“英文大写字符”。当然小t好想好想为民除害，但是小t从来不打没有准备的战争。知己知彼，百战不殆，小t首先要做的是知道这个病毒网站特征：包含多少不同的病毒，每种病毒出现了多少次。大家能再帮帮他吗？
Input
第一行，一个整数N（1<=N<=1000），表示病毒特征码的个数。
接下来N行，每行表示一个病毒特征码，特征码字符串长度在1—50之间，并且只包含“英文大写字符”。任意两个病毒特征码，不会完全相同。
在这之后一行，表示“万恶之源”网站源码，源码字符串长度在2000000之内。字符串中字符都是ASCII码可见字符（不包括回车）。
Output
按以下格式每行一个，输出每个病毒出现次数。未出现的病毒不需要输出。
病毒特征码: 出现次数
冒号后有一个空格，按病毒特征码的输入顺序进行输出。
Sample Input
3
AA
BB
CC
ooxxCC%dAAAoen….END
Sample Output
AA: 2
CC: 1

Hint
Hit：
题目描述中没有被提及的所有情况都应该进行考虑。比如两个病毒特征码可能有相互包含或者有重叠的特征码段。
计数策略也可一定程度上从Sample中推测。

分析

还是一个比较裸的AC自动机，用一个ans数组记录每一个单词出现的次数即可。
单词的字母最大26种，但主串的字母最多有128种，所以树分支有128。

AC代码

//202ms 17.8MB
#include <cstdio>
#include <cstring>
#include <algorithm>
#include <set>
#include <queue>
using namespace std;
const int maxn=50*1000+100; //单词长度*单词数
const int maxlen=2e6+100; //主串长度
const int maxc=128; //树分支数
int n;
int ans[1005];
//在trie上添加后缀结点
//后缀结点：x的后缀结点为非x结点y，y的路径字符串为x的路径字符串的最长后缀且y的路径字符串是单词
struct ac_au
{
    int child[maxn][maxc];//trie
    int fail[maxn];//后缀结点
    int sta[maxn];//该路径字符串单词有无，有存编号
    int cnt;//结点的个数
    int root;
    int newnode()
    {
        for(int i=0;i<maxc;i++)
            child[cnt][i]=-1;
        sta[cnt++]=0;
        return cnt-1;
    }
    void init()
    {
        cnt=0;
        root=newnode();
    }
    void insert(char *s,int id)//trie构建
    {
        int len=strlen(s),p=0;
        for(int i=0;i<len;i++)
        {
            int c=s[i];
            int &num=child[p][c];
            if(num==-1)
                num=newnode();
            p=num;
        }
        sta[p]=id;
    }
    void build()//后缀结点的添加
    {
        int p=0;
        queue<int> que;
        fail[p]=p;
        for(int i=0;i<maxc;i++)
        {
            int &num=child[p][i];
            if(num==-1)
                num=p;//为下文铺垫，和树意义无关
            else
            {
                fail[num]=p;//有点多余，初始化都是0，理解算法过程有用
                que.push(num);
            }
        }
        while(!que.empty())
        {
            p=que.front();que.pop();
            for(int i=0;i<maxc;i++)
            {
                int &num=child[p][i];
                if(num==-1)
                    num=child[fail[p]][i];
                else
                {
                    fail[num]=child[fail[p]][i];
                    que.push(num);
                }
            }
        }
    }
    void query(char *s)//解决多模式串匹配单主串问题
    {
        int len=strlen(s),p=0;
        memset(ans,0,sizeof(ans));
        for(int i=0;i<len;i++)
        {
            int c=s[i];
            p=child[p][c];
            int tmp=p;
            while(tmp)
            {
                if(sta[tmp])
                    ans[sta[tmp]]++;
                tmp=fail[tmp];
            }
        }

    }
}ac;
char str[maxlen];
char s[1005][55];

int main()
{

    int m;
    while(~scanf("%d",&n))
    {
        ac.init();
        for(int i=1; i<=n; i++)
        {
            scanf("%s",s[i]);
            ac.insert(s[i],i);
        }
        ac.build();
        scanf("%s",str);
        ac.query(str);
        for(int i=1; i<=n; i++)
            if(ans[i])
                printf("%s: %d\n",s[i],ans[i]);
    }
    return 0;
}