最长公共子串和最长公共序列

1. 最长公共子串

注意子串是连续的。有下列动态转移方程

c[i][j] = c[i-1][j-1] + 1   when X[i] = Y[j]

c[i][j] = 0   when X[i] != Y[j]

int c[100][100];

int LCS(char x[], int len_x, char y[], int len_y){

       int max_len = 0;

       for(int i =0; i < len_x ; i++){
           for(int j = 0; j < len_y;  j++){
                 if(x[i] == y[j]){
                      if(i && j)
                          c[i][j] = c[i-1][j-1] + 1;
                      if(i == 0 || j == 0)
                          c[i][j] = 1;
                      if(c[i][j] > max_len)
                           max_len = c[i][j];
                }
           }
      }

      return max_len;
}    

2. 最长公共子序列

子序列不要求连续。

有下列状态转移方程

c[i][j] = c[i-1][j-1] + 1 if x[i] == y[j]

c[i][j] = max{c[i-1][j], c[i][j-1]} if x[i] != y[j]

// x = "AB"   y = "CD"

int c[100][100] = 0;

int lcs(char x[], char y[]){
     int m = strlen(x);
     int n = strlen(y);

     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;
             else
                   c[i][j] = max(c[i][j-1], c[i-1][j]);
       }

    return c[m][n];
}