10.25T2 二维线段树

Description

　　为了准备校庆庆典，学校招募了一些学生组成了一个方阵，准备在庆典上演出。
　　这个方阵是一个n*m的矩形，第i行第j列有一名学生，他有一个能力值Ai,j。
　　校长会定期检查一个p*q的方阵，询问这个方阵的学生能力值之和，或是学生能力值的最大值，或是学生能力值的最小值。由于校长不喜欢一个方阵长宽之比差太多，他每次询问的方阵的长不会超过宽的两倍。作为校庆筹办组组长的你，应该迅速并准确的回答校长所问的问题。

Input

　　第一行包含两个整数n,m，表示这个方阵的两条边的长度。
　　接下来n行，每行m个数，表示每个学生的能力值。
　　接下来一行包含一个整数q，表示校长的询问数。
　　接下来q行，每行先一个字符串s，接下来4个整数x1，y1，x2，y2，保证x1≤x2，y1≤y2，设以第x1行y1列为左上角，第x2行y2列为右下角的方阵为P。(本题为0下标)
　　若字符串内容为“SUM”，请求出P中所有学生的能力值之和。
　　若字符串内容为“MAX”，请求出P中所有学生的能力值的最大值。
　　若字符串内容为“MIN”，请求出P中所有学生的能力值的最小值。

Output

　　输出总共q行，第i行的数为第i组询问对应的答案ansi。

Sample Input

3 3 1 2 3 4 5 6 7 8 9 3 SUM 0 0 1 1 MAX 0 0 2 2 MIN 0 1 1 1

Sample Output

12 9 2

Hint

【样例解释】
　　对于第一组询问，能力值之和为1+2+4+5=12。
　　对于第二组询问，能力值最大的位置为第2行第2列。
　　对于第三组询问，能力值最小的位置为第0行第1列。
【数据规模和约定】
　　对于40%的数据，n,m≤200,q≤200
　　对于60%的数据，n,m≤300,q≤100000
　　对于80%的数据，n,m≤500,q≤500000
　　对于100%的数据，n,m≤800,q≤500000,0≤Aij≤3000,每个询问的方阵的长不超过宽的两倍。

 

 

 

 

二维线段树的板子。。。emmmm（好难写啊）

code：

#include<iostream>
#include<cstdio>
#include<string>
#define lc (p<<1)
#define rc (p<<1|1)
#define N 3301
#define max(i,j) i>j?i:j
#define min(i,j) i<j?i:j
using namespace std;
const int inf=0x3f3f3f3f;
int MAX=0,MIN=9999999;
int n,m;
int g[801][801],sum[801][801];
struct segY {
    struct node {
        int l,r,min,max;
        node() {
            min=9999999;
        }
    } t[N];
    int up,down;
    void pushup(int p) {
        t[p].max=max(t[lc].max,t[rc].max);
        t[p].min=min(t[lc].min,t[rc].min);
    }
    void build(int p,int l,int r) {
        t[p].l=l,t[p].r=r;
        if(l==r) {
            for(int i=down; i<=up; i++)t[p].max=max(t[p].max,g[i][l]),t[p].min=min(t[p].min,g[i][l]);
            return;
        }
        int mid=l+r>>1;
        build(lc,l,mid);
        build(rc,mid+1,r);
        pushup(p);
    }
    int query_max(int p,int ql,int qr) {
        if(ql<=t[p].l&&t[p].r<=qr) {
            return t[p].max;
        }
        int a1=0,a2=0;
        int mid=t[p].l+t[p].r>>1;
        if(ql<=mid)a1=query_max(lc,ql,qr);
        if(qr>mid)a2=query_max(rc,ql,qr);
        return a1>a2?a1:a2;
    }
    int query_min(int p,int ql,int qr) {
        if(ql<=t[p].l&&t[p].r<=qr) {
            return t[p].min;
        }
        int a1=inf,a2=inf;
        int mid=t[p].l+t[p].r>>1;
        if(ql<=mid)a1=query_min(lc,ql,qr);
        if(qr>mid)a2=query_min(rc,ql,qr);
        return a1>a2?a2:a1;
    }
};
struct segX {
    int l,r;
    segY Y;
} t[800*4];
void Build(int p,int l,int r) {
    t[p].l=l,t[p].r=r;
    t[p].Y.down=l,t[p].Y.up=r;
    t[p].Y.build(1,1,m);
    if(l==r)return;
    int mid=l+r>>1;
    Build(lc,l,mid);
    Build(rc,mid+1,r);
}
int Query_min(int p,int x1,int y1,int x2,int y2) {
    if(x1<=t[p].l&&t[p].r<=x2) {
        return t[p].Y.query_min(1,y1,y2);
    }
    int a1=inf,a2=inf;
    int mid=t[p].l+t[p].r>>1;
    if(x1<=mid)a1=Query_min(lc,x1,y1,x2,y2);
    if(x2>mid)a2=Query_min(rc,x1,y1,x2,y2);
    return a1>a2?a2:a1;
}
int Query_max(int p,int x1,int y1,int x2,int y2) {
    if(x1<=t[p].l&&t[p].r<=x2) {
        return t[p].Y.query_max(1,y1,y2);
    }
    int a1=0,a2=0;
    int mid=t[p].l+t[p].r>>1;
    if(x1<=mid)a1=Query_max(lc,x1,y1,x2,y2);
    if(x2>mid)a2=Query_max(rc,x1,y1,x2,y2);
    return a1<a2?a2:a1;
}
long long read() {
    long long x=0,f=1;
    char c=getchar();
    while(!isdigit(c)) {
        if(c=='-')f=-1;
        c=getchar();
    }
    while(isdigit(c)) {
        x=(x<<3)+(x<<1)+c-'0';
        c=getchar();
    }
    return x*f;
}
int main() {
    n=read(),m=read();
    for(int i=1; i<=n; i++) {
        for(int j=1; j<=m; j++) {
            g[i][j]=read();
            sum[i][j]=sum[i-1][j]+sum[i][j-1]-sum[i-1][j-1]+g[i][j];
        }
    }
    Build(1,1,n);
    int Q;
    Q=read();
    char opt[5];
    while(Q--) {
        scanf("%s",opt);
        if(opt[1]=='U') {
            int x1,x2,y1,y2;
            x1=read(),y1=read(),x2=read(),y2=read();
            x1++,x2++,y1++,y2++;
            cout<<sum[x2][y2]-sum[x1-1][y2]-sum[x2][y1-1]+sum[x1-1][y1-1]<<'
';
        }
        if(opt[1]=='A') {
            int x1,x2,y1,y2;
            x1=read(),y1=read(),x2=read(),y2=read();
            x1++,x2++,y1++,y2++;
            cout<<Query_max(1,x1,y1,x2,y2)<<'
';
        }
        if(opt[1]=='I') {
            int x1,x2,y1,y2;
            x1=read(),y1=read(),x2=read(),y2=read();
            x1++,x2++,y1++,y2++;
            cout<<Query_min(1,x1,y1,x2,y2)<<'
';
        }
    }
    return 0;
}

over