一.题目描述
在字母表上分别给出两个长度为m和n的字符串x和y,确定在x和y中的最长公共子序列。例如,x为“ABCBDAB”,y为“BDCABA”,则两者的最长公共子序列为BDAB。
二.思路分析
假设两个字符串为x1,x2…xi和y1,y2…yj,再令L【i,j】表示这两个的最长公共子序列。思路如下:如果xi==yj(即最后一个字母相等),则,L[i,j]=L【i-1,j-1】+1.如果xi!=yj,则,L【i,j】=max{L【i-1,j】,L【i,j-1】}.
说白了就是,如果最后一个字母相等,则最长子序列等于前面的子序列再加上最后这个相等的。如果不想等,就在两者间取个最大的(或者理解为,平移)。
其实最根本的思想就是递归。当然这是动态规划的思想。动态规划与递归的区别在于,为了节约重复求相同子问题的时间,引入一个数组,不管它们是否对最终解有用,把所有子问题的解存于该数组中,这就是动态规划法所采用的基本方法。
三.代码实现
#include<iostream>;
#include<cstring>;
using namespace std;
#define MAXLEN 100
void LCSLength(char *x,char*y,int m,int n,int c[][MAXLEN],int b[][MAXLEN]){
//二维数组c,存储的是LCS(最长公共子序列)的中间值,b存储的是三种情况,即xi==yj
//和不等时的两种情况。x,y代表两个字符串,m和n是他们的长度。
for(int i=0;i<=m;i++){
c[i][0]=0;//根据前面讲的数组c的定义,很明显,c[i][0]全部为0.
}
for(int i=1;i<=n;i++){
c[0][i]=0;
}
for(int i=1;i<=m;i++){
for(int j=1;j<=n;j++){
if(x[i-1]==y[j-1]){
c[i][j]=c[i-1][j-1]+1;
b[i][j]=0;
}
else if(c[i][j-1]>=c[i-1][j])
{c[i][j]=c[i][j-1];
b[i][j]=1;
}
else {
c[i][j]=c[i-1][j];
b[i][j]=-1;
}
}
}
}
void printLCS(char*x,int m,int n,int b[][MAXLEN]){
if(m==0||n==0)
return;
if(b[m][n]==0)
{
printLCS(x,m-1,n-1,b);
cout<<x[m-1];
}
else if(b[m][n]==1){
printLCS(x,m,n-1,b);
}
else printLCS(x,m-1,n,b);
}
int main(){
char x[MAXLEN]="ABCBDAB";
char y[MAXLEN]="BDCABA";
int b[MAXLEN][MAXLEN];
int c[MAXLEN][MAXLEN];
int m=strlen(x);
int n=strlen(y);
LCSLength(x,y,m,n,c,b);
printLCS(x,m,n,b);
return 0;
}