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

相关文章

怎样通过分析GC日志来定位Java进程的内存问题

《怎样通过分析GC日志来定位Java进程的内存问题》:本文主要介绍怎样通过分析GC日志来定位Java进程的内存问题,具有很好的参考价值,希望对大家有所帮助,如有错误或未考虑完全的地方,望不吝赐教... 目录一、GC 日志基础配置1. 启用详细 GC 日志2. 不同收集器的日志格式二、关键指标与分析维度1.

MySQL查询JSON数组字段包含特定字符串的方法

《MySQL查询JSON数组字段包含特定字符串的方法》在MySQL数据库中,当某个字段存储的是JSON数组,需要查询数组中包含特定字符串的记录时传统的LIKE语句无法直接使用,下面小编就为大家介绍两种... 目录问题背景解决方案对比1. 精确匹配方案(推荐)2. 模糊匹配方案参数化查询示例使用场景建议性能优

Java 线程安全与 volatile与单例模式问题及解决方案

《Java线程安全与volatile与单例模式问题及解决方案》文章主要讲解线程安全问题的五个成因(调度随机、变量修改、非原子操作、内存可见性、指令重排序)及解决方案,强调使用volatile关键字... 目录什么是线程安全线程安全问题的产生与解决方案线程的调度是随机的多个线程对同一个变量进行修改线程的修改操

关于集合与数组转换实现方法

《关于集合与数组转换实现方法》:本文主要介绍关于集合与数组转换实现方法,具有很好的参考价值,希望对大家有所帮助,如有错误或未考虑完全的地方,望不吝赐教... 目录1、Arrays.asList()1.1、方法作用1.2、内部实现1.3、修改元素的影响1.4、注意事项2、list.toArray()2.1、方

Redis出现中文乱码的问题及解决

《Redis出现中文乱码的问题及解决》:本文主要介绍Redis出现中文乱码的问题及解决,具有很好的参考价值,希望对大家有所帮助,如有错误或未考虑完全的地方,望不吝赐教... 目录1. 问题的产生2China编程. 问题的解决redihttp://www.chinasem.cns数据进制问题的解决中文乱码问题解决总结

全面解析MySQL索引长度限制问题与解决方案

《全面解析MySQL索引长度限制问题与解决方案》MySQL对索引长度设限是为了保持高效的数据检索性能,这个限制不是MySQL的缺陷,而是数据库设计中的权衡结果,下面我们就来看看如何解决这一问题吧... 目录引言:为什么会有索引键长度问题?一、问题根源深度解析mysql索引长度限制原理实际场景示例二、五大解决

Springboot如何正确使用AOP问题

《Springboot如何正确使用AOP问题》:本文主要介绍Springboot如何正确使用AOP问题,具有很好的参考价值,希望对大家有所帮助,如有错误或未考虑完全的地方,望不吝赐教... 目录​一、AOP概念二、切点表达式​execution表达式案例三、AOP通知四、springboot中使用AOP导出

Python中Tensorflow无法调用GPU问题的解决方法

《Python中Tensorflow无法调用GPU问题的解决方法》文章详解如何解决TensorFlow在Windows无法识别GPU的问题,需降级至2.10版本,安装匹配CUDA11.2和cuDNN... 当用以下代码查看GPU数量时,gpuspython返回的是一个空列表,说明tensorflow没有找到

解决未解析的依赖项:‘net.sf.json-lib:json-lib:jar:2.4‘问题

《解决未解析的依赖项:‘net.sf.json-lib:json-lib:jar:2.4‘问题》:本文主要介绍解决未解析的依赖项:‘net.sf.json-lib:json-lib:jar:2.4... 目录未解析的依赖项:‘net.sf.json-lib:json-lib:jar:2.4‘打开pom.XM

IDEA Maven提示:未解析的依赖项的问题及解决

《IDEAMaven提示:未解析的依赖项的问题及解决》:本文主要介绍IDEAMaven提示:未解析的依赖项的问题及解决,具有很好的参考价值,希望对大家有所帮助,如有错误或未考虑完全的地方,望不吝... 目录IDEA Maven提示:未解析的依编程赖项例如总结IDEA Maven提示:未解析的依赖项例如