CH5101 LCIS 简单dp

题目及思路来源：李煜东 《算法竞赛进阶指南》

题目链接：LCIS

如何优化dp是一个值得思考的问题  dp与搜索的区别在于 dp避免了大多的重复运算 以此来达到高效 （记忆化搜索）  不过有时候写的dp思路会比较直接 会算很多重复的  这时候要适当优化  比如这题  一开始 写一个n^3的dp dp[i][j]代表  a1到ai 与b1到bj 以bj结尾的最长公共上升子序列的长度

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
int dp[3002][3002],a[3002],b[3002]; 
int main(){
	int n;
	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++){
		for(int j = 1; j <= n; j++){
			if(a[i]==b[j]){
				for(int k = 0; k < j; k++)
				if(b[k]<a[i]) 
				dp[i][j]=max(dp[i][j],dp[i-1][k]+1);
			}else dp[i][j]=dp[i-1][j];
		}
	}
	int ans = 0;
	for(int i = 1; i <= n; i++) ans=max(ans,dp[n][i]);
	printf("%d\n",ans); 
}

这样复杂度显然不行（不过oj上还是过了  可能是数据不强） 我们发现第三重循环似乎做了很多没必要的工作  对于每个j 我们都从前往后去寻找它能转移的最大dp[i-1][k]  （专业的说是叫啥  决策集合） 而对于 j+1来说  决策集合几乎没变  也是 dp[i-1][k] k从0到j-1或0到j（我们需要判断）   那我们可以想到 我们其实只要维护最大的dp[i-1][k]就行  其他的都不需要  那我们设置一个val 来动态更新就行了  这样 直接就把复杂度降到n^2了 是质跃

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
int dp[3002][3002],a[3002],b[3002]; 
int main(){
	int n;
	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++){
		int val = 0;
		for(int j = 1; j <= n; j++){
			if(a[i]==b[j])
				dp[i][j]=val+1;
			else dp[i][j]=dp[i-1][j];
			if(b[j]<a[i]) val=max(val,dp[i-1][j]);
		}
	}
	int ans = 0;
	for(int i = 1; i <= n; i++) ans=max(ans,dp[n][i]);
	printf("%d\n",ans); 
}