最长公共子序列 过程图解

最长公共子序列 过程图解

顾名思义，就是判断两个序列之间的最长的一样的子序列，这里需要注意最长公共子序列（longest commom sequence LCS)和最长公共子串(longest commom substring)是两个概念。

给定字符串中任意个连续的字符组成的子序列称为该串的子串，而子序列是将给定序列中零个或多个元素去掉之后的结果

求解LCS问题，可以使用动态规划来做，使用从c[i][j]来记录序列$X_i$和$Y_i$的最长公共子序列的长度。由最优字结构可建立递归关系如下

$c[i][j]=\{\begin{array}{l}0,i>0;j=0\\ c[i-1][j-1]+1,i,j>0;x_i=y_i\\ max(c[i][j-1],c[i-1][j]),i,j>0;x_i!=y_i\end{array}$

计算最优值算法伪代码

void LCSLength(int m ,int n, char *x,char *y,int **c,int **b){
    int i,j;
    for(i = 1; i < m; i++){
        c[i][0] = 0;
    }
    for(i = 1; i < n; i++){
        c[0][i] = 0;
    }
    for(i = 1; i <= m; i++){
        for(j = 1; j <= n; j++){
            if(x[i] == y[j]){
                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;
            }
        }
    }
}

构造最长公共子序列

由算法LCSLength计算得到数组b可用于快速构造最长公共子序列。首先从b[m][n]开始，依其值在数组b中搜索。

• 当b[i][j] = 1时，表示$X_i$和$Y_i$的最长公共子序列是由$X_{i-1}$和$Y_{j-1}$的最长公共子序列在尾部加上$x_i$所得到的子序列
• 当b[i][j] = 2时，表示$X_i$和$Y_i$的最长公共子序列与$X_{i-1}$和$Y_j$的最长公共子序列相同
• 当b[i][j] = 3时，表示$X_i$和$Y_i$的最长公共子序列与$X_i$和$Y_{j-1}$的最长公共子序列相同

伪代码如下

void LCS(int i,int j,char *x,char **b){
    if(i == 0 || j == 0)
    return;
    if(b[i][j] == 1){
        LCS(i-1,j-1,x,b);
        cout<<x[i];
    }
    else if(b[i][j] == 2){
        LCS(i-1,j,x,b);
    }
    else
        LCS(i,j-1,x,b);
}

案例学习

求$X=ABACDABCBA$ ,$Y=BCABDCAB$的最长公共子序列