hdu4819 Mosaic（二维线段树）

题目链接

分析：
二维线段树的模板
线段树都很熟悉了
二维只不过是在每一个线段树结点中再开一棵线段树

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;
}