BZOJ 1030 ac自动机+dp

13 篇文章 0 订阅
6 篇文章 0 订阅

题目:文本生成器
头一次做这类题目,果然智商又不够用了
我认为主要思想是将一个虚拟的串在ac自动机上进行匹配,再把匹配的状态记录下来,而不用生成所有的串一个个进行匹配,因为前缀相同的串在ac自动机上匹配的结果是相同的,只需要在相同前缀的分支后边进行26种不同的选择,每一种再分别对应到ac自动机上的一个状态;
而dp[i][j]就是记录的长度为i,能够匹配到j的不合法串
(答案就是总的可能串-不合法串嘛~)
这里再详细一点解释就是,如果有一个长度为i的串,他有26^i种不同的可能,其中正常在ac自动机上走能走到j这个节点的可能数;(这里的可能数是指不合法串的数量)
这里需要注意,长度为i不一定之前的前缀一定相同(好像肯定不同。。。)但是如果符合条件能按ac自动机的匹配规则走到j点就可以;
明确了状态之后就好做了(我太弱了这个状态搞了我很久。。。)
那么考虑从每一个已知的状态向其他状态转移(也就是向i+1的状态转移)
枚举接下来的26种可能,设选择的字符为k
首先看接下来是不是可能包含一个已知串(也就是不合法的串)
由于任何不合法的点都不会被转移到,所以当前点沿fail指针一路跳上去的所有点肯定都是合法的,但是接下来的k这个节点就不一定了,即使当前fail指针跳到的这个位置没有k这个子节点,也不能保证他继续往上跳会不会出现不符合的情况;
所以要检查包括根节点的所有fail指针能够跳到的节点,检查其如果有k这个子节点,看看这个子节点是否有end标记(标记是一个单词的结尾)
当然也可以一开始就从一个单词的结尾沿fail一路染色,好像这样会更优?
当所有检查到的点都没有end标记时,就可以模拟ac自动机的过程进行转移,看如果接下来添加的这个节点真的是k,那么沿着ac自动机会跳到哪里去;
将跳到的节点的dp+=dp[i][j]即可
当然有可能直接跳到0节点,这就代表失配的过程;
完了;
(感觉我好啰嗦。。。如果不多写点以后忘了就惨了,希望和我有同样脑回路的人能看懂吧qAq)
(还有我要提醒自己,ac自动机的板子不能再打错了!!!for循环找到fail指针不能再忘了break!!!)

#include<iostream>
#include<algorithm>
#include<string>
#include<cstring>
#include<cstdio>
#define LL long long
#define random(a,b) (a+rand()%((b)-(a)+1))
const int mod=10007;
const int maxn=6005;
int que[maxn];
int dp[105][maxn];
struct qwe
{
    int next[26];
    int fail;
    bool end;
}node[maxn]; 
struct asd{
    int cnt;
    asd()
    {
        cnt=0;
        node[cnt].fail=-1;
    }
    void insert(char *s)
    {
        int len=strlen(s);
        int now=0;
        for(int i=0;i<len;i++)
        {
            int nv=s[i]-'A';
            if(!node[now].next[nv])node[now].next[nv]=++cnt;
            now=node[now].next[nv];
        }
        node[now].end=true;
    }
    void bfs()
    {
        int tail=0,head=0;
        que[tail++]=0;
        while(head<tail)
        {
            int now=que[head++];
            for(int i=0;i<26;i++)
                if(node[now].next[i])
                {
                    que[tail++]=node[now].next[i];
                    for(int k=node[now].fail;k!=-1;k=node[k].fail)
                        if(node[k].next[i])
                        {
                            node[node[now].next[i]].fail=node[k].next[i];
                            break;//break!!!
                        }
                }
        }   
    } 
}acto; 
char s[105];
int n,m;
int main()
{
    scanf("%d%d",&n,&m);
    for(int i=1;i<=n;i++){
        scanf("%s",s);
        acto.insert(s);
    }
    acto.bfs();
    dp[0][0]=1;
    for(int i=0;i<m;i++)
    {
        for(int k=0;k<=acto.cnt;k++)
            if(dp[i][k])
                for(int c=0;c<26;c++)
                {
                    int flag=0;
                    for(int tn=k;tn!=-1;tn=node[tn].fail)//检查沿fail指针的每一个子节点
                        if(node[node[tn].next[c]].end){
                            flag=1;
                            break;  
                        }
                    if(flag)continue;
                    int now;
                    for(now=k;now;now=node[now].fail)if(node[now].next[c])break;
                    now=node[now].next[c];//如果for循环到0节点有下一个节点就转移,没有就直接到0节点
                    dp[i+1][now]=(dp[i+1][now]+dp[i][k])%mod;
                }
    }
    int ans=1;
    for(int i=1;i<=m;i++)
        ans=(ans*26)%mod; 
    for(int i=0;i<=acto.cnt;i++)
        ans=(ans-dp[m][i]+mod)%mod; 
    printf("%d\n",ans);
    return 0;
} 
  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值