BZOJ 1030-文本生成器（DP+AC自动机）

1030: [JSOI2007]文本生成器

Time Limit: 1 Sec   Memory Limit: 162 MB
Submit: 4888   Solved: 2030
[ Submit][ Status][ Discuss]

Description

　　JSOI交给队员ZYX一个任务，编制一个称之为“文本生成器”的电脑软件：该软件的使用者是一些低幼人群，
他们现在使用的是GW文本生成器v6版。该软件可以随机生成一些文章―――总是生成一篇长度固定且完全随机的文
章—— 也就是说，生成的文章中每个字节都是完全随机的。如果一篇文章中至少包含使用者们了解的一个单词，
那么我们说这篇文章是可读的（我们称文章a包含单词b，当且仅当单词b是文章a的子串）。但是，即使按照这样的
标准，使用者现在使用的GW文本生成器v6版所生成的文章也是几乎完全不可读的?。ZYX需要指出GW文本生成器 v6
生成的所有文本中可读文本的数量，以便能够成功获得v7更新版。你能帮助他吗？

Input

　　输入文件的第一行包含两个正整数，分别是使用者了解的单词总数N (<= 60)，GW文本生成器 v6生成的文本固
定长度M；以下N行，每一行包含一个使用者了解的单词。这里所有单词及文本的长度不会超过100，并且只可能包
含英文大写字母A..Z

Output

　　一个整数，表示可能的文章总数。只需要知道结果模10007的值。

Sample Input

2 2
A
B

Sample Output

100

HINT

距离学AC自动机差不多满一年了，想当年学的AC自动机基本上没怎么理解，于是最近有花时间搞了搞

在加深对AC自动机的理解的同时，也学会了不同的姿势，AC自动机思想及工作原理见：大牛的博客

（虽然觉得讲解的有一点点问题），可是还是理解不足吧，还是要多搞搞。

言归正传，说说这题的思路：

题目是让你求一篇文章含有给定单词（任意数量）的方案数，求方案数当然是DP了，然而难点在于求一

篇文章里是否含有该单词。咦，这不就是AC自动机要处理的吗（好吧，模板早忘了，以下模板借鉴了别人，

以后应该都会用这个了。。。。）

首先将建一棵树将可读单词插入，然后构造失败指针。。。。这些都是AC自动机常规步骤就不多说了，

剩下的就是如何DP的问题了，简单想想（其实一点都不简单！！），发现正着求方案数很麻烦，怎么个

麻烦法，比如说你要求的文章中可能含若干个数量的可读单词，单词与单词之间也会影响最终的答案，光

想想都很头疼，呢该怎么求呢，答案很简单，倒过来求一发即可。。

dp[i][j]:  长度为i串的在第j个节点没有经过任何一个end节点的方案数 （真的很机智啊啊啊）

#include<map>
#include<stack>
#include<queue>
#include<vector>
#include<math.h>
#include<stdio.h>
#include<iostream>
#include<string.h>
#include<stdlib.h>
#include<algorithm>
using namespace std;
typedef long long  ll;
#define  maxn  6999
#define  mod 10007
struct node
{
	int next[30];
	int fail;
	int end;
}tr[maxn];
queue<int>q;
int dp[105][maxn],cnt;
char s[105];
void insert()
{
	int p=0,i;
	for(i=1;s[i]!='\0';i++)
	{
		if(tr[p].next[s[i]-'A'+1]==0)
			tr[p].next[s[i]-'A'+1]=++cnt;
		p=tr[p].next[s[i]-'A'+1];
	}
	tr[p].end=1;
}
void build()
{
	int tmp,i,j;
	q.push(0);
	while(q.empty()==0)
	{
		tmp=q.front();
		q.pop();
		for(i=1;i<=26;i++)
		{
			if(tr[tmp].next[i]==0)
				tr[tmp].next[i]=tr[tr[tmp].fail].next[i];
			else
			{
				if(tmp)
					tr[tr[tmp].next[i]].fail=tr[tr[tmp].fail].next[i];
				tr[tr[tmp].next[i]].end+=tr[tr[tr[tmp].fail].next[i]].end;
				q.push(tr[tmp].next[i]);
			}
		}
	}
}
int main(void)
{
	int n,m,i,j,ans=1,k;
	scanf("%d%d",&n,&m);
	for(i=1;i<=n;i++)
	{
		scanf("%s",s+1);
		insert();
	}
	build();dp[0][0]=1;
	for(i=1;i<=m;i++)
	{
		ans*=26;
		ans%=mod;
		for(j=0;j<=cnt;j++)
			for(k=1;k<=26;k++)
				if(tr[j].end==0)
				{
					dp[i][tr[j].next[k]]+=dp[i-1][j];
					dp[i][tr[j].next[k]]%=mod;
				}
	}   
	for(i=0;i<=cnt;i++)
		if(tr[i].end==0)
			ans=(ans-dp[m][i])%mod;
	printf("%d\n",(ans+mod)%mod);
	return 0;
}