注意事项:
本题为"线性dp—最长上升子序列的长度"和"线性dp—最长公共子序列"的扩展题,请确保已经清晰理解代码和思路!
题目:
熊大妈的奶牛在小沐沐的熏陶下开始研究信息题目。
小沐沐先让奶牛研究了最长上升子序列,再让他们研究了最长公共子序列,现在又让他们研究最长公共上升子序列了。
小沐沐说,对于两个数列 A和 B,如果它们都包含一段位置不一定连续的数,且数值是严格递增的,那么称这一段数是两个数列的公共上升子序列,而所有的公共上升子序列中最长的就是最长公共上升子序列了。
奶牛半懂不懂,小沐沐要你来告诉奶牛什么是最长公共上升子序列。
不过,只要告诉奶牛它的长度就可以了。
数列 A和 B的长度均不超过 3000。
输入格式
第一行包含一个整数 N,表示数列 A,B的长度。
第二行包含 N个整数,表示数列 A。
第三行包含 N个整数,表示数列 B。
输出格式
输出一个整数,表示最长公共上升子序列的长度。
数据范围
1≤N≤3000,
序列中的数字均不超过 2^31−1。
输入:
4
2 2 1 3
2 1 2 3
输出:
2
#include <cmath>
#include <cstring>
#include <iostream>
#include <algorithm>
using namespace std;
const int N = 3010;
int n;
int a[N], b[N], f[N][N];
//朴素版本,O(n^3)
//完全按照dp思路来就可以,线性dp
void base() {
for (int i = 1; i<=n; i++) {
for (int j = 1; j<=n; j++) {
f[i][j] = f[i-1][j]; //先更新不包括当前a[i]的情况
if (a[i] == b[j]) { //只有当a[i]=b[j],才能计算包括a[i]的情况
// maxv代表f[a[i],b[1~j-1]]的最长公共上升子序列长度
int maxv = 1;
for (int k = 1; k < j; k++) {
if (b[k] < b[j]) maxv = max(f[i][j], f[i][k] + 1);
}
f[i][j] = max(maxv, f[i][j]); //把maxv更新给f[i][j];
}
}
}
}
// 优化版本,O(n^2)
//通过将maxv看作前缀最大值的方式,来进行优化,省去每次都需要遍历b[1~j-1]的循环
//注意,优化是通过代码层面来分析的,不要想意义层面,会想疯的(
void op() {
for (int i = 1; i<=n; i++) {
int maxv = 1;
for (int j = 1; j<=n; j++) {
f[i][j] = f[i-1][j];
if (a[i] == b[j]) f[i][j] = max(f[i][j], maxv);
//这里是优化b[k]<b[j],因a[i]=b[j], 有b[k]<a[i], 而maxv使得我们不需要重复计算b[k],直接和b[j]比即可
if (b[j] < a[i]) maxv = max(maxv, f[i][j] + 1);
}
}
}
int main()
{
cin >> n;
for (int i = 1; i<=n; i++) cin >> a[i];
for (int i = 1; i<=n; i++) cin >> b[i];
// base();
op();
//找到在前a[n]和前b[i]中的,以b[i]为结尾的最长公共上升子序列,全找一遍输出最大即可
int res = 0;
for (int i = 1; i<=n; i++) res = max(res, f[n][i]);
cout << res;
return 0;
}
思路:
经典的y式dp法
1.状态表示
f[i][j]
:
所有从a[1-i]
和b[1-j]
中以b[j]
为结尾的公共上升子序列的集合,
属性为Max。
2.状态计算
根据公共上升子序列中是否包含a[i]来进行划分集合,不重不漏。
不包含a[i]的集合:
直接从前一项转移:f[i-1][j]
。
包含a[i]的集合(通过 遍历/前缀最大值 找到):
继续划分为 a[i]和b[1到j-1],找到其中的最大值,也就是:max(f[i][1~j-1] + 1)
。
最后:优化思路
优化就是直接从代码层面上来看,
发现每次循环求得的maxv是满足a[i] > b[k]的f[i][k] + 1的前缀最大值。
因此可以直接将maxv提到第一层循环外面,减少重复计算,此时只剩下两重循环。
脑壳痛的题,有状态压缩那味了呜呜呜。
声明:
算法思路来源为y总,详细请见https://www.acwing.com/
本文仅用作学习记录和交流