最长公共子序列动态规划c语言,《算法导论》动态规划—最长公共子序列(不连续)--c语言实现...

本文介绍了最长公共子序列的概念,以及如何使用动态规划方法在C语言中解决这个问题。通过一个递归解法和动态规划的自底向上计算,详细展示了算法的实现过程,包括构建二维数组来存储子序列长度和构造最长公共子序列的过程。最后,给出了C语言代码实现和算法的时间、空间复杂度分析。
摘要由CSDN通过智能技术生成

1、基本概念

一个给定序列的子序列就是该给定序列中去掉零个或者多个元素的序列。形式化来讲就是:给定一个序列X={x1,x2,……,xm},另外一个序列Z={z1、z2、……,zk},如果存在X的一个严格递增小标序列1,i2……,ik>,使得对所有j=1,2,……k,有xij =zj,则Z是X的子序列。例如:Z={B,C,D,B}是X={A,B,C,B,D,A,B}的一个子序列,相应的小标为<2,3,5,7>。从定义可以看出子序列直接的元素不一定是相邻的。

公共子序列:给定两个序列X和Y,如果Z既是X的一个子序列又是Y的一个子序列,则称序列Z是X和Y的公共子序列。例如:X={A,B,C,B,D,A,B},Y={B,D,C,A,B,A},则序列{B,C,A}是X和Y的一个公共子序列,但不不是最长公共子序列。因为它的长度等于3,而子序列{B,C,A,B}其长度等于4,所以序列{B,C,B,A}才是X和Y的一个最长公共子序列。

最长公共子序列(LCS)问题描述:给定两个序列X={x1,x2,……,xm}和Y={y1,y2,……,yn},找出X和Y的最长公共子序列。

2、动态规划解决过程

1)描述一个最长公共子序列

如果序列比较短,可以采用蛮力法枚举出X的所有子序列,然后检查是否是Y的子序列,并记录所发现的最长子序列。如果序列比较长,这种方法需要指数级时间,不切实际。

LCS的最优子结构定理:设X={x1,x2,……,xm}和Y={y1,y2,……,yn}为两个序列,并设Z={z1、z2、……,zk}为X和Y的任意一个LCS,则:

(1)如果xm=yn,那么zk=xm=yn,而且Zk-1是Xm-1和Yn-1的一个LCS。

(2)如果xm≠yn,那么zk≠xm蕴含Z是是Xm-1和Yn的一个LCS。

(3)如果xm≠yn,那么zk≠yn蕴含Z是是Xm和Yn-1的一个LCS。

定理说明两个序列的一个LCS也包含两个序列的前缀的一个LCS,即LCS问题具有最优子结构性质。

2)一个递归解

根据LCS的子结构可知,要找序列X和Y的LCS,根据xm与yn是否相等进行判断的,如果xm=yn则产生一个子问题,否则产生两个子问题。设C[i,j]为序列Xi和Yj的一个LCS的长度。如果i=0或者j=0,即一个序列的长度为0,则LCS的长度为0。LCS问题的最优子结构的递归式如下所示:

4ccdcb05c9836a789e54ff2ce491d1be.png

3)计算LCS的长度

采用动态规划自底向上计算解。书中给出了求解过程LCS_

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
最长公共子序列是指在两个序列中,找到一个最长的子序列,使得它在两个序列中的相对位置都相同。下面是使用动态规划算法求解最长公共子序列C语言代码: ```c #include <stdio.h> #include <string.h> #define MAXLEN 100 //最大序列长度 int dp[MAXLEN][MAXLEN]; //记录最长公共子序列的长度 int LCS(char *str1, char *str2, int len1, int len2) { int i, j; //初始化第一行和第一列 for (i = 0; i <= len1; i++) dp[i][0] = 0; for (j = 0; j <= len2; j++) dp[0][j] = 0; //动态规划求解 for (i = 1; i <= len1; i++) { for (j = 1; j <= len2; j++) { if (str1[i-1] == str2[j-1]) dp[i][j] = dp[i-1][j-1] + 1; else dp[i][j] = (dp[i][j-1] > dp[i-1][j]) ? dp[i][j-1] : dp[i-1][j]; } } //返回最长公共子序列的长度 return dp[len1][len2]; } int main() { char str1[MAXLEN], str2[MAXLEN]; int len1, len2, lcs_len; //读入两个字符串 printf("请输入第一个字符串:"); gets(str1); printf("请输入第二个字符串:"); gets(str2); //计算最长公共子序列的长度 len1 = strlen(str1); len2 = strlen(str2); lcs_len = LCS(str1, str2, len1, len2); //输出结果 printf("最长公共子序列的长度为:%d\n", lcs_len); return 0; } ``` 其中,dp[i][j]表示str1的前i个字符和str2的前j个字符的最长公共子序列的长度。在动态规划求解时,如果str1[i-1]等于str2[j-1],则当前字符可以加入最长公共子序列中,此时dp[i][j]等于dp[i-1][j-1]+1;否则,当前字符不能加入最长公共子序列中,此时dp[i][j]等于dp[i][j-1]和dp[i-1][j]中的较大值。最终,dp[len1][len2]即为最长公共子序列的长度。
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值