HDU 4819 Mosaic

二维线段树。先找到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;
}