【bzoj1030】[JSOI2007]文本生成器

1030: [JSOI2007]文本生成器

Time Limit: 1 Sec Memory Limit: 162 MB
Submit: 3797 Solved: 1553
[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自动机，然后在AC自动机上做dp，f[i][j]表示当前长度为i，匹配到AC自动机上第j个结点时的方案数。

则f[i][temp]=f[i][temp]+f[i-1][j] (temp是j结点往后匹配到的下一个结点)

 1 #include<iostream>
 2 #include<cstdio>
 3 #include<cstdlib>
 4 #include<cstring>
 5 #include<ctime>
 6 #include<cmath>
 7 #include<algorithm>
 8 using namespace std;
 9 #define MAXN 110
10 #define mod 10007
11 struct node{int v,fail,Link[26];}tr[MAXN*60];
12 int n,m,len,ans,mul(1),q[MAXN*60],f[MAXN][MAXN*60];
13 char ch[MAXN];
14 inline int read()
15 {
16     int x=0,f=1;  char ch=getchar();
17     while(!isdigit(ch))  {if(ch=='-')  f=-1;  ch=getchar();}
18     while(isdigit(ch))  {x=x*10+ch-'0';  ch=getchar();}
19     return x*f;
20 }
21 void insert(int k,int x)
22 {
23     int dex=ch[k]-'A';
24     if(!tr[x].Link[dex]) tr[x].Link[dex]=++len;
25     int y=tr[x].Link[dex];
26     if(k==strlen(ch+1))  {tr[y].v++;  return;}
27     insert(k+1,y);
28 }
29 void build()
30 {
31     int head=0,tail=1;
32     while(++head<=tail)
33     {
34         int x=q[head];
35         for(int i=0;i<26;i++)
36             if(tr[x].Link[i])
37             {
38                 int y=tr[x].Link[i];
39                 q[++tail]=y;
40                 if(!x) continue;
41                 int temp=tr[x].fail;
42                 while(temp&&!tr[temp].Link[i])  temp=tr[temp].fail;
43                 tr[y].fail=tr[temp].Link[i];
44                 if(tr[tr[temp].Link[i]].v) tr[y].v=1;
45             }
46     }
47 }
48 void dp()
49 {
50     f[0][0]=1;
51     for(int i=1;i<=m;i++)
52         for(int j=0;j<=len;j++)
53             if(f[i-1][j]&&(!tr[j].v))
54                 for(int k=0;k<26;k++)
55                 {
56                     int temp=j;
57                     while(temp&&!tr[temp].Link[k]) temp=tr[temp].fail;
58                     temp=tr[temp].Link[k];
59                     if(!tr[temp].v)  f[i][temp]=(f[i][temp]+f[i-1][j])%mod;
60                 }
61 }
62 int main()
63 {
64     //freopen("cin.in","r",stdin);
65     //freopen("cout.out","w",stdout);
66     n=read();  m=read();
67     for(int i=1;i<=n;i++){scanf("%s",ch+1);insert(1,0);}
68     build();
69     dp();
70     for(int i=0;i<=len;i++)  ans=(ans+f[m][i])%mod;
71     for(int i=1;i<=m;i++)  mul=(mul*26)%mod;
72     printf("%d
",(mul-ans+mod)%mod);
73     return 0;
74 }