bzoj1030 [JSOI2007]文本生成器

http://www.elijahqi.win/archives/2878
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
因为求>=1的不好求 所以转化成求有多少个不会出现 然后用总数减一下即可
设dp[i][j]表示当前匹配到文本的第几位 然后我现在处在AC自动机上第j 个节点上不会匹配形成模式串的方案数 直接转移即可 设trans[x][i]表示第x号节点转移到第i个字符应该是什么 这样的好处是避免了每次跳fail指针的复杂度

#include<queue>
#include<cstdio>
#include<cctype>
#include<cstring>
#include<algorithm>
#define mod 10007
#define N 110
using namespace std;
inline int read(){
    int x=0;char ch=getchar();
    while(!isdigit(ch)) ch=getchar();
    while(isdigit(ch)) x=x*10+ch-'0',ch=getchar();
    return x;
}
int cnt=1,trans[N*N][26],dp[N][N*N],fail[N*N];
bool end[N*100];
inline void insert1(char s[]){
    int len=strlen(s+1);int now=1;
    for (int i=1,nxt;i<=len;++i) {
        if (!trans[now][s[i]-'A']) trans[now][s[i]-'A']=nxt=++cnt;
        else nxt=trans[now][s[i]-'A'];now=nxt;
    }end[now]=1;
}
inline void buildAC(){
    queue<int>q;q.push(1);
    while(!q.empty()){
        int x=q.front();q.pop();
        for (int i=0;i<26;++i){
            if (trans[x][i]) fail[trans[x][i]]=trans[fail[x]][i],q.push(trans[x][i]);
            else {trans[x][i]=trans[fail[x]][i];continue;}end[trans[x][i]]|=end[fail[trans[x][i]]];
        }
    }
}
inline int ksm(int b,int t){
    int tmp=1;for (;t;b=b*b%mod,t>>=1) if (t&1) tmp=tmp*b%mod;return tmp;
}
char s[N];int n,m; 
int main(){
    freopen("bzoj1030.in","r",stdin);
    n=read();m=read();for (int i=0;i<26;++i) trans[0][i]=1;dp[0][1]=1;
    for (int i=1;i<=n;++i) scanf("%s",s+1),insert1(s);buildAC();
    for (int i=0;i<=m;++i)
        for (int j=1;j<=cnt;++j){
            if (!dp[i][j]||end[j]) continue;
            for (int ii=0;ii<26;++ii)
                (dp[i+1][trans[j][ii]]+=dp[i][j])%=mod;
        }
    int ans=ksm(26,m);
    for (int i=1;i<=cnt;++i) if (!end[i]) ans=ans-dp[m][i]<0?ans-dp[m][i]+mod:ans-dp[m][i];
    printf("%d\n",ans);
    return 0;
}
  • 0
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值