JSOI交给队员ZYX一个任务，编制一个称之为“文本生成器”的电脑软件：该软件的使用者是一些低幼人群，
他们现在使用的是GW文本生成器v6版。该软件可以随机生成一些文章―――总是生成一篇长度固定且完全随机的文
章—— 也就是说，生成的文章中每个字节都是完全随机的。如果一篇文章中至少包含使用者们了解的一个单词，
那么我们说这篇文章是可读的（我们称文章a包含单词b，当且仅当单词b是文章a的子串）。但是，即使按照这样的
标准，使用者现在使用的GW文本生成器v6版所生成的文章也是几乎完全不可读的?。ZYX需要指出GW文本生成器 v6
生成的所有文本中可读文本的数量，以便能够成功获得v7更新版。你能帮助他吗？
输入文件的第一行包含两个正整数，分别是使用者了解的单词总数N (<= 60)，GW文本生成器 v6生成的文本固
定长度M；以下N行，每一行包含一个使用者了解的单词。这里所有单词及文本的长度不会超过100，并且只可能包
含英文大写字母A..Z
一个整数，表示可能的文章总数。只需要知道结果模10007的值。
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 }
原文：http://www.cnblogs.com/chty/p/5998439.html