BZOJ 1030: [JSOI2007]文本生成器 AC自动机

1030: [JSOI2007]文本生成器

Time Limit: 1 Sec  Memory Limit: 162 MB
Submit: 4972  Solved: 2058
[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

AC自动机上dp

最终答案=全部方案-不可行方案

这样就好搞一些

这个题建AC自动机时要注意，如果一个节点的fail指向danger节点，那么这个节点也是danger的

dp[i][j]表示前i个字符，且后缀为trie树上由根到j的不可行方案数

#include<cmath>
#include<ctime>
#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<complex>
#include<iostream>
#include<algorithm>
#include<iomanip>
#include<vector>
#include<string>
#include<bitset>
#include<queue>
#include<set>
#include<map>
using namespace std;

inline int read()
{
	int x=0,f=1;char ch=getchar();
	while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
	while(ch<='9'&&ch>='0'){x=(x<<1)+(x<<3)+ch-'0';ch=getchar();}
	return x*f;
	
}

const int M=110,N=6010,mod=10007;

bool danger[N];

int sz=1,tr[N][27],fail[N],dp[M][N],n,m;

inline void insert(char s[])
{
	int len=strlen(s+1);
	register int i,k,now=1;
	for(i=1;i<=len;++i)
	{
		k=s[i]-'A'+1;
		if(!tr[now][k])tr[now][k]=++sz;now=tr[now][k];
		if(i==len)danger[now]=1;
	}
}

int q[N];

void build()
{
	register int head=0,tail=1,i,k,u;
	q[head]=1;
	while(head<tail)
	{
		u=q[head++];
		for(i=1;i<=26;++i)
		{
			if(!tr[u][i])tr[u][i]=tr[fail[u]][i];
			else
			{
				k=fail[u];
				fail[tr[u][i]]=tr[k][i];
				if(danger[tr[k][i]])danger[tr[u][i]]=1;
				q[tail++]=tr[u][i];
			}
		}
	}
}

void solve(int x)
{
	register int i,j;
	for(i=1;i<=sz;++i)
	{
		if(danger[i]||!dp[x-1][i])continue;
		for(j=1;j<=26;++j)(dp[x][tr[i][j]]+=dp[x-1][i])%=mod;
	}
}

char s[M];

int main()
{
	n=read();m=read();
	register int i,j;
	for(i=1;i<=26;++i)tr[0][i]=1;
	while(n--){scanf("%s",s+1);insert(s);}
	build();int ans=1;dp[0][1]=1;
	for(i=1;i<=m;++i)ans=ans*26%mod;
	for(i=1;i<=m;++i)solve(i);
	for(j=1;j<=sz;++j)if(!danger[j])ans-=dp[m][j];
	ans%=mod;ans+=mod;printf("%d\n",ans%mod);return 0;
}
/*
2 2
A
B

100
*/