问题描述:
给定一个数组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);
}