P1439 【模板】最长公共子序列

书接上文~我来补坑~本章讲讲上文中未能AC的一道题

之前使用的的代码，完美超时，如下

#include<bits/stdc++.h>
using namespace std;
int dp[1001][1001],a1[2001],a2[2001],n,m;
int main()
{
   cin>>n>>m;
   for(int i=1;i<=n;i++)scanf("%d",&a1[i]);
   for(int i=1;i<=n;i++)scanf("%d",&a2[i]);
   for(int i=1;i<=n;i++)
    for(int j=1;j<=n;j++)
     {
     	if(a1[i]==a2[j]) dp[i][j]=max(dp[i][j],dp[i-1][j-1]+1); 
     	else dp[i][j]=max(dp[i-1][j],dp[i][j-1]);
     }
   cout<<dp[n][n];
}

耗时，写不了一点

是时候想想怎么优化了...

因为两个序列都是1 到 n的全排列，那么两个序列元素互异且相同，也就是说只是位置不同罢了，那么我们通过一个map数组将A序列的数字在B序列中的位置表示出来——

因为最长公共子序列是按位向后比对的，所以a序列每个元素在b序列中的位置如果递增，就说明b中的这个数在a中的这个数整体位置偏后，可以考虑纳入LCS——那么就可以转变成nlogn用来记录新的位置的map数组中的LIS。

#include<iostream>
#include<cstdio>
using namespace std;
int a[100001],b[100001],map[100001],f[100001];
int main()
{
	int n;
	cin>>n;
	for(int i=1;i<=n;i++){
        scanf("%d",&a[i]);
        map[a[i]]=i;
    }
	for(int i=1;i<=n;i++){
        scanf("%d",&b[i]);
        f[i]=0x7fffffff;
    }
	int len=0;
	f[0]=0;
	for(int i=1;i<=n;i++)
	{
		int l=0,r=len,mid;
		if(map[b[i]]>f[len])f[++len]=map[b[i]];
		else 
		{
		while(l<r)
		{	
		    mid=(l+r)/2;
		    if(f[mid]>map[b[i]])r=mid;
			else l=mid+1; 
		}
		f[l]=min(map[b[i]],f[l]);
     	}
    }
    cout<<len;
    return 0
}