吃西瓜 最大子矩阵 三维的。 rqnoj93

#include<cstdio>
#include<cctype>
#include<cstring>
#define rep(i,j,k) for(int i = j; i <= k; i++)
#define maxn 60
using namespace std;

int key[maxn][maxn][maxn], f[maxn][maxn] = {0}, dp[maxn];
int h, x, y;

inline int read() //重要的输入优化
{
    int s = 0, t = 1; char c = getchar();
    while( !isdigit(c) ){
        if( c == ‘-‘ ) t = -1; c = getchar();
    }
    while( isdigit(c) ){
        s = s * 10 + c - ‘0‘; c = getchar();
    }
    return s * t;
}

int maxsum3()
{
    int ans = -0xfffffff;
    rep(i,1,h){
        rep(j,i,h){
            rep(k,1,x){
                rep(l,1,y){
                    f[k][l] = key[j][k][l] - key[i-1][k][l];
                }
            }
            //重要优化所在，不要枚举k和l之间的所有值，把他们一次性算出来后，清零，下一次再重算，而是采用这种不断累加的方式。
            rep(k,1,x){
                memset(dp,0,sizeof(dp));
                rep(l,k,x){
                    rep(t,1,y) dp[t] += f[l][t];
                    int sum = 0;
                    rep(t,1,y){
                        if( sum > 0 ) sum += dp[t];
                        else sum = dp[t];
                        if( ans < sum ) ans = sum;
                    }
                }
            }
        }
    }
    return ans;
}

int main()
{
    h = read(), x = read(), y = read();
    rep(i,1,h){
        rep(j,1,x){
            rep(k,1,y){
                key[i][j][k] = read();
                key[i][j][k] += key[i-1][j][k];
            }
        }
    }
    int ans = maxsum3();
    printf("%d\n", ans);
    return 0;
}