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

相关文章

好题——hdu2522(小数问题:求1/n的第一个循环节)

好喜欢这题,第一次做小数问题,一开始真心没思路,然后参考了网上的一些资料。 知识点***********************************无限不循环小数即无理数,不能写作两整数之比*****************************(一开始没想到,小学没学好) 此题1/n肯定是一个有限循环小数,了解这些后就能做此题了。 按照除法的机制,用一个函数表示出来就可以了,代码如下

hdu1043(八数码问题,广搜 + hash(实现状态压缩) )

利用康拓展开将一个排列映射成一个自然数,然后就变成了普通的广搜题。 #include<iostream>#include<algorithm>#include<string>#include<stack>#include<queue>#include<map>#include<stdio.h>#include<stdlib.h>#include<ctype.h>#inclu

csu 1446 Problem J Modified LCS (扩展欧几里得算法的简单应用)

这是一道扩展欧几里得算法的简单应用题,这题是在湖南多校训练赛中队友ac的一道题,在比赛之后请教了队友,然后自己把它a掉 这也是自己独自做扩展欧几里得算法的题目 题意:把题意转变下就变成了:求d1*x - d2*y = f2 - f1的解,很明显用exgcd来解 下面介绍一下exgcd的一些知识点:求ax + by = c的解 一、首先求ax + by = gcd(a,b)的解 这个

hdu2241(二分+合并数组)

题意:判断是否存在a+b+c = x,a,b,c分别属于集合A,B,C 如果用暴力会超时,所以这里用到了数组合并,将b,c数组合并成d,d数组存的是b,c数组元素的和,然后对d数组进行二分就可以了 代码如下(附注释): #include<iostream>#include<algorithm>#include<cstring>#include<stack>#include<que

购买磨轮平衡机时应该注意什么问题和技巧

在购买磨轮平衡机时,您应该注意以下几个关键点: 平衡精度 平衡精度是衡量平衡机性能的核心指标,直接影响到不平衡量的检测与校准的准确性,从而决定磨轮的振动和噪声水平。高精度的平衡机能显著减少振动和噪声,提高磨削加工的精度。 转速范围 宽广的转速范围意味着平衡机能够处理更多种类的磨轮,适应不同的工作条件和规格要求。 振动监测能力 振动监测能力是评估平衡机性能的重要因素。通过传感器实时监

hdu 1166 敌兵布阵(树状数组 or 线段树)

题意是求一个线段的和,在线段上可以进行加减的修改。 树状数组的模板题。 代码: #include <stdio.h>#include <string.h>const int maxn = 50000 + 1;int c[maxn];int n;int lowbit(int x){return x & -x;}void add(int x, int num){while

uva 10025 The ? 1 ? 2 ? ... ? n = k problem(数学)

题意是    ?  1  ?  2  ?  ...  ?  n = k 式子中给k,? 处可以填 + 也可以填 - ,问最小满足条件的n。 e.g k = 12  - 1 + 2 + 3 + 4 + 5 + 6 - 7 = 12 with n = 7。 先给证明,令 S(n) = 1 + 2 + 3 + 4 + 5 + .... + n 暴搜n,搜出当 S(n) >=

缓存雪崩问题

缓存雪崩是缓存中大量key失效后当高并发到来时导致大量请求到数据库,瞬间耗尽数据库资源,导致数据库无法使用。 解决方案: 1、使用锁进行控制 2、对同一类型信息的key设置不同的过期时间 3、缓存预热 1. 什么是缓存雪崩 缓存雪崩是指在短时间内,大量缓存数据同时失效,导致所有请求直接涌向数据库,瞬间增加数据库的负载压力,可能导致数据库性能下降甚至崩溃。这种情况往往发生在缓存中大量 k

6.1.数据结构-c/c++堆详解下篇(堆排序,TopK问题)

上篇:6.1.数据结构-c/c++模拟实现堆上篇(向下,上调整算法,建堆,增删数据)-CSDN博客 本章重点 1.使用堆来完成堆排序 2.使用堆解决TopK问题 目录 一.堆排序 1.1 思路 1.2 代码 1.3 简单测试 二.TopK问题 2.1 思路(求最小): 2.2 C语言代码(手写堆) 2.3 C++代码(使用优先级队列 priority_queue)

poj 3723 kruscal,反边取最大生成树。

题意: 需要征募女兵N人,男兵M人。 每征募一个人需要花费10000美元,但是如果已经招募的人中有一些关系亲密的人,那么可以少花一些钱。 给出若干的男女之间的1~9999之间的亲密关系度,征募某个人的费用是10000 - (已经征募的人中和自己的亲密度的最大值)。 要求通过适当的招募顺序使得征募所有人的费用最小。 解析: 先设想无向图,在征募某个人a时,如果使用了a和b之间的关系