PTA---01-复杂度1 最大子列和问题

给定K个整数组成的序列{ N₁,N₂,…,N_K}，“连续子列”被定义为{N_i,N_i+1,…,N_j}，其中 1≤i≤j≤K。“最大子列和”则被定义为所有连续子列元素的和中最大者。例如给定序列{ -2, 11, -4, 13, -5, -2 }，其连续子列{ 11, -4, 13 }有最大的和20。现要求你编写程序，计算给定整数序列的最大子列和。
本题旨在测试各种不同的算法在各种数据情况下的表现。各组测试数据特点如下：

• 数据1：与样例等价，测试基本正确性；
• 数据2：102个随机整数；
• 数据3：103个随机整数；
• 数据4：104个随机整数；
• 数据5：105个随机整数；

输入格式:

输入第1行给出正整数K (≤100000)；第2行给出K个整数，其间以空格分隔。

输出格式:

在一行中输出最大子列和。如果序列中所有整数皆为负数，则输出0。

输入样例:

6
-2 11 -4 13 -5 -2

输出样例:

20

---

根据浙江大学的数据结构课程，陈越姥姥的讲解，这里给出她在课堂上四个算法的实现代码。

算法1：复杂度T(N) = O(n^3)

int MaxSubseqSum1(int a[], int n)
{
    int max = 0;
    int i,j,k;
    int thisSum;
    for (i = 0; i < n; i++)//i是子列左端位置
    {
        for(j = i; j < n; j++)//j是子列右端位置
        {
            thisSum = 0;//thisSum是从a[i]到a[j]的子列和
            for(k = i; k < j; k++)
                thisSum += a[k];
            if(thisSum > max)//如果刚得到的这个子列和更大
                max = thisSum;//则更新结果
        }//j循环结束
    }//i循环结束
    return max;
}

算法2：复杂度T(N) = O(n^2)

int MaxSubseqSum2(int a[], int n)
{
    int max = 0;
    int i,j;
    int thisSum;
    for (i = 0; i < n; i++)//i是子列左端位置
    {
        thisSum = 0;//thisSum是从a[i]到a[j]的子列和
        for(j = i; j < n; j++)//j是子列右端位置
        {
            thisSum += a[j];//对于相同的i，不同的j，只要在j-1次循环的基础上累加1项即可
            if(max < thisSum)//如果刚得到的这个子列和更大
                max = thisSum;//则更新结果
        }//j循环结束
    }//i循环结束
    return max;
}

算法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，在线测试，复杂度T(N) = O(N)

int MaxSubseqSum4(int A[], int N )  
{   int ThisSum, MaxSum;
    int i;
    ThisSum = MaxSum = 0;
    for( i = 0; i < N; i++ ) {
        ThisSum += A[i]; /* 向右累加 */
            if( ThisSum > MaxSum )
                MaxSum = ThisSum; /* 发现更大和则更新当前结果 */
            else if( ThisSum < 0 ) /* 如果当前子列和为负 */
                ThisSum = 0; /* 则不可能使后面的部分和增大，抛弃之 */
    }
    return MaxSum;  
}

完整运行代码

#include "stdio.h"
#define size 100000

int MaxSubseqSum1(int a[], int n)
{
    int max = 0;
    int i,j,k;
    int thisSum;
    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 > max)
                max = thisSum;
        }
    }
    return max;
}

int MaxSubseqSum2(int a[], int n)
{
    int max = 0;
    int i,j;
    int thisSum;
    for (i = 0; i < n; i++)
    {
        thisSum = 0;
        for(j = i; j < n; j++)
        {
            thisSum += a[j];
            if(max < thisSum)
                max = thisSum;
        }
    }
    return max;
}


// 第三个算法
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 );
}

int MaxSubseqSum4(int A[], int N )  
{   int ThisSum, MaxSum;
    int i;
    ThisSum = MaxSum = 0;
    for( i = 0; i < N; i++ ) {
        ThisSum += A[i]; /* 向右累加 */
            if( ThisSum > MaxSum )
                MaxSum = ThisSum; /* 发现更大和则更新当前结果 */
            else if( ThisSum < 0 ) /* 如果当前子列和为负 */
                ThisSum = 0; /* 则不可能使后面的部分和增大，抛弃之 */
    }
    return MaxSum;  
}
int main(int argc, char const *argv[])
{
    int i, n, a[size];
    while(scanf("%d", &n) != EOF)
    {
        for (i = 0; i < n; ++i)
        {
            scanf("%d", &a[i]);
        }
        printf("%d\n", MaxSubseqSum4(a, n));
    }
    return 0;
}