最长公共子序列问题LCS

      一个给定序列的子序列是指在原序列顺序不变的基础上删去若干元素后得到的序列。给定两个序列X和Y，当一个序列Z既是X的子序列又是Y的子序列时，称Z序列为X和Y 的公共子序列。
      例如，X=【A,B,C,B,D,A,B】,Y= 【B,D,C,A,B,A】则序列【B,C,A】是X和Y的一个公共子序列，但不是X和Y的最长子序列。序列【B，C，B，A】是X和Y的最长公共子序列。
      最长公共子序列问题（LCS）：给定两个子序列X=【x1,x2,x3....xm】,Y =【y1,y2,y3,.....yn】，找出X和Y的最长公共子序列。
      1.最优子结构
       设序列X=【x1,x2,x3....xm】,Y =【y1,y2,y3,.....yn】的最长公共子序列为Z=【z1,z2,z3......zk】
则
     （1）若Xm = Yn则Zk = Xm = Yn,且Zk-1是Xm-1和Yn-1的最长公共子序列；
     （2）若Xm != Yn,且Zk != Xm,则Z是Xm-1和Y的最长公共子序列；
       (3)  若 Xm != Yn,且Zk != Yn,则Z是X和Yn-1的最长公共子序列；
      2.子问题递归结构
       当Xm ==yn时，找出Xm-1和Yn-1的最长公共子序列，当Xm !=Yn时，必须解决两个子问题，及找出Xm-1和Y的一个最长公共子序列或者X和Yn-1的一个最长公共子序列。
       求结果过程中会有子问题重叠。
   
        首先建立子问题最优解的递归关系。用C[i][j] 记录X和Y的最长公共子序列的长度。其中Xi ={x1,x2,x3......xm},Yj = {y1,y2,y3,.....yn},当i = 0或j =0 时，空序列是Xi和Yj的最长公共子序列。此时C[i][j] = 0.
        可建立如下递归关系：
                  0                                     i=0,j=0
    C[i][j] =   C[i-1][j-1] +1                   i,j>0;Xi == Yj
                  max{C[i-1][j],C[i][j-1]}      i,j >0,Xi!=Yj

        X,Y两个数组作为输入，输出两个数组c,b.c[i][j]储存Xi和Yj的最长公共子序列的长度，b[i][j]记录c[i][j]的值是由哪一个子问题求解得到的。
void  LCSLength(int m,int n,char *x,char *y,char **c,char **b)
{
       int i,j;
       for(i =0; i <= m; i++)
       {
             c[i][0] = 0;
       }
       for(j = 0; j <= n; j++)
       {
             c[0][i] = 0;
       }
       for(i = 1; i <= m; i++)
       {
              for(j = 1; j <= n;j++)
              {
                      if(x[i] == y[i])
                      {
                              c[i][j] = c[i-1][j-1] + 1;
                              b[i][j] = 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;
                      }
              }
       }
}