【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[],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=;
     memset(t[x].son,,sizeof(t[x].son));
 }

 void read_trie(int x)
 {
     scanf("%s",ss[x]+);
     int len=strlen(ss[x]+);
     int now=;
     for(int j=;j<=len;j++)
     {
         int ind=ss[x][j]-'A'+;
         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();
     while(!q.empty())
     {
         x=q.front();
         y=t[x].fail;
         for(j=;j<=;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+);
     int len=strlen(s+);
     int now=;
     for(int j=;j<=len;j++)
     {
         if(s[j]<'A'||s[j]>'Z') {now=;continue;}
         int ind=s[j]-'A'+;
         now=t[now].son[ind];
         c[t[now].num]++;
     }
 }

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

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

[HDU3065]

2016-06-17 13:46:45