最长公共子序列(Longest Common Subsequence)与最长公共子串(Longest Common Substring)的区别: 子序列问题只需保持字符相对顺序一致,并不要求连续,而子串问题要求字符在原字符串中是连续的。
一、最长公共子序列(Longest Common Subsequence)(动态规划)
根据下列公式来编写程序:
程序步骤如下:
1、构造dp数组(二维数组多出一行一列来保存dp初始状态的0值);
2、正序构造dp数组:根据公式给dp数组赋值(注意下标值(str1[i-1]
与str2[j-1]
)与dp[i][j]
对应));
3、结果输出:倒序遍历dp数组,如果当前位置
的str1[i-1]==str2[j-1]
, 则将该点对应的值插入的结果字符串的头部并且i--; j--;
,如果当前位置
的dp值
与上方
或左方
的dp值
相同,则向上或向左移动。如此循环,直到i < 0 || j < 0
。 (因为最长公共子序列不一定唯一,这样只是输出一种方案,所有的LCS输出见PrintAllLCS()
)。
4、输出LCSubsequence;
#include <iostream>
#include <string>
#include <vector>
#include <set>
using namespace std;
set<string> allRes;
int Max(int a, int b) {
return a > b ? a : b;
}
void PrintAllLCSubsequence(string &str1, string &str2, int i, int j,
vector<vector<int>> &dpvec, string res) {
//print all LCS
if (i == 0 || j == 0) {
if (allRes.find(res) == allRes.end()) {
cout << res << endl;
allRes.insert(res);
}
return ;
}
if (i > 0 && j > 0) {
if (str1[i-1] == str2[j-1]) {
string tmp = str1[i-1] + res;
PrintAllLCSubsequence(str1, str2, i-1, j-1, dpvec, tmp);
} else {
if (dpvec[i][j] == dpvec[i-1][j])
PrintAllLCSubsequence(str1, str2, i-1, j, dpvec, res);
if (dpvec[i][j] == dpvec[i][j-1])
PrintAllLCSubsequence(str1, str2, i, j-1, dpvec, res);
}
}
}
void PrintOneLCSubsequence(string &str1, string &str2, vector<vector<int>> &dpvec) {
//print one LCS
string res = "";
int i = str1.size();
int j = str2.size();
while (i > 0 && j > 0) {
if (str1[i-1] == str2[j-1]) {
res = str1[i-1] + res;
i--;
j--;
} else {
if (dpvec[i][j] == dpvec[i-1][j])
i--;
else if (dpvec[i][j] == dpvec[i][j-1])
j--;
}
}
cout << res << endl;
}
int LCSubsequence(string &str1, string &str2) {
int len1 = str1.size();
int len2 = str2.size();
if (len1 == 0 || len2 == 0)
return 0;
vector<vector<int>> dp(len1+1, vector<int>(len2+1));
for (int i = 0; i <= len1; i++) {
for (int j = 0; j <=len2; j++) {
if (i == 0 || j == 0)
dp[i][j] = 0;
else if (str1[i-1] == str2[j-1])
dp[i][j] = dp[i-1][j-1] + 1;
else
dp[i][j] = Max(dp[i-1][j], dp[i][j-1]);
}
}
PrintOneLCSubsequence(str1, str2, dp);
//PrintAllLCSubsequence(str1, str2, len1, len2, dp, "");
return dp[len1][len2];
}
二、最长公共子串(Longest Common Substring)(动态规划)
这里的最大公共字串要求的字串是连续的。
求字串的方法和求子序列方法类似:
当str1[i] == str2[j]
时,子序列长度dp[i][j] = dp[i - 1][j - 1] + 1
;只是当str1[i] != str2[j]
时,dp[i][j] = 0
,而不是max{dp[i-1][j], dp[i][j-1]}
。
#include <iostream>
#include <string>
#include <vector>
using namespace std;
void PrintOneLCSubstring(string &str1, string &str2, vector<vector<int>> &dpvec, int maxLen) {
//print LCS
if (maxLen != 0) {
for (int i = str1.size(); i > 0; i--) {
for (int j = str2.size(); j > 0; j--) {
if (dpvec[i][j] == maxLen) {
string res = "";
int a = i;
int b = j;
while (a > 0 && b > 0) {
res = str1[a-1] + res;
--a;
--b;
}
cout << res << endl;
return ;
}
}
}
}
return ;
}
int LCSubstring(string &str1, string &str2) {
int len1 = str1.size();
int len2 = str2.size();
if (len1 == 0 || len2 == 0)
return 0;
int maxLen = 0;
vector<vector<int>> dp(len1+1, vector<int>(len2+1));
for (int i = 0; i <= len1; i++) {
for (int j = 0; j <=len2; j++) {
if (i == 0 || j == 0) {
dp[i][j] = 0;
} else if (str1[i-1] == str2[j-1]) {
dp[i][j] = dp[i-1][j-1] + 1;
if (dp[i][j] > maxLen)
maxLen = dp[i][j];
} else
dp[i][j] = 0;
}
}
PrintOneLCSubstring(str1, str2, dp, maxLen);
return maxLen;
}
int main() {
string X = "asdfas";
string Y = "werasdfaswer";
cout << "The length of LCS is " << LCSubstring(X, Y);
cout << endl;
return 0;
}