二维线段树。先找到X的对应,再找Y的对应。
注意UPDATE要写好 一开始没写错了 单点更新的UPDATE要注意对准范围更新,否则TLE。。。
#include <algorithm>
#include <iostream>
#include <cstdio>
using namespace std;
int a[805][805], n;
int maxv[805<<2][805<<2], minv[805<<2][805<<2];
int x, y, x1, y1, x2, y2, max_ans, min_ans, val;
void update(int xrt,int l,int r,int rt)
{
int ls, rs;
ls=xrt<<1, rs=xrt<<1|1;
maxv[xrt][rt]=max(maxv[ls][rt],maxv[rs][rt]);
minv[xrt][rt]=min(minv[ls][rt],minv[rs][rt]);
if(l==r) return;
int mid=(l+r)>>1;
update(xrt,l,mid,rt<<1);
update(xrt,mid+1,r,rt<<1|1);
}
void one_update(int xrt,int l,int r,int rt)
{
int ls, rs;
ls=xrt<<1, rs=xrt<<1|1;
maxv[xrt][rt]=max(maxv[ls][rt],maxv[rs][rt]);
minv[xrt][rt]=min(minv[ls][rt],minv[rs][rt]);
if(l==r) return;
int mid=(l+r)>>1;
if(y<=mid)
one_update(xrt,l,mid,rt<<1);
else
one_update(xrt,mid+1,r,rt<<1|1);
}
void build_1D(int xrt,int x,int l,int r,int rt)
{
if(l==r){
maxv[xrt][rt]=minv[xrt][rt]=a[x][l];
return ;
}
int mid=(l+r)>>1;
build_1D(xrt,x,l,mid,rt<<1);
build_1D(xrt,x,mid+1,r,rt<<1|1);
maxv[xrt][rt]=max(maxv[xrt][rt<<1],maxv[xrt][rt<<1|1]);
minv[xrt][rt]=min(minv[xrt][rt<<1],minv[xrt][rt<<1|1]);
}
void build_2D(int l,int r,int rt)
{
if(l==r){
build_1D(rt,l,1,n,1);
return;
}
int mid=(l+r)>>1;
build_2D(l,mid,rt<<1);
build_2D(mid+1,r,rt<<1|1);
update(rt,1,n,1);
}
void query_1D(int xrt,int l,int r,int rt)
{
if(y1<=l&&r<=y2){
max_ans=max(max_ans,maxv[xrt][rt]);
min_ans=min(min_ans,minv[xrt][rt]);
return;
}
int mid=(l+r)>>1;
if(y1<=mid) query_1D(xrt,l,mid,rt<<1);
if(y2>mid) query_1D(xrt,mid+1,r,rt<<1|1);
}
void query_2D(int l,int r,int rt)
{
if(x1<=l&&r<=x2){
query_1D(rt,1,n,1);
return;
}
int mid=(l+r)>>1;
if(x1<=mid) query_2D(l,mid,rt<<1);
if(x2>mid) query_2D(mid+1,r,rt<<1|1);
}
void modify_1D(int xrt,int l,int r,int rt)
{
if(l==r){
maxv[xrt][rt]=minv[xrt][rt]=val;
return;
}
int mid=(l+r)>>1;
if(y<=mid) modify_1D(xrt,l,mid,rt<<1);
else modify_1D(xrt,mid+1,r,rt<<1|1);
maxv[xrt][rt]=max(maxv[xrt][rt<<1],maxv[xrt][rt<<1|1]);
minv[xrt][rt]=min(minv[xrt][rt<<1],minv[xrt][rt<<1|1]);
}
void modify_2D(int l,int r,int rt)
{
if(l==r){
modify_1D(rt,1,n,1);
return;
}
int mid=(l+r)>>1;
if(x<=mid) modify_2D(l,mid,rt<<1);
else modify_2D(mid+1,r,rt<<1|1);
one_update(rt,1,n,1);
}
int main()
{
int i, j, m, t, l, cs=0;
scanf("%d",&t);
while(t--){
printf("Case #%d:\n",++cs);
scanf("%d",&n);
for(i=1;i<=n;i++) for(j=1;j<=n;j++) scanf("%d",&a[i][j]);
build_2D(1,n,1);
scanf("%d",&m);
while(m--){
scanf("%d%d%d",&x,&y,&l);
l/=2;
x1=max(1,x-l), x2=min(n,x+l);
y1=max(1,y-l), y2=min(n,y+l);
min_ans=1<<30, max_ans=-(1<<30);
query_2D(1,n,1);
val=(max_ans+min_ans)/2;
printf("%d\n",val);
modify_2D(1,n,1);
}
}
return 0;
}