求二维矩阵最大子矩阵和

//暴力枚举的话需要On6的时间复杂度，用矩阵前缀和+最大子序列则可以将时间复杂度减少到On3，同时使用矩阵数组空间复杂度只要On2

//输入一个矩阵，求和最大的子矩阵
#include <iostream>
#include <cstdio>
#include <algorithm>

#define MAX_SIZE 100000
int r,c;
int martix[100][100];
//输入矩阵，值可正可负
void input()
{
	for(int i=0; i<r; i++)
		for(int j=0; j<c; j++)
			scanf("%d",&martix[i][j]);
}
//存储0,0为起点，右下角坐标为x,y的子矩形和，第一行第一列为sum[1][1]
int sum[100][100];
//获得宽度为1，长度为l,右下角坐标为x,y的子矩阵的和
#define COM(x,y,l) (sum[(x)][(y)]+sum[(x)-(l)][(y)-1]-sum[(x)][(y)-1]-sum[(x)-(l)][(y)])

int solve()
{
	int ans=0,s,temp;
	//计算矩阵前缀和sum,耗费平方复杂度
	for(int i=0; i<r; i++)
		for(int j=0; j<c; j++)
			sum[i+1][j+1]=sum[i][j+1]+sum[i+1][j]-sum[i][j]+martix[i][j];
	//枚举矩阵的行坐标，然后将y-x+1长度的单列求和看成一个数字用最大连续和算法求出结果，耗费立方复杂度
	for(int y=1; y<=r; y++)
		for(int x=y; x<=r; x++)
		{
			//求最大连续子序列算法
			temp=0;s=0;
			for(int i=1;i<=c;i++)
			{
				if(s<=0)
					s = COM(x,i, x-y+1);
				else
					s+= COM(x,i, x-y+1);
				if(s>temp)
					temp=s;
			}
			if(temp>ans) 
				ans=temp;
		}
		return ans;
}

int main() 
{
    scanf("%d %d",&r,&c);
    input();
    printf("%d",solve());
    return 0;
}

//别人blog拷来的，最大连续和算法On复杂,貌似在某本书上看到过这个属于什么在线联机算法
int maxsum(int a[], int len)
{
    int maxsum, maxhere;
    maxsum = maxhere = a[0];   //初始化最大和为a【0】
    for (int i=1; i<len; i++) {
        if (maxhere <= 0)
            maxhere = a[i];  //如果前面位置最大连续子序列和小于等于0，则以当前位置i结尾的最大连续子序列和为a[i]
        else
            maxhere += a[i]; //如果前面位置最大连续子序列和大于0，则以当前位置i结尾的最大连续子序列和为它们两者之和
        if (maxhere > maxsum) {
            maxsum = maxhere;  //更新最大连续子序列和
        }
    }
    return maxsum;
}

//别人blog拷来的，三维矩阵的情况，表示看不太懂，先放在这

#include <cstdlib>
#include <cstdio>
#include <iostream>
#include <algorithm>
#include <cstring>
using namespace std;
typedef long long LL;

const LL INF = 1LL<<60;
int A, B, C;
LL sumz[25][25][25]; // 平面上[1,1]到[i][j]在高度为k时的区域前缀和

LL get(int x1, int y1, int x2, int y2, int z) {
    return sumz[x2][y2][z] - sumz[x1-1][y2][z] - sumz[x2][y1-1][z] + sumz[x1-1][y1-1][z];
}

void solve() {
    // 需要枚举平面上的二维组合
    LL ret = -INF;
    LL Min, s;
    for (int i = 1; i <= A; ++i) {
        for (int j = i; j <= A; ++j) {
            for (int k = 1; k <= B; ++k) {
                for (int h = k; h <= B; ++h) {
                    Min = 0;
                    for (int p = 1; p <= C; ++p) {
                        s = get(i, k, j, h, p);
                        ret = max(ret, s - Min);
                        if (s < Min) Min = s;
                    }
                }
            }
        }
    }
    printf("%lld\n", ret);
}

int main() {
    int T;
    scanf("%d", &T);
    while (T--) {
        scanf("%d %d %d", &A, &B, &C);
        LL tot, x;
        for (int i = 1; i <= A; ++i) {
            for (int j = 1; j <= B; ++j) {
                tot = 0;
                for (int k = 1; k <= C; ++k) {
                    scanf("%lld", &x);
                    tot += x;
                    sumz[i][j][k] = sumz[i-1][j][k] + sumz[i][j-1][k] - sumz[i-1][j-1][k] + tot;
                }
            }
        }
        solve();
        if (T) puts("");
    }
    return 0;
}