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;}