最长公共子序列问题

输入：两个字符串序列

输出：两个字符串序列的最长公共子序列的长度和子序列

例如：

输入：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项的最长公共子序列的长度，则会有

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]
;                                                                      //返回最长子序列的长度
}

public static void lcsprint(char B[][],String str1,int m,int n){
if(m==0||n==0)return;
if(B[m]
=='X'){lcsprint(B,str1,m-1,n-1);System.out.print(str1.charAt(m-1));}
else if(B[m]
=='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);

}

}

View Code
输出为：

第一个字符串str1为：BDCABA
第二个字符串str2为：ABCBDAB
最长公共子序列的长度为：4
最长公共子序列为：BDAB