二维前缀和——模板

最大子矩阵

描述
给你一个m×n的整数矩阵，在上面找一个x×y的子矩阵，使子矩阵中所有元素的和最大。

输入
输入数据的第一行为一个正整数 T，表示有 T 组测试数据。每一组测试数据的第一行为四个正整数 n,m,x,y（0<m,n<1000 AND 0<x<=m AND 0<y<=n），表示给定的矩形有 n 行 m 列。接下来这个矩阵，有 n 行，每行有 m 个不大于 1000 的正整数。
输出
对于每组数据，输出一个整数，表示子矩阵的最大和。

输入样例 1
1
4 5 2 2
3 361 649 676 588
992 762 156 993 169
662 34 638 89 543
525 165 254 809 280

输出样例 1
2474

代码

#include <iostream>
#include <stdio.h>
#include <string.h>
#include <math.h>
#include <algorithm>
#include <string>
#include <string>
#include <map>
#include <bitset>
#include <vector>
void fre() { system("clear"), freopen("A.txt", "r", stdin); freopen("Ans.txt","w",stdout); }
void Fre() { system("clear"), freopen("A.txt", "r", stdin);}
void Run(int x = 0) {     
#ifdef ACM  //宏定义免注释 freopen
    if(! x) fre(); else Fre();
#endif
}
#define ios ios::sync_with_stdio(false)
#define Pi acos(-1)
#define pb push_back
#define fi first
#define se second
#define db double
#define ll long long
#define ull unsigned long long
#define Pir pair<ll, ll>
#define m_p make_pair
#define for_(i, s, e) for(ll i = (ll)(s); i <= (ll)(e); i ++)
#define rep_(i, e, s) for(ll i = (ll)(e); i >= (ll)(s); i --)
#define memset(a, b, c) memset(a, (int)b, c);
#define size() size() * 1LL
#define sc scanf
#define pr printf
#define sd(a) scanf("%lld", &a)
#define ss(a) scanf("%s", a)
#define inf 0x3f3f3f3f
#define INF 0x3f3f3f3f3f3f3f3f
#define esp 1e-7
#define mod (ll)(1e9 + 7)
using namespace std;
/*=========================ACMer===========================*/
const ll mxn = 1e3 + 10;
ll mz[mxn][mxn];


int main()
{
    Run();
    ll T; sd(T);
    while(T --)
    {
        ll n, m, x, y;
        sc("%lld %lld %lld %lld", &n, &m, &x, &y);
        for_(i, 1, n)
        {
            for_(j, 1, m)
            {
                sd(mz[i][j]);
                mz[i][j] += mz[i - 1][j] + mz[i][j - 1] - mz[i - 1][j - 1];
            }
        }

        ll ans = 0;
        for_(i, 1, n - x + 1)
        {
            for_(j, 1, m - y + 1)
            {
                ans = max(ans, mz[i + x - 1][j + y - 1] - mz[i + x - 1][j - 1] - mz[i - 1][j + y - 1] + mz[i - 1][j - 1]);
            }
        }

        pr("%lld\n", ans);
    }

    return 0;
}