动态规划解决方案最长公共子序列问题(开启)

动态规划解决方案最长公共子序列问题(开启)

大家好,又见面了,我是全栈君,今天给大家准备了Idea注册码。



动态规划

常常会遇到复杂问题不能简单地分解成几个子问题,而会分解出一系列的子问题。简单地採用把大问题分解成子问题。并综合子问题的解导出大问题的解的方法,问题求解耗时会按问题规模呈幂级数添加。

为了节约反复求同样子问题的时间,引入一个数组。无论它们是否对终于解实用。把全部子问题的解存于该数组中,这就是动态规划法所採用的基本方法。

【问题】 求两字符序列的最长公共字符子序列

问题描写叙述:字符序列的子序列是指从给定字符序列中任意地(不一定连续)去掉若干个字符(可能一个也不去掉)后所形成的字符序列。

令给定的字符序列X=“x0,x1,…。xm-1,序列Y=“y0,y1。…。yk-1X的子序列。存在X的一个严格递增下标序列<i0,i1,…。ik-1>,使得对全部的j=0,1,…,k-1,有xij=yj

比如,X=“ABCBDAB”,Y=“BCDB”是X的一个子序列。

考虑最长公共子序列问题怎样分解成子问题,设A=“a0。a1。…,am-1B=“b0,b1,…。bm-1,并Z=“z0。z1,…,zk-1为它们的最长公共子序列。不难证明有下面性质:

(1假设am-1=bn-1,则zk-1=am-1=bn-1,且“z0。z1,…。zk-2是“a0,a1。…,am-2和“b0,b1。…,bn-2的一个最长公共子序列;

(2假设am-1!=bn-1,则若zk-1!=am-1,蕴涵“z0。z1,…。zk-1是“a0。a1,…。am-2和“b0,b1,…,bn-1的一个最长公共子序列;

(3假设am-1!=bn-1,则若zk-1!=bn-1,蕴涵“z0,z1。…,zk-1是“a0。a1,…,am-1和“b0。b1。…,bn-2的一个最长公共子序列。

这样,在找A和B的公共子序列时,如有am-1=bn-1。则进一步解决一个子问题,找“a0,a1。…,am-2和“b0。b1,…,bm-2的一个最长公共子序列。假设am-1!=bn-1,则要解决两个子问题,找出“a0。a1。…,am-2和“b0,b1,…,bn-1的一个最长公共子序列和找出“a0,a1,…,am-1和“b0,b1。…。bn-2的一个最长公共子序列,再取两者中较长者作为A和B的最长公共子序列。

 

 

求解:

引进一个二维数组c[][],用c[i][j]记录X[i]与Y[j] 的LCS 的长度。b[i][j]记录c[i][j]是通过哪一个子问题的值求得的,以决定搜索的方向。

我们是自底向上进行递推计算,那么在计算c[i,j]之前,c[i-1][j-1],c[i-1][j]与c[i][j-1]均已计算出来。

此时我们依据X[i] = Y[j]还是X[i] != Y[j]。就能够计算出c[i][j]。

问题的递归式写成:

动态规划解决方案最长公共子序列问题(开启)

回溯输出最长公共子序列过程:

动态规划解决方案最长公共子序列问题(开启)

 

算法分析:
因为每次调用至少向上或向左(或向上向左同一时候)移动一步,故最多调用(m + n)次就会遇到i = 0或j = 0的情况,此时開始返回。

返回时与递归调用时方向相反,步数同样,故算法时间复杂度为Θ(m + n)。

 

 

代码:

 

#include <stdio.h>
#include <string.h>
#define MAXLEN 100

void LCSLength(char *x, char *y, int m, int n, int c[][MAXLEN], int b[][MAXLEN])

{
    int i, j;
    
    for(i = 0; i <= m; i++)
        c[i][0] = 0;
    for(j = 1; j <= n; j++)
        c[0][j] = 0;
    for(i = 1; i<= m; i++)
    {
        for(j = 1; j <= n; j++)
        {
            if(x[i-1] == y[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;
            }
        }
    }
}

void PrintLCS(int b[][MAXLEN], char *x, int i, int j)
{
    if(i == 0 || j == 0)
        return;
    if(b[i][j] == 0)
    {
        PrintLCS(b, x, i-1, j-1);
        printf("%c ", x[i-1]);
    }
    else if(b[i][j] == 1)
        PrintLCS(b, x, i-1, j);
    else
        PrintLCS(b, x, i, j-1);
}

int main(int argc, char **argv)
{
    char x[MAXLEN] = ...{"ABCBDAB"};
    char y[MAXLEN] = ...{"BDCABA"};
    int b[MAXLEN][MAXLEN];
    int c[MAXLEN][MAXLEN];
    int m, n;
    
    m = strlen(x);
    n = strlen(y);
    
    LCSLength(x, y, m, n, c, b);
    PrintLCS(b, x, m, n);
    
    return 0;
}


认为思路很清晰,所以转载。

感谢作者:http://blog.csdn.net/yysdsyl/article/details/4226630

版权声明:欢迎follow我的开源项目:https://github.com/skyqinsc

版权声明:本文内容由互联网用户自发贡献,该文观点仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 举报,一经查实,本站将立刻删除。

发布者:全栈程序员-用户IM,转载请注明出处:https://javaforall.cn/116852.html原文链接:https://javaforall.cn

【正版授权,激活自己账号】: Jetbrains全家桶Ide使用,1年售后保障,每天仅需1毛

【官方授权 正版激活】: 官方授权 正版激活 支持Jetbrains家族下所有IDE 使用个人JB账号...

(0)
blank

相关推荐

  • 第三届“SLAM技术及应用”暑期学校视频分享(完整版)

    第三届“SLAM技术及应用”暑期学校视频分享(完整版)B站录播视频,up主上传资源:视频链接:https://www.bilibili.com/video/BV1Nk4y127qX?p=15所有课件PPT下载:http://www.cad.zju.edu.cn/home/gfzhang/download/2020-SLAM-Summer-School-slides.zip

  • 令人心动的WEB HOSTING 虚拟主机服务[通俗易懂]

    令人心动的WEB HOSTING 虚拟主机服务[通俗易懂]你正在为网站寻找虚拟主机(WEBHOSTING)吗?哪里能找到功能卓越,价格低廉的ISP服务商呢?跟着我,你便会找到惊喜.一个偶然的机会使我有机会访问到国外的网站服务公司HostExcellence. ,它的服务套餐有以下几种LinuxHostingPlans:StarterPlan-1Domain&50mb/$2.95InstantPlan-1Domain&

  • 思科交换机AAA认证

    思科交换机AAA认证aaanew-modelusernameadminprivilege15passwordadminaaaauthenticationlogindefaultgrouptacacs+localaaaauthorizationexecdefaultgrouptacacs+localaaaaccountingexecdefaultstar…

  • 二叉树层序遍历(C语言)[通俗易懂]

    二叉树层序遍历(C语言)[通俗易懂]二叉树的层序遍历即从上到下,在每一层从左到右依次打印数据。如下:层序遍历结果:ABCDEFG基本思路即将根节点入队后,之后每次都将队首元素出队,打印队首元素数据,并将队首元素左右子树入队,一直重复上述过程。自然,本题还可以用数组来实现。代码:#include<stdio.h>#include<stdlib.h>#defineQueueMax100typedefstructNode{chardata;structNode*

  • 扫描主机漏洞的工具_漏洞扫描工具有哪些

    扫描主机漏洞的工具_漏洞扫描工具有哪些0x00说明:这是一款基于主机的漏洞扫描工具,采用多线程确保可以快速的请求数据,采用线程锁可以在向sqlite数据库中写入数据避免databaseislocked的错误,采用md5哈希算法确保数据不重复插入。本工具查找是否有公开exp的网站为shodan,该网站限制网络发包的速度,因而采用了单线程的方式,且耗时较长。功能:查找主机上具有的CVE查找具有公开EXP的CVE0x…

  • TerminateThread可能导致程序崩溃

    TerminateThread可能导致程序崩溃DWORD__stdcallmythread(void*){   while(true)   {      char*p=newchar[1024];      deletep;   }}int_tmain(intargc,_TCHAR*argv[]){   HANDLEh=CreateThread(NULL,0,mythread,N…

    2022年10月30日

发表回复

您的电子邮箱地址不会被公开。

关注全栈程序员社区公众号