hdu4819 & poj2019 二维线段树

树套树~~
x线段树的每一个点都对应一个y线段树
这样就形成了二维线段树

二维线段树最大的理解难点在于对点的push_up，即由右往左更新和由下往上更新这两个地方，还有就是x线段的每一个点都对应一个y线段树

hdu 4819

#include <map>
#include <set>
#include <queue>
#include <stack>
#include <math.h>
#include <vector>
#include <cstdio>
#include <string>
#include<string.h>
#include <fstream>
#include <iostream>
#include <algorithm>
using namespace std;
#define exp 1e-8
#define INF 0x3f3f3f3f
#define ll long long
#define set(a,b) memset(a,b,sizeof(a));
#define set(a,b) memset(a,b,sizeof(a));
#define for1(a,b) for(int a=1;a<=b;a++)//1---(b)
#define for0(a,b) for(int a=0;a<=b;a++)//0---(b)
void bug(string st="bug")
{cout<<st<<endl;}
template<typename __ll>
inline void READ(__ll &m){
    __ll x=0,f=1;char ch=getchar();
    while(!(ch>='0'&&ch<='9')){if(ch=='-')f=-1;ch=getchar();}
    while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}
    m=x*f;
}
template<typename __ll>
inline void read(__ll &m){READ(m);}
template<typename __ll>
inline void read(__ll &m,__ll &a){READ(m);READ(a);}
template<typename __ll>
inline void read(__ll &m,__ll &a,__ll &b){READ(m);READ(a);READ(b);}
#define lson l,mid,rt<<1
#define rson mid+1,r,rt<<1|1
#define MAXN 900
int n,m;  //n*m的矩阵
int maxn[MAXN<<2][MAXN<<2];
int minn[MAXN<<2][MAXN<<2];
/*
更新方式有两种，一种是往上压，还有是往左压。
往左压：x是一个区间，y是一个值，往左压
往上压：y是一个区间，x是一个值或者说是一个固定区间
一般执行先后是：先往左压，后往上压

一般来说更新x，都是从下往上更新，所以一般没有固定的push_up_x的函数
如果硬是写的话也可以，只不过需要花费额外的时间
所以最好的方式是根据更新y的函数来更新x
如build_y,,,modify_x  就有更新x的功能。
因为省去了调用一些多余的区间，所以时间大大节省了
*/
void push_up_y(int rt,int f)  
{
    maxn[f][rt]=max(maxn[f][rt<<1],maxn[f][rt<<1|1]);
    minn[f][rt]=min(minn[f][rt<<1],minn[f][rt<<1|1]);
}
void build_y(int l,int r,int rt,int f,bool flag)
{                //flag决定往不往左压，不往的情况就是x是一个值
    if(l==r)
    {
        if(flag)  
        {              //左压
            maxn[f][rt]=max(maxn[f<<1][rt],maxn[f<<1|1][rt]);
            minn[f][rt]=min(minn[f<<1][rt],minn[f<<1|1][rt]);
            return ;
        }
        read(maxn[f][rt]);  
        minn[f][rt]=maxn[f][rt];
        return ;
    }
    int mid=l+r>>1;
    build_y(lson,f,flag);
    build_y(rson,f,flag);
    push_up_y(rt,f);  //上压
}
void build_x(int l,int r,int rt)
{
    if(l==r)
    {
        build_y(1,m,1,rt,0);
        return ;
    }
    int mid=l+r>>1;
    build_x(lson);
    build_x(rson);
    build_y(1,m,1,rt,1);  //每个点都对应一个线段树~~所以建立并更新
}
int query_maxn_y(int l,int r,int rt,int f,int y1,int y2)
{
    if(y1<=l&&r<=y2)
        return maxn[f][rt];
    int mid=r+l>>1;
    int tmaxn=-INF;
    if(y1<=mid)
        tmaxn=max(tmaxn,query_maxn_y(lson,f,y1,y2));
    if(y2>mid)
        tmaxn=max(tmaxn,query_maxn_y(rson,f,y1,y2));
    return tmaxn;
}
int query_maxn_x(int l,int r,int rt,int x1,int x2,int y1,int y2)
{
    if(x1<=l&&r<=x2)
        return query_maxn_y(1,m,1,rt,y1,y2);
    int mid=l+r>>1;
    int tmaxn=-INF;
    if(x1<=mid)
        tmaxn=max(tmaxn,query_maxn_x(lson,x1,x2,y1,y2));
    if(x2>mid)
        tmaxn=max(tmaxn,query_maxn_x(rson,x1,x2,y1,y2));
    return tmaxn;
}
int query_minn_y(int l,int r,int rt,int f,int y1,int y2)
{
    if(y1<=l&&r<=y2)
        return minn[f][rt];
    int mid=r+l>>1;
    int tminn=INF;
    if(y1<=mid)
        tminn=min(tminn,query_minn_y(lson,f,y1,y2));
    if(y2>mid)
        tminn=min(tminn,query_minn_y(rson,f,y1,y2));
    return tminn;
}
int query_minn_x(int l,int r,int rt,int x1,int x2,int y1,int y2)
{
    if(x1<=l&&r<=x2)
        return query_minn_y(1,m,1,rt,y1,y2);
    int mid=l+r>>1;
    int tminn=INF;
    if(x1<=mid)
        tminn=min(tminn,query_minn_x(lson,x1,x2,y1,y2));
    if(x2>mid)
        tminn=min(tminn,query_minn_x(rson,x1,x2,y1,y2));
    return tminn;
}
void modify_y(int l,int r,int rt,int f,int y,int val,bool flag)
{
    if(l==r)
    {
        if(flag)  //左压
        {
            maxn[f][rt]=max(maxn[f<<1][rt],maxn[f<<1|1][rt]);
            minn[f][rt]=min(minn[f<<1][rt],minn[f<<1|1][rt]);
        }
        else  maxn[f][rt]=minn[f][rt]=val;
        return ;
    }
    int mid=l+r>>1;
    if(y<=mid)
        modify_y(lson,f,y,val,flag);
    else modify_y(rson,f,y,val,flag);
    push_up_y(rt,f);  //上压
}
void modify_x(int l,int r,int rt,int x,int y,int val)
{
    if(l==r)
    {
        modify_y(1,m,1,rt,y,val,0);
        return ;
    }
    int mid=l+r>>1;
    if(x<=mid)
        modify_x(lson,x,y,val);
    else modify_x(rson,x,y,val);
    modify_y(1,m,1,rt,y,val,1);  //这个区间也要更新
}
int main()
{
    int cas,op,x,y,l;
    read(cas);
    for(int t=1;t<=cas;t++)
    {
        read(n);m=n;
        build_x(1,n,1);
        read(op);
        printf("Case #%d:\n",t);
        while(op--)
        {
            read(x,y,l);
            int x1=max(1,x-l/2),y1=max(1,y-l/2);
            int x2=min(n,x+l/2),y2=min(m,y+l/2);
            int smax=query_maxn_x(1,n,1,x1,x2,y1,y2);
            int smin=query_minn_x(1,n,1,x1,x2,y1,y2);
            printf("%d\n",(smax+smin)>>1);
            modify_x(1,n,1,x,y,(smax+smin)>>1);
        }
    }
    return 0;
}