动态规划算法学习四:最大上升子序列问题(LIS:Longest Increasing Subsequence)

news2024/9/20 20:49:57

文章目录

  • 前言
  • 一、问题描述
  • 二、DP步骤
    • 1、最优子结构
      • a、限界上升子序列
      • b、最优子结构性质
    • 2、状态表示和递推方程
    • 3、计算最优值
    • 4、算法实现
  • 三、优化:非DP /二分法
    • 1、新问题
    • 2、算法实现

前言

一、问题描述

在这里插入图片描述
在这里插入图片描述

二、DP步骤

1、最优子结构

  • 给定序列𝑆=[𝑠1,𝑠2,⋯,𝑠𝑛],如果子序列[𝑆(𝑖1 ),𝑆(𝑖2 ), ⋯,𝑆(𝑖𝑘 )]是其最大上升子序列,则[𝑆(𝑖1 ),𝑆(𝑖2 ), ⋯,𝑆(𝑖(𝑘−1) )]是子问题𝑆=[𝑠1,𝑠2,⋯,𝑆(𝑖(𝑘−1) )]的最大上升子序列吗?
  • 例:给定 S = [1, 3, 4, 2, 7, 9, 6, 8],最大上升子序列可为 [1, 3, 4, 6, 8],最大长度为5。
  • [1, 3, 4, 6]显然不是[1, 3, 4, 2, 7, 9, 6]的最大上升子序列,因为[1, 3, 4, 7, 9]是其最大上升子序列。
    注意:可以同时有多个最大上升子序列!每个的最大值都不相同
  • 如果直接把最大上升子序列的长度作为规划目标,那么该问题不具备最优子结构性质(因为可能同时有多个最大上升子序列,使得最优子结构不成立)。只能采用间接的办法,引入一些中间目标作为动态规划的对象。

a、限界上升子序列

在这里插入图片描述
在这里插入图片描述

b、最优子结构性质

在这里插入图片描述
在这里插入图片描述

2、状态表示和递推方程

在这里插入图片描述
在这里插入图片描述
公式讲解:只求长度,不求具体的子序列

  • 逐一找出1—(m-1)个元素中所有小于Sm的元素Si及其对应的Len(Si)值,然后从中取最大的一个Len(.)值,然后加1。【 1–i个元素中所有大于Sm的元素的Len(.)值全都不予考虑】

  • 如果精简为max{Len(i)+1 | 1≤i<m, Sm>Si},可以吗?
    sure

3、计算最优值

  • 𝐿𝑒𝑛(𝑚)的计算同样按照自底向上的顺序进行:所有子问题的总数为n个(即,分别以每个元素为限界);m越小,其子问题𝐿𝑒𝑛(𝑚)求解的复杂度就越低
  • 当所有的𝐿𝑒𝑛(𝑚)(1≤𝑚≤𝑛)都计算完毕后,统计其中最大值,即可得到输入序列 S 的最大上升子序列
  • 注意:最大上升子序列的长度不一定等于𝑳𝒆𝒏(𝒏)
    在这里插入图片描述
    答案详解:
  1. Len(1) = 1
  2. Len(3) = 2:找到3之前小于3的数字中,最大的Len(.)值加1,即:len(1)+1=1。
  3. Len(4)=3:找到4之前小于4的数字:1和3,最大Len(.)值加1,即:len(3)+1=3。
  4. Len(2)=2:找到2之前小于2的数字:1,最大Len(.)值加1,即:Len(1)+1=2。
  5. Len(7)=4:找到7之前小于7的数字:1,3,4,2,最大Len(.)值加1,即:Len(4)+1=4。
  6. Len(6)=4:找到6之前小于6的数字:1,3,4,2,最大Len(.)值加1,即:Len(4)+1=4。
  7. Len(8)=5:找到8之前小于8的数字:1,2,3,2,7,6,最大Len(.)值加1,即:Len(7)+1=5 或者 Len(6)+1=5

4、算法实现

时间复杂度:n2

/**
 * DP算法之最大上升子序列问题
 */
public class Main4 {
    public static int MAXN = 100;
    public static void main(String[] args) {
        int[] seqSrc = {1, 3, 4, 2, 7, 6, 8};
        int i = LISLength(7, seqSrc);
        System.out.println(i);
    }

    public static int LISLength(int num, int[] seqSrc) {
        int[] Len = new int[MAXN];
        int res = 1;
        //设第m个数的值为上界
        for (int m = 0; m < num; m++) {  
        	//每个新m为上界时,Len[m]总是从1开始
            Len[m] = 1;     
            for (int i = 0; i < m; i++)
            	//遍历所有以第i个数为上界的长度,从中选出符合max公式条件的最大值加1,就是Len[m]。
                if (seqSrc[i] < seqSrc[m] && Len[i] + 1 > Len[m])   
                    Len[m] = Len[i] + 1;
            //记录从len[1]到len[m]中的最大值,res为最大公共子序列长度,最后的解
            res = (res > Len[m] ? res : Len[m]);  
        }
        return res;
    }
}

在这里插入图片描述
代码解说:

  1. 第一个for遍历:有几个元素就遍历几个,计算每个的Len(.)值。
  2. 第二个for遍历:遍历当前元素的最大Len(.)值。
    • seqSrc[i] < seqSrc[m] :过滤小于当前元素的值,
    • Len[i] + 1 > Len[m]:找到最大值

三、优化:非DP /二分法

时间复杂度 nlog2n。

1、新问题

在这里插入图片描述
上图中的公式中的 i ,为数组的下标
例:

  1. k = 4 的值有两个,Len(5) = 4,Len(6) = 4
  2. S[5] = 7, S[6] = 6
  3. 所以 T[k=4] = 6。
    在这里插入图片描述
    在这里插入图片描述
    在这里插入图片描述

2、算法实现

/**
 * DP之最大上升子序列:时间复杂度  nlog2n
 */
public class Main5 {

    public static void main(String[] args) {
        int[] seqSrc = {1, 3, 4, 2, 7, 6, 8};
        int i = lengthOfLIS(seqSrc);
        System.out.println(i);
    }

    public static int lengthOfLIS(int[] arr) {
        if (arr == null || arr.length == 0) {
            return 0;
        }
        int[] ends = new int[arr.length];
        ends[0] = arr[0];
        int right = 0;
        int l = 0;
        int r = 0;
        int m = 0;
        int max = 1;
        for (int i = 1; i < arr.length; i++) {
            l = 0;
            r = right;
            while (l <= r) {
                m = (l + r) / 2;
                if (arr[i] > ends[m]) {
                    l = m + 1;
                } else {
                    r = m - 1;
                }
            }
            right = Math.max(right, l);
            ends[l] = arr[i];
            max = Math.max(max, l + 1);
        }
        return max;
    }
}

在这里插入图片描述

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.coloradmin.cn/o/31107.html

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈,一经查实,立即删除!

相关文章

“300万”只是新起点,比亚迪将开启下一个 “黄金周期”

比亚迪再次创造全球新能源汽车市场新标杆。 11月16日&#xff0c;比亚迪第300万辆新能源汽车正式下线。成为首个达成这一里程碑的中国品牌。 正如比亚迪股份有限公司董事长兼总裁王传福说&#xff0c;从“第1辆新能源汽车到第100万辆新能源汽车”用时13年、从“100万辆到200万…

MCE | 磁珠 Protocol,如何快速捕获您心仪的蛋白~

磁珠的优势 ◎ 蛋白荷载量高 ◎ 特异性强、非特异性结合性低 ◎ 样品损失小 ◎ 操作方便 如何操作 ■ 磁珠预处理 将磁珠充分混悬&#xff0c;取 25-50 μL 磁珠&#xff0c;置于 1.5 mL EP 管中&#xff0c;加入 400 μL 结合/洗涤缓冲液&#xff0c;充分混悬&#xff0c;置…

Word处理控件Aspose.Words功能演示:使用 Java 将文本转换为 PDF

TXT格式的文本文档包含行形式的纯文本。TXT 文件是存储没有任何格式的纯文本的最简单和最简单的方法。我们可以在任何文本编辑器或文字处理应用程序中轻松创建、打开和编辑 TXT 文件。在某些情况下&#xff0c;我们可能需要将文本转换为只读格式&#xff0c;例如PDF。在本文中&…

flutter 怎么消除按钮事件的点击溅射背景

flutter 怎么消除按钮事件的点击溅射背景前言一、设置 ThemeData二、Theme 设置三、单独设置总结前言 在flutter 中&#xff0c;大部分事件组件都有一个溅射背影&#xff0c;但是假如某天需求让我们取消点击溅射效果&#xff0c;我们该怎么办呢&#xff1f;本篇文章将记录怎么…

画法几何及机械制图复习题及答案

机 械 制 图复习题及参考答案 一、填空题 1&#xff0e;比例的种类有 、 、 。 2&#xff0e;图样中的可见轮廓线用 绘制&#xff1b;图样中尺寸线和尺寸界线用 绘制。 3&#xff0e;正投影的基本性质包括 、 、 。 4&#xff0e;三视图的投影关系表现为&#xff1a;主、俯…

【ML特征工程】第 1 章 :机器学习管道

&#x1f50e;大家好&#xff0c;我是Sonhhxg_柒&#xff0c;希望你看完之后&#xff0c;能对你有所帮助&#xff0c;不足请指正&#xff01;共同学习交流&#x1f50e; &#x1f4dd;个人主页&#xff0d;Sonhhxg_柒的博客_CSDN博客 &#x1f4c3; &#x1f381;欢迎各位→点赞…

数字逻辑·时序线路分析【常见的时序线路】

这一篇和之前那一篇讲的是时序线路 之前学过的是组合线路 寄存器 有3个D触发器控制 C1 − C3 用来寄存二进制代码。 下面的与或非门用来接收要寄存的二进制代码。 上面的与非门用来发送寄存的二进制代码。 输入信号&#xff1a; RD&#xff1a;清除信号。 WAC&#xff1a;直送…

spring-cloud-dubbo基本使用

创建模块 api模块使用mave quick-start 构建: &#xff0c;provider模块使用 下面方式创建&#xff1a; 点击下一步&#xff0c;会看到一些基于阿里的cloud的依赖&#xff1a; 上面这个是基于阿里云的&#xff0c;下面的Spring Cloud Alibaba是开源的组件依赖&#xff1a; …

torch.as_tensor()、torch.Tensor() 、 torch.tensor() 、transforms.ToTensor()的区别

1&#xff09;torch.as_tensor(data, dtypeNone,deviceNone)->Tensor : 为data生成tensor&#xff0c;保留 autograd 历史记录并尽量避免复制&#xff08;dtype和devices相同&#xff0c;尽量浅拷贝&#xff09;。 如果data已经是tensor&#xff0c;且dtype和device与参数相…

基于复杂环境下的雷达目标检测技术(Matlab代码实现)

&#x1f352;&#x1f352;&#x1f352;欢迎关注&#x1f308;&#x1f308;&#x1f308; &#x1f4dd;个人主页&#xff1a;我爱Matlab &#x1f44d;点赞➕评论➕收藏 养成习惯&#xff08;一键三连&#xff09;&#x1f33b;&#x1f33b;&#x1f33b; &#x1f34c;希…

轻量级模型设计与部署总结

前言一些关键字定义及理解 计算量 FLOPs内存访问代价 MACGPU 内存带宽Latency and Throughput英伟达 GPU 架构 CNN 架构的理解手动设计高效 CNN 架构建议 一些结论&#xff1a; 一些建议轻量级网络模型部署总结轻量级网络论文解析文章 参考资料 文章同步发于 github 仓库 和 知…

论文阅读11——《Mutual Boost Network for Attributed Graph Clustering》

原文地址&#xff1a; 论文阅读11——《Mutual Boost Network for Attributed Graph Clustering》 作者&#xff1a;Xiaoqiang Yan, Xiangyu Yu, Shizhe Hu, Yangdong Ye 发表时间&#xff1a;预印本 论文地址&#xff1a;https://papers.ssrn.com/sol3/papers.cfm?abstract_i…

Camtasia Studio2023喀秋莎免费实用的屏幕录像工具

CamtasiaStudio2023是一款非常不错的软件。总的来说CamtasiaStudio的功能从专业度来说&#xff0c;分别有&#xff1a;录制桌面&#xff0c;录制视频教程&#xff0c;录制音频&#xff1b;剪截视频&#xff0c;拼接合成视频&#xff0c;制作小视频Camtasia Studio是TechSmith的…

web网页设计期末课程大作业:环境保护主题网站设计——农业三级带表单带js(14页)HTML+CSS+JavaScript

&#x1f380; 精彩专栏推荐&#x1f447;&#x1f3fb;&#x1f447;&#x1f3fb;&#x1f447;&#x1f3fb; ✍️ 作者简介: 一个热爱把逻辑思维转变为代码的技术博主 &#x1f482; 作者主页: 【主页——&#x1f680;获取更多优质源码】 &#x1f393; web前端期末大作业…

自顶向下计算机网络学习 传输层

自顶向下计算机网络学习 传输层一、概述和传输层服务1.1 传输层提供的服务1.2传输层和网络层的关系1.3 传输层协议可靠的、保序的传输&#xff1a;TCP不可靠、不保序的传输&#xff1a;UDP二、多路复用与解复用2.1 什么是复用与解复用2.2 面向连接(TCP)的多路复用与分解2.3 无连…

Verilog 实现无毛刺时钟切换电路,RTL代码设计+testbench代码测试,波形前仿真

Verilog 实现无毛刺时钟切换电路 1,原理2,无毛刺时钟切换3,RTL代码设计4,testbench测试代码5,RTL+testbench综合的 Netlist6,前仿真波形验证参考文献1 1,原理 想要切换时钟电路,最简单的方法肯定是使用一个MUX,control作为控制信号; control = 1, clk_output = clk_…

JVM 双亲委派模型

一言以蔽之&#xff0c;向上委托&#xff0c;向下委派。 向上委托&#xff1a;如果一个类加载器收到了类加载请求&#xff0c;它并不会自己先去加载&#xff0c;而是把这个请求委托给父类的加载器去执行&#xff0c;如果父类加载器还存在其父类加载器&#xff0c;则进一步向上…

家居建材如何在线管理订单?数商云采购系统实现订单发收货、退换货流程化管控

新时代的流程化业务&#xff0c;利用信息互通的时效性&#xff0c;提高员工执行力&#xff0c;提升市场竞争力&#xff0c;在悄无声息地促进企业发展。订单管理作为企业客户关系管理的有效延伸&#xff0c;能更好的把个性化、差异化服务有机的融入到客户管理中去&#xff0c;能…

220kV降压变电所电气部分初步设计33号

目录 1 绪论 5 &#xff11;.&#xff11; 概述 5 &#xff11;.&#xff12; 本次设计内容 5 &#xff11;.&#xff13; 本次设计任务 5 2 变电站总体分析 6 2.&#xff11;变电所总体分析 6 2.2主变压器选择 7 3 电气主接线选择 10 3.1 电气主接线的概念 10 3.2 电气主接…

Linux —— 进程控制

1.进程控制的四个概念 进程控制分为四类&#xff0c;分别是&#xff1a; 进程创建进程终止进程等待进程替换 2.进程创建 2.1初识fork fork的作用是通过拷贝当前进程创建一个子进程&#xff0c;这两个进程的区别在于PID不同(还有一些资源、统计量也不同&#xff0c;但PID是我…