3.7bzoj1030 [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

题解：这道题有多个字串，为了在后面的处理中能一起处理，故很明显使用AC自动机。

使用了时候，我们发现，如果直接处理这个问题的话难度不小，所以不放采用取补集的方式，即算出不含的有多少种，然后拿总方案数减去即可

处理“算出不含的有多少种”这个问题时，明显是个DP

我们记f[i][j]表示当前走到i点,总共j步的方案数
那么如果一个点可以到达的话（也就是这个点的标记tag值为0，即非末尾节点），我们就可以把它的方案数加给它的子节点，这样总方案数=sigma（f[i][m]）
初始状态为第一个点（根节点）走了1步，方案数为1种

#include <iostream>
#include <stdio.h>
#include <string.h>
#include <stdlib.h>
#include <memory.h>
#include <math.h>
#include <queue>
#include <stack>
#include <map>
#include <vector>
#include <limits.h>
#include <malloc.h>
#include <ctype.h>
#include <float.h>
using namespace std;

int i,trilen,j,k,n,m;
char str[105];
int tri[10005][30];
int fail[10005];
int tag[10005];
int f[10005][105];

void add(char *str){
    int p=1;
    int len=strlen(str);
    for(int i=0;i<len;i++){
        if(tri[p][str[i]-'A']==0){
            trilen++;
            tri[p][str[i]-'A']=trilen;
        }
        p=tri[p][str[i]-'A'];
    }
    tag[p]++;
}

void ACcreat(){
    queue <int> q;
    fail[1]=1;
    for(i=0;i<26;i++){
        if(tri[1][i]){
            fail[tri[1][i]]=1;
            q.push(tri[1][i]);
        }
        else
            tri[1][i]=1;
    }
    int u;
    while(!q.empty()){
        u=q.front();
        q.pop();
        for(i=0;i<26;i++){
            if(tri[u][i]){
                fail[tri[u][i]]=tri[fail[u]][i];
                q.push(tri[u][i]);
            }
            else
                tri[u][i]=tri[fail[u]][i];
        }
        tag[u]+=tag[fail[u]];
    }
}

int pow(int a,int x){
    if(x==1)
        return a;
    int tmp=pow(a,x/2);
    if(x%2==0)
        return tmp*tmp%10007;
    else
        return (tmp*tmp%10007)*a%10007;
}

void print(){
    puts("**********");
    printf("%d:::\n",trilen);
    for(int i=1;i<=trilen;i++){
        for(int j=0;j<26;j++)
            printf("%d:%d ",j,tri[i][j]);
        puts("");
    }
    puts("**********");
}

//f[i][j]表示当前走到i点,总共j步 
int main(){
    freopen("test.in","r",stdin);
    freopen("test.out","w",stdout);
    scanf("%d %d",&n,&m);
    trilen++;
    for(i=1;i<=n;i++){
        scanf("\n%s",str);
        add(str);
    }
    ACcreat();
    //print();
    f[1][0]=1;
    //for(i=1;i<=trilen;i++)
    //    printf("%d:%d\n",i,tag[i]);
    for(j=1;j<=m;j++)
        for(i=1;i<=trilen;i++)
            for(k=0;k<26;k++)
                if(!tag[tri[i][k]])
                    f[tri[i][k]][j]=(f[tri[i][k]][j]+f[i][j-1])%10007;
    /*for(i=1;i<=trilen;i++){
        for(j=1;j<=m;j++)
            printf("%d ",f[i][j]);
        puts("");
    }*/
    for(i=1;i<=trilen;i++)
        f[0][0]=(f[0][0]+f[i][m])%10007;
    printf("%d\n",(pow(26,m)-f[0][0]+10007)%10007);
    return 0;
}