最大子矩阵:前缀和、动态规划

最近在学习动态规划,在牛客上刷题时碰到了这一题。其实最初的想法是暴力和前缀和,但是时间复杂度极高,需要套4层循环。后来去网上搜了一下相关的题解和做法,进而了解到了前缀和+线性动态规划的做法。但是在成功做出这题之前,个人感觉所搜到的博主讲解偏向于代码的编写,对于我这种初入算法的小白来说还是蛮费力气的,所以本节内容我将和大家一起从算法原理到代码一一剖析,争取写出清晰且容易理解的算法思路。

题目链接:最大子矩阵_牛客题霸_牛客网 (nowcoder.com)

 

首先我们来明确一下题目要求:本题要求的是我们在一个给定的N x N 的“矩阵”中找到 “元素和”最大的“子矩阵”。这个N x N 的矩阵其实就是一个二维数组,我们把它理解为一个矩阵。那么它的子矩阵就是矩阵中的某一片矩阵型的区域,那所要求的“最大子矩阵”的自然是找到在所有子矩阵中,矩阵元素和最大的那个矩阵。

我们先来看示例1来帮助我们理解一下题意:

我们可以发现,除该子矩阵外,在其他任意地方随意圈出一个子矩阵中的元素和均比图示矩阵的元素和小,此时我们就找到了该矩阵的最大子矩阵。

以下是两种解法:

 

 

说实话,第一种纯暴力解法就不用看了,数据量稍微大一点,就超时了。

第二种使用二维数组前缀和对原始数组进行了预处理,得到了二维前缀和数组,之后在我们求圈定范围的子矩阵元素和时会带来不小便利,但是4层循环O(n^4)的时间复杂度仍然不可小觑。 

那么我们有什么方法能减少时间复杂度呢?

我们先来了解一维数组的前缀和:

一维数组的前缀和是指将数组中从起始位置到当前位置的所有元素相加的结果,即上图中的prefix数组。

我们再来复习一下求一维数组最大子序列所使用到的线性dp算法:

算法原理:

状态表示:dp[i]所表示的是以i结尾的所有子数组中元素的最大和。 

状态转移方程:用来更新动态规划数组dp[i]:

  • nums[i - 1]表示当前位置i对应的原始数组 nums 中的元素值。

  • dp[i - 1] + nums[i - 1] 表示从当前位置向前延伸的子序列的和,即以 nums[i - 1] 结尾的子序列和加上当前位置的元素值。这个值表示了当前位置开始的新子序列的和。

  • max(nums[i - 1], dp[i - 1] + nums[i - 1]) 选择了两种情况中的较大值:第一种情况是只包含当前元素 nums[i - 1],即以当前位置 i 结尾的子序列。第二种情况是将当前位置的元素加入到之前的子序列中,即从当前位置开始新的子序列。

  • 将较大值更新到 dp[i],表示以 nums[i - 1] 结尾的最大子序列和。

这样,通过动态规划数组 dp 的更新,每个位置 i 都计算出了以该位置结尾的最大子序列和,最终找到整个数组的最大子序列和。

由于最后要输出的是dp数组中的最大值,我们可以使用一个变量来记录以 i 位置为结尾的最大子序列和,这样空间复杂度就从O(N)减少到了O(1)。

 其实,对于我们这一题,也可以使用一维数组前缀和与线性dp来解决,我们只需要将二维数组“压缩”为一维数组就好了。那么压缩方式自然是使用前缀和了。

怎么压缩呢?压缩后如何使用前缀和搭配线性dp解题呢?我们接着往下看:

进行完这一步操作后,第 i 行中第 j 列的元素即为:第 j 列从起始位置到当前位置 i 的所有元素相加的结果。 

其实通过上面的例子,我们不难发现进行预处理后的前缀和数组,仍然可以表示原数组的元素,更方便的是对于求原数组圈定矩阵元素和,在处理后的数组中仅通过使用末行数组元素减去起始行前一行的数组元素就可以得到所求矩阵的元素和。

通过不同的末行对起始行的减法操作,我们最后可以得到如下序列:

 通过前后两张图的解析,其实我们不难发现在最终生成的任意一个子序列中,随机取一段连续的数字即可表示原数组中的子矩阵。即原矩阵中所有的子矩阵均可由生成的子序列得到。

大家一定要好好理解并验证上面的两张图和两段话,当理解通透了才便于进行后续代码的编写。

既然 原矩阵中所有的子矩阵均可由生成的子序列得到,那么最大子矩阵必定是所有子序列数组中的最大连续子序列(注意:是上面通过两层循环得到的10个子序列数组中最大连续子序列元素和)。

那么我们接下来的目标很简单,就是对10个子序列数组中的每一个进行对最大连续子序列元素和的求解。所以我们需要再加上一层循环,用来遍历子序列数组。并且我们需要一个变量 ans 来记录每个子序列数组中的最大子序列元素和,需要注意的是 ans 在每次进入循环之前要更新为0,防止对后续最大子序列的求解造成影响。同时需要一个变量maxi 来记录所有子序列数组中最大的连续子序列元素和,而这个值就是我们的最大子矩阵元素和。 

最终代码:

通常我们会对二维数组多申请一行和一列的空间并初始化为0,是为了dp的状态转移方程在使用时不需要对边界情况进行特殊处理,并且不对dp数组元素的结构造成影响,提高代码的简洁性。

#include <iostream>
#include <vector>
using namespace std;const int N = 101; // 数组的最大大小int main() {int n = 0;cin >> n; // 初始化大小为 (n+1) x (n+1) 的二维数组,所有元素为 0vector<vector<int> > arr(n + 1, vector<int>(n + 1, 0)); // 初始化大小为 (n+1) x (n+1) 的二维向量 dp,所有元素为 0vector<vector<int> > dp(n + 1, vector<int>(n + 1, 0)); // 输入矩阵的元素,并计算每列的前缀和for(int i = 1; i <= n; i++) {for(int j = 1; j <= n; j++) {cin >> arr[i][j];dp[i][j] = dp[i - 1][j] + arr[i][j]; // 计算每列的前缀和}}int maxi = -128 * 1E4; // 初始化最大值为一个该题中最小的数int ans = 0;// 遍历所有的行for(int i = 0; i < n; i++)  // 从第 0 行到第 n-1 行{for(int j = i + 1; j <= n; j++) // 从第 1 行到第 n 行{ ans = 0;// 遍历每一列,并计算当前行对的最大子序列和for(int k = 1; k <= n; k++)  // 遍历每一列的元素{// 计算当前行对的最大子序列和,并更新 ansans = max(dp[j][k] - dp[i][k], ans + dp[j][k] - dp[i][k]);// 更新目前为止找到的最大子序列和(maxi)maxi = max(ans, maxi);}}}cout << maxi << endl; // 输出最大子序列和return 0;
}

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

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

相关文章

JVM 类的加载过程详解

文章目录 1. 哪些类需要加载2. 类加载步骤2.1 装载2.1.1 这个过程都做了什么事2.1.2 类的模板对象2.1.3 二进制流获取方式2.1.4 Class 实例的位置2.1.5 数组类的加载有什么不同 2.2 链接2.2.1 验证2.2.2 准备2.2.3 解析 2.3 初始化 1. 哪些类需要加载 在 Java 中数据类型分为 …

最好的网校在线教育系统,怎样才能有效的提高听课质量效率?

课堂是学生获得知识的主要渠道&#xff0c;课堂听课效率是会影响到学习效果的&#xff0c;听课质量是关系着学生学习成绩好和坏的因素之一&#xff0c;那怎样提高听课效率&#xff1f; 第一、养成课前预习。学生课前不预习&#xff0c;缺乏知识准备&#xff0c;容易导致上课听不…

Python 整数类型(int)详解:无限范围与多种进制

引言 在编程中&#xff0c;整数是最基本的数据类型之一。不同编程语言对整数的处理方式各不相同&#xff0c;这往往影响到程序的性能和开发者的选择。本文将深入探讨 Python 中的整数类型&#xff08;int&#xff09;&#xff0c;其独特的处理方式&#xff0c;以及它在日常编程…

Ubuntu24 文件目录结构——用户——权限 详解

目录 权限 用户 文件目录结构 一个目录可以有程序&#xff0c;目录&#xff0c;文件&#xff0c;以及这三者的链接。可以看到还分别有使用者和权限信息。 每个文件和目录都有与之关联的三个主要属性&#xff1a;所有者&#xff08;owner&#xff09;、组&#xff08;group&a…

软件工程期末复习(3)

软件生命周期 一般问题的解决过程&#xff1a; 问题的阐述&#xff1a;界定问题&#xff0c;用较宽的范围而不是细节来定义和描述待解问题&#xff1b; 问题的分析&#xff1a;问题定义的提炼&#xff0c;把问题分成可以理解和处理的子问题&#xff0c;进而提供基本细节&…

小区物业管理系统

文章目录 小区物业管理系统一、项目演示二、项目介绍三、部分功能截图四、部分代码展示五、底部获取项目源码&#xff08;9.9&#xffe5;带走&#xff09; 小区物业管理系统 一、项目演示 小区物业管理系统 二、项目介绍 基于springbootvue的前后端分离物业管理系统 系统角…

Java锁事

1.1 大厂面试题 一、Synchronized相关问题 1.Synchronized 用过吗&#xff0c;其原理是什么? 2.你刚才提到获取对象的锁&#xff0c;这个“锁”到底是什么?如何确定对象的领? 3.什么是可重入性&#xff0c;为什么说Synchronized是可重入锁? 4.JVM对Java的原生锁做了哪些优化…

到底什么是ASIC和FPGA?

ASIC&#xff08;Application Specific Integrated Circuit&#xff0c;应用特定集成电路&#xff09;和FPGA&#xff08;Field Programmable Gate Array&#xff0c;现场可编程门阵列&#xff09;都是集成电路&#xff08;IC&#xff09;的重要类型&#xff0c;它们在电子设计…

上海市计算机学会竞赛平台2022年3月月赛丙组正负数判断

题目描述 给定一个整数 &#x1d45b;n&#xff0c;若 &#x1d45b;n 为正数&#xff0c;输出 Positive&#xff0c;若 &#x1d45b;n 为负数&#xff0c;输出 Negative&#xff0c;若 &#x1d45b;n 恰好为零&#xff0c;输出 Zero。 输入格式 单个整数&#xff1a;表示…

Ubuntu 24 换国内源及原理 (阿里源 清华源 中科大源 网易源)

备份原文件 sudo cp /etc/apt/sources.list.d/ubuntu.sources /etc/apt/sources.list.d/ubuntu.sources.bak 编辑源文件 sudo gedit /etc/apt/sources.list.d/ubuntu.sources 粘贴到文本&#xff08;其中一个即可&#xff09;&#xff1a; &#xff08;阿里源&#xff09…

数据结构与算法学习笔记之线性表三---顺序表的动态分配存储表示和实现(C++)

目录 前言 1.顺序表的动态分配存储表示 2.顺序表的静态存储表示 1.初始化 2.销毁 3.清空 4.判空 5.表长 6.数据元素 7.获取下标 8.前驱节点 9.后继节点 10.插入 11.删除 12.遍历 13.测试代码 前言 这篇文章讲的是线性表的动态分配存储表示。 1.顺序表的…

【JAVA进阶篇教学】第十三篇:Java中volatile关键字讲解

博主打算从0-1讲解下java进阶篇教学&#xff0c;今天教学第十三篇&#xff1a;volatile关键字讲解。 在 Java 中&#xff0c;volatile关键字是一种轻量级的同步机制&#xff0c;用于确保变量的可见性和禁止指令重排序。本文将详细解释volatile关键字的工作原理、可见性保证以及…

买卖股票的最佳时机 II(LeetCode 122)

❤️❤️❤️ 欢迎来到我的博客。希望您能在这里找到既有价值又有趣的内容&#xff0c;和我一起探索、学习和成长。欢迎评论区畅所欲言、享受知识的乐趣&#xff01; 推荐&#xff1a;数据分析螺丝钉的首页 格物致知 终身学习 期待您的关注 导航&#xff1a; LeetCode解锁100…

react组件渲染性能优化之函数组件-memo使用

函数组件赋值相同的值视图不会渲染&#xff0c;类组件会&#xff0c;因为函数组件默认就阻止了 import {useState} from react export default function App() {const [counter, setCounter] useState(1)console.log(App组件渲染了);return(<div><h1>{counter}&l…

已解决java.lang.AbstractMethodError: 抽象方法错误的正确解决方法,亲测有效!!!

已解决java.lang.AbstractMethodError: 抽象方法错误的正确解决方法&#xff0c;亲测有效&#xff01;&#xff01;&#xff01; 目录 问题分析 报错原因 解决思路 解决方法 更新和重新编译依赖 确认类和接口的版本一致性 类加载器配置检查 总结 问题分析 java.lang.…

防城港知识付费系统,教学活动在设计中需要注意什么?如何进行教学设计?

老师们的教学能否让学生学到知识&#xff0c;让学生懂得书本上的道理&#xff0c;那么教学活动很重要&#xff0c;因此在设计的时候&#xff0c;一定要根据教学的目的以及孩子的特点来进行设计。 假设老师为学生进行美术绘画的教学活动设计&#xff0c;那么需要注意以下三个方面…

实现字符串复制(C语言)

一、N-S流程图&#xff1b; 二、运行结果&#xff1b; 三、源代码&#xff1b; # define _CRT_SECURE_NO_WARNINGS # include <stdio.h>int main() {//初始化变量值&#xff1b;int i 0;char a[100], b[100];//获取字符串&#xff1b;printf("请为数组a输入字符串…

使用模拟SPI接口驱动串行接口的LCD( STM32F4)

目录 概述 1. 硬件介绍 1.1 ST7796-LCD 1.2 MCU IO与LCD PIN对应关系 2 代码实现 2.1 STM32CubeMX 6.11生成工程 2.2 IO模拟SPI接口 2.3 实现LCD的驱动 3 测试 测试代码下载地址&#xff1a; stm32-f407-lcd-ft6336-proj资源-CSDN文库 gitee下载地址&#xff1a; h…

「PHP系列」PHP AJAX运用之数据库实例/XML实例

文章目录 一、AJAX 数据库实例1. 创建数据库和表2. HTML页面 (index.html)3. PHP脚本 (search.php)4. 配置和运行 二、AJAX XML实例1. PHP 脚本 (get_data.php)2. HTML 页面 (index.html)3. 配置和运行 三、相关链接 一、AJAX 数据库实例 一个PHP和AJAX结合使用来从数据库获取…

【Spring】验证 @ServerEndpoint 的类成员变量线程安全

文章目录 前言猜想来源验证方法Controller 的情况ServerEndpoint 的情况 后记 前言 最近有 websocket 的需求。探索 ServerEndpoint 的类成员变量特点。 这里类比 Controller 讨论 ServerEndpoint 类成员变量是否线程安全。 猜想来源 网上的教程大多数都这么展示程序&#…