最长公共子序列 动态规划法（C++实现）

问题：对给定序列X，Y，最长子序列问题（Longest Public Subsequence Problem）就是在序列X和Y中查找长度最长的公共子序列。

算法：递归m（X的长度）乘 n（Y的长度），得到每个点（e.g.X的前两个数据和Y的前4个数据的子序列）的最大子序列长度，并记录下状态，回溯。

代码如下：

same[] 存储公共子序列，length[][] 存储最长公共子序列的长度，state[][]存储相应的状态

int CommonOrder(char x[],int m,char y[],int n,char same[]){
	int i,j,k,length[m+1][n+1],state[m+1][n+1];//length公共长度，state状态
	for(j=0;j<=n;j++)//初始化第0行
		length[0][j] = 0;
	for(i=0;i<=m;i++)//初始化第0列
		length[i][0] = 0;
	for(i=1;i<=m;i++)
		for(j=1;j<=n;j++)
			if(x[i] == y[j]){
				length[i][j] = length[i-1][j-1] + 1;
				state[i][j] = 1;
			}
			else if(length[i][j-1] >= length[i-1][j]){
				length[i][j] = length[i][j-1];
				state[i][j] = 2;
			}
			else{
				length[i][j] = length[i-1][j];
				state[i][j] = 3;
			}
	i = m;j = n;k = length[m][n];//回溯
	while(i>0 && j>0){
		if(state[i][j] == 1){
			same[--k] = x[i];
			i--;j--;
		}
		else if(state[i][j] == 2)	j--;
		else i--;
	}
	return length[m][n];
}

测试如下：

int main(){
	char x[7] = {'0','a','b','c','b','d','b'};//从下标1开始存储
	char y[10] = {'0','a','c','b','b','a','b','d','b','b'};
	char same[10];
	int length = CommonOrder(x,7,y,10,same);
	cout<<"公共子序列长度为 "<<length<<endl;
	cout<<"公共子序列为 ";
	for(int i=0;i<length;i++)
		cout<<same[i]<<" ";
	return 0;
}

输出如下：

公共子序列长度为 5
公共子序列为 a c b d b