BZOJ 1030 文本生成器

1030: [JSOI2007]文本生成器
Time Limit: 1 Sec Memory Limit: 162 MB
Submit: 4777 Solved: 1986
[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

#include<iostream>
#include<cstring>
#include<cstdio>
#include<queue>
using namespace std;
#define MAXN 10010
#define Mod 10007
int n,m,ans,sz;
bool end[MAXN]; char s[MAXN];
int son[MAXN][27],fail[MAXN],f[110][MAXN];
queue<int> q;
void Insert(char *s){
    int cur=1;
    for(int i=0;i<strlen(s);i++){
        if(!son[cur][s[i]-'A'+1]) son[cur][s[i]-'A'+1]= ++sz,cur=sz;
        else cur=son[cur][s[i]-'A'+1];
    }
    end[cur]=1;
}
void Build_Fail(){
    q.push(1);fail[1]=0;int cur=1;
    while(!q.empty()){
        int now=q.front();q.pop();
        for(int i=1;i<=26;i++){
            if(!son[now][i]) continue;
            int ff=fail[now];
            while(!son[ff][i]) ff=fail[ff];
            fail[son[now][i]]=son[ff][i];
            if(end[son[ff][i]]) end[son[now][i]]=1;
            q.push(son[now][i]);
        }
    }
}
#define LL long long
int Fast_Pow(int a,int b){
    LL ret=1;
    for(int i=b;i;i>>=1,a=a*a%Mod)
        if(i&1) ret=ret*a%Mod;
    return (int) ret;
}
int main(){
    sz=1;
    for(int i=1;i<=26;i++) son[0][i]=1;
    scanf("%d%d",&n,&m);
    for(int i=1;i<=n;i++) scanf("%s",s),Insert(s);
    Build_Fail();f[0][1]=1;
    for(int i=1;i<=m;i++)
        for(int j=1;j<=sz;j++){
            if(!end[j]&&f[i-1][j]){
                for(int k=1;k<=26;k++){
                    int tmp=j;
                    while(!son[tmp][k]) tmp=fail[tmp];
                    f[i][son[tmp][k]]+=f[i-1][j];
                    f[i][son[tmp][k]]%=Mod;
                }
            }
        }
    for(int i=1;i<=sz;i++) 
        if(!end[i]) ans+=f[m][i],ans%=Mod;
    int tot=Fast_Pow(26,m);
    printf("%d\n",(tot-ans+Mod)%Mod);
    return 0;
}
//f[i][j] 表示匹配到模式串上第i个位置，在AC自动机上第j
//个位置，合法的个数 
//统计吧不合法的 用总数减去