hdu-4819-线段树套线段树

http://acm.hdu.edu.cn/showproblem.php?pid=4819

 给出一个N*N的矩阵，每次询问一个m*m的子矩阵里的floor((maxv+minv)/2)并把中间的元素修改为这个值。

 线段树套线段树，第一层X表示对行建立的线段树，内层表示对Y也就是列建立的线段树。

 分别对X和Y建立相应的函数来完成操作，当更新X树的节点时，再更新完当前X的节点的左右儿子之后，回头对X的这个节点对应的Y树进行更新，相当于X的左右儿子里的Y树来更新X的Y树，能理解这一点就很简单了。

　　

#include<bits/stdc++.h>
using namespace std;
#define pii pair<int,int>
#define mp make_pair
#define LL long long 
#define lc (id<<1)
#define rc (id<<1|1)
#define flc (fid<<1)
#define frc (fid<<1|1)
#define mid ((L+R)>>1)
const int maxn=810;
int maxv[maxn<<2][maxn<<2];
int minv[maxn<<2][maxn<<2]; 
int a[maxn][maxn],N,x1,x2,y1,y2;
void buildY(int,int,int,int,int,int);
void buildX(int id,int L,int R){
    if(L==R){
        buildY(id,L,R,1,1,N);
    }
    else{
        buildX(lc,L,mid);
        buildX(rc,mid+1,R);
        buildY(id,L,R,1,1,N);
    }
}
void buildY(int fid,int fl,int fr,int id,int L,int R){
    if(L==R){
        if(fl==fr)maxv[fid][id]=minv[fid][id]=a[fl][L];
        else{
            maxv[fid][id]=max(maxv[flc][id],maxv[frc][id]);
            minv[fid][id]=min(minv[flc][id],minv[frc][id]);
        }
    }
    else{
        buildY(fid,fl,fr,lc,L,mid);
        buildY(fid,fl,fr,rc,mid+1,R);
        maxv[fid][id]=max(maxv[fid][lc],maxv[fid][rc]);
        minv[fid][id]=min(minv[fid][lc],minv[fid][rc]);
    } 
}
pii askY(int fid,int id,int L,int R){
    if(L>=y1&&R<=y2){
        return mp(maxv[fid][id],minv[fid][id]);
    }
    else{
        if(y2<=mid) return askY(fid,lc,L,mid);
        else if(y1>mid) return askY(fid,rc,mid+1,R);
        else{
            pii pl=askY(fid,lc,L,mid),
                pr=askY(fid,rc,mid+1,R);
            return mp(max(pl.first,pr.first),min(pl.second,pr.second));
        }
    }
}
pii askX(int id,int L,int R){
    if(L>=x1&&R<=x2){
        return askY(id,1,1,N);
    }
    else{
        if(x2<=mid) return askX(lc,L,mid);
        else if(x1>mid) return askX(rc,mid+1,R);
        else{
            pii pl=askX(lc,L,mid),
                pr=askX(rc,mid+1,R);
            return mp(max(pl.first,pr.first),min(pl.second,pr.second));
        }
    }
}
void updateY(int fid,int fl,int fr,int id,int L,int R,int x,int y,int v){
    if(L==R){
        if(fl==fr){
            maxv[fid][id]=minv[fid][id]=v;
        }
        else{
            maxv[fid][id]=max(maxv[flc][id],maxv[frc][id]);
            minv[fid][id]=min(minv[flc][id],minv[frc][id]);
        }
    } 
    else{
        if(y<=mid){
            updateY(fid,fl,fr,lc,L,mid,x,y,v);
        }    
        else{
            updateY(fid,fl,fr,rc,mid+1,R,x,y,v);
        }
        maxv[fid][id]=max(maxv[fid][lc],maxv[fid][rc]);
        minv[fid][id]=min(minv[fid][lc],minv[fid][rc]);
    }
}
void updateX(int id,int L,int R,int x,int y,int v){
    if(L==R){
        updateY(id,L,R,1,1,N,x,y,v);
    }
    else{
        if(x<=mid) updateX(lc,L,mid,x,y,v);
        else updateX(rc,mid+1,R,x,y,v);
        updateY(id,L,R,1,1,N,x,y,v);
    }
}
int main()
{
    int T,M,i,j,x,y,d;
    scanf("%d",&T);
    for(int cas=1;cas<=T;++cas){
        scanf("%d",&N);
        for(i=1;i<=N;++i){
            for(j=1;j<=N;++j){
                scanf("%d",&a[i][j]);
            }
        }
        buildX(1,1,N);
        scanf("%d",&M);
        printf("Case #%d:
",cas);
        while(M--){
            scanf("%d%d%d",&x,&y,&d);
            x1=x-d/2,x2=x+d/2,y1=y-d/2,y2=y+d/2;
            x1=max(x1,1),x2=min(x2,N);
            y1=max(y1,1),y2=min(y2,N);
            pii ans=askX(1,1,N);
            updateX(1,1,N,x,y,(ans.first+ans.second)/2);
            printf("%d
",(ans.first+ans.second)/2);
        }
    }
    return 0;
}