贪心-最大子矩阵

最大子矩阵
总时间限制: 1000ms 内存限制: 65536kB
描述
已知矩阵的大小定义为矩阵中所有元素的和。给定一个矩阵，你的任务是找到最大的非空(大小至少是1 * 1)子矩阵。

比如，如下4 * 4的矩阵

0 -2 -7 0
9 2 -6 2
-4 1 -4 1
-1 8 0 -2

的最大子矩阵是

9 2
-4 1
-1 8

这个子矩阵的大小是15。

输入
输入是一个N * N的矩阵。输入的第一行给出N (0 < N <= 100)。再后面的若干行中，依次（首先从左到右给出第一行的N个整数，再从左到右给出第二行的N个整数……）给出矩阵中的N2个整数，整数之间由空白字符分隔（空格或者空行）。已知矩阵中整数的范围都在[-127, 127]。
输出
输出最大子矩阵的大小。

样例输入
4
0 -2 -7 0
9 2 -6 2
-4 1 -4 1
-1 8 0 -2
样例输出
15

贪心只是一种算法，而一道题往往要用到多种算法，比如这道题还需要用到枚举，因此，不要局限于找最大的某一个点，要先枚举子矩阵，再来判断子矩阵的最大值。
对于这道题而言，首先是枚举的子矩阵左上角，然后枚举子矩阵的高度，最后通过贪心来枚举子矩阵左上角的列数和子矩阵的宽度。
这道题的贪心策略实际上为：求一个数组中连续元素的和的最大值

参考代码

#include <cstdio>
#include <cstdlib>
#include <cmath>
#include <cstring>
#include <iostream>
#include <algorithm>
using namespace std;

int a;
int rect[105][105];
int arrays[105]; //保存列的和的数组

//贪心：求数组中连续的元素的和的最大值：如果已经累加的和大于0，肯定有助于和变大，就直接加上；如果小于等于0，那不如不加。
int sum()
{
    int max1 = arrays[0], temp = 0; //max1等于第一个元素而不能等于0，因为有可能是负数
    for (int i = 0; i < a; i++)
    {
        if (temp > 0) temp += arrays[i];
        else temp = arrays[i];
        max1 = max1 > temp ? max1 : temp; //考虑到会有当前数为负数但是相减后和仍为正的情况，所以每一步都要判断最大值
    }
    return max1;
}

int main()
{
#ifdef _DEBUG
    freopen("r.txt", "r", stdin);
#endif
    scanf("%d", &a);
    for (int i = 0; i < a; i++)
    {
        for (int j = 0; j < a; j++)
        {
            scanf("%d", &rect[i][j]);
        }
    }

    int res = rect[0][0]; //用于保存结果的变量，因为有可能是负数，所以初始化为可能的结果之一，即某一个子矩阵的值
    //以下过程实际上是枚举
    //i - 用于枚举的变量，不是当前行数，而是起始行数, j - 当前行数, k - 列数
    for (int i = 0; i < a; i++) //相当于在枚举子矩阵左上角的行数
    {
        memset(arrays, 0, sizeof(arrays)); //初始化用于保存和的变量
        for (int j = i; j < a; j++) //从第i行开始，相当于在枚举矩阵的高度
        {
            for (int k = 0; k < a; k++) //处理一行
            {
                arrays[k] += rect[j][k];
            }
            int temp = sum(); //判断当前的数据 - 相当于在枚举子矩阵左上角的列数及子矩阵的宽度
            if (temp > res)
            {
                res = temp;
            }
        }
    }

    printf("%d", res);
    return 0;
}