HDU-3065 病毒侵袭持续中 AC自动机

病毒侵袭持续中

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


[align=left]Problem Description[/align]
小t非常感谢大家帮忙解决了他的上一个问题。然而病毒侵袭持续中。在小t的不懈努力下，他发现了网路中的“万恶之源”。这是一个庞大的病毒网站，他有着好多好多的病毒，但是这个网站包含的病毒很奇怪，这些病毒的特征码很短，而且只包含“英文大写字符”。当然小t好想好想为民除害，但是小t从来不打没有准备的战争。知己知彼，百战不殆，小t首先要做的是知道这个病毒网站特征：包含多少不同的病毒，每种病毒出现了多少次。大家能再帮帮他吗？

[align=left]Input[/align]
第一行，一个整数N（1<=N<=1000），表示病毒特征码的个数。
接下来N行，每行表示一个病毒特征码，特征码字符串长度在1—50之间，并且只包含“英文大写字符”。任意两个病毒特征码，不会完全相同。
在这之后一行，表示“万恶之源”网站源码，源码字符串长度在2000000之内。字符串中字符都是ASCII码可见字符（不包括回车）。

[align=left]Output[/align]
按以下格式每行一个，输出每个病毒出现次数。未出现的病毒不需要输出。
病毒特征码: 出现次数
冒号后有一个空格，按病毒特征码的输入顺序进行输出。

[align=left]Sample Input[/align]

3
AA
BB
CC
ooxxCC%dAAAoen....END

[align=left]Sample Output[/align]

AA: 2
CC: 1

Hint

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

　　其实该题没有想象中的那么复杂，仔细一想就知道，AC自动机自身不是有一个重置操作吗，即找的的子串曾经被我们删除过，该题只要不进行删除操作就行了，这都多亏了在该算法中，本身的fail指针是不停的回溯的，例如AAA匹配AA时，前面的AA计算一次，到达第三个A时，由于后面的AA只有两个字符，算法将自动跳到AA的第二个A来匹配AAA中的第三个A，就是这样。
　　代码如下：

#include <cstdio>
#include <cstring>
#include <cstdlib>
#include <queue>
using namespace std;

char ss[1005][55], T[2000005];

int cnt[1005];

struct Node
{
int flag;
Node *ch[26], *fail;
};

void getstr( char *s )
{
int p = 0;
char c;
while( c = getchar(), c == '\n' ) ;
s[p++] = c;
while( c = getchar(), c != '\n' )
s[p++] = c;
s[p] = '\0';
}

Node *init(  )
{
Node *p = new Node;
p->flag = -1;
p->fail = NULL;
memset( p->ch, 0, sizeof( p->ch ) );
return p;
}

void BuildTree( Node *p, char *in, int No )
{
int dx = *in - 'A';
if( *in )
{
if( !p->ch[dx] )
p->ch[dx] = init();
BuildTree( p->ch[dx], in + 1, No );
}
else
p->flag = No;
}

void BuildFailPoint( Node *r )
{
queue<Node *>q;
q.push( r );
while( !q.empty() )
{
Node *pos = q.front();
q.pop();
for( int i = 0; i <26; ++i )
{
Node *f = pos->fail;
if( pos->ch[i] )
{
while( f ) // f已经退了步，所以后面不许要让步
{
if( f->ch[i] )
{
pos->ch[i]->fail = f->ch[i];
break;
}
f = f->fail;
}
if( !f )
{
pos->ch[i]->fail = r;
}
q.push( pos->ch[i] );
}
}
}
}

void AcAuto( Node *r, char *text )
{
int len = strlen( text );
Node *f = r;
for( int i = 0; i < len; ++i )
{
int dx = text[i] - 'A';
if( dx < 0 || dx >= 26 )
{
f = r;
continue;
}
while( f )
{
if( f->ch[dx] )
{
f = f->ch[dx]; // 匹配成功，向后移一位
break;
}
f = f->fail;
}
if( !f )
f = r;
else
{
Node *t = f;
while( t )
{
if( t->flag!= -1 )
{
cnt[t->flag]++;
}
t = t->fail;
}
}
}
}

void _free( Node *p )
{
for( int i = 0; i < 26; ++i )
{
if( p->ch[i] )
_free( p->ch[i] );
}
free( p );
}

int main()
{
int N;
while( scanf( "%d", &N ) != EOF )
{
Node *r = init();
memset( cnt, 0, sizeof( cnt ) );
for( int i = 0; i < N; ++i )
{
getstr( ss[i] );
BuildTree( r, ss[i], i );
}
BuildFailPoint( r );
getstr( T );
AcAuto( r, T );
for( int i = 0; i < N; ++i )
{
if( cnt[i] )
printf( "%s: %d\n", ss[i], cnt[i] );
}
_free( r );
}
return 0;
}