【分治】最大子数组问题

学习目标

分治
  • 分:将大问题分成小问题
  • 治:通过递归的办法解决一系列小问题
  • 归:最后再将小问题合起来
学习递归式的求解办法
  • 代入法
  • 递归树法
  • 主方法
  • 证明主定理

最大子数组问题
  • 输入:数组A[1…n]

  • 输出:(i,j,s)使得A[i…j]为最大数组,且和为s

    分:取数组中点mid,最大子数组要么完全位于左侧或者右侧,要么跨越中点

    治:完全位于一侧可递归解决,需要2T(n/2);跨越中点则从中点向两侧扩张遍历解决,需要O(n)

    合:三种情况取和最大的即可,需O(1)

    代价T(n) = 2T(n/2)+O(n)=nlogn

例题

给定一个整数数组nums,找到一个具有最大和的连续子数组(子数组最少包含一个元素),返回其最大和

示例:

输入:[-2,1,-3,4,-1,2,1,-5,4]
输出:6
解释:连续子数组[4,-1,2,1]的和最大,为6
#define	MAX(a,b,c) ((a)>(b)>(c)?(b):(c))?(a):((b))>(c)?(b):(c)))
int maxSubArray(int* nums, int numsSize) {
	int i = 0;
	int	iTmp = 0;
	int Max_L = 0;
	int Max_L = 0;
	int Max_R = 0;
	int Max_R = 0;
	int Max_l = 0 ;
	int Max_m = 0;
	if ((numsSize == 0) || (numsSize == 1)) {
		return nums[0];
	}
	else{
		Max_L = maxSubArray(&nums[0], (numsSize - 1) / 2);
		Max_R = maxSubArray(&nums[(numsSize + 1)/2] , numsSize / 2);
	}
	iTmp = 0;
	Max_l = nums[(numsSize - 1) / 2];//前半段右边界
	for (i = (numsSize - 1) / 2; i >= 0; i--) {
		iTmp += nums[i];
		if (iTmp > Max_l) {
			Max_l = iTmp; //Max_l为左侧部分最大和
		}
	}
	/*
	* 处理中间部分
	*/
	iTmp = Max_l;
	Max_m = Max_l;
	for (i = (numsSize + 1) / 2; i < numsSize; i++) {
		iTmp += nums[i];
		if (iTmp > Max_m)
		{
			Max_m = iTmp;
		}
		return MAX(Max_L, Max_R, Max_m);
	}
	
}
### 使用分治算法与递归解决最大子数组问题 #### 1. 基本概念 分治算法的核心思想是将一个问题划分为若干个规模较小的相同问题,分别求解这些子问题,最后将它们的结果合并得到原问题的解。在最大子数组问题中,目标是从给定的一维数组 `A` 中找到连续的一个子数组,使得该子数组的和达到最大。 #### 2. 解决方案描述 为了利用分治法解决问题,可以按照以下逻辑划分: - **分解阶段**:假设当前处理的是数组的一部分 `[low, high]`,将其从中点 `mid` 划分成两部分 `[low, mid]` 和 `[mid+1, high]`。 - **递归阶段**:对这两部分分别调用相同的函数来寻找各自的最大子数组。 - **合并阶段**:除了考虑左半边和右半边的最大子数组外,还需要额外计算跨越中间位置的最大子数组,并最终返回三者的最大值。 具体来说,任意一个子数组要么位于左侧区间 `[low, mid]`,要么位于右侧区间 `[mid+1, high]`,或者跨越这两个区间的边界形成新的组合[^2]。 以下是基于此理论的具体实现方法及其代码展示。 #### 3. 实现细节 下面是完整的 Python 函数定义以及辅助说明: ```python def find_max_crossing_subarray(A, low, mid, high): """ 找到跨越中间区域的最大子数组 """ left_sum = float('-inf') sum_ = 0 max_left = mid for i in range(mid, low - 1, -1): # 向左遍历 sum_ += A[i] if sum_ > left_sum: left_sum = sum_ max_left = i right_sum = float('-inf') sum_ = 0 max_right = mid + 1 for j in range(mid + 1, high + 1): # 向右遍历 sum_ += A[j] if sum_ > right_sum: right_sum = sum_ max_right = j return (max_left, max_right, left_sum + right_sum) def find_maximum_subarray_recursive(A, low, high): """ 主体递归函数用于查找整个范围内的最大子数组 """ if high == low: # 只有一个元素的情况 return (low, high, A[low]) mid = (low + high) // 2 left_low, left_high, left_sum = find_maximum_subarray_recursive(A, low, mid) right_low, right_high, right_sum = find_maximum_subarray_recursive(A, mid + 1, high) cross_low, cross_high, cross_sum = find_max_crossing_subarray(A, low, mid, high) if left_sum >= right_sum and left_sum >= cross_sum: return (left_low, left_high, left_sum) elif right_sum >= left_sum and right_sum >= cross_sum: return (right_low, right_high, right_sum) else: return (cross_low, cross_high, cross_sum) # 测试样例 if __name__ == "__main__": test_array = [-2, 1, -3, 4, -1, 2, 1, -5, 4] result = find_maximum_subarray_recursive(test_array, 0, len(test_array)-1) print(f"Maximum Subarray Indices: {result[:2]}, Sum: {result[-1]}") ``` 以上代码实现了三个主要功能模块: 1. 计算跨越两个分区之间的最大子数组; 2. 对于单侧数据集执行递归操作; 3. 将两者结合起来得出全局最优解。 #### 4. 时间复杂度分析 由于每次都将输入序列一分为二并独立处理每一段,再加上线性的扫描过程以确定跨界的最佳选项,总体的时间消耗大致相当于 O(n log n)[^4]。 ---
评论
成就一亿技术人!
拼手气红包6.0元
还能输入1000个字符
 
红包 添加红包
表情包 插入表情
 条评论被折叠 查看
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值