【HDU3065】 病毒侵袭持续中（AC自动机）

Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 32768/32768 K (Java/Others)
Total Submission(s): 10504    Accepted Submission(s): 3701

Problem Description

小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中推测。

 

【分析】

　　做了BZOJ2434，这题就显得很简单了。反正我的做法跟上一题其实差不多，就是先建AC自动机，然后求DFS序，再走一遍目标串，打好标记，最后求一个区间的和即可。不过目测其他慢一点的方法也是可以的。

 

代码如下：

#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<iostream>
#include<algorithm>
#include<queue>
using namespace std;
#define Maxn 100100
#define Maxl 2000010
#define Ml 60
#define Mn 1010

char ss[Mn][Ml];
char s[Maxl];
int p[Maxn],c[Maxn];
int n;

struct hp
{
    int x,y,next;
}qy[Maxn];int al;

struct node
{
    int son[30],cnt,fail;
    int num,rt;
}t[Maxn];int tot;

int first[Maxn];

void ins(int x,int y)
{
    qy[++al].x=x;qy[al].y=y;
    qy[al].next=first[x];first[x]=al;
}

void upd(int x)
{
    t[x].cnt=0;
    memset(t[x].son,0,sizeof(t[x].son));
}

void read_trie(int x)
{
    scanf("%s",ss[x]+1);
    int len=strlen(ss[x]+1);
    int now=0;
    for(int j=1;j<=len;j++)
    {
        int ind=ss[x][j]-'A'+1;
        if(!t[now].son[ind])
         t[now].son[ind]=++tot,upd(tot); 
        now=t[now].son[ind];
        if(j==len) t[now].cnt++,p[x]=now;
    }
}

queue<int > q;
void build_AC()
{
    int i,j,x,y;
    while(!q.empty()) q.pop();
    q.push(0);
    while(!q.empty())
    {
        x=q.front();
        y=t[x].fail;
        for(j=1;j<=26;j++)
        {
            if(t[x].son[j])
            {
                q.push(t[x].son[j]);
                t[t[x].son[j]].fail=x?t[y].son[j]:x;
                ins(t[t[x].son[j]].fail,t[x].son[j]);
            }
            else t[x].son[j]=t[y].son[j];
        }
        q.pop();
    }
}

void dfs(int x)
{
    t[x].num=++al;t[x].rt=t[x].num;
    for(int i=first[x];i;i=qy[i].next)
      dfs(qy[i].y),t[x].rt=t[qy[i].y].rt;
}

void dfs2(int x)
{
    scanf("%s",s+1);
    int len=strlen(s+1);
    int now=0;
    for(int j=1;j<=len;j++)
    {
        if(s[j]<'A'||s[j]>'Z') {now=0;continue;}
        int ind=s[j]-'A'+1;
        now=t[now].son[ind];
        c[t[now].num]++;
    }
}

void init()
{
    memset(first,0,sizeof(first));
    tot=0;upd(0);
    for(int i=1;i<=n;i++)
    {
        read_trie(i);
    }
    al=0;build_AC();al=0;
    dfs(0);
    memset(c,0,sizeof(c));
    dfs2(0);c[0]=0;
    for(int i=2;i<=al;i++) c[i]=c[i-1]+c[i];
    for(int i=1;i<=n;i++)
    {
        int ans=0;
        ans=c[t[p[i]].rt]-c[t[p[i]].num-1];
        if(ans!=0) printf("%s: %d
",ss[i]+1,ans);
    }
}

int main()
{
    while(scanf("%d",&n)!=EOF)
    {
        init();
    }
    return 0;
}

2016-06-17 13:46:45