LCS-入门

最长公共子序列问题（longest-common-subsequence problem）

问题解释：对两个序列X,Y，求他们的长度最长的公共子序列。子序列：例如，X = {5, 6,  7},那么{5} {6} {7} {5,6} {5,7} {6,7}  {5,6,7}, 空集都是X的子序列。

暴力求解：用暴力求解，那么就要穷举X的所有子序列，然后一个个检查X的子序列是否也是Y的子序列，找到最长的公共子序列，m长度的X序列有2^m个，所以暴力方法的时间复杂度为指数阶，并不实用。

1.符号：X[i],Y[j]:X,Y序列;    c[i][j]:表示X序列前i个和Y序列前j个的LCS长度，   b[i][j]:记录i, j的走向

2.递归方程：

两个序列的LCS包含两个序列前缀的LCS，因此，LCS问题具有最优子结构性质。可用动态规划求解。

解释：举个例子，X = {A,B,C,B},Y = {A,B,A,B}，LCS = {A, B, B}, 若设定LCS最后一个B为去掉，它对应的是XY最后一个B，此时X1 = {A, B, C} 是X的前缀， Y2 = {A, B, A}是Y的前缀， 这时，X2和Y2的LCS为{A, B}是X和Y的前缀。

递推公式

c[i][j] = 0;   i = 0 或 j = 0   //当有序列长度为0时是没有公共子序列的

         = c[i - 1][j - 1] + 1;  X[i] = Y[j]   //X序列第i个和Y序列第j个相等， 则加一

         = max{c[i][j - 1],   c[i - 1][j]}    //在X[i]、Y[j - 1]和X[i - 1]、Y[j]中找最长的LCS赋给c[i][j]

用次递归方法人工查找：

X = A B C B D A B

Y = B D C A B A

1.X[1]是否等于Y[1]  ，否

2.找X[2]是否等于Y[1]  或者X[1]是否等于Y[2]， X[2]等于Y[1]，LCS第一个为X[2]

3.X[2 + 1]是否等于Y[1 + 1], 否

。。。重复找LCS

代码：

#include<iostream>
using namespace std;
//定义c[i,j]表示x链前i个和y链前j个的LCS的长度
char X[100]; //x链
char Y[100]; //y链
int c[100][100]; //LCS个数
int b[100][100]; //记录走势,1表示i, j都减一，2表示i-1，3表示j-1 
void LCS(char X[], char Y[], int n, int m)  {  //n，m为X,Y链的长度 
	//第一步先设定临界值 
	for(int i = 0; i <= m; i++)
	    c[0][i] = 0;
	for(int i = 1; i <= n; i++)
	    c[i][0] = 0;
	//主体部分 
	for(int i = 1; i <= n; i++)
	    for(int j = 1; j <= m; 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;
			}
} 
void print(char X[], char Y[], int i, int j)  {
	if(i == 0 || j == 0)  
		return;
	if(b[i][j] == 1)  {
	    print(X, Y, i - 1, j - 1);//从后往前搜索所有的字母后再打印字母 
	    cout<<X[i];
	}
	else if(b[i][j] == 2)
	    print(X, Y, i - 1, j);
	else
	    print(X, Y, i, j - 1);
}
int main()  {
	int n, m;
	X[0] = '0'; Y[0] = '0';  //链的第0个不用,可写可不写 
	cin>>n>>m;
	for(int i = 1; i <= n; i++)
	    cin>>X[i];
	for(int i = 1; i <= m; i++)
	    cin>>Y[i];
	LCS(X,Y,n,m);
	cout<<c[n][m]<<endl;
	print(X, Y, n, m);
} 
/*
7 6
A B C B D A B
B D C A B A
*/

自学算法导论上的LCS，欢迎指错，菜鸟一枚，大神勿喷