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

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

大家好,又见面了,我是全栈君,今天给大家准备了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)


相关推荐

  • 码农盖房记,纯图片

    码农盖房记,纯图片

  • pycharm如何配置编译环境_python不配置环境变量会怎样

    pycharm如何配置编译环境_python不配置环境变量会怎样随便打开一个.py文件时,右上角三角形运行按钮不能选中,需要配置编译环境配置编译环境有两个部分:1、添加编译器(interpreter)位置;2、添加脚本(所打开的.py文件)位置,具体步骤如下:①pycharm左上角依次点击File→Settings…→ProjectInterpreter点击右上的小齿轮→Add,选择下面的Existingenvironment选项,再在目录中寻找带python.exe的文件夹(本文最下面贴了如何寻找这个文件位置),然后点OK这里可以看到所有已安装的第

  • 15款android设备上的代码编辑器

    15款android设备上的代码编辑器如果你希望你的Android设备,如智能手机和平板电脑,在任何时间和任何地方都能够编写代码,那么,不妨看看下面我将介绍的15款Android代码编辑器,它们必将成为你的理想工具。1.DeuterIDE一款基于移动开发的文本编辑器,在线编译器,内置支持Android的集成开发环境。它的虚拟按键可提供常用的编程符号和语法高亮功能,并支持超过40种编程语言,从而能让你在旅途中写代码变得从未如此简单。

  • Android应用程序启动过程源代码分析(4)

    Android应用程序启动过程源代码分析(4)

  • redisson读写锁使用场景_Redisson酒店

    redisson读写锁使用场景_Redisson酒店读写锁一次只有一个线程可以占有写模式的读写锁,但是可以有多个线程同时占有读模式的读写锁.正是因为这个特性,当读写锁是写加锁状态时,在这个锁被解锁之前,所有试图对这个锁加锁的线程都会被阻塞.当读写锁在读加锁状态时,所有试图以读模式对它进行加锁的线程都可以得到访问权,但是如果线程希望以写模式对此锁进行加锁,它必须直到所有的线程释放锁.通常,当读写锁处于读模式锁住状态时,如果有另外线程试图以写模式加锁,读写锁通常会阻塞随后的读模式锁请求,这样可以避免读模式锁长期占用,而等待的写模式

  • Java哈希表以及哈希冲突

    Java哈希表以及哈希冲突文章目录Java哈希表概念冲突避免冲突哈希函数的设计方法常见哈希函数负载因子调节解决哈希冲突两种常见的方法是:闭散列和开散列哈希表和java类集的关系Java哈希表概念顺序结构以及平衡树中,元素关键码与其存储位置之间没有对应的关系,因此在查找一个元素时,必须要经过关键码的多次比较。顺序查找时间复杂度为O(N),平衡树中为树的高度,即O(logN),搜索的效率取决于搜索过程中元素的比较次…

发表回复

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

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