最长公共子序列输出路径

本文详细探讨了如何找到两个序列的最长公共子序列(LCS),并介绍了通过动态规划求解LCS的过程。内容包括输入输出说明、要点分析,特别是关于如何打印LCS的路径。此外,还指出在程序实现中可能出现的问题,如局部变量的使用,以及如何调整以正确地打印所有可能的LCS路径。
摘要由CSDN通过智能技术生成
输入输出
输入样例
2
7 6
A B C B D A B
B D C A B A
8 9
b a a b a b a b
a b a b b a b b a
输出
Case 1
4 LCS(X,Y):B C B A
Case 2
6 LCS(X,Y):a b a b a b
要点分析

设序列X = {x1,x2,…,xm}和Y = {y1,y2,…,yn}的最长公共子序列为Z = {z1,z2,…,zk} ,则

  1. 若xm = yn,则zk = xm = yn,且zk-1是xm-1和yn-1的最长公共子序列
  2. 若xm≠yn且zk≠xm,则Z是xm-1和Y的最长公共子序列
  3. 若xm≠yn且zk≠yn,则Z是X和yn-1的最长公共子序列
      由此可见,2个序列的最长公共子序列包含了这2个序列的前缀的最长公共子序列。用C[i][j]记录序列和的最长公共子序列的长度,Xi = {x1,x2,…,xi};Yj = {y1,y2,…,yj}。当i = 0或j = 0时,空序列是Xi和Yj的最长公共子序列,此时C[i][j] = 0。
      由最优子结构性质可建立关系如下:
    公式
打印路径
而在打印所有的路径的时候注意:

  • 0
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
以下是C语言求解最长公共子序列的伪代码: 1. 首先定义一个函数Lcs_length,该函数接受两个字符串str1和str2以及两个指向指针的指针c和b作为参数,其中c和b用于存储最长公共子序列的长度和路径信息。 2. 在Lcs_length函数中,首先获取str1和str2的长度,然后创建一个二维数组c,用于存储最长公共子序列的长度,以及一个二维数组b,用于存储路径信息。 3. 然后使用两个for循环遍历str1和str2的所有字符,如果当前字符相同,则将c[i][j]设置为c[i-1][j-1]+1,并将b[i][j]设置为左上角。 4. 如果当前字符不同,则将c[i][j]设置为max(c[i-1][j], c[i][j-1]),并将b[i][j]设置为c[i-1][j] > c[i][j-1] ? 上方 : 左方。 5. 最后返回c和b。 6. 定义一个函数Find_Lcs,该函数接受两个字符串str1和str2作为参数。 7. 在Find_Lcs函数中,首先调用Lcs_length函数获取最长公共子序列的长度和路径信息。 8. 然后根据路径信息,倒序遍历c数组,将最长公共子序列存储到一个字符数组中。 9. 最后输出最长公共子序列的长度和内容。 ```c // 定义Lcs_length函数 int Lcs_length(char *str1, char *str2, int **c, int **b) { int len1 = strlen(str1); int len2 = strlen(str2); int i, j; for (i = 0; i <= len1; i++) { c[i][0] = 0; } for (j = 0; j <= len2; j++) { c[0][j] = 0; } for (i = 1; i <= len1; i++) { for (j = 1; j <= len2; j++) { if (str1[i-1] == str2[j-1]) { c[i][j] = c[i-1][j-1] + 1; b[i][j] = 0; } else if (c[i-1][j] >= c[i][j-1]) { c[i][j] = c[i-1][j]; b[i][j] = 1; } else { c[i][j] = c[i][j-1]; b[i][j] = -1; } } } return c[len1][len2]; } // 定义Find_Lcs函数 void Find_Lcs(char *str1, char *str2) { int len1 = strlen(str1); int len2 = strlen(str2); int **c = (int **)malloc((len1+1) * sizeof(int *)); int **b = (int **)malloc((len1+1) * sizeof(int *)); int i, j; for (i = 0; i <= len1; i++) { c[i] = (int *)malloc((len2+1) * sizeof(int)); b[i] = (int *)malloc((len2+1) * sizeof(int)); } int len = Lcs_length(str1, str2, c, b); char *lcs = (char *)malloc((len+1) * sizeof(char)); lcs[len] = '\0'; i = len1; j = len2; while (i > 0 && j > 0) { if (b[i][j] == 0) { lcs[--len] = str1[i-1]; i--; j--; } else if (b[i][j] == 1) { i--; } else { j--; } } printf("最长公共子序列长度为:%d\n", len); printf("最长公共子序列为:%s\n", lcs); free(lcs); for (i = 0; i <= len1; i++) { free(c[i]); free(b[i]); } free(c); free(b); } ```
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值