272. 最长公共上升子序列
题意
熊大妈的奶牛在小沐沐的熏陶下开始研究信息题目。
小沐沐先让奶牛研究了最长上升子序列,再让他们研究了最长公共子序列,现在又让他们研究最长公共上升子序列了。
小沐沐说,对于两个数列 A 和 B,如果它们都包含一段位置不一定连续的数,且数值是严格递增的,那么称这一段数是两个数列的公共上升子序列,而所有的公共上升子序列中最长的就是最长公共上升子序列了。
奶牛半懂不懂,小沐沐要你来告诉奶牛什么是最长公共上升子序列。
不过,只要告诉奶牛它的长度就可以了。
数列 A 和 B 的长度均不超过 3000。
思路
状态表示: f [ i ] [ j ] f[i][j] f[i][j] 表示 第一个数组(a)前 i i i 个数字 和 第二个数组(b)前 j j j 个数字组成的 并且以 $ b[j]$ 结尾的 最长公共上升子序列的最大长度
把 f [ i ] [ j ] f[i][j] f[i][j] 分成两种状态 : { 包 含 a [ i ] 不 包 含 a [ i ] \begin{cases}包含a[i] \\不包含a[i] \end{cases} {包含a[i]不包含a[i]
对于不包含 a [ i ] a[i] a[i] 即第一个数组 前 i − 1 i - 1 i−1 数字 和 第二个数组前 j j j 个数字组成的 并且以 b [ j ] b[j] b[j] 结尾的 最长公共上升子序列的最大长度 可以用 f [ i − 1 ] [ j ] f[i - 1][j] f[i−1][j] 直接表示
f[i][j] = f[i - 1][j];
对于包含 a [ i ] a[i] a[i] 的 首先要满足公共子序列的性质 有 a [ i ] = = b [ j ] a[i] == b[j] a[i]==b[j] 但其无法直接用 f [ x ] [ y ] f[x][y] f[x][y] 这种形式表示 所以我们根据 b b b 数组倒数第二个数字将上述集合进行划分:
即分成 以 a a a 数组前 i i i 个数字和 b b b 数组前 k k k ( k ∈ [ 1 , j − 1 ] ) (k \in [1,j-1]) (k∈[1,j−1])个数字组成 并且以 b [ k ] b[k] b[k] 结尾的 最长公共子序列
但是因为是最长公共上升子序列 所以要满足 b [ k ] < b [ j ] b[k] < b[j] b[k]<b[j]
所以可以得到以下的代码
#include<bits/stdc++.h>
#define INF 0x3f3f3f3f
#define mod 998244353
#define endl '\n'
using namespace std;
typedef long long LL;
typedef pair<int, int>PII;
const int N = 3010;
int n;
int a[N], b[N];
int f[N][N];
// f[i][j] 表示 第一个数组前 i 个数字 和 第二个数组前 j 个数字 并且以b[j] 结尾的
// 最长公共上升子序列长度的最大值
int main() {
cin >> 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) {
f[i][j] = f[i - 1][j]; // 不包含 a[i]
if (a[i] == b[j]) { // 包含 a[i] 的条件为 a[i] == b[j]
f[i][j] = max(f[i][j], 1); // 至少为1 因为当前已经有 a[i] == b[j]
for (int k = 1; k < j; ++k) { // 从前 j - 1 个中找最大值
if(b[k] < b[j]) //因为是上升子序列 所以要保持严格递增关系
f[i][j] = max(f[i][j], f[i][k] + 1);
}
}
}
}
int res = -1;
for (int i = 1; i <= n; ++i) {
res = max(res, f[n][i]);
}
cout << res << endl;
return 0;
}
但是这种算法是 O ( n 3 ) O(n^3) O(n3) 的 显然不合题意
所以我们对原本的代码进行优化
注意到:
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] 的条件为 a[i] == b[j]
f[i][j] = max(f[i][j], 1); // 至少为1 因为当前已经有 a[i] == b[j]
for (int k = 1; k < j; ++k) { // 从前 j - 1 个中找最大值
if(b[k] < b[j]) //因为是上升子序列 所以要保持严格递增关系
f[i][j] = max(f[i][j], f[i][k] + 1);
}
}
}
}
在这个循环中 当 a [ i ] = = b [ j ] a[i] == b[j] a[i]==b[j] 且 b [ k ] < b [ j ] b[k] < b[j] b[k]<b[j]时 循环找到 m a x ( f [ i ] [ k ] + 1 ) max(f[i][k] + 1) max(f[i][k]+1) 因为 a [ i ] = = b [ j ] a[i] == b[j] a[i]==b[j] 所以 b [ k ] < a [ i ] b[k] < a[i] b[k]<a[i]
那么这个循环就跟 j j j 没有任何关系了 可以放在 f o r j for\ \ \ j for j 循环的外面
for (int i = 1; i <= n; ++i) {
int maxn = 1;
for (int j = 1; j <= n; ++j) {
f[i][j] = f[i - 1][j]; // 不包含 a[i]
if (a[i] == b[j]) f[i][j] = max(f[i][j], maxn); //包含 a[i]
// 将下方的判断放在 f[i][j] 更新之后 可以保证f[i][j] 更新时 用到的 maxn 是前 j - 1 个数的
if (b[j] < a[i]) maxn = max(maxn, f[i][j] + 1);
}
}
完整代码
#include<bits/stdc++.h>
#define INF 0x3f3f3f3f
#define mod 998244353
#define endl '\n'
using namespace std;
typedef long long LL;
typedef pair<int, int>PII;
const int N = 3010;
int n;
int a[N], b[N];
int f[N][N];
// f[i][j] 表示 第一个数组前 i 个数字 和 第二个数组前 j 个数字 并且以b[j] 结尾的
// 最长公共上升子序列长度的最大值
int main() {
cin >> 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 maxn = 1;
for (int j = 1; j <= n; ++j) {
f[i][j] = f[i - 1][j]; // 不包含 a[i]
if (a[i] == b[j]) f[i][j] = max(f[i][j], maxn);
if (b[j] < a[i]) maxn = max(maxn, f[i][j] + 1);
}
}
int res = 0;
for (int i = 1; i <= n; ++i) {
res = max(res, f[n][i]);
}
cout << res << endl;
return 0;
}