浅谈LCS

danxiangfa

已于 2023-01-15 11:15:45 修改

阅读量96

点赞数

文章标签：算法 c++

于 2023-01-15 11:14:48 首次发布

本文链接：https://blog.csdn.net/danxiangfa/article/details/128692725

版权

LCS （最长公共子序列）

$O(n^2)$ dp。

求 $a 1$ 和 $a 2$ 的最长公共子序列， $f [i] [j]$ 表示 $a 1 前 i 项所组成的序列$ 与 $a 2 前 j 项所组成的序列$ 的 LCS长度。

如何转移？显然，如果 $a 1 [i]$ 与 $a 2 [j]$ 相同，那么他们就一定是LCS的一部分。即使不相同，那么这个 $f [i] [j]$ 也可以从 $f [i - 1] [j]$ 或 $f [i] [j - 1]$ 继承而来。

最终转移方程：

$\begin{cases} max(f[i-1][j-1]+1,f[i-1][j],f[i][j-1])& a1[i]=a2[j] \\ max(f[i-1][j],f[i][j-1])& a1[i] \neq a2[j] \end{cases}$

可以转化成 LIS

关于为什么可以转化成LIS问题，这里提供一个解释。
A:3 2 1 4 5
B:1 2 3 4 5
我们不妨给它们重新标个号：把3标成a,把2标成b，把1标成c……于是变成：
A: a b c d e
B: c b a d e
这样标号之后，LCS长度显然不会改变。但是出现了一个性质：
两个序列的子序列，一定是A的子序列。而A本身就是单调递增的。
因此这个子序列是单调递增的。
换句话说，只要这个子序列在B中单调递增，它就是A的子序列。
哪个最长呢？当然是B的LIS最长。
自此完成转化。——阮行止

模板题P1439

#include <bits/stdc++.h>
using namespace std;
int ans, n, a1[100001], a2[100001], id[100001], f[100001];
int main() {
	ios::sync_with_stdio(false), cin.tie(0), cout.tie(0);
	cin >> n;
	for(int i = 1; i <= n; i++) cin >> a1[i], id[a1[i]] = i;
	for(int j = 1; j <= n; j++) cin >> a2[j], a2[j] = id[a2[j]];
	ans = 1, f[1] = a2[1];
	for(int i = 2; i <= n; i++)
		if(a2[i] > f[ans])
			f[++ans] = a2[i];
		else *lower_bound(f + 1, f + ans + 1, a2[i]) = a2[i];
	cout << ans << endl;
	return 0;
}