[Luogu P1439] LCS升级版

badge：

普通O(n^2)很容易想到，但是TLE。

#include<cstdio>
using namespace std;
int a[1010],b[1010],dp[1010][1010];
int maxf(int x,int y){return x>y?x:y;}
int minf(int x,int y){return x<y?x:y;}
int main()
{
    int n,m,i,j,k;
    scanf("%d",&n);
    for(i=1;i<=n;i++)scanf("%d",&a[i]);
    for(i=1;i<=n;i++)scanf("%d",&b[i]);
    for(i=1;i<=n;i++)
    {
        for(j=1;j<=n;j++)
        {
            dp[i][j]=maxf(dp[i][j],maxf(dp[i-1][j],dp[i][j-1]));
            if(a[i]==b[j])
            dp[i][j]=maxf(dp[i][j],dp[i-1][j-1]+1);
        }
    }
    printf("%d\n",dp[n][n]);
    return 0;
}

我们思考这题的关键点：两个序列都是1~n的排列。

排列？这意味着两个序列里数都一样，只是位置不同。

有点类似于友好城市，只需要求一堆不相交的配对即可。

画出图来，把右边的序列按照左边每个数所在的位置重新填好，求一个最长上升子序列即可。

即：记录a序列里每个数填的位置f[a[i]]=i，然后把b[i]赋成f[b[i]]。在b序列里O(nlogn)求一下LIS即可。

#include<cstdio>
using namespace std;
int a[1000010],b[1000010];
int c[1000010],f[1000010],tot;
int pos(int l,int r,int val)
{
    int mid,ans;
    while(l<=r)
    {
        mid=(l+r)>>1;
        if(c[mid]>val)r=mid-1,ans=mid;
        else l=mid+1;
    }
    return ans;
}
int main()
{
    int n,m,i,j,k;
    scanf("%d",&n);
    for(i=1;i<=n;i++)scanf("%d",&a[i]);
    for(i=1;i<=n;i++)scanf("%d",&b[i]);
    for(i=1;i<=n;i++)f[a[i]]=i;
    for(i=1;i<=n;i++)b[i]=f[b[i]];
    //for(i=1;i<=n;i++)printf("%d ",f[i]);printf("\n");
    //for(i=1;i<=n;i++)printf("%d ",b[i]);printf("\n");
    c[++tot]=b[1];
    for(i=2;i<=n;i++)
    {
        if(b[i]>c[tot])c[++tot]=b[i];
        else
        {
            k=pos(1,tot,b[i]);
            //printf("found pos %d\n",k);
            c[k]=b[i];
        }
    }
    printf("%d\n",tot);
    return 0;
}
/*
5 
3 2 1 4 5
1 2 3 4 5
*/