输入:两个字符串序列
输出:两个字符串序列的最长公共子序列的长度和子序列
例如:
输入:str1="BDCABA", str2="ABCBDAB"
输出:4 “BDAB”
若此问题采用暴力法求解:枚举str1的每个子序列str,检查str是否为序列str2的子序列,记录当前找到的长度的最大的公共子序列。可以看出暴力法的效率十分低下,并且每次检查str是否为str2的子序列的时间复杂度为O(n),然而str有2m 个子序列,因此采用暴力法的最坏时间复杂度为O(n2m)。
针对此问题可以采用动态规划去做。
令C[i,j]表示str1的前i项和str2的前j项的最长公共子序列的长度,则会有
计算长度算法为LcsLength
打印最长公共子串的算法为LcsPrint
用Java语言的具体实现代码如下:
public class LcsLength {
public static int lcslength(String str1,String str2,char b[][]){
int m=str1.length();
int n=str2.length();
int c[][]=new int[m+1][n+1]; //统计长度
for(int i=0;i<=m;i++)c[i][0]=0;
for(int j=0;j<=n;j++)c[0][j]=0;
for(int i=1;i<=m;i++)
for(int j=1;j<=n;j++)
if(str1.charAt(i-1)==str2.charAt(j-1)){c[i][j]=c[i-1][j-1]+1;b[i][j]='X';} //X表示左上箭头
else if(c[i-1][j]>=c[i][j-1]){c[i][j]=c[i-1][j];b[i][j]='U';} //U表示向上箭头
else {c[i][j]=c[i][j-1];b[i][j]='L';} //L表示向左箭头
return c[m][n]; //返回最长子序列的长度
}
public static void lcsprint(char B[][],String str1,int m,int n){
if(m==0||n==0)return;
if(B[m][n]=='X'){lcsprint(B,str1,m-1,n-1);System.out.print(str1.charAt(m-1));}
else if(B[m][n]=='U')lcsprint(B,str1,m-1,n);
else lcsprint(B,str1,m,n-1);
}
public static void main(String[] args) {
// TODO 自动生成的方法存根
String str1 = new String("BDCABA");
String str2 = new String("ABCBDAB");
System.out.println("第一个字符串str1为:"+str1);
System.out.println("第二个字符串str2为:"+str2);
int m=str1.length();
int n=str2.length();
char B[][]=new char[m+1][n+1];
System.out.println("最长公共子序列的长度为:"+lcslength(str1,str2,B));
System.out.print("最长公共子序列为:");
lcsprint(B,str1,m,n);
}
}
输出为:
第一个字符串str1为:BDCABA
第二个字符串str2为:ABCBDAB
最长公共子序列的长度为:4
最长公共子序列为:BDAB