最大子矩阵和

给你一个m×n的整数矩阵，在上面找一个x×y的子矩阵，使子矩阵中所有元素的和最大。
输入描述：
输入数据的第一行为一个正整数T，表示有T组测试数据。每一组测试数据的第一行为四个正整数m,n,x,y（0<m,n<1000 AND 0<x<=m AND 0<y<=n），表示给定的矩形有m行n列。接下来这个矩阵，有m行，每行有n个不大于1000的正整数。
输出描述：
对于每组数据，输出一个整数，表示子矩阵的最大和。
代码：

#include<bits/stdc++.h>
using namespace std;
int a[1001][1001];
int main()
{
	int t;
	cin>>t;
	while(t--)
	{
		int m,n,x,y,i,j,k;
		cin>>m>>n>>x>>y;
		long long sumz,sumy,sum=0,maxsum=0;
		for(i=0; i<m; i++)
			for(j=0; j<n; j++)
				cin>>a[i][j];
		for(i=0; i<m-x; i++)
		{
			sum=0;
			for(j=0; j<n-y; j++)
			{
				if(j==0)
				{
					for(k=i; k<i+x; k++)
					{
						for(int z=0; z<y; z++)
							sum+=a[k][z];
					}
				}
				else
				{
					sumy=sumz=0;
					for(k=i; k<i+x; k++)
					{
						sumy+=a[k][j+y-1];
						sumz+=a[k][j-1];
					}
					sum=sum-sumz+sumy;
				}
				maxsum=maxsum>sum?maxsum:sum;
			}
		}
		cout<<maxsum<<endl;
	}
	return 0;
}