最长公共子序列

最新推荐文章于 2024-09-19 08:38:15 发布

acm_JL

最新推荐文章于 2024-09-19 08:38:15 发布

阅读量1k

点赞数

分类专栏：动态规划

本文链接：https://blog.csdn.net/acm_JL/article/details/50926606

版权

动态规划专栏收录该内容

9 篇文章 0 订阅

订阅专栏

#include <stdio.h>
#define NUM 100
int c[NUM][NUM];
int b[NUM][NUM];
void LCSLength(int m, int n, const char x[],char y[])//求解最长公共子序列的长度的最优解值c[m][n] 
{  
	int i,j;
	for (i = 1; i <= m; i++) c[i][0] = 0;//将第一行第一列都设置成0 
	for (i = 1; i <= n; i++) c[0][i] = 0;
	
	for (i = 1; i <= m; i++)
	for (j = 1; j <= n; j++) 
	{
		if (x[i]==y[j]) //如果字符相同就是斜向上方向加1作为结果 
			{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 LCS(int i,int j,char x[])//求解最长公共子序列长度的最优解 
{
	if (i ==0 || j==0) return;
	if (b[i][j]== 1){ LCS(i-1,j-1,x); printf("%c",x[i]); }
	else if (b[i][j]== 2) LCS(i-1,j,x);
	else LCS(i,j-1,x);
}

int main() 
{
	char x[NUM];
	char y[NUM];
	int m, n;
	scanf("%d\n", &m);
	for (int i=1; i<=m; i++)
		scanf("%c", &x[i]);
	scanf("%d\n", &n);
	for (int i=1; i<=n; i++)
		scanf("%c", &y[i]);
	LCSLength(m, n, x, y);
	printf("%d\n", c[m][n]);
	LCS(m,n,x);
	return 0;
}