合工大 程序设计艺术 实验四 最长公共子序列

题目

求最长公共子序列
X 的一个子序列是相应于 X 下标序列{1, 2, …, m}的一个子序列，求解两个序列的所有 子序列中长度最大的，例如输入：pear, peach 输出：pea。

题意梳理

这里要求的最长公共子序列与字符串匹配不同，这里的字符可以不相连，只要相对顺序是一样的即可。
例如：hfutxcxq与hdtxqb的最长公共子序列就是htxq。
最长公共子序列的个数可能不是唯一的，但是长度一定是唯一的。

解题思路

暴力枚举是我们最容易想到的方法，找出两个字符串的所有子序列，然后比较出二者最长的子序列即可，但是这个时间复杂度太高了。
因此我这里采用动态规划法。
我用Ax表示序列A的连续前x项构成的子序列，即Ax= a1,a2,……ax, By= b1,b2,……by, 用dp(x, y)表示它们的最长公共子序列长度，那原问题等价于求dp(m,n)。
若A的第i个字符和B的第j个字符相等时，dp[i][j] = dp[i-1][j-1] + 1;
若不相等时，dp[i][j] = max{dp[i-1][j],dp[i][j-1]};
当都为空时dp[i][j]=0;
可得状态转移方程：
//状态转移方程
for (int i = 1; i <= len1; i++) {
for (int j = 1; j <= len2; j++) {
if (s1[i - 1] == s2[j - 1]) {//等于斜上方的值+1
dp[i][j] = dp[i - 1][j - 1] + 1;
}
else {//等于上方或者左方的最大值
dp[i][j] = max(dp[i - 1][j], dp[i][j - 1]);
}
}
}
到这里只是获得了长度，我们还需要得到最长公共子序列。
由这张图我们可以看出，灰色的痕迹就是我们要求的序列。
递归输出公共子序列即可。递归结束条件为i0||j0;
运行结果：

源代码

#include <iostream>
#include <algorithm>
using namespace std;
#define MAXN 1010
char s1[MAXN], s2[MAXN];//存储1和2字符串
int dp[MAXN][MAXN];//状态数组
//找公共子序列
void LCS(int i, int j) {
	if (i == 0 || j == 0) {
		return;
	}
	if (s1[i - 1] == s2[j - 1]) {
		LCS(i - 1, j - 1);//斜右上方递归
		cout << s1[i - 1];
	}
	else if (dp[i][j - 1] > dp[i - 1][j]) {//左边大
		LCS(i, j - 1);
	}
	else {
		LCS(i - 1, j);
	}
}
int main() {
	//读入字符串
	cin >> s1;
	cin >> s2;
	//获取两个字符串长度
	int len1, len2;
	len1 = strlen(s1);
	len2 = strlen(s2);
	//状态转移方程
	for (int i = 1; i <= len1; i++) {
		for (int j = 1; j <= len2; j++) {
			if (s1[i - 1] == s2[j - 1]) {//等于斜上方的值+1
				dp[i][j] = dp[i - 1][j - 1] + 1;
			}
			else {//等于上方或者左方的最大值
				dp[i][j] = max(dp[i - 1][j], dp[i][j - 1]);
			}
		}
	}
	cout << "最大长度:" << dp[len1][len2] << endl;
	cout << "最长公共子序列为:";
	LCS(len1, len2);
}
//时间复杂度O（mn）+O（m+n）