1,问题
求两个已知序列的最长公共子序列,一个序列的子序列满足下标严格递增。
2,解析
如果Zk 是 Xi 和 Yj 的最长公共子序列(下标表示该序列的长度)
则有以下三点规则:
(1) 如果X[i] = Y[j],那么Z[k] = X[i] =Y[j],Zk-1 是 Xi-1 和 Yj-1 的最长公共子序列;
(2) 如果X[i] ≠ Y[j],Z[k] ≠ X[i] ,那么Zk 是 Xi-1 和 Yj 的最长公共子序列;
(3) 如果X[i] ≠ Y[j],Z[k] ≠ Y[j] ,那么Zk 是 Xi 和 Yj-1 的最长公共子序列;
得出算法递推关系:
C[i][j]表示X序列前i个元素和Y序列前j个元素的最长公共子序列长度(int), i < Xlen, j < Ylen
C[i][j] = C[i-1][j-1] + 1 , X[i] = Y[j]
C[i][j] = MAX{ C[i][j-1], C[i-1][j] }, X[i] != Y[j]
3,设计
C[i][j] 存放两个序列的最长公共子序列的长度
B[i][j] 存放两个序列与[i-1],[j-1]两个序列的关系
stack,存放子序列的栈
初始化:C[0][j] 和C[i][0]都为0
for i < Xlen
for j < Ylen
if(X[i]==Y[j]) C[i][j] = C[i-1][j-1] +1 B[i][j] = ‘2’;
else C[i][j] = C[i-1][j] >= C[i][j-1] ? C[i-1][j] : C[i][j-1] B[i][j] =‘X’ || ‘Y’
for i = Xlen, j = Ylen; i>0 && j>0
if (B[i][j] == ‘X’) i–;
else if (B[i][j] == ‘Y’) j–;
else stack.push(X[i]) i-- j–
逐个pop stack;
4,分析
T = O(Xlen * Ylen)
5,源码
#include<stdio.h>
#include<stdlib.h>
#include<string.h>
#include<stack>
using namespace std;
#define Xlen 5
#define Ylen 4
void getLCS(char X[],char Y[]) {
int C[Xlen + 1][Ylen + 1];
char B[Xlen + 1][Ylen + 1];
stack<char> s;
for (int i = 0; i < Xlen; i++) {
C[i][0] = 0;
}
for (int j = 0; j < Ylen; j++) {
C[0][j] = 0;
}
for (int k = 1; k <= Xlen; k++) {
for (int l = 1; l <= Ylen; l++) {
if (X[k] == Y[l]) {
C[k][l] = C[k - 1][l - 1] + 1;
B[k][l] = '2';
}
else if (C[k - 1][l] >= C[k][l - 1]) {
C[k][l] = C[k - 1][l];
B[k][l] = 'X';
}
else {
C[k][l] = C[k][l - 1];
B[k][l] = 'Y';
}
}
}
printf("%d\n", C[Xlen][Ylen]);
for (int i = Xlen,j = Ylen;j > 0 && i > 0;) {
if (B[i][j] == 'X') i--;
else if (B[i][j] == 'Y') j--;
else {
s.push(X[i]);
i--;
j--;
}
}
while (!s.empty()) {
printf("%c ", s.top());
s.pop();
}
}
int main() {
char X[Xlen+2];
char Y[Ylen+2];
for (int i = 1;; i++) {
char ch = getchar();
if (i == Xlen + 1 || i == '\n') {
X[i] = '\0';
break;
}
X[i] = ch;
}for (int i = 1;; i++) {
char ch = getchar();
if (i == Ylen + 1 || i == '\n') {
Y[i] = '\0';
break;
}
Y[i] = ch;
}
getLCS(X, Y);
return 0;
}