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

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

题目传送门

分析：

如果大家看它是道最长公共子序列的模板题，就直接打个$O(n^2)$的模板的最长公共子序列上去，你会发现你的程序愉快地RE了，只拿到了50分

RE瞬间

50分的CODE：

#include<iostream>
#include<cstdio>
#include<cstring>
using namespace std;
int d[10001][10001],n,a[10001],b[10001];
int main()
{
	cin>>n;
	for (int i=0;i<n;i++)
		scanf("%d",&a[i]);
	for (int i=0;i<n;i++)
		scanf("%d",&b[i]);
	for (int i=0;i<n;i++)
	{
		for (int j=0;j<n;j++)
		{
			if (a[i]==b[j])
			{
				d[i+1][j+1]=d[i][j]+1;
			}
			else d[i+1][j+1]=max(d[i][j+1],d[i+1][j]);
		}
	}
	cout<<d[n][n];
	return 0;
}

拜托这可是一道绿题会让你那么容易过的吗

我们看一下题目的数据范围

对于 $100\%$ 的数据， $n\le 10^5$

如果d数组开$10^5\ast 10^5$的话妥妥的爆内存

我们尝试换一种做法

题目输入：

5
3 2 1 4 5
1 2 3 4 5

我们将A序列的内容和下标联系起来

3对应下标1，2对应下标2，1对应下标3，4对应下标4，5对应下标5

由此我们可以对照这个顺序将B序列转换

转换成 3，2，1，4，5

这样转换之后，最长公共子序列的长度是不会改变。但是出现了一个性质：

两个序列的子序列，一定是A的子序列。而A本身就是单调递增的。

因此这个子序列是单调递增的。

换句话说，只要这个子序列在B中单调递增，它就是A的子序列

也就把这个求最长公共子序列的问题转换成了一个求最长不下降子序列的问题

但是呢有些人可能要问了啊，最长不下降子序列的问题不还是$O(n^2)$的时间复杂度吗，只能将空间占用减少吗

其实我们还可以用单调队列对其进行优化

放一张图片让大家理解一下

这里的a数组表示的是输入的数

d数组则是单调不降的队列

每次$a_i$就只要和d数组的最后一位数作比较就行了

如果队列的最后一位大，我们就可以直接放，再把队列的长度加一就行了

如果放不下的话，我们就只要在d数组里二分查找一个比当前数大但又是最小的一个数，只要替换掉就行了

这样的话时间复杂度就从$O(n^2)$降到了$O(nlogn)$

是不是very easy

100分的CODE:

#include<iostream>
#include<cstdio>
#include<cstring> 
using namespace std;
const int MAXN=100001;
int n,k;
int a[MAXN];
int b[MAXN];
int c[MAXN];
int d[MAXN];
void input()
{
	scanf("%d",&n);
	for (int i=1;i<=n;i++)
	{
		scanf("%d",&a[i]);
		c[a[i]]=i;
		a[i]=0x7ffffff;
	}
	for (int i=1;i<=n;i++)
	{
		scanf("%d",&b[i]);
		b[i]=c[b[i]];//转换成最长不下降子序列的问题 
	}
}
int search(int l,int r,int x)//二分查找 
{
	while (l<=r)
	{
		int mid=(l+r)/2;
		if (x<=a[mid])
			r=mid-1;
		else
			l=mid+1;
	}
	return l;
}
void work()
{
	k=0;
	for (int i=1;i<=n;i++)
	{
		if (b[i]>a[k])
			a[++k]=b[i];//如果比队列的最后一位大，直接放，再把队列的长度加一
		else {
			a[search(1,k,b[i])]=b[i];//否则二分查找一个比它大又是最小的一个，将其替换掉 
		}
	}
}
void output()
{
	cout<<k<<endl;//输出队列的长度，也就是最长不下降子序列的长度 
}
int main()
{
	input();
	work();
	output();
	return 0;
}

AC瞬间

注意：队列里的数并不就是序列的最长不下降子序列，也就是说这个优化只能求出数量并不能够求出序列

总结：

dp博大精深不要以为学到了一点皮毛就觉地自己学完了dp，毕竟dp可是在每一本编程书里都是独占一个单元的呢

另外学dp挺废脑子的