分析:
二维线段树的模板
线段树都很熟悉了
二维只不过是在每一个线段树结点中再开一棵线段树
tip
一定要头脑清楚,很容易写错(可能是因为第一次写)
注意:x坐标的update放在了y坐标的构建中
因此传入的时候,不能妨碍到正常的区间边界判断
数组版,比较优美(就是变量名不是很高大上)
#include<cstdio>
#include<cstring>
#include<iostream>
using namespace std;
const int INF=1e9;
const int N=805;
int t[N<<2][N<<2],mn[N<<2][N<<2],mx[N<<2][N<<2];
int n,m,maxx,minn,a[N][N];
void update(int x,int y)
{
mx[x][y]=max(mx[x][y<<1],mx[x][y<<1|1]);
mn[x][y]=min(mn[x][y<<1],mn[x][y<<1|1]);
}
void build_y(int bh,int l,int r,int x,int xrt)
{
if (l==r)
{
if (x!=-1) mx[xrt][bh]=mn[xrt][bh]=a[x][l];
else
{
mn[xrt][bh]=min(mn[xrt<<1][bh],mn[xrt<<1|1][bh]);
mx[xrt][bh]=max(mx[xrt<<1][bh],mx[xrt<<1|1][bh]);
}
return;
}
int mid=(l+r)>>1;
build_y(bh<<1,l,mid,x,xrt);
build_y(bh<<1|1,mid+1,r,x,xrt);
update(xrt,bh);
}
void build_x(int bh,int l,int r)
{
if (l==r)
{
build_y(1,1,n,l,bh);
return;
}
int mid=(l+r)>>1;
build_x(bh<<1,l,mid);
build_x(bh<<1|1,mid+1,r);
build_y(1,1,n,-1,bh); //
}
void ask_y(int bh,int l,int r,int L,int R,int xrt)
{
if (l>=L&&r<=R)
{
maxx=max(maxx,mx[xrt][bh]);
minn=min(minn,mn[xrt][bh]);
return;
}
int mid=(l+r)>>1;
if (L<=mid) ask_y(bh<<1,l,mid,L,R,xrt);
if (R>mid) ask_y(bh<<1|1,mid+1,r,L,R,xrt);
}
void ask_x(int bh,int l,int r,int xl,int xr,int yl,int yr)
{
if (l>=xl&&r<=xr)
{
ask_y(1,1,n,yl,yr,bh);
return;
}
int mid=(l+r)>>1;
if (xl<=mid) ask_x(bh<<1,l,mid,xl,xr,yl,yr);
if (xr>mid) ask_x(bh<<1|1,mid+1,r,xl,xr,yl,yr);
}
void change_y(int bh,int l,int r,int y,int z,int xrt)
{
if (l==r)
{
if (z!=-1) mx[xrt][bh]=mn[xrt][bh]=z;
else
{
mx[xrt][bh]=max(mx[xrt<<1][bh],mx[xrt<<1|1][bh]);
mn[xrt][bh]=min(mn[xrt<<1][bh],mn[xrt<<1|1][bh]);
}
return;
}
int mid=(l+r)>>1;
if (y<=mid) change_y(bh<<1,l,mid,y,z,xrt);
else change_y(bh<<1|1,mid+1,r,y,z,xrt);
update(xrt,bh);
}
void change_x(int bh,int l,int r,int x,int y,int z)
{
if (l==r)
{
change_y(1,1,n,y,z,bh);
return;
}
int mid=(l+r)>>1;
if (x<=mid) change_x(bh<<1,l,mid,x,y,z);
else change_x(bh<<1|1,mid+1,r,x,y,z);
change_y(1,1,n,y,-1,bh); //
}
int main()
{
int T;
scanf("%d",&T);
for (int cas=1;cas<=T;cas++)
{
scanf("%d",&n);
int x,y,z;
for (int i=1;i<=n;i++)
for (int j=1;j<=n;j++)
scanf("%d",&a[i][j]);
build_x(1,1,n);
printf("Case #%d:\n",cas);
scanf("%d",&m);
for (int i=1;i<=m;i++)
{
scanf("%d%d%d",&x,&y,&z);
z/=2; maxx=-INF,minn=INF;
ask_x(1,1,n,max(1,x-z),min(x+z,n),max(1,y-z),min(y+z,n));
printf("%d\n",(maxx+minn)/2);
change_x(1,1,n,x,y,(maxx+minn)/2);
}
}
return 0;
}