最大子段和问题

问题描述：

       给定一个数组a[n]（其中可能有负数），求下标i和j，使a[i]+a[i+1]+…+a[j]最大。规定：如果所有数均是负数时，其最大子段和为0。

编程思路：

       这里介绍几种不同的思路来求解，分别如下：

（1）常规思路

       对所有可能的情况进行遍历，找到最优解，参考代码如下：

#include <stdio.h>
#define n 6
void main()
{
	int a[n]={-2,11,-4,13,-5,-2};
	int i,j,besti, bestj;
	int sum=0, thissum;
	for(i=0;i<n;i++){
		for(j=i;j<n;j++){
			thissum=0;
			for(int k=i;k<=j;k++){
				thissum+=a[k];
				if (thissum>sum){
					sum=thissum; besti=i; bestj=j;
}
}
}
}
printf(“the best i, j and sum is %d, %d and %d.\n”, i, j, sum);
}

（2）对常规思路的改进

      上述程序对所有可能的下标进行了遍历，因而可以得到最大子段和。但在计算的过程中，存在许多重复的情况，因而造成了算法的效率较为低下。考虑到

                                                       \sum_{k=1}^j a[k]=\sum_{k=1}^{j-1}a[k]+a[j]

因而可以根据该公式对算法进行改进，参考程序如下：

#include <stdio.h>
#define n 6
void main()
{
	int a[n]={-2,11,-4,13,-5,-2};
	int i,j,besti, bestj;
	int sum=0, thissum;
	for(i=0;i<n;i++){
		thissum=0;	//计算以i为起点的子段和的最大值
		for(j=i;j<n;j++){
				thissum+=a[j];
				if (thissum>sum){
					sum=thissum; besti=i; bestj=j;
}
}
}
}
printf(“the best i, j and sum is %d, %d and %d.\n”, i, j, sum);
}

（3）分治法

       基于分治法的策略，可以将给定的数组一分为二，分别计算以下三种情况的最大子段和：

       a. 前半段的某一子段的和是最大子段和；

       b. 后半段的某一子段的和是最大子段和；

       c. 中间的某一子段的和是最大子段和，该子段的起点一定在前半段，终点一定在后半段。而中间这一段的子段和可以分为两部分，以中点为终点的最大子段和和以中点为起点的最大子段和，这样才可以保证得到的子段和为最大的子段和。

       参考程序如下：    

   #include <stdio.h>
       int maxSum(int n, int *a);
       int maxSubSum(int *a, int left,int right);
       int main()
       {
              int a[6]={ -2,11,-4,13,-5,-2};
              printf(“the max sum is%d.\n”, maxSum(6, a));
}
int maxSum(int n, int *a)
{
       returnmaxSubSum(a, 0,n-1);
}
int maxSubSum(int *a, int left,int right)
{
       intsum;
if (left= =right)sum=a[left]>0?a[left]:0;
else{
       int mid=left+(right-left)/2;
       int leftSum=maxSubSum(a,0,mid);
       int rightSum=maxSumSum(a,mid+1,right);
       int s1=0, leftSum=0, s2=0, rightSum=0;
       for(int i=mid;i>=left;i--)
       {
              leftSum+=a[i];
if (leftSum>s1)
       s1=leftSum;
}
for(int i=mid+1;i<=right;i++)
{
       rightSum+=a[i];
       if(rightSum>s2)
              s2=rightSum;
}
sum=s1+s2;
if (leftSum>sum) sum=leftSum;
if (rightSum>sum)sum=rightSum;
}
return sum;
}

（4）动态规划法

       前面分别用了常规算法、改进的常规算法以及分治法对最大子段和进行了求解，下面采用动态规划法求解。

       将以j下标为终点的最大子段和记为b[j]，而在计算b[j]时，首先考虑b[j-1]的符号，如果b[j-1]>0，则b[j]=b[j-1]+a[j]；否则b[j]=a[j]。基于这一点，可以设计出求解最大子段和的动态规划算法，参考代码如下所示。

void main()
{
	int a[n]={-2,11,-4,13,-5,-2};
int sum=0, b=0;
for(int i=0;i<=n;i++){
	if (b>0) 
b+=a[i];
	else 
b=a[i];
	if (b>sum) 
sum=b;
}
printf(“the best sum is %d.\n”, sum);
}