【刷题】动态规划——线性DP：最长公共子序列【LCS、相同元素转化成最长上升nlogn】

题目链接

f[i][j]表示所有在第一个序列的前i个字母中出现，并在第二个序列前j个字母中出现的公共子序列长度的 最大值。

考虑a[i]和b[j]是否包含在子序列当中。
1、不含a[i]、不含b[j]：f[i][j]=f[i-1][j-1]
2、含a[i]、含b[j]：f[i][j]=f[i-1][j-1] + 1 且 a[i]=b[j]
3、不含a[i]、含b[j]：f[i][j]=f[i-1][j]
4、含a[i]、不含b[j]：f[i][j]=f[i][j-1]

于是可以得到
当a[i]=b[j]时f[i][j] = max(f[i-1][j], f[i][j-1], f[i-1][j-1]+1)
当a[i]!=b[j]时f[i][j] = max(f[i-1][j], f[i][j-1], f[i-1][j-1])

---

注意f[i-1][j]是前j个字母出现在b中的最优解，不只包含了含b[j]的最优解，还包含了不含b[j]的最优解（也就是f[i-1][j-1]）。所以可用f[i-1][j]表示所有不含a[i]、含b[j]的情况。

同理，可用f[i][j-1]表示所有含a[i]、不含b[j]的情况（当然还包括了不含a[i]的情况f[i-1][j-1]）。

---

综上，由于f[i][j]=f[i-1][j] 和f[i][j]=f[i][j-1] 已经包括了f[i][j]=f[i-1][j-1]
所以：
当a[i]=b[j]时f[i][j] = max(f[i-1][j], f[i][j-1], f[i-1][j-1]+1)
当a[i]!=b[j]时f[i][j] = max(f[i-1][j], f[i][j-1])

这样时间复杂度是O(n^2)

#include <iostream>
#include <algorithm>
using namespace std;
const int N = 1005;
int n, f[N][N];
int a[N], b[N];
int main() {
	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++) {
			f[i][j] = max(f[i][j-1], f[i-1][j]);
			if (a[i] == b[j]) {
				f[i][j] = max(f[i][j], f[i-1][j-1] + 1);
			}
		}
	}
	printf("%d\n", f[n][n]);
	return 0;
}

---

下面介绍O(nlogn)的做法

注意到这两个序列长度相同，并且数字都是1~n，以题目为例。
5
3 2 1 4 5
1 2 3 4 5
假设其中一个序列是从小到大顺序（这边的序列B），因为序列B是1~n的上升序列，任意上升子序列都包含在B中，所以只要在A中找最长上升子序列就是AB的最长公共子序列，这边是1 4 5。最长上升子序列时间复杂度是O(nlogn)的。

所以对于俩序列，只要将一个序列A通过映射P处理成1~n的上升序列，序列B也通过这个映射P，对映射后的结果查上升子序列。
同样以样例为例，序列A映射后是1 2 3 4 5

原来数字	映射后数字
3	1
2	2
1	3
4	4
5	5
序列B映射后是3 2 1 4 5
数组p[i]表示把i映射成p[i]

#include <iostream>
#include <algorithm>
using namespace std;
const int N = 100005;
int n, f[N], lf;
int a[N], b[N], p[N];
int main() {
	scanf("%d", &n);
	for (int i = 1; i <= n; i++) scanf("%d", &a[i]), p[a[i]] = i;
	for (int i = 1; i <= n; i++) scanf("%d", &b[i]);
	for (int i = 1; i <= n; i++) {
		b[i] = p[b[i]];
	}
	
	lf = 0;
	for (int i = 1; i <= n ; i++) {
		int k = lower_bound(f + 1, f + lf + 1, b[i]) - f;
		if (k > lf) {
			f[++lf] = b[i];
		}
		else {
			f[k] = b[i];
		}
	}
	printf("%d\n", lf);
	return 0;
}