最大子段和及其变形小结

最近做刷51nod的时候，看到了这类题目，挺有意思的，小结一下

51nod1049 最大子段和

题意：O（-1）

思路：经典问题

复杂度：O（n）

代码：

#include <iostream>  
#include <cstdio>  
using namespace std;  
const int MAXN = 50000+5;  
int a[MAXN];  
int main()  
{  
    int n;  
    while(~scanf("%d",&n))  
    {  
        for(int i=1; i<=n; i++)  
            scanf("%d",&a[i]);  
        long long Max = -1, sum = 0;  
        for(int i=1; i<=n; i++)  
        {  
            if(sum >= 0)  
                sum += a[i];  
            else  
                sum = a[i];  
            if(sum > Max)  
                Max = sum;  
        }  
        printf("%lld\n",Max);  
    }  
    return 0;  
}  

51nod1065 最小正子段和

题意：O（-1）

思路：先预处理一下前缀和，再对其排序，再用计算相邻两个前缀和贪心策略做，详解见Here

复杂度：O（nlogn）

代码：

#include <bits/stdc++.h>
using  namespace  std;
#define ll __int64
const int N=5e4+10;

struct node
{
    ll val;
    int pos;
    bool operator < (const node& tmp)const
    {
        if(val == tmp.val)return pos>tmp.pos;
        return val<tmp.val;
    }
}p[N];


int  main(){
    int n, mx;
    scanf("%d", &n);

    p[0].val=0;
    p[0].pos=0;

    for(int i=1; i<=n; i++){
        scanf("%lld", &p[i].val);
        p[i].val += p[i-1].val;
        p[i].pos=i;
    }
    sort(p, p+n+1);  

    ll ans=0;
    for(int i=1; i<=n; i++){
        if(p[i].pos>p[i-1].pos && p[i].val>p[i-1].val){
            if(ans == 0)ans = p[i].val - p[i-1].val;
            else ans = min(ans, p[i].val - p[i-1].val);
        }
    }
    printf("%lld\n", ans);
    return 0;
}

51nod1050 循环数组最大子段和

题意：给定一个长度为50000的数组，求它的循环数组的最大子段和。

思路：

本题与普通的最大子段和问题不同的是，最大子段和可以是首尾相接的情况，即可以循环。那么这个题目的最大子段和有两种情况

（1）正常数组中间的某一段和最大。这个可以通过普通的最大子段和问题求出。
（2）此数组首尾相接的某一段和最大。这种情况是由于数组中间某段和为负值，且绝对值很大导致的，那么我们只需要把中间的和为负值且绝对值最大的这一段序列求出，用总的和减去它就行了。

即，先对原数组求最大子段和，得到ans，然后把数组中所有元素符号取反，再求最大子段和，得到ans1，
原数组的所有元素和为sum，那么最终答案就是ans=max(ans, sum + ans1)。（这里的加号稍微想一下）

复杂度：O（n）

代码：

#include <iostream>
#include <cstdio>
#include <cstring>
#include <cstdlib>
#include <cmath>
#include <algorithm>
#include <set>
#define pr(x) cout << #x << "= " << x << "  " ；
#define pl(x) cout << #x << "= " << x << endl;
using namespace std;
typedef long long ll;
const int N = 50005;

int a[N], n;

ll maxsum(){
    ll sum=0, mx=0;
    for(int i=1; i<=n; i++){
        sum += a[i];
        if(sum<0)sum=0;
        if(sum > mx)mx=sum;
    }
    return mx;
}

int main(){
    ll sum = 0, ans;
    scanf("%d", &n);
    for(int i=1; i<=n; i++){
        scanf("%d", &a[i]);
        sum += a[i];
    }
    ans = maxsum();
    for(int i=1; i<=n; i++)a[i]=-a[i];
    ll ans1 = maxsum();
    ans = max(ans, sum+ans1);
    printf("%lld\n", ans);
    return 0;
}

51nod1051 最大矩阵和

题意：O（-1）

思路：

最大矩阵和其实就是二维的而已，我们可以转化为最大子段和求解。 首先我们枚举行数，在将第i行到第j行的每一列加起来，看成一个数，最后枚举按照最大字段和的思路枚举列数。 再记录下枚举过程中记录下最大值即可。

复杂度：O（m*n^2）

代码：

#include <cstdio>
#include <cstring>
#include <cstdlib>
#include <cmath>
#include <iostream>
#include <algorithm>
#include <set>
using namespace std;
typedef long long ll;

ll mat[505][505], ans=0;
int n, m;

int main(){
    scanf("%d%d", &m, &n);
    for(int i=1; i<=n; i++)
        for(int j=1; j<=m; j++){
            ll x;
            scanf("%lld", &x);
            mat[i][j] = mat[i-1][j]+x;
        }
    for(int i=1; i<=n; i++){
        for(int j=i; j<=n; j++){
            ll sum=0;
            for(int k=1; k<=m; k++){
                sum += mat[j][k] - mat[i-1][k];
                if(sum<0)sum=0;
                if(sum>ans)ans=sum;
            }
        }
    }
    printf("%lld\n", ans);
    return 0;
}

51nod1052最大M子段和

题意：O（-1）

思路：

用两个数组，pre[MAXN]和dp[MAXN]。
首先m次循环，第x次循环代表的是把整个序列分成x个子段所能得到的最大x子段和。
pre[i]数组记录的是，从第1个数到第i个数被分成x个子段所能得到的最大子段和。

dp[i]数组记录的是，从第1个数到第i个数（包括第i个数在内）被分成x个子段所能得到的最大子段和

假设当前已分成了x个子段(即最外层循环执行了x次)，然后需要执行第x+1次时：
则，dp[i]的转移是从下列两种情况取max值。

1.前i-1个数分为x个子段得到的最大和pre[i-1]加上单独把input[i]作为第x+1个子段的开头；
2.从前i-1个数中已经分出了x+1个子段，且第x+1个子段的尾部需要加上input[i]，即dp[i-1]+input[i]

复杂度：O（n*m）

代码：

#include <iostream>
#include <cstdio>
#include <algorithm>
#include <cstring>
#include <cstdlib>
#include <cmath>
#include <set>
#define pr(x) cout << #x << "= " << x << "  " ；
#define pl(x) cout << #x << "= " << x << endl;
using namespace std;
typedef long long ll;
const int N=5005;
const int inf=0x3f3f3f3f;

ll dp[N], pre[N], input[N];
int n, m;

int main(){
    scanf("%d%d", &n, &m);
    dp[0]=-inf;
    pre[0]=0;
    for(int i=1; i<=n; i++){
        scanf("%lld", &input[i]);
        pre[i]=0;
        dp[i]=-inf;
    }
    ll ans;
    while(m--){
        ans=-inf;
        for(int i=1; i<=n; i++){
            if(i==1)dp[i]=pre[i-1]+ input[i];
            if(dp[i-1]>pre[i-1])
                dp[i]=dp[i-1]+input[i];//第i个数与它的上一个数在同一子段内
            else 
                dp[i]=pre[i-1]+input[i];//从第i个分出一个新的子段
             //dp[i] = max(dp[i - 1], pre[i - 1]) + input[i]; 或者直接这么取值也可
            pre[i-1]=ans; //先用分成x个子段的pre，然后再在这里更新pre
            ans=max(ans, dp[i]);
        }
        pre[n]=ans;
    }
    printf("%lld\n", pre[n]);
    return 0;
}

推荐习题：

POJ2479 Maximum sum  （最大m字段和特例，令m=2）

51nod1254 最大子段和 V2

51nod1053 最大M子段和 V2

51nod1115 最大M子段和 V3