作业9——动态规划之LCS问题

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;
}