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

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

 

　　其实该题没有想象中的那么复杂，仔细一想就知道，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;
}