算法设计——最长公共子序列（动态规划）

问题

给定2个序列X和Y，当另一序列Z既是X的子序列又是Y的子序列时，称Z是 序列X和Y的公共子序列。
子序列的长度：子序列的元素个数。
问题描述：给定2个序列X={x1,x2,…,xm}和Y={y1,y2,…,yn}，求X和Y的最长公共子序列

分析

由此可见，2个序列的LCS包含了这2个序列前缀的LCS。因此， LCS问题具有最优子结构性质和子问题重叠性质。

因此，可是用动态规划算法
①当两序列（m、n个）最后一个元素不相等时，那么最长公共子序列为前m-1个和n-1个序列的LCS。
②当两序列（m、n个）最后一个元素相等时，那么最长公共子序列为前m-1个和n-1个序列的LCS+1。
状态转移方程为：

有了状态转移方程后，我们通过标记函数记录状态方程方向
可通过标记函数输出公共子序列，当b[i][j]=1时打印，递归后打印

时间复杂度为 O（mn）(m、n为序列元素个数)
空间复杂度O（mn）

代码实现

#include<stdio.h>
#include<string.h>
void fun(int i,int j,char str[],int flag[][101]){
	if(i==0||j==0) return ;
	if(flag[i][j]==1){
		fun(i-1,j-1,str,flag);
		printf("%c",str[i]);
	} 
	else if(flag[i][j]==2) fun(i-1,j,str,flag);
	else if(flag[i][j]==3) fun(i,j-1,str,flag);
}
int main()
{
	char str1[101],str2[101],ch;//默认str1长于str2 
	int flag[101][101],a[101][101];
	int len1,len2,i,j,k=1;
	while(scanf("%c",&ch)&&ch!='\n'){//下标从一开始 
		str1[k++]=ch;
	}
	len1=k;
	printf("%d",len1); 
	k=1;
	while(scanf("%c",&ch)&&ch!='\n'){//下标从一开始 
		str2[k++]=ch;
	}
	len2=k;
	for(i=0;i<len1;i++) a[i][0]=0;//0行0列置0 
	for(i=0;i<len2;i++) a[0][i]=0;
	for(i=1;i<len1;i++){
		for(j=1;j<len2;j++){//找标记和分情况 
			if(str1[i]==str2[j]){
				a[i][j]=a[i-1][j-1]+1;
				flag[i][j]=1;
			}else if(a[i-1][j]>=a[i][j-1]){
				a[i][j]=a[i-1][j];
				flag[i][j]=2;
			}else{
				a[i][j]=a[i][j-1];
				flag[i][j]=3;
			}
		}
	}
	fun(len1-1,len2-1,str1,flag);//这里长度要减一,因为第0行0列未使用 
	return 0;
}