【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\n",(mul-ans+mod)%mod);
73     return 0;
74 }

 

 
posted @ 2016-10-25 21:57  chty  阅读(137)  评论(0编辑  收藏