题目描述:
题目链接:
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;
}