[JZOJ5424]【NOIP2017提高A组集训10.25】凤凰院凶真

Description

给出两个序列A，B
求出这两个序列的LCIS，即最长公共上升子序列
输出任意一组解

Solution

先离散化

$N^3$的做法容易想到，$F_{i,j,k}$表示A到了i，B到了j，当前末尾为j，直接转移。

考虑优化
可以用二维表示状态
$F_{i,j}$表示A到了i，B到了j，且分别以i,j结尾，那么末尾最大值这一维可以省略

先枚举i，然后在枚举j时我们发现，假设$f_{i,j}$要从末尾为p的某个$f_{i',j'}$转移而来
那么$i'$一定是最靠近$i$的$a[i']=p$的位置，即i’应该在小于i的同时最大

这个可以预处理。

考虑$j'$
既然$a[i]$是已知的，只需要在枚举$j$的同时顺便记录一下$b[j']<a[i]$且$F_{i',j'}$最大的位置即可。

复杂度$O(N^2)$

Code

#include <cstdio>
#include <algorithm>
#include <cstdlib>
#include <iostream>
#include <cstring>
#include <cmath>
#define fo(i,a,b) for(int i=a;i<=b;i++)
#define fod(i,a,b) for(int i=a;i>=b;i--)
#define N 5005
using namespace std;
int n,a[N],m,b[N],f[N][N],ls[N][N],fr[N][N][2],sx[N],vl,t[N];
struct node
{
    int v,w,p;
    friend bool operator<(node x,node y)
    {
        return x.v<y.v;
    }
}a1[2*N];
int main()
{
    cin>>n;
    fo(i,1,n) scanf("%d",&a[i]),a1[i].v=a[i],a1[i].w=i,a1[i].p=0;
    cin>>m;
    fo(i,1,m) scanf("%d",&b[i]),a1[i+n].v=b[i],a1[i+n].w=i,a1[i+n].p=1;
    sort(a1+1,a1+n+m+1);
    fo(i,1,n) a[i]=-1;
    fo(i,1,m) b[i]=-1;
    bool bx=0,by=0;
    vl=0;
    fo(i,1,n+m+1)
    {
        if(a1[i].v!=a1[i-1].v)
        {
            if(bx&&by)
            {
                int j=i-1;
                vl++;
                t[vl]=a1[i-1].v;
                while(a1[j].v==a1[i-1].v) 
                {
                    if(a1[j].p==0) a[a1[j].w]=vl;
                    else b[a1[j].w]=vl;
                    j--;
                }
            }
            bx=0,by=0;
        }
        if(a1[i].p==0)bx=1;
        else by=1;
    } 
    int n1=0,m1=0;
    fo(i,1,n) if(a[i]!=-1) a[++n1]=a[i];
    fo(i,1,m) if(b[i]!=-1) b[++m1]=b[i];
    fo(i,1,n1)
    {
        fo(j,1,vl) ls[i][j]=ls[i-1][j];
        ls[i][a[i]]=i;
    }
    int ans=0,sl=0,sr=0;
    fo(i,1,n1)  
    {   
        int mx=0,ml=0,mr=0;
        fo(j,1,m1)
        {
            if(a[i]==b[j]) 
            {
                f[i][j]=mx+1,fr[i][j][0]=ml,fr[i][j][1]=mr;
            }
            if(b[j]<a[i]&&f[ls[i-1][b[j]]][j]>mx) mx=f[ls[i-1][b[j]]][j],ml=ls[i-1][b[j]],mr=j;
            if(f[i][j]>ans) ans=f[i][j],sl=i,sr=j;
        }
    }
    printf("%d\n",ans);
    fod(i,ans,1) 
    {
        sx[i]=t[a[sl]];
        int p=sl;
        sl=fr[p][sr][0];
        sr=fr[p][sr][1];
    }
    fo(i,1,ans) printf("%d ",sx[i]);
}