bzoj 1030: [JSOI2007]文本生成器

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自动机的题目,好激动啊。

首先这题要求长度为m的包含给定串中至少一个的字符串的个数。

那么不妨求它的反面,即长度为m的不包含任意一个给定串的字符串个数。

考虑dp+ac自动机。

一个长度为m的满足条件的字符串,

相当于从根节点开始,在AC自动机上走m步,且未经过不合法节点(有给定串在当前节点结束)。

f[i][j]表示字符串前i位已确定,当前走到j号点,的方案数。

ok了。

 

代码

 1 #include <cstdio>
 2 #include <cstring>
 3 #include <iostream>
 4 using namespace std;
 5 #define ref(i,x,y)for(int i=x;i<=y;++i)
 6 const int N=6010;
 7 const int mod=10007;
 8 int n,m,ans,q[N],dp[101][N],tot=1;
 9 char s[101];
10 struct node{
11     bool ok;
12     int fail,ch[26];
13 }a[N];
14 int mi(int a,int b)
15 {
16     int res=1;
17     ref(i,1,b)res=res*a%mod;
18     return res;
19 }
20 void insert()
21 {
22     int len=strlen(s),now=1;
23     ref(i,0,len-1)
24     {
25         int &t=a[now].ch[s[i]-'A'];
26         if(!t)t=++tot; now=t;
27     }
28     a[now].ok=1;
29 }
30 void getfail()
31 {
32     ref(i,0,25)a[0].ch[i]=1;
33     q[1]=1;a[1].fail=0;
34     for(int h=1,t=1,x=q[h];h<=t;++h,x=q[h])
35     ref(i,0,25)
36     {
37         if(!a[x].ch[i])continue;
38         int p=a[x].fail,s=q[++t]=a[x].ch[i];
39         while(!a[p].ch[i])p=a[p].fail;
40         a[s].fail=a[p].ch[i];
41         if(a[a[p].ch[i]].ok)a[s].ok=1;
42     }
43 }
44 void work()
45 {
46     dp[0][1]=1;
47     ref(i,1,m)ref(j,1,tot)if(!a[j].ok&&dp[i-1][j])
48     ref(k,0,25)
49     {
50         int p=j;
51         while(!a[p].ch[k])p=a[p].fail;
52         (dp[i][a[p].ch[k]]+=dp[i-1][j])%=mod;
53     }
54 }
55 int main()
56 {
57     cin>>n>>m;
58     ref(i,1,n)
59     {
60         cin>>s;
61         insert();
62     }
63     getfail();
64     work();
65     ref(i,1,tot)if(!a[i].ok&&dp[m][i])
66         (ans+=dp[m][i])%=mod;
67     ans=(mod+mi(26,m)-ans)%mod;
68     cout<<ans<<endl;
69 }

 

posted @ 2017-04-02 20:43  I_m_Eden  阅读(203)  评论(0编辑  收藏  举报