最大子段和问题

问题提出

有一由n个整数组成的序列A={a1,a2,…an,}A={a1,a2,…an,}，截取其中从i→j的子段并计算字段和，那么最大的字段和为多少？

举例说明

例如：

i	1	2	3	4	5	6	7	8
a[i]	1	3	-5	3	-2	5	-4	3

则其最大字段和为：3+(−2)+5=6

---

1. 动态规划算法

假设最大字段和我们设为 M
我们设 1→j中包括a[j]最大字段和为 b[j]
即：
                             b[j]=max{0,b[j−1]}+a[j]

我们求出包括a[i]→a[n]的所有最大字段和后，只需要求出这些最大字段和中最大的那个就可以了！
即：
                             M=max b[j]    1≤j≤n

i	1	2	3	4	5	6	7	8
a[i]	1	3	-5	3	-2	5	-4	3
b[i]	1	4	-1	3	1	6	-2	3

故可知 M=6

算法描述

#include<stdio.h>
int getMaxSum(int a[], int n, int &begin, int &end); 
int main(){
	int n, begin, end, maxSum, i;
	scanf("%d", &n);
	int arr[n];
	for(i = 0; i < n; i++)
	scanf("%d", &arr[i]);
	maxSum = getMaxSum(arr, n-1, begin, end);
	printf("maxSum = %d  begin:%d end:%d", maxSum, begin, end);
	
} 
int getMaxSum(int a[], int n, int &begin, int &end)
{
	int i;
	int temp = 0;
	int max = 0;
	for(i = 0; i < n; i++)
	{
		if(temp > 0)
		{
			end = i - 1;//当之前temp<0当前的temp>0时，i的值就是当之前temp<0时数组的元素所对应的i，而end应该是之前的那个元素所对应的i 
			temp += a[i];
		}
			
		else
		{
			temp = a[i];
			begin = i;
		}
			
		if(temp > max)
			max = temp;	
	}
	
	return max;
}

2.分治递归的算法

• 思想：将序列划分为左右两部分，则最大子段和可能在三处出现：左半部、右半部以及跨越左右边界的部分。递归的终止条件是：left == right。

算法描述

#include <iostream>
 
using namespace std;
int maxSubSum(int a[] , int left , int right , int &begin , int &end)
{
	int center , left_sum , right_sum , sum , i ;
	int sLeft = 0 , sRight = 0 ;//sLeft是以数组中央为界的指针向左移动时数字相加之和， sRight是以数组中央为界的指针向右移动时数字相加之和
	if(left == right)
		return a[left] > 0 ? a[left] : 0;
	else
	{
		center = (left + right) / 2;
		left_sum = maxSubSum(a , left , center , begin , end);//处理情况1
		right_sum = maxSubSum(a , center + 1 , right , begin , end); //处理情况2
		/*处理情况3*/
		int s1 = 0;//s1是以数组中央为界的左边最大的字段和 
		for(i = center ; i >= left ; i --)
		{
			sLeft += a[i] ;
			if(sLeft > s1)
			{
				s1 = sLeft;
				begin = i;	
			} 
		}
		
		int s2 = 0;//s2是以数组中央为界的右边最大的字段和 
		for(i = center + 1 ; i <= right ; i ++)
		{
			sRight += a[i] ;
			if(sRight > s2)
			{
				s2 = sRight;
				end = i;	
			} 
		}
		
		sum=s1+s2;
        if(sum < left_sum)
        	sum = left_sum;
        if(sum < right_sum)
        	sum = right_sum;
	} 
	
	return sum;
}	

int maxSum(int *a, int n, int &start, int &end)
{
    return maxSubSum(a,0,n,start,end);
}

int main()
{
    int n,start,end,maxsum;
    cin>>n;
    int *a=new int[n];
    for(int i=0;i<n;i++)
        cin>>a[i];
        
    maxsum = maxSum(a,n-1,start,end);
    cout<<maxsum<<" "<<++start<<" "<<++end;

    return 0;
}