HDU 1559 最大子矩阵

#include <iostream>
#include <cstdio>
#include <cstdlib>
#include <algorithm>
#include <queue>
#include <stack>
#include <map>
#include <cstring>
#include <climits>
#include <cmath>
#include <cctype>

const int inf = 0x3f3f3f3f;//1061109567
typedef long long ll;
using namespace std;

int a[1010][1010];
int b[1010];

int main()
{
int t;
int m,n,x,y;
int i,j,k;
scanf("%d",&t);
while(t--)
{
scanf("%d%d%d%d",&m,&n,&x,&y);
int max1 = -1;
for(i=0; i<m; i++)
{
for(int j=0; j<n; j++)
{
scanf("%d",&a[i][j]);
}
}
for(i=0; i<=m-x; i++)//枚举起始行
{
memset(b,0,sizeof(b));
for(j=i; j<=i+x-1; j++)
{
for(k=0; k<n; k++)
{
b[k] += a[j][k];
}
}
for(j=0; j<=n-y; j++)
{
int sum = 0;
for(k=j; k<=j+y-1; k++)
{
sum += b[k];
}
if(sum > max1)
{
max1 = sum;
}
}
}
printf("%d\n",max1);
}
return 0;
}