题目
求最长公共子序列
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)