最大子数组(线性时间算法思路)
从头系统学习算法个人随想
习题来源:算法导论4.1-5
在题干中写道,已知[1…j]的最大子数组,则[1…j+1]的最大子数组或者来源于[1…j]的最大子数组,或者来源于[i…j+1]的最大子数组(1<=i<=j+1)。
一开始,我对这段话的第一反应是动态规划,f(i,j)表示[i…j]之间的最大子数组,但是想了很久,没有想出来一个正常的状态转移方程(可能有吧,我太菜了)。后来,看了大佬对这种算法的解释,我就如醍醐灌顶一般,明白了线性算法的思路。
先放代码,再解释:
#include<iostream>
#include<cstdio>
using namespace std;
int sum_max(int num[], int lenth);
int Max(int num1, int num2);
int main(){
int num[] = {13, -3, -25, 20, -3, -16, -23, 18, 20, -7, 12, -5, -22, 15, -4, 7};
int length = sizeof(num) / sizeof(num[0]);
int maxx = sum_max(num, length);
cout << maxx;
return 0;
}
int sum_max(int num[], int lenth){
int num_max[lenth];
num_max[0] = num[0];
for(int i=1; i<lenth; i++){
if(num_max[i-1] + num[i] > num[i]){
num_max[i] = num_max[i-1] + num[i];
}
else{
num_max[i] = num[i];
}
}
int maxx = num_max[0];
for(int i=1; i<lenth; i++){
maxx = Max(maxx, num_max[i]);
}
return maxx;
}
int Max(int num1, int num2){
return num1>num2 ? num1 : num2;
}
这个算法还是和动态规划同源,f(i)表示以第i个数结尾的最大子数组,那么,就有状态转移方程: f(i) = max( f(i-1)+num[i], num[i] )。
意思就是以第i个数结尾的最大子数组有两种可能:
- 以第i-1个数结尾的最大子数组+第i个数。
- 第i个数
在这两种情况中取得最大值,就是f(i)。
实现代码如下:
for(int i=1; i<lenth; i++){
if(num_max[i-1] + num[i] > num[i]){
num_max[i] = num_max[i-1] + num[i];
}
else{
num_max[i] = num[i];
}
}
之后,我们只需要去寻找f(i)中的最大值,就是我们要找的最大子数组啦!