最长公共上升子序列:题解

1前言

本文是acwing提高课动态规划专题最长公共上升子序列的题解
会先介绍最长上升子序列和最长公共子序列，因为优化和本题无关，请勿提交

2最长上升子序列

作者认为最简单的dp
思考步骤如下
（1）可以设置状态为dp[i]为含有*a[i]项的最长上升子序列
（2）状态转移，可以发现，含有a[i]项，dp[i]可初始化为1，并且对于所有j(j<i且a[j]<a[i])，a[i]可以接在a[j]后面
由此得状态转移方程 dp[i] = max(dp[j]+1,dp[i]) (j<i&&a[j]<a[i])
时间复杂度O(nn)

3最长公共子序列

这个是典型的二维dp
状态表示为dp[i][j]为a串到i，b串到j的子串的最长公共子序列
如果a[i]=b[j]
则dp[i][j] = max(dp[i][j],dp[i-1][j-1]+1)
否则dp[i][j] = max(dp[i-1][j],dp[i][j-1])
时间复杂度也是O(nn)*

4最长公共上升子序列

首先考虑暴力

可以看到，由于公共性质，上升的话只需要考虑一个b数组就行了
因为要考虑上升，要比大小，这个子序列最后一位要是定值，那么，dp[i][j]状态要含有b[i]
因为这个性质，要含有b跟上一个有点不同，如果a[i]不等于b[j]，不能改j
则有dp[i][j] = dp[i-1][j]，也就是说直接复制上一个i
那如果相等呢，就得扫一遍dp[i]k,因为公共的性质，不用考虑i了，只要考虑所有令b[k]<b[j]中dp[i][k]的最大值
则有dp[i][j] = max(dp[i][k])+1 不要忘了+1,这代表着j
一看时间复杂度，O(nnn)，上一个三次方的还是floyd，这个不超时就怪了

优化

前两维循环不可省，但是那个扫一遍的多余了
我们可以开一个max来记录max(dp[i][k])
如果a[i]不等于b[j]，未来可能会有一个b[x]和a[i]相等，所以如果b[j]<a[i]，更新max
因为i固定，k在变，甚至不用开数组
这下变成了O(nn)*，不会超时了
附代码（c++)

#include<bits/stdc++.h>
using namespace std;
int a[114514];
int b[114514];
int dp[3333][3333];
int n,ans = 0;
int main(){
	cin>>n;
	for(int i = 1;i<=n;i++){
		cin>>a[i];
	}
	for(int i = 1;i<=n;i++){
		cin>>b[i];
	}//输入 
	for(int i = 1;i<=n;i++){
		int maxs = 1;//max(dp[i][k]) 
		for(int j = 1;j<=n;j++){
			dp[i][j] = dp[i-1][j];//不相等，不用if了，相等的话直接覆盖掉了 
			if(a[i]==b[j]){//相等 
				dp[i][j] = max(dp[i][j],maxs);
			}else if(a[i]>b[j]){//特判一下小于a[i],因为公共所以和a[i]比较即可 
				maxs = max(maxs,dp[i][j]+1);
			}
		}
	}
	for(int i = 1;i<=n;i++){//因为i大的dp[i][j]值永远最大，遍历j值即可 
		ans = max(ans,dp[n][i]);
	}
	cout<<ans;
	return 0;
}

# 5后记

暴力的做法作者并未在网上发现，如有错误请各位神犇指点
**森林古猿出品，必属精品，请认准CSDN森林古猿1**