HDU 1559 最大子矩阵

#include<stdio.h>
#include<string.h>
#include<algorithm>
using namespace std;
#define N 1005
int add

;
int lower_bit(int i)
{
return  i&(-i);
}
void updata(int L,int R,int c)
{
for(int i=L;i<=N;i+=lower_bit(i))
{
for(int j=R;j<=N;j+=lower_bit(j))
{
add[i][j]+=c;
}
}
}
int sum(int x,int y)
{
int s=0;
for(int i=x;i>0;i-=lower_bit(i))
{
for(int j=y;j>0;j-=lower_bit(j))
{
s+=add[i][j];
}
}
return s;
}
int main()
{
int t;
scanf("%d",&t);
while(t--)
{
int m,n,x,y;
scanf("%d%d%d%d",&m,&n,&x,&y);
int v;
int i,j;
memset(add,0,sizeof(add));
for(i=1;i<=m;i++)
{
for(j=1;j<=n;j++)
{
scanf("%d",&v);
updata(i,j,v);
}
}
int M=0;
for(i=1;i<=m-x+1;i++)
{
for(j=1;j<=n-y+1;j++)
{
int x1,y1,x2,y2;
x1=i;
y1=j;
x2=i+x-1;
y2=j+y-1;
int s=(sum(x2,y2)-sum(x2,y1-1)-sum(x1-1,y2)+sum(x1-1,y1-1));
M=max(M,s);
}
}
printf("%d\n",M);
}
return 0;
}