http://www.elijahqi.win/archives/2878
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
因为求>=1的不好求 所以转化成求有多少个不会出现 然后用总数减一下即可
设dp[i][j]表示当前匹配到文本的第几位 然后我现在处在AC自动机上第j 个节点上不会匹配形成模式串的方案数 直接转移即可 设trans[x][i]表示第x号节点转移到第i个字符应该是什么 这样的好处是避免了每次跳fail指针的复杂度
#include<queue>
#include<cstdio>
#include<cctype>
#include<cstring>
#include<algorithm>
#define mod 10007
#define N 110
using namespace std;
inline int read(){
int x=0;char ch=getchar();
while(!isdigit(ch)) ch=getchar();
while(isdigit(ch)) x=x*10+ch-'0',ch=getchar();
return x;
}
int cnt=1,trans[N*N][26],dp[N][N*N],fail[N*N];
bool end[N*100];
inline void insert1(char s[]){
int len=strlen(s+1);int now=1;
for (int i=1,nxt;i<=len;++i) {
if (!trans[now][s[i]-'A']) trans[now][s[i]-'A']=nxt=++cnt;
else nxt=trans[now][s[i]-'A'];now=nxt;
}end[now]=1;
}
inline void buildAC(){
queue<int>q;q.push(1);
while(!q.empty()){
int x=q.front();q.pop();
for (int i=0;i<26;++i){
if (trans[x][i]) fail[trans[x][i]]=trans[fail[x]][i],q.push(trans[x][i]);
else {trans[x][i]=trans[fail[x]][i];continue;}end[trans[x][i]]|=end[fail[trans[x][i]]];
}
}
}
inline int ksm(int b,int t){
int tmp=1;for (;t;b=b*b%mod,t>>=1) if (t&1) tmp=tmp*b%mod;return tmp;
}
char s[N];int n,m;
int main(){
freopen("bzoj1030.in","r",stdin);
n=read();m=read();for (int i=0;i<26;++i) trans[0][i]=1;dp[0][1]=1;
for (int i=1;i<=n;++i) scanf("%s",s+1),insert1(s);buildAC();
for (int i=0;i<=m;++i)
for (int j=1;j<=cnt;++j){
if (!dp[i][j]||end[j]) continue;
for (int ii=0;ii<26;++ii)
(dp[i+1][trans[j][ii]]+=dp[i][j])%=mod;
}
int ans=ksm(26,m);
for (int i=1;i<=cnt;++i) if (!end[i]) ans=ans-dp[m][i]<0?ans-dp[m][i]+mod:ans-dp[m][i];
printf("%d\n",ans);
return 0;
}