问题
给定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;
}