LCS（最长公共子序列）介绍

LCS（最长公共子序列）

LCS是Longest Common Subsequence的缩写，即最长公共子序列。一个序列，如果是两个或多个已知序列的子序列，且是所有子序列中最长的，则为最长公共子序列。

解决方法

对于一般的LCS问题，都属于NP问题。当数列的量为一定的时，都可以采用动态规划去解决

算法

设两个字符串x，y，c[i][j]代表x的i位和y的j位之前的最长公共子序列的长度。

则如果x[i]==y[j]，那么c[i][j]=c[i-1]c[j-1]+1。如果x[i]!=y[j]，那么c[i][j]一定等于c[i-1][j]或者c[i][j-1]。

如果我们需要把最长公共子序列找出，我们可以通过记录每次最长公共子序列长度的情况来找出子序列。

以下是代码实现：

void lcs(char *x, int m, char *y, int n)
{
    int i, j;
    memset(c, 0, sizeof(c));
    memset(b, 0, sizeof(b));
    for (i=1; i<=m; i++) {
        for (j=1; j<=n; j++) {
            if (x[i] == y[j]) {
                c[i][j] = c[i-i][j-1] + 1;
                b[i][j] = 1;
            } else if (c[i][j-1] < c[i-1][j]) {
                c[i][j] = c[i-1][j];
                b[i][j] = 2;
            } else {
                c[i][j] = c[i][j-1];
                b[i][j] = 3;
            }
        }
    }
}

void show(int i, int j, char *c) {
    if (i == 0 || j == 0) {
        return;
    }
    if (b[i][j] == 1) {
        printf("%c", c[i-1]);
        show(i-1, j-1, c);
    } else if (b[i][j] == 2) {
        show(i-1, j, c);
    } else {
        show(i, j-1, c);
    }
}

int main(int argc, const char * argv[])
{
    
    char x[] = "aabcd";
    char y[] = "12abcabcd";
    int m = (int)strlen(x);
    int n = (int)strlen(y);
    lcs(x, m, y, n);
    show(m, n, x);
    printf("\n");
    return 0;
}