问题描述:
最大子段和问题是将一个n个整数的序列a[1],a[2]….a[n]中字段a[first]….a[last]之和,(1<=first<=last<=n)求这些子段和中最大的。
例如(a[1],a[2],a[3],a[4],a[5],a[6])=(-2,11,-4,13,-5,-2)时,最大子段和为20,子段为a[2],a[3],a[4]。
求解方法:
如果不会算法,那就用时间复杂度为O(n^3)的枚举,i为从1到n的起点,j为从i到n的终点,k为从i到j的子段之和。
还是枚举,改进一下,得到O(n^2)的枚举算法,就是将k去掉,在找其终点j的时候就将子段和记录下来,因为从i到j的子段和就是从i到j-1的子段和加上a[j]。
再改进一下,将这个序列分成1到(1+n)/2的序列与(1+n)/2到n的序列。那么最大的子段有可能出现在:
1.左侧序列。2.右侧序列。3.跨越中间点的序列。
我们从中间点两侧找最大子段,再找越过中间点的最大子段,就形成了我们所说的分治算法,得到复杂度为O(nlogn)的算法。
其实,我们在选择一个元素a[j]的时候,只有两种情况,将a[i]至a[j-1]加上,或者从a[j]以j为起点开始。我们用一个数组dp[i]表示以i为结束的最大子段和,对于每一个a[i],加上dp[i-1]成为子段,或以a[i]开始成为新段的起点。因为我们只需要记录dp值,所以复杂度是O(n)。
这就是最大子段和的动态规划算法。
我们甚至不需要dp数组,只需要定义一个dp变量,因为最后要求的dp值也是最大的,所以我们可以在求dp的时候更新为最大的。
暴力法:
#include<iostream>
using namespace std;
int MaxSum(int n,int *a,int &besti,int &bestj)
{
int sum=0;
for(int i=0;i<n;i++)
{
for(int j=i;j<n;j++)
{
int tempSum=0;
for(int k=i;k<j;k++)
{
tempSum+=a[k];
if(tempSum>sum)
{
sum=tempSum;
besti=i;
bestj=j;
}
}
}
}
return sum;
}
int main()
{
int a[100]={-2,11,-4,13,-5,-2};
int n=6;
int besti=0;
int bestj=0;
MaxSum(n,a,besti,bestj);
cout<<"输出最大子段和:";
for(int i=besti;i<bestj;i++)
{
cout<<a[i]<<" ";
}
}
暴力法改进:
#include<iostream>
using namespace std;
int MaxSum(int n,int *a,int &besti,int &bestj)
{
int sum=0;
for(int i=0;i<n;i++)
{
int tempSum=0;
for(int j=i;j<n;j++)
{
tempSum+=a[j];
if(tempSum>sum)
{
sum=tempSum;
besti=i;
bestj=j;
}
}
}
return sum;
}
int main()
{
int a[100]={-2,11,-4,13,-5,-2};
int n=6;
int besti=0;
int bestj=0;
MaxSum(n,a,besti,bestj);
for(int i=besti;i<=bestj;i++)
{
cout<<a[i]<<" ";
}
}
动态规划:
#include<iostream>
using namespace std;
int MaxSum(int n,int *a,int &besti,int &bestj)
{
int sum=0;
int tempSum=0;
for(int i=0;i<n;i++)
{
if(tempSum>0)
{
tempSum+=a[i];
}
else
{
tempSum=a[i];
}
if(tempSum>sum)
{
sum=tempSum;
}
}
return sum;
}
int main()
{
int a[100]={-2,11,-4,13,-5,-2};
int n=6;
int besti=0;
int bestj=0;
cout<<MaxSum(n,a,besti,bestj);
}