最大子序列和以及最大矩阵和

做题的时候遇到了一个最大矩阵和，搜了一下是从最大子序列和演变而来的，所以先做一下最大子序列和。

最大子序列和

这个有一个需要注意的事情是，如果矩阵里面所有的元素都是负值，可以返回0，也可以返回最大负数，看具体题目要求

一、暴力枚举 O(n^3)

对数组来说，一个最大的子序列可以通过遍历从每一个位置开始的序列中得到，所以我们遍历如a1,a2,…,an; a2,a3,…,an等不同长度不同起始位置的序列，不过时间复杂度太高O(n^3)，优化一下（保存加和）也不过是O(n^2)。

二、分治法 O(nlgn)

在这里，分治法的思想主要体现在，首先将原始序列分成左右两个部分，然后将最终的最大子序列的情况分成了三种。

1. 最大的子序列完全在原始序列左侧。
2. 最大的子序列完全在原始序列右侧。
3. 最大的子序列横跨左右两侧。
   其中可能让你迷惑的是第三个情况，如果横跨的话，那不就又相当于最原始的，不划分左右的情况吗？其实不是，我们在找左右最大两个子序列的时候，左侧是从middle遍历到left，右侧是从middle+1遍历到right，也就是说，我们得到的左右最大的子序列是连续的，横跨左右两侧的情况，也只是一个合并的过程，就是left_max+right_max。

#include <iostream>
using namespace std;
//求三个元素的最大值
int max_of_three(int i, int j, int k){
    if (i >= j && i >= k) {
        return i;
    }
    return max_of_three(j, k, i);
}

int maxsequence(int a[], int left, int right){
    if (left > right) return 0;
    if (left == right) return a[left];//只有一个元素，递归的基本情况
    int middle = (left + right) / 2;
    //先找左侧的最大元素
    int left_max = a[middle], left_sum = 0;
    for (int i = middle; i >= left; i--) {
        left_sum += a[i];
        if (left_max < left_sum) {
            left_max = left_sum;
        }
    }
    //再找右侧的最大元素
    int right_max = a[middle+1], right_sum = 0;
    for (int i = middle+1; i <= right; i++) {
        right_sum += a[i];
        if (right_max < right_sum) {
            right_max = right_sum;
        }
    }
    //返回中比较的对象分别是，横跨的，左侧的，右侧的最大子序列。
    return max_of_three(left_max+right_max, maxsequence(a, left, middle), maxsequence(a, middle+1, right));
}

三、线性查找

我们再来理性分析一下。

1. 最大子序列不可能从负数开始。如果从负数开始，那么假设a1是负数，a1+a2+a3

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

最大子矩阵

引用自这一篇文章
我实在是不会动态规划……我明天就去学动态规划呜呜呜

#include <iostream>  
#include <cstring>  
using namespace std;  
#define inf 0x3f3f3f3f;  
int n,m;  
int a[505][505],b[505],dp[505];  
int res=-inf;  
int main()  
{  
    cin>>n>>m;  
    for(int i=0;i<n;i++)  
       for(int j=0;j<m;j++)  
          cin>>a[i][j];  
    for(int x=0;x<n;x++){//设置起始行   
        memset(b,0,sizeof(b));//设置结束行   
        for(int y=x;y<n;y++){  
            for(int j=0;j<m;j++)  
                b[j]+=a[y][j];//b存储第起始行x到y行到中每一列的值   
            dp[0]=b[0];  

            if(dp[0]>res)//不要忘了这个   
                res=dp[0];  

            for(int i=1;i<m;i++){  
                if(dp[i-1]<0)  
                    dp[i]=b[i];  
                else   
                   dp[i]=dp[i-1]+b[i];  

                if(dp[i]>res)          //维护最后的结果   
                   res=dp[i];  
            }  

        }  
    }  
    cout<<res<<endl;  
    return 0;  
}