ZOJ - 2859 Matrix Searching (二维线段树)

ZOJ - 2859 

这个题目以前学ST表的时候写过，但是ST表对于一些带修改的矩阵最值查询却没有办法。所以就整了一下二维线段树。

 

我对于二维线段树的一些理解：

1.二维线段树像是下面的结构。

像这样。。。

对于X坐标来说，每一个节点都是一颗线段树。所以总共的空间应该是16*n*n；

这样建树的话我们儿子节点于父节点之间的关系依旧保持着2*x 2*x+1的关系。

所以对于x来说  横坐标在（3-4）纵坐标在（1-n）的节点的信息可以有 （3-3）（1-n） 和 （4-4）（1-n）得到。

在矩阵上表现为

也就是后两列的最大值可以有第三列和第四列的最大值得到。

#include<bits/stdc++.h>
using namespace std;
#define LL long long
#define sca(x) scanf("%d",&x)
#define rep(i,j,k) for(int i=j;i<=k;i++)
#define pb(x) push_back(x)
#define N 310
#define inf 0x3f3f3f3f
#define mp(x,y) make_pair(x,y)


int sum[N*4][N*4];
int n;

void build_y(int rt,int l,int r,int id,int op)
{
    if(l==r)
    {
        if(op==1) scanf("%d",&sum[id][rt]);//找到y的叶子节点
        else sum[id][rt]=min(sum[id<<1][rt],sum[id<<1|1][rt]); //合并x的区间信息
        return ;
    }
    int m=(l+r)>>1;
    build_y(rt<<1,l,m,id,op);
    build_y(rt<<1|1,m+1,r,id,op);
    if(op==1)
    sum[id][rt]=min(sum[id][rt<<1],sum[id][rt<<1|1]); //op=1的时候我们合并y的区间信息
    else
    sum[id][rt]=min(sum[id<<1][rt],sum[id<<1|1][rt]);//op=2的时候合并该x区间所对应的y区间的信息。
}

void build_x(int rt,int l,int r)
{
    if(l==r)      //找到x的叶子节点。
    {
        build_y(1,1,n,rt,1);
        return ;
    }
    int m=(l+r)>>1;
    build_x(rt<<1,l,m);
    build_x(rt<<1|1,m+1,r);
    build_y(1,1,n,rt,2);  //合并x区间的信息
}


int query_y(int rt,int l,int r,int ql,int qr,int id)
{
    if(l>=ql && r<=qr)
    {
        return sum[id][rt];
    }
    int m=(l+r)>>1;
    int q1=inf,q2=inf;
    if(ql<=m) q1=query_y(rt<<1,l,m,ql,qr,id);
    if(qr>m) q2=query_y(rt<<1|1,m+1,r,ql,qr,id);
    return min(q1,q2);
}

int query_x(int rt,int ql,int qr,int l,int r,int y1,int y2)
{
    if(l>=ql && r<=qr)
    {
        return query_y(1,1,n,y1,y2,rt);
    }
    int m=(l+r)>>1;
    int q1=inf,q2=inf;
    if(ql<=m) q1=query_x(rt<<1,ql,qr,l,m,y1,y2);
    if(qr>m) q2=query_x(rt<<1|1,ql,qr,m+1,r,y1,y2);
    return min(q1,q2);
}
int main()
{
    int t;
    sca(t);
    while(t--)
    {
        sca(n);
        build_x(1,1,n);
        int m;
        sca(m);
        while(m--)
        {
            int r1,c1,r2,c2;
            scanf("%d%d%d%d",&r1,&c1,&r2,&c2);
            printf("%d\n",query_x(1,r1,r2,1,n,c1,c2));
        }
    }
}