最大子数组问题:
在一个数列当中寻找一个子数列,使得这个子数列的元素之和最大。为了使代码简单化,在这里我们只求出最大值。
①暴力搜索,即在C(n,2)种可能的方法中寻找符合题意得一种,显然时间复杂度为θ(n²);
#include <iostream>
#include <climits>
using namespace std;
int a[16] = {13,-3,-25,20,-3,-16,-23, 18,20,-7,12 ,-5,-22,15,-4,7};
int main()
{
int sum = -INT_MAX;
int cur_sum;
for(int i=0;i<16;++i)
{
cur_sum = 0;
for(int j=i;j<16;++j)
{
cur_sum += a[j];
if(cur_sum > sum)
sum = cur_sum;
}
}
cout<<sum;
return 0;
}
②分治法
分:将原数列拆成两个大小相同的子数列(对半分);
治:对两个子数列递归求解,当子数列只有一个元素时,返回这个元素的值;
合并:返回左数列,右数列及中间数列(必须贯穿左右数列)的最大值
根据分析,此方法的时间复杂度为Θ(n*lgn)
#include <iostream>
#include <climits>
using namespace std;
int a[16] = {13,-3,-25,20,-3,-16,-23, 18,20,-7,12 ,-5,-22,15,-4,7};
int findCrossbigSubarray(int low,int mid,int high) //和一定要包含a[mid]和a[mid+1]
{
int sum = 0;
int left_big = -INT_big;
for(int i=mid; i>=low; --i) //mid属于左数组
{
sum += a[i];
if(sum > left_big) left_big = sum;
}
sum = 0;
int right_big = -INT_big;
for(int i=mid+1; i<=high; ++i) //mid+1属于右数组
{
sum += a[i];
if(sum > right_big) right_big = sum;
}
return right_big + left_big;
}
int findbigSubarray(int low, int high)
{
if(low == high) return a[low];
int mid = (low + high) / 2;
int left_big = findbigSubarray(low,mid); //左数组的最大值
int right_big = findbigSubarray(mid+1,high); //右数组的最大值
int cross_big = findCrossbigSubarray(low,mid,high); //穿越数组的最大值
int big = big(left_big,right_big);
if(cross_big > big) big = cross_big;
return big; //返回三者的最大值
}
int main()
{
cout<<findbigSubarray(0,15);
return 0;
}
③这题最快的解法还是要属DP(动态规划),其时间复杂度能够降低至θ(n)
我们需要一个辅助数组sum,其中sum[i]为以a[i]为尾元素的最大子数列【精髓】,待sum数组全部完成后,找出sum数组里元素的最大值,即为答案。
动态转移方程:
sum[i+1] = max(sum[i]+a[i+1], a[i+1])
#include <iostream>
using namespace std;
int a[16] = {13,-3,-25,20,-3,-16,-23, 18,20,-7,12 ,-5,-22,15,-4,7};
//sum = {13,10,-15,20,17, 1, -22, 18,38,31,43 ,38, 16,31,27,34};
int sum[16] = {a[0]};
int main()
{
for(int i=0; i<16; ++i)
sum[i+1] = max(sum[i]+a[i+1], a[i+1]);
//取sum中的最大值
int Max = sum[0];
for(int i=1; i<16; ++i)
if(sum[i] > Max)
Max = sum[i];
cout<<Max;
return 0;
}
以上就是最大子数组问题的三大解法。