NYOJ 104 最大和（求最大子矩阵）

OJ题目：猛戳~~

描述

给定一个由整数组成二维矩阵（r*c），现在需要找出它的一个子矩阵，使得这个子矩阵内的所有元素之和最大，并把这个子矩阵称为最大子矩阵。 
例子：
0 -2 -7 0 
9 2 -6 2 
-4 1 -4 1 
-1 8 0 -2 
其最大子矩阵为：

9 2 
-4 1 
-1 8 
其元素总和为15。 

输入

第一行输入一个整数n（0<n<=100）,表示有n组测试数据；
每组测试数据：
第一行有两个的整数r，c（0<r,c<=100），r、c分别代表矩阵的行和列；
随后有r行，每行有c个整数；

输出

输出矩阵的最大子矩阵的元素之和。

样例输入

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

样例输出

15

//参考在数组中求最大子段和的方法
//枚举行与行 ， 将行之间的列看作数组中的数
#include <iostream>
#include <string>
#include <string.h>
#include <map>
#include <stdio.h>
#include <algorithm>
#include <queue>
#include <vector>
#include <math.h>
#include <set>
#define Max(a,b) ((a)>(b)?(a):(b))
#define Min(a,b) ((a)<(b)?(a):(b))
#define SWAP(a , b){ int temp = a; a = b; b = temp;}
#define inf 100000000

using namespace std;
typedef long long LL;

int n, m;
int x[102][102];
LL ps[102][102];

LL BC(int i , int j , int c)//返回i行到j行的第c列的值
{
    return ps[j][c] - ps[j][c - 1] - ps[i - 1][c] + ps[i-1][c - 1];
}

LL MaxSum()
{
    LL maxsum = -inf;
    LL nend = 0 , nall;
    for(int i = 1;i <= n;i++)
        for(int j = i;j <= n;j++)
    {
        nend = BC(i , j , 1);//初始值为i行到j行的第1列
        nall = BC(i , j , 1);
        for(int k = 2;k <= m;k++)
        {
            if(nend < 0)
                nend = 0;
            nend += BC(i , j , k);
            nall = Max(nall , nend);
        }
        maxsum = Max(maxsum , nall);
    }
    return maxsum;
}

int main()
{
    int t;
    cin >> t;
    while(t--)
    {
        scanf("%d%d" , &n , &m);
        for(int i = 1;i <= n;i++)
            for(int j = 1;j <= m;j++)
            {
                scanf("%d" , &x[i][j]);
            }
        //下面求部分和
        for(int i = 0;i <= n;i++)
            ps[i][0] = 0;//边界值
        for(int j = 0;j <= m;j++)
            ps[0][j] = 0;//边界值
        for(int i = 1;i <= n;i++)
            for(int j = 1;j <= m;j++)
            ps[i][j] = ps[i - 1][j] + ps[i][j - 1] - ps[i - 1][j - 1] + x[i][j];//求部分和
        printf("%d\n" , MaxSum());
    }
    return 0;
}