算法设计与分析--最长公共子序列

最长公共子序列

问题描述：
给定两个字符串，求解这两个字符串的最长公共子序列（Longest Common Sequence）。比如字符串1：BDCABA；字符串2：ABCBDAB
则这两个字符串的最长公共子序列长度为4，最长公共子序列是：BCBA
问题解决：
1>.穷举法：列举出所有的可能，但其时间复杂度较高不可取。
2>.递归法：
设X=<x1,x2,x3,x4…,xm>，Y=<y1,y2,y3,y4…,yn>为两个序列，Z=<z1,z2,z3,z4…,zk>是他们的任意公共子序列

经过分析，我们可以知道：

1、如果xm = yn，则zk = xm = yn 且 Zk-1是Xm-1和Yn-1的一个LCS

2、如果xm != yn 且 zk != xm，则Z是Xm-1和Y的一个LCS

3、如果xm != yn 且 zk != yn，则Z是X和Yn-1的一个LCS

所以如果用一个二维数组c表示字符串X和Y中对应的前i，前j个字符的LCS的长度话，可以得到以下公式：
代码如下：

// 最长公共子序列
    public static int num=0;
    public static int LcsLength(String x ,String y,int i,int j,int[][] c,int[][] s){
        //记录程序运行的次数
        num++;
        if(i==0||j==0)return 0;
        if(c[i][j]>0) return c[i][j];
      if(x.charAt(i)==y.charAt(j)){
          c[i][j]= LcsLength(x,y,i-1,j-1,c,s)+1;
          s[i][j]=1;
      }
      else {
          int max1=LcsLength(x,y,i-1,j,c,s);
          int max2=LcsLength(x,y,i,j-1,c,s);
         // c[i][j]= max1>max2?max1:max2;
          if(max1>max2) {
              c[i][j]=max1;
              s[i][j]=2;
          }
          else {
              c[i][j]=max2;
              s[i][j]=3;
          }
      }
      return c[i][j];
    }
    
 public static void TrackBack(int[][]s,String x,int i,int j){
        if(i==0||j==0) return;
        if(s[i][j]==1){
            TrackBack(s,x,i-1,j-1);
            System.out.println(x.charAt(i));
        }
        //行-1
        else if(s[i][j]==2){
            TrackBack(s,x,i-1,j);
        }
        else TrackBack(s,x,i,j-1);
 }
public static void main(String[] args){
String x=new String("#ABCBDAB");
        String y=new String("#BDCABA");
        int[][] c=new int[x.length()][y.length()];
        int[][] s=new int[x.length()][y.length()];
        int maxlen=LcsLength(x,y,x.length()-1,y.length()-1,c,s);
        System.out.println(maxlen);
        System.out.println(num);
        //打印c
        for (int i = 0; i < c.length; i++) {
            for (int j = 0; j < c[i].length; j++) {
                System.out.print(c[i][j]+" ");
            }
            System.out.println();
        }
        //打印s
        for (int i = 0; i < s.length; i++) {
            for (int j = 0; j < s[i].length; j++) {
                System.out.print(s[i][j]+" ");
            }
            System.out.println();
        }
        TrackBack(s,x,x.length()-1,y.length()-1);
    }