最长公共上升子序列(AcWing, 算法竞赛进阶指南)

题目描述: 

 题目链接:

https://www.acwing.com/problem/content/description/274/

 思路: 非常经典的动态规划问题，值得思考。

首先不难想到以下dp数组的定义:

f[i][j] 表示a[1~i] 和b[1~j] 且以b[j]结尾的最长公共上升子序列的长度。

考虑f[i][j]的更新，根据a[i] 能否和b[j] 相匹配，可得:

1. 若a[i] != b[j] , f[i][j]=max(f[i][j],f[i-1][j])

2. 若a[i]==b[j] , 则需要在序列b中找到上一个匹配位置，和a[1~(i-1)] 相匹配

i.e. f[i][j]=max(f[i][j],max_pre)， 其中max_pre=max(f[i-1][k]), 

DP的优化:

 不难看出，上面的解法是正确的，但是时间复杂度为O(n3), 需要对dp进行优化。

yxc老师说过,dp 优化一般与dp的思路无关，只是对dp更新式子做等价变形。

我们注意到: 在max_pre的计算中，实际上b[j] 一定是等于a[i] 的，而a[i] 的值与j 的取值无关， 也就是说，考虑以下循环:

for(int k=1;k<j;k++){
    if(b[k]<const)
        // 更新max_pre
}

注意到，上面的东西实际上就是一个“前缀最大值”， i.e. 考虑1~(j-1)中那些满足b[k]<a[i] 的位置k的f[i-1][k] 的最大值。

那我们就可以采用“动态更新”的策略，边求值，边更新. i.e. 先利用前缀最大值max_pre去更新f[i][j], 若这个位置j满足b[j] < a[i], 则用f[i-1][j] 去更新max_pre 。

通过这样的做法，就可以去掉一层循环，从而使得整个算法的复杂度变为O(n2)

代码：

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
const int N=3010;
int a[N],b[N];
int n;
int f[N][N];

int main(void){
	
	scanf("%d",&n);
	for(int i=1;i<=n;i++){
		scanf("%d",&a[i]);
	}
	for(int i=1;i<=n;i++){
		scanf("%d",&b[i]);
	}
	for(int i=1;i<=n;i++) f[i][0]=f[0][i]=0;
	for(int i=1;i<=n;i++){
		int pre=0;
		for(int j=1;j<=n;j++){
			if(a[i]!=b[j]) f[i][j]=max(f[i][j],f[i-1][j]);
			else{
				f[i][j]=max(f[i][j],pre+1);
			}
			if(b[j]<a[i]) pre=max(pre,f[i-1][j]);
		}
	}
	int ans=0;
	for(int i=1;i<=n;i++) ans=max(ans,f[n][i]);
	printf("%d",ans);
	return 0;
}