1. 描述最大字段和的分治算法
题目
思路
判断最大子段和,可以用分治的思想,每次将序列一分为二,选择两个序列的最大子段和。
但是这里还有一种可能,就是子段可以横跨两个子序列,所以我们的最大子段和就是:
MAX(左边序列最大字段和,横跨两序列的最大子段和,右边序列的最大子段和)。
对于左右两边的最大子段和,可以用分治递归的方法来做,临界条件就是序列中只剩一个数了,这时候最大子段和就是这个数,而递归函数就是对左右两边分别求最大子段和(调用自身),而且还得求跨序列的最大子段和,取三者的最大值来返回。
那么怎么求跨序列的最大子段和呢?其实很简单,首先要对原来的大序列添加几个指针,开头的是指针l,最右边的是指针r,因为要分治,所以再设置一个中间的指针mid,此时序列就可以分为两个部分,分别是(l,mid)和(mid+1,r),这时候的跨序列子段,必须包含mid和mid+1这两个地方,当然也可以向左或向右延申,所以,我们只需要求出从mid开始向左延申的最大字段和,还有从mid+1开始向右延申的最大子段和,将两者相加,就能得到跨序列的最大子段和了。
思路很好理解,照着上面的描述画出图来就一目了然了。下面来看看代码实现吧。
代码
#include<bits/stdc++.h> using namespace std; const int N = 1e5; int n, a[N]; int maxSum (int left, int right) { if (left == right) return a[left]; int mid = left + right >> 1; int lmax = maxSum (left, mid); int rmax = maxSum (mid + 1, right); int sum = a[mid]; int clmax = a[mid]; for (int i = mid - 1; i >= left; i--) { sum += a[i]; if (sum > clmax) clmax = sum; } sum = a[mid + 1]; int crmax = a[mid + 1]; for (int i = mid + 2; i <= right; i++) { sum += a[i]; if (sum > crmax) crmax = sum; } int cmax = clmax + crmax; int maxsum = max (cmax, max (lmax, rmax)); if (maxsum < 0) maxsum = 0; return maxsum; } int main () { cin >> n; for (int i = 0; i < n; i++) cin >> a[i]; cout << maxSum (0, n - 1); return 0; }
2. 分析该算法的时间复杂度
分解子问题:O(1)
求解子问题:2T(n/2)
合并子问题:O(n)
故时间复杂度为T(n)=2T(n/2)+O(n)=nlogn
3. 对分治法的体会和思考
分治法的设计思想是,将一个难以直接解决的大问题,分割成一些规模较小的相同问题,以便各个击破,分而治之。
其中的划分再击破,和递归的分解再解决异曲同工,其实同样用到了递归的思想,只不过分治法先分再治,最后还得合并。