算法——最大子列和问题
问题:给定N个整数的序列{A1,A2,…AN},求函数f(i,j)=max{0,∑Ak}的最大值。
四种算法:
1.暴力求解(时间复杂度O(n^3))
把所有连续子列和全部算出,选最小。
`int MaxSubseqSum1(int A[], int N)
{
int ThisSum, MaxSum = 0;
int i, j, k;
for(i = 0;i < N;i++){
for(j = i;j < N;j++){
ThisSum = 0;
for(k = i;k <= j;k++)
ThisSum += A[k];
if(ThisSum > MaxSum)
MaxSum = ThisSum;
}
}
return MaxSum;
}
2.在当前已算出i到j的和的基础上,加新的元素值(时间复杂度O(n^2))
int MaxSunseqSum2(int A[], int N)
{
int ThisSum,MaxSum = 0;
int i, j;
for(i = 0;i < N;i++){
ThisSum=0; //ThisSum是A[i]到A[j]的和
for(j = i;j < N;j++){
ThisSum += A[j];
if(ThisSum > MaxSum)
MaxSum=ThisSum;
}
}
return MaxSum;
}
3.分治法(时间复杂度O(NlogN))
利用二分法
int Max3( int A, int B, int C )
{ /* 返回3个整数中的最大值 */
return A > B ? A > C ? A : C : B > C ? B : C;
}
int DivideAndConquer( int List[], int left, int right )
{ /* 分治法求List[left]到List[right]的最大子列和 */
int MaxLeftSum, MaxRightSum; /* 存放左右子问题的解 */
int MaxLeftBorderSum, MaxRightBorderSum; /*存放跨分界线的结果*/
int LeftBorderSum, RightBorderSum;
int center, i;
if( left == right ) { /* 递归的终止条件,子列只有1个数字 */
if( List[left] > 0 ) return List[left];
else return 0;
}
/* 下面是"分"的过程 */
center = ( left + right ) / 2; /* 找到中分点 */
/* 递归求得两边子列的最大和 */
MaxLeftSum = DivideAndConquer( List, left, center );
MaxRightSum = DivideAndConquer( List, center+1, right );
/* 下面求跨分界线的最大子列和 */
MaxLeftBorderSum = 0; LeftBorderSum = 0;
for( i=center; i>=left; i-- ) { /* 从中线向左扫描 */
LeftBorderSum += List[i];
if( LeftBorderSum > MaxLeftBorderSum )
MaxLeftBorderSum = LeftBorderSum;
} /* 左边扫描结束 */
MaxRightBorderSum = 0; RightBorderSum = 0;
for( i=center+1; i<=right; i++ ) { /* 从中线向右扫描 */
RightBorderSum += List[i];
if( RightBorderSum > MaxRightBorderSum )
MaxRightBorderSum = RightBorderSum;
} /* 右边扫描结束 */
/* 下面返回"治"的结果 */
return Max3( MaxLeftSum, MaxRightSum, MaxLeftBorderSum + MaxRightBorderSum );
}
int MaxSubseqSum3( int List[], int N )
{ /* 保持与前2种算法相同的函数接口 */
return DivideAndConquer( List, 0, N-1 );
}
4.在线处理(时间复杂度O(N))
int MaxSubseqSum4(int List[], int N){
int MaxSum = 0,ThisSum = 0;
for(int i = 0;i < N;i++){
ThisSum += List[i];
if(ThisSum > MaxSum){
MaxSum = ThisSum;
}
else if(thismax < 0){
ThisSum = 0;
}
}
return MaxSum;
}