bzoj 5454: Subsequence

题意

给出两个由小于等于k的正整数构成的数列A与B,长度为n与m
现在需要一个由小于等于k的正整数构成的数列C,使得C既不是A的子序列,也不是B
的子序列,请求出C的最小长度

题解

继续开发第45版。。

有两个做法,都说一说把

看到序列自然想到序列自动机
建立自动机
然后bfs一下
f [ x ] [ y ] f[x][y] f[x][y]表示第一个串到 x x x,第二个串到 y y y的最短步数
然后 O ( k ) O(k) O(k)转移?
最后 f [ 0 ] [ 0 ] f[0][0] f[0][0]就是答案
看起来复杂度是 O ( n k 2 ) O(nk^2) O(nk2)
但是实际上, f f f有值的不多,所以跑得飞快

PP说序列自动机的题都可以用DP来解决
那么我们就来考虑DP
考虑我给你一个串
我们怎么判断答案
自然是DP
一位一位插入
f i , j f_{i,j} fi,j表示第一个串匹配到 i i i,第二个匹配到 j j j是否可能
然后如果最后有任何一个位置OK就是存在了
考虑用DP模拟这个过程
f i , j f_{i,j} fi,j表示第一个串匹配到 i i i,第二个匹配到 j j j的最小步数是什么
然后每一次也是 O ( k ) O(k) O(k)转移
转移到什么呢?
当然是转移到 i i i后面的第一个 k k k j j j后面的第一个 k k k
如果不存在这个k,就可以输出答案了
看起来也是 O ( n m k ) O(nmk) O(nmk)
但是实际上跑得飞快

或许可以证明是更优的复杂度?
但是我不会证
当然我也没想证
我认为是没有办法卡掉的
欢迎hack(大雾

CODE(序列自动机):

 
#include<cstdio>
#include<algorithm>
#include<iostream>
#include<cstring>
#include<queue>
using namespace std;
const int N=4005;
int n,m,k,ans;
struct qq
{
    short int son[N][N];
    short int par[N];
    short int last[N];
    short int tot;
    void Ins (int x)
    {
        int np=++tot;
        for (int i=1;i<=k;i++)
        for (int j=last[i];j!=0&&son[j][x]==0;j=par[j])
        son[j][x]=np;
        par[np]=last[x];    last[x]=np;
    }
}a,b;
struct qt
{
    short int rt1,rt2,dep;
    qt () {};
    qt (short int _rt1,short int _rt2,short int _dep)   {rt1=_rt1;rt2=_rt2;dep=_dep;}
};
bool vis[N][N];
void bfs ()
{
    queue<qt> q;
    q.push(qt(1,1,0));
    while (!q.empty())
    {
        qt x=q.front();q.pop();
        for (int u=1;u<=k;u++)
        {
            int y=a.son[x.rt1][u],z=b.son[x.rt2][u];
            if (vis[y][z]) continue;
            vis[y][z]=true;
            if (y==0&&z==0)
            {
                printf("%d\n",x.dep+1);
                return ;
            }
            q.push(qt(y,z,x.dep+1));
        }
    }
}
int main()
{
    memset(vis,false,sizeof(vis));
    scanf("%d%d%d",&n,&m,&k);
    ans=n+m+1;
    for (int u=1;u<=k;u++)   b.last[u]=a.last[u]=1;
    for (int u=1;u<=n;u++)
    {
        int x;
        scanf("%d",&x);
        a.Ins(x);
    }
    for (int u=1;u<=m;u++)
    {
        int x;
        scanf("%d",&x);
        b.Ins(x);
    }
    bfs();
    return 0;
}

  • 1
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

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

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

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值