最长公共子序列

问题叙述：

子序列形式化定义：
给定一个序列X=<x1,x2,x3,x4...,xm>，另一个序列Z=<z1,z2,z3,z4...,zk>，若存在一个严格递增的X的下标序列<i1,i2,i3,...,ik>对所有的1,2,3,...,k，都满足x(ik)=zk，则称Z是X的子序列

比如Z=<B,C,D,B>是X=<A,B,C,B,D,A,B>的子序列

公共子序列定义：
如果Z既是X的子序列，又是Y的子序列，则称Z为X和Y的公共子序列

最长公共子序列（以下简称LCS）：
2个序列的子序列中长度最长的那个
 

动态规划求解最长公共子序列：
分析规律：
设X=<x1,x2,x3,x4...,xm>，Y=<y1,y2,y3,y4...,yn>为两个序列，Z=<z1,z2,z3,z4...,zk>是他们的任意公共子序列C[i][j]来代表目前最长公共序列的长度

经过分析，我们可以知道：

1、如果xm = yn，则zk = xm = yn 且 Zk-1是Xm-1和Yn-1的一个LCS

2、如果xm != yn 且 zk != xm，则Z是Xm-1和Y的一个LCS

3、如果xm != yn 且 zk != yn，则Z是X和Yn-1的一个LCS

所以如果用一个二维数组c表示字符串X和Y中对应的前i，前j个字符的LCS的长度话，可以得到以下公式：

 我们可以根据图来理解（求ABCBDAB和BDCABA的LCS:）

我们采用好理解的方式去想：

如果箭头是↖，则说明c[i][j]的公共长度也是c[i-1][j-1]的公共长度+1，条件为X[i]=Y[j]

如果箭头是←，则说明c[i][j]的公共长度也是c[i][j-1]的公共长度，但此时X[i]!=Y[j]，条件为

c[i-1][j]<c[i][j-1],如i=3，j=4时;

如果箭头是↑，则说明c[i][j]的公共长度也是c[i-1][j]的公共长度，但此时X[i]!=Y[j]，条件为

c[i-1][j]>c[i][j-1];如i=3，j=4时;

c[i-1][j]与c[i][j-1]的比较时为了检验目前的c[i][j]由谁传过来的值，即谁目前的最长公共子序列大

/**
题目就是让我们求两个字符串的 LCS 长度：

输入: str1 = "abcde", str2 = "ace"

输出: 3

解释: 最长公共子序列是 "ace"，它的长度是 3
*/
#include<iostream>
#include<cstring>
using namespace std;

int len=0;
void LCSLength(int m,int n,char *x,char *y,int **c,int **b){
	for(int i=1;i<=m;i++){
		c[i][0]=0;
	}
	for(int j=1;j<=n;j++){
		c[0][j]=0;
	}
	for(int i=1;i<=m;i++){
		for(int j=1;j<=n;j++){
			if(x[i]==y[j]){//两个相等，则他们x[i-1][j-1]也相等，所以c[i][j]=c[i-1][j-1]+1
				c[i][j]=c[i-1][j-1]+1;
				b[i][j]=1;
			} 
			else if(c[i-1][j]>c[i][j-1]){//x的元素多，y的元素少则c[i][j]从上方继承 
				c[i][j]=c[i-1][j];
				b[i][j]=2;
			}
			else {      //y的元素多，x的元素少则c[i][j]从左方继承 
				c[i][j]=c[i][j-1];
				b[i][j]=3;
			}
		} 
	}
}
void LCS(int i,int j,char *x,int **b){
	if(i==0||j==0){
		return;
	}
	//b[i][j]==1,说明这个值是从左上角来的即c[i-1][j-1];	
	if(b[i][j]==1) {
		LCS(i-1, j-1, x, b);
		cout<<x[i];
        //若输出长度直接在后边定义一个len，输出一个len++；
        //len++;
	}
	//b[i][j] == 2说明这个值是从上方来的即c[i-1][j];
	else if(b[i][j] == 2) {
		LCS(i-1, j, x, b);	
	}
	//b[i][j] == 3明这个值是从左方来的即c[i][j-1];
	else 
	LCS(i, j-1, x, b);
} 

int main()
{
	int m = 0, n = 0;
	char str1[10000], str2[10000];
	int *b[10000], *c[10000];//**b匹配的是指针数组 而不是b[10000][10000]
	for (int i = 0; i < 10000; i++)//划重点
	{
		b[i] = new int[10000];
		c[i] = new int[10000];
	}

	cout << "请输入第一串数字/字母的长度：";
	cin >> m;
	cout << "请输入第二串数字/字母的长度：";
	cin >> n;


	cout << "请输入第一串数字/字母:";
	for (int i = 1; i <= m; i++)//这里用1~m 而不是 0~m-1 因为在上面的函数中初始化为0的那两行是没意义的
		cin >> str1[i];
	cout << "请输入第二串数字/字母:";
	for (int i = 1; i <= n; i++)
		cin >> str2[i];
	LCSLength(m,n,str1,str2,c,b);
	LCS(m,n,str1,b);
	//cout<<"  "<<"最长公共子序列的长度为："<<len;
	return 0;
}