luogu P4052 [JSOI2007]文本生成器

背景:

暑假集训七月份最后一天,写写原来没有写过的 blog \text{blog} blog吧。

题目传送门:

https://www.luogu.org/problemnew/show/P4052

题意:

给出 n n n个单词,要求这 n n n个单词至少有一个出现在长度为 m m m的文章中的方案数。

思路:

如果改成没有出现过的呢?
首先可以发现答案就是所有串的个数减去不包含可读串的串的个数。
前半部分是 2 6 m 26^m 26m。后半部分 dp \text{dp} dp求解。
d p i , j dp_{i,j} dpi,j表示匹配到第 i i i个单词,在自动机第 j j j个位置失配的方案数。
若当前不是 end \text{end} end,则可以更新状态。
最后匹配完 m m m位即可。

代码:

#include<cstdio>
#include<cstring>
#include<queue>
#include<algorithm>
#define mod 10007
using namespace std;
    int n,m,len=0;
    char s[10010];
    int dp[110][10100];
queue<int> f;
struct node
{
    int fail,tot;
    int son[30];
    node()
    {
        fail=tot=0;
        memset(son,0,sizeof(son));
    }
} tr[10010];
void build(char *s)
{
    int l=strlen(s),x=0;
    for(int i=0;i<l;i++)
    {
        if(!tr[x].son[s[i]-'A']) tr[x].son[s[i]-'A']=++len;
        x=tr[x].son[s[i]-'A'];
    }
    tr[x].tot=1;
}
void get_fail()
{
    while(!f.empty()) f.pop();
    for(int i=0;i<26;i++)
        if(tr[0].son[i])
        {
            tr[tr[0].son[i]].fail=0;
            f.push(tr[0].son[i]);
        }
    while(!f.empty())
    {
        int x=f.front();
        f.pop();
        for(int i=0;i<26;i++)
        {
            int y=tr[x].son[i];
            if(y)
            {
                tr[y].fail=tr[tr[x].fail].son[i];
                tr[y].tot|=tr[tr[y].fail].tot;
                f.push(y);
            }
            else tr[x].son[i]=tr[tr[x].fail].son[i];
        }
    }
}
int dg(int x,int k)
{
    if(!k) return 1;
    int op=dg(x,k>>1);
    if(k&1) return op*op%mod*x%mod; else return op*op%mod;
}
int main()
{
    scanf("%d %d",&n,&m);
    for(int i=1;i<=n;i++)
    {
        scanf("%s",s);
        build(s);
    }
    tr[0].fail=0;
    get_fail();
    dp[0][0]=1;
    for(int i=1;i<=m;i++)
        for(int j=0;j<=len;j++)
            for(int k=0;k<26;k++)
                if(!tr[tr[j].son[k]].tot) dp[i][tr[j].son[k]]=(dp[i][tr[j].son[k]]+dp[i-1][j])%mod;
    int ans=dg(26,m);
    for(int i=0;i<=len;i++)
        ans=(ans-dp[m][i]+mod)%mod;
    printf("%d",ans);
}
  • 1
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值