打印最长的公共子序列

时间:2014-10-12 14:07:00

标签: c++ algorithm

我在最长的子序列中遇到的问题: 对于Ex:

“ABCDGH” and “AEDFHR” is “ADH” of length 3

代码:

void lcs( char *X, char *Y, int m, int n )
{
   int L[m+1][n+1];

   /* Following steps build L[m+1][n+1] in bottom up fashion. Note
      that L[i][j] contains length of LCS of X[0..i-1] and Y[0..j-1] */
   for (int i=0; i<=m; i++)
   {
     for (int j=0; j<=n; j++)
     {
       if (i == 0 || j == 0)
         L[i][j] = 0;
       else if (X[i-1] == Y[j-1])
         L[i][j] = L[i-1][j-1] + 1;
       else
         L[i][j] = max(L[i-1][j], L[i][j-1]);
     }
   }

我不明白为什么这行代码:

L[i][j] = max(L[i-1][j], L[i][j-1]);

如果我想打印序列怎么办,即ADH我该怎么做?

1 个答案:

答案 0 :(得分:1)

  

我不明白为什么这行代码:

L[i][j] = max(L[i-1][j], L[i][j-1]);

如果X[0..i-1]的最后一个字符与Y[0..j-1]的最后一个字符不匹配,则对于每个公共子序列,这些字符中至少有一个不属于。因此,答案由X[0..i-2]Y[0..j-1]的最大长度或X[0..i-1]Y[0..j-2]的最大长度给出。

要恢复实际的子序列,我们必须追溯这些决定。

char lcs[min(m, n) + 1];
char *end = lcs;
int i = m;
int j = n;
while (i > 0 && j > 0) {
    if (X[i - 1] == Y[j - 1]) {
        *end = X[i - 1];
        end++;
        i--;
        j--;
    } else if (L[i - 1][j] >= L[i][j - 1]) {
        i--;
    } else {
        j--;
    }
}
reverse(lcs, end);
*end = '\0';
puts(lcs);