首先明白两个区别:
1、最长公共子串:子串是串的一个连续的部分,在原字符串中位置是连续的
2、最长公共子序:不改变序列的顺序,从序列中去掉任意的元素而获得新的序列,也就是说子序在原字符串中位置可以不连续。
递推方程:
步骤:1序列str1和序列str2
·长度分别为m和n;
·创建1个二维数组L[m.n];
·初始化L数组内容为0
·m和n分别从0开始,m++,n++循环:
- 如果str1[m] == str2[n],则L[m,n] = L[m - 1, n -1] + 1;
- 如果str1[m] != str2[n],则L[m,n] = max{L[m,n - 1],L[m - 1, n]}
·最后L[m,n]中的数字一定是最大的,且这个数字就是序列str1和序列str2最长公共子序列的长度
步骤:2从数组L中查找一个最长的公共子序列
i和j分别从m,n开始,递减循环直到i = 0,j = 0。其中,m和n分别为两个串的长度。
·如果str1[i] == str2[j],则将str[i]字符插入到子序列内,i--,j--;
·如果str1[i] != str[j],则比较L[i,j-1]与L[i-1,j],L[i,j-1]大,则j--,否则i--;(如果相等,则任选一个)
#include<stdio.h>
#include<iostream>
#include<string>
using namespace std;
#define max(a,b) ((a)>(b)?(a):(b));
void LCS(string str1, string str2);
int main()
{
string str1 = "ABCBDAB";
string str2 = "BDCABA";
LCS(str1,str2);
return 0;
}
void LCS(string str1, string str2)
{
int i,j;
//申请空间
int ** martix = new int *[str1.length()+1];
for(i = 0;i<=str1.length();i++)
martix[i] = new int[str2.length()+1];
//初始化边界,过滤掉0的情况
for (i = 0; i <= str1.length(); i++)
martix[i][0] = 0;
for (j = 0; j <= str2.length(); j++)
martix[0][j] = 0;
//填充矩阵
for (i = 1; i <= str1.length(); i++)
{
for (j = 1; j <= str2.length(); j++)
{
//相等的情况
if (str1[i - 1] == str2[j - 1])
{
martix[i][j] = martix[i-1][j-1] + 1;
}
else
{
martix[i][j]=max(martix[i-1][j],martix[i][j-1]);
}
}
}
cout<<"最长公共子序列长度为:"<<martix[str1.length()][str2.length()]<<endl;
cout<<"输出矩阵为:"<<endl;
for(i = 0 ; i <=str1.length(); i++)
{
for( j = 0; j <=str2.length(); j++)
{
cout << martix[i][j] << " ";
}
cout << endl;
}
cout<<"最长公共子序列长度为:"<<endl;
for(i=str1.length(), j=str2.length(); i >= 1 && j >= 1;)
{
if(str1[i - 1] == str2[j - 1])
{
cout << str1[i - 1] << " ";//倒序打印的
i--;
j--;
}
else
{
if(martix[i][j -1] > martix[i - 1][j])
{
j--;
}
else
{
i--;
}
}
}
//释放空间
for(i = 0;i < str1.length();i++)
{
delete martix[i];
martix[i] = NULL;
}
delete [str1.length()]martix;
martix = NULL;
}