算法设计北航童咏昕老师-最长公共子序列问题(动态规划)C语言实现

zyjyyds113

已于 2024-10-13 23:58:33 修改

阅读量257

点赞数 8

CC 4.0 BY-SA版权

分类专栏： c语言学习算法设计与分析学习文章标签：动态规划算法 c语言

于 2024-10-13 17:52:15 首次发布

本文链接：https://blog.csdn.net/m0_69550216/article/details/142901219

c语言学习同时被 2 个专栏收录

11 篇文章

订阅专栏

算法设计与分析学习

9 篇文章

订阅专栏

课程链接：算法设计与分析_中国大学MOOC(慕课) (icourse163.org)

课程链接：算法设计与分析_中国大学MOOC(慕课) (icourse163.org)

最长公共子序列-动态规划函数

/// <summary>
/// 最长公共子序列问题
/// </summary>
/// <param name="arr1"></序列arr1>
/// <param name="arr2"></序列arr2>
/// <returns></returns>
int LongestCommonSubsequence(char* arr1, char* arr2) {
	int arr1len = strlen(arr1);
	int arr2len = strlen(arr2);
	int C[MAX][MAX] = {0};//初始化二维数组
	for (int i = 0; i < arr1len + 1; i++)
	{
		for (int j = 0; j < arr2len + 1; j++)
		{
			C[i][j] = 0;
			rec[i][j] = 0;
		}
	}
	int i = 1;
	int j = 1;
	//依次计算子问题
	for (; i < arr1len + 1; i++)
	{
		for (j = 1; j < arr2len + 1; j++)
		{
			int ii = i - 1, jj = j - 1;//arr1、arr2的下标要从0开始
			if (arr1[ii] == arr2[jj])
			{
				C[i][j] = C[i - 1][j - 1] + 1;//记录长度和决策
				rec[i][j] = 1;
			}
			else if (C[i - 1][j] >= C[i][j - 1])
			{
				C[i][j] = C[i - 1][j];
				rec[i][j] = 2;
			}
			else
			{
				C[i][j] = C[i][j - 1];
				rec[i][j] = 3;
			}
		}
	}
	return C[i - 1][j - 1];
}

倒序追踪方案-递归函数

/// <summary>
/// 倒序追踪方案
/// </summary>
/// <param name="X"></序列arr1/arr2>
/// <param name="i"></当前位置i>
/// <param name="j"></当前位置j>
void PrintLCS(char* X, int i, int j)
{
	if (i == 0 || j == 0)
	{
		return;
	}
	if (rec[i][j] == 1)
	{
		PrintLCS(X, i-1, j-1);
		printf("%c", X[i-1]);//正序输出打印
	}
	if (rec[i][j] == 2)
	{
		PrintLCS(X, i - 1, j);
	}
	if (rec[i][j] == 3)
	{
		PrintLCS(X, i, j - 1);
	}
}

测试main函数

int main()
{
	char arr1[] = "ABCDBAB";
	char arr2[] = "DBCABA";
	int len=LongestCommonSubsequence(arr1,arr2);
	int arr1len = strlen(arr1);
	int arr2len = strlen(arr2);	
	printf("最长公共子序列长度是:%d\n", len);
	printf("序列是:");
	PrintLCS(arr1, arr1len, arr2len);
	return 0;
}