最大子矩阵

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

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

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

Sample Output

2474

#include <cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
typedef long long ll;
ll dp[1010][1010];

int main()
{
    int T;
    scanf("%d",&T);
    while(T--){
        memset(dp,0,sizeof(dp));
        int m,n,x,y;
        scanf("%d %d %d %d",&m,&n,&x,&y);
        for(int i=1;i<=m;i++)
            for(int j=1;j<=n;j++)
            scanf("%I64d",&dp[i][j]);
        ll ans=0;
        for(int i=1;i<=m;i++){
            for(int j=1;j<=n;j++){
                dp[i][j]+=dp[i-1][j]+dp[i][j-1]-dp[i-1][j-1];
                if(i>=x&&j>=y){
                    ll temp_ans=dp[i][j]-dp[i-x][j]-dp[i][j-y]+dp[i-x][j-y];
                    ans=max(ans,temp_ans);
                }
            }
        }
        printf("%I64d\n",ans);
    }
    return 0;
}