动态规划（4）

最长公共子序列问题（LCS问题）

1.问题描述：给定两个序列X<x1,x2,x3....xn>与Y<y1,y2,y3....yn>，求X,Y长度最长的公共子序列。

注意：子序列不一定是连续的，例如：

X=<A,B,C,B,D,A,B>,Y=<B,D,C,A,B,A>两个序列的最长公共子序列为<B,D,A,B>。

下面利用动态规划的思想解决LCS问题。

步骤一：刻画最长公共子序列的特征

 步骤二：构造一个递归解

通过最长公共子序列的特征我们可以知道，在构造递归的时候有两种情况。

在求解X<x1,x2,x3....xn>与Y<y1,y2,y3....yn>的一个LCS时。

1.如果xm=yn，我们就应该求解X（n-1）和Y（n-1）的LCS,然后将xm与yn追加到结尾即可。

2.如果xm！=yn，我们必须求解两个子问题求X（n-1）和Y的LCS与X和Y（n-1）的LCS，选择其中较长的为X和Y的LCS。

因此得到了如下的计算公式：

 步骤三：计算LCS的长度

通过公式我们利用动态规划自底向上的计算。

代码：

 public static void LCS_LENGTH(char X[],char Y[],int b[][],int c[][]){
        int m=X.length;
        int n=Y.length;
        //int b[][]=new int [m+1][n+1];//记录最优解
       // 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(X[i-1]==Y[j-1]){
                    c[i][j]=c[i-1][j-1]+1;
                    b[i][j]=1;//1只是一个标识，没有特殊含义。
                }else if(c[i-1][j]>=c[i][j-1]){
                    c[i][j]=c[i-1][j];
                    b[i][j]=2;
                }else{
                    c[i][j]=c[i][j-1];
                    b[i][j]=3;
                }
            }
        }
    }

运行过程图解：

 步骤四：构造LCS

通过b[][]将最终的LCS打印出来。

代码：

public static void PRINT_LCS(int b[][],char X[],int m,int n){
        if(m==0||n==0){
            return;
        }
        if(b[m][n]==1){
            PRINT_LCS(b,X,m-1,n-1);
            System.out.print(X[m-1]+"->");
        }else if(b[m][n]==2){
            PRINT_LCS(b,X,m-1,n);
        }else{
            PRINT_LCS(b,X,m,n-1);
        }

    }

最终代码：

public static void LCS_LENGTH(char X[],char Y[],int b[][],int c[][]){
        int m=X.length;
        int n=Y.length;
        //int b[][]=new int [m+1][n+1];//记录最优解
       // 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(X[i-1]==Y[j-1]){
                    c[i][j]=c[i-1][j-1]+1;
                    b[i][j]=1;//1只是一个标识，没有特殊含义。
                }else if(c[i-1][j]>=c[i][j-1]){
                    c[i][j]=c[i-1][j];
                    b[i][j]=2;
                }else{
                    c[i][j]=c[i][j-1];
                    b[i][j]=3;
                }
            }
        }
    }

    public static void PRINT_LCS(int b[][],char X[],int m,int n){
        if(m==0||n==0){
            return;
        }
        if(b[m][n]==1){
            PRINT_LCS(b,X,m-1,n-1);
            System.out.print(X[m-1]+"->");
        }else if(b[m][n]==2){
            PRINT_LCS(b,X,m-1,n);
        }else{
            PRINT_LCS(b,X,m,n-1);
        }

    }

    public static void main(String[] args) {
        char X[]={'A','B','C','B','D','A','B'};
        char Y[]={'B','D','C','A','B','A'};
        int m=X.length;
        int n=Y.length;
        int b[][]=new int[m+1][n+1];
        int c[][]=new int[m+1][n+1];
        LCS_LENGTH(X,Y,b,c);
        System.out.println(c[m][n]);
        PRINT_LCS(b,X,m,n);
    }

结果：