Maximum subarray problem(“最大和子数组”问题)与Kadane’s algorithm

本文主要是介绍Maximum subarray problem(“最大和子数组”问题)与Kadane’s algorithm,希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

最大和子串:在计算机科学中,最大子数组问题就是在给定一维数组中寻找和最大的连续的子数组,并确定起点索引i,终点索引j和最大和。(该定义翻译自wikipedia)

                                                          

最直接的方法就是穷举每一个子串计算和:暴力搜索1(Brute force)

def findMaxSum(arr):n=len(arr)maxSum=arr[0]beginIndex=0endIndex=0for i in range(n):for j in range(i,n):tempSum=0for k in range(i,j+1):tempSum += arr[k]if maxSum<tempSum:maxSum=tempSumbeginIndex=iendIndex=jreturn maxSum,beginIndex,endIndex

该算法复杂度为O(n^3)很容易发现,在确定子数组起点索引i后,j在增大的过程中求和运算存在大量冗余计算,故改进为:采用一个临时变量将从i到j的和保存下来,每次j增大时,这个临时变量只要加上一个数组值就行了:

def findMaxSum1(arr):n=len(arr)maxSum=arr[0]beginIndex=0endIndex=0for i in range(n-1):tempSum=0  #这里用临时变量存储i到j的和for j in range(i,n):tempSum += arr[j]if maxSum<tempSum:maxSum=tempSumbeginIndex=iendIndex=jreturn maxSum,beginIndex,endIndex

该算法时间复杂度为O(n^2),还能不能继续优化?

 

分治(Divide & Conquer):递归地将原问题二分得到子问题,直到子问题为只有一个元素的子数组为止。然后再把子问题归并,归并两个子数组时需要将子数组A的“最大和子数组”的和子数组B的“最大和子数组”的和以及“AB合并数组”的“最大和子数组”的和三者进行比较得到最大值。(回忆一下归并排序中的分而治之的归并思维

def findMaxSum2(arr,start,end):if start>=end:return arr[start],start,startmid=start+(end-start)//2maxSumLeft,beginLeftIndex,endLeftIndex=findMaxSum2(arr,start,mid)maxSumRight,beginRightIndex,endRightIndex=findMaxSum2(arr,mid+1,end)leftExtendMaxSum=leftExtendTempSum=0tempIndex1=0for i in range(mid,start-1,-1):leftExtendTempSum+=arr[i]if leftExtendTempSum>leftExtendMaxSum:leftExtendMaxSum=leftExtendTempSumtempIndex1=irightExtendMaxSum=rightExtendTempSum=0tempIndex2=0for i in range(mid+1,end+1):rightExtendTempSum+=arr[i]if rightExtendTempSum>rightExtendMaxSum:rightExtendMaxSum=rightExtendTempSumtempIndex2=imaxSum=leftExtendMaxSum+rightExtendMaxSumbeginIndex=tempIndex1endIndex=tempIndex2if maxSum<maxSumLeft:beginIndex=beginLeftIndexendIndex=endLeftIndexmaxSum=maxSumLeftif maxSum<maxSumRight:beginIndex=beginRightIndexendIndex=endRightIndexmaxSum=maxSumRightreturn maxSum,beginIndex,endIndex

该算法的时间复杂度为O(nlogn)

下面来膜拜一下线性时间复杂的解决此问题的算法!!!

Kadane’s algorithm算法的intuition竟然是Dynamic Programme!!!

如果我们知道数组中以索引位置对应元素结尾的“最大和子数组”的和为,那么以索引位置对应元素结尾的“最大和子数组”的和为多少?或者说以结尾的“最大和子数组”将结尾的“最大数组和”作为前缀。即:

                                                                     

因为该算法使用了最优子结构(在以每个位置结尾的最大和子数组是以一个更小且有所重叠的子问题来计算),可以看作是Dynamic Programme的一个应用,给大神跪了。

def findMaxSum3(arr):maxSum=-float("inf")beginIndex=0endIndex=0max_ending_here=maxSumcurrentStartIndex=currentEndIndex=0for currentEndIndex in range(0,len(arr)):max_ending_here+=arr[currentEndIndex]if max_ending_here>maxSum:maxSum,beginIndex,endIndex=max_ending_here,currentStartIndex,currentEndIndexif max_ending_here<0:max_ending_here=0currentStartIndex=currentEndIndex+1return maxSum,beginIndex,endIndex

只返回最大值的版本更简洁,直接从维基百科上copy过来:

def max_subarray(A):max_ending_here = max_so_far = A[0]for x in A[1:]:max_ending_here = max(x, max_ending_here + x)max_so_far = max(max_so_far, max_ending_here)return max_so_far

下面测试代码:

import numpy as np
arr=np.random.randint(-100,100,1000)
print(findMaxSum(arr))
print(findMaxSum1(arr))
print(findMaxSum2(arr,0,len(arr)-1))
print(findMaxSum3(arr))

输出一致,时间消耗差别明显!!!

(1520, 219, 554)
(1520, 219, 554)
(1520, 219, 554)
(1520, 219, 554)

这篇关于Maximum subarray problem(“最大和子数组”问题)与Kadane’s algorithm的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/920403

相关文章

IDEA和GIT关于文件中LF和CRLF问题及解决

《IDEA和GIT关于文件中LF和CRLF问题及解决》文章总结:因IDEA默认使用CRLF换行符导致Shell脚本在Linux运行报错,需在编辑器和Git中统一为LF,通过调整Git的core.aut... 目录问题描述问题思考解决过程总结问题描述项目软件安装shell脚本上git仓库管理,但拉取后,上l

JavaScript对象转数组的三种方法实现

《JavaScript对象转数组的三种方法实现》本文介绍了在JavaScript中将对象转换为数组的三种实用方法,文中通过示例代码介绍的非常详细,对大家的学习或者工作具有一定的参考学习价值,需要的朋友... 目录方法1:使用Object.keys()和Array.map()方法2:使用Object.entr

idea npm install很慢问题及解决(nodejs)

《ideanpminstall很慢问题及解决(nodejs)》npm安装速度慢可通过配置国内镜像源(如淘宝)、清理缓存及切换工具解决,建议设置全局镜像(npmconfigsetregistryht... 目录idea npm install很慢(nodejs)配置国内镜像源清理缓存总结idea npm in

pycharm跑python项目易出错的问题总结

《pycharm跑python项目易出错的问题总结》:本文主要介绍pycharm跑python项目易出错问题的相关资料,当你在PyCharm中运行Python程序时遇到报错,可以按照以下步骤进行排... 1. 一定不要在pycharm终端里面创建环境安装别人的项目子模块等,有可能出现的问题就是你不报错都安装

idea突然报错Malformed \uxxxx encoding问题及解决

《idea突然报错Malformeduxxxxencoding问题及解决》Maven项目在切换Git分支时报错,提示project元素为描述符根元素,解决方法:删除Maven仓库中的resolv... 目www.chinasem.cn录问题解决方式总结问题idea 上的 maven China编程项目突然报错,是

Python爬虫HTTPS使用requests,httpx,aiohttp实战中的证书异步等问题

《Python爬虫HTTPS使用requests,httpx,aiohttp实战中的证书异步等问题》在爬虫工程里,“HTTPS”是绕不开的话题,HTTPS为传输加密提供保护,同时也给爬虫带来证书校验、... 目录一、核心问题与优先级检查(先问三件事)二、基础示例:requests 与证书处理三、高并发选型:

前端导出Excel文件出现乱码或文件损坏问题的解决办法

《前端导出Excel文件出现乱码或文件损坏问题的解决办法》在现代网页应用程序中,前端有时需要与后端进行数据交互,包括下载文件,:本文主要介绍前端导出Excel文件出现乱码或文件损坏问题的解决办法,... 目录1. 检查后端返回的数据格式2. 前端正确处理二进制数据方案 1:直接下载(推荐)方案 2:手动构造

Python绘制TSP、VRP问题求解结果图全过程

《Python绘制TSP、VRP问题求解结果图全过程》本文介绍用Python绘制TSP和VRP问题的静态与动态结果图,静态图展示路径,动态图通过matplotlib.animation模块实现动画效果... 目录一、静态图二、动态图总结【代码】python绘制TSP、VRP问题求解结果图(包含静态图与动态图

MyBatis/MyBatis-Plus同事务循环调用存储过程获取主键重复问题分析及解决

《MyBatis/MyBatis-Plus同事务循环调用存储过程获取主键重复问题分析及解决》MyBatis默认开启一级缓存,同一事务中循环调用查询方法时会重复使用缓存数据,导致获取的序列主键值均为1,... 目录问题原因解决办法如果是存储过程总结问题myBATis有如下代码获取序列作为主键IdMappe

k8s容器放开锁内存限制问题

《k8s容器放开锁内存限制问题》nccl-test容器运行mpirun时因NCCL_BUFFSIZE过大导致OOM,需通过修改docker服务配置文件,将LimitMEMLOCK设为infinity并... 目录问题问题确认放开容器max locked memory限制总结参考:https://Access