一、子序列的定义
给定一个序列X = <x1,x2,x3…xm>,对于Z = <z1, z2, z3…zk>,满足如下条件时,就成为X的子序列:存在一个严格递增的X的下标序列,对所有j = 1,2…k满足xij = zj,比如X = <A,B,C,B,D,A,B>,Y = <B,D,C,A,B,A>,那么<B,C,A>就是X和Y的公共子序列,但最长的公共子序列是<B,C,B,A>。
如何解决问题
步骤1:刻画最长公共子序列的特征
步骤2:一个递归解
1、上图意味着,在求X=<x1,x2,…,xm〉和Y=<y1,y2,…,yn,>的一个LCS时,我们需要求解一个或两个子问题。如果xm=yn,我们应该求解Xm-1和Yn-1的一个LCS。将xm=yn追加到这个LCS的末尾,就得到X和Y的一个LCS.如果xm≠y,我们必须求解两个子问题:求Xm-1和Y的一个LCS与X和Yn-1的一个LCS。两个LCS较长者即为X和Y的一个LCS。
2、我们很容易就看出其中的重叠子问题性质。我们都要求Xm-1和Yn-1的LCS(可以依据情况在后面进行追加或者是不追加)。
3、用c[i,j]表示Xi和Yj的LCS的长度。
步骤3:计算LCS的长度
LCS-LENGTH(X,Y)
m = X.length
n = Y.length
let b[1...m,1...n] and c[0...m,0...n] be new tables
for i = 1 to m
c[i,0] = 0
for j = 0 to n
c[0,j] = 0
for i = 1 to m
for j = 1 to n
if xi = yj
c[i,j] = c[i-1,j-1] + 1
b[i,j] = "↖"
elseif c[i-1,j]>=c[i,j-1]
c[i,j] = c[i-1,j]
b[i,j] = "↑"
else c[i,j]=c[i,j-1]
b[i,j] = "←"
return c and b
构造LCS
按照上图给出的箭头,我们就可以追踪到所有的元素。
PRINT-LCS(b,X,i,j)
if i==0 or j==0
return
if b[i,j]=="↖"
PRINT-LCS(b,X,i-1,j-1)
print xi
elseif b[i,j]=="↑"
PRINT-LCS(b,X,i-1,j)
else PRINT-LCS(b,X,i,j-1)
C++代码
#include<iostream>
#include<string>
#define NUM 20
using namespace std;
int c[NUM][NUM];
char b[NUM][NUM];
void LCS_LENGTH(string X, string Y) {
int m = X.length();
int n = Y.length();
for (int i = 1; i <= m; i++) {
c[i][0] = 0;
}
for (int j = 0; j <= n; j++) {
c[0][j] = 0;
}
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;
b[i][j] = 'a';
}
else if (c[i - 1][j] >= c[i][j - 1]) {
c[i][j] = c[i - 1][j];
b[i][j] = 'b';
}
else {
c[i][j] = c[i][j - 1];
b[i][j] = 'c';
}
}
}
}
void PRINT_LCS(string X, int m, int n)
{
if (m == 0 || n == 0)
return;
if (b[m][n] == 'a') {
PRINT_LCS(X, m - 1, n - 1);
cout << X[m - 1];
}
else if (b[m][n] == 'b') {
PRINT_LCS(X, m - 1, n);
}
else
PRINT_LCS(X, m, n - 1);
}
int main()
{
string X, Y;
cin >> X >> Y;
LCS_LENGTH(X, Y);
PRINT_LCS(X, X.length(), Y.length());
return 0;
}