求两个字符串的最长公共子序列

 

具体的算法思想参考以下文章：

http://blog.csdn.net/lisonglisonglisong/article/details/41548557

http://blog.csdn.net/zhongkeli/article/details/8847694

https://www.kancloud.cn/digest/pieces-algorithm/163624

 

 

#include "iostream"
#include "string"
#include "vector"
using namespace std;

std::vector<char> common;
int lcs(string A, string B) {
	std::vector<vector<int> > len;
	len.resize(A.size() + 1);
	for (int i = 0; i <= static_cast<int>(A.size()); i++) {
		len[i].resize(B.size() + 1, 0);
	}
	for (int i = 1; i <= static_cast<int>(A.size()); ++i)
	{
		for (int j = 1; j <= static_cast<int>(B.size()); ++j)
		{
			if (A[i - 1] == B[j - 1])
			{
				len[i][j] = len[i - 1][j - 1] + 1;
			}
			else if (len[i - 1][j] >= len[i][j - 1])
			{
				len[i][j] = len[i - 1][j];
			}
			else {
				len[i][j] = len[i][j - 1];
			}
		}
	}
	int apos = A.size();
	int bpos = B.size();
	int commonlen = len[apos][bpos];
	int k = commonlen;
	common.resize(commonlen);
	while (apos && bpos) {
		if (len[apos][bpos] == len[apos - 1][bpos-1] + 1) {
			common[--k] = A[--apos];
			--bpos;
		}
		else if (len[apos - 1][bpos] >= len[apos][bpos - 1])
		{
			--apos;
		}
		else {
			--bpos;
		}
	}
	for (int i = 0; i < commonlen; i++) {
		cout << common[i];
	}
	cout << endl;
	return commonlen;
}


int main()
{
	string A = "aQ11";
	string B = "a11df";
	cout << lcs(A, B);

	system("pause");
	return 0;
}

动态规划解决LCS问题的时间复杂度为O(mn)，这比简单的递归实现要快多了。空间复杂度是O(mn)，因为使用了一个动态规划表。

 

要输出所有LCS的内容

两个字符串对应的最长公共子序列不一定唯一，这个程序输出所有的LCS内容。

基本思想是：

 

具体参考文章：http://blog.csdn.net/lisonglisonglisong/article/details/41596309

#include <vector>
#include <iomanip>
#include <set>
#include <string>
#include <map>
#include <iostream>
#include <string.h>
#include <sstream>
using namespace std;

set<string> all_lcs; //注意这里要用set去除重复的LCS
//最长公共子串（LCS）
//二维数组veca[i][j]记录的是两个字符串Xi和Yj的LCS长度
int LCS_length(const string &str1, const string &str2, vector<vector<int> > &veca) {
	int i, j;
	int biggest = 0;
	if (str1 == "" || str2 == "")
		return 0;

	for (i = 0; i <= str1.length(); i++) {
		veca[i][0] = 0;
	}
	for (j = 0; j <= str2.length(); j++) {
		veca[0][j] = 0;
	}
	for (i = 1; i <= str1.length(); i++) {
		for (j = 1; j <= str2.length(); j++) {
			if (str1[i - 1] == str2[j - 1]) {
				veca[i][j] = veca[i - 1][j - 1] + 1;
			}
			else {
				if (veca[i - 1][j] >= veca[i][j - 1])
					veca[i][j] = veca[i - 1][j];
				else
					veca[i][j] = veca[i][j-1];
			}
		}
	}
	return veca[str1.length()][str2.length()];
}

//该函数找出所有的LCS的序列，并将其存在vector中
void PrintAllLCS(string &str1, string &str2, int i, int j, 
	             vector<vector<int> > &veca, string lcs_str) {
//注意这里形参lcs_str不可以为引用，这里需要每次调用lcs_str都重新生成一个对象
	while (i > 0 && j > 0) {
		if (str1[i - 1] == str2[j - 1]) {
			lcs_str = str1[i - 1] + lcs_str; //逆向存放
			--i;
			--j;
		}
		else {
			if (veca[i - 1][j] > veca[i][j - 1]) //向左走
				--i;
			else if (veca[i - 1][j] < veca[i][j - 1]) //向上走
				--j;
			else { //此时向上向右均为LCS的元素
				PrintAllLCS(str1, str2, i - 1, j, veca, lcs_str);
				PrintAllLCS(str1, str2, i, j - 1, veca, lcs_str);
				return;
			}
		}
	}
	cout << "   " << lcs_str << endl;
	all_lcs.insert(lcs_str);
}
int main() {
	string input;
	getline(cin, input);
	stringstream ss(input);
	string str1, str2;
	ss >> str1;
	ss >> str2;
	//将veca初始化为一个二维数组,其行列值分别为str1和str2的长度加1
	//二维数组veca记录的是两个字符串Xi和Yj的LCS长度
	vector<vector<int> > veca(str1.length() + 1, vector<int>(str2.length() + 1));
	cout << LCS_length(str1, str2, veca) << endl;

	string lcs_str;
	PrintAllLCS(str1, str2, str1.length(), str2.length(), veca, lcs_str);
	set<string>::iterator iter = all_lcs.begin();
	while (iter != all_lcs.end()) {
		cout << *iter++ << endl;
	}
	return 0;
}

如图所示的两个字符串共有三个LCS。