JZOJ5417. 【NOIP2017提高A组集训10.24】 二维线段树

时间:2022-03-17 20:47:57

题意:给一个矩阵,q次询问子矩阵,查询最大值最小值,区间和, 每个询问的方阵的长不超过宽的两倍。

屁股后面那个两倍什么玩意儿的你想用就用,我不想用,懒得打什么鬼畜的rmq,去他的,二维线段树多爽,直接过。
然后我比赛把前缀和打错了啊啊啊啊啊!!!
/喷血 犯了史上最傻逼错误没有之一,更让我绝望的是,我打拍的时候直接把前缀和那个地方黏过去用所以根本没拍出来/喷血。
血的教训,希望我能铭记。

#include<cstdio>
#include<algorithm>
#include<cstring>
#define fo(i,a,b) for(int i=a;i<=b;i++)
#define fd(i,a,b) for(int i=a;i>=b;i--)
using namespace std;
const int N=805;
const int M=1e5+5;
const int inf=0x3f3f3f3f;
int n,m,q;
int mx[N<<2][N<<2],mn[N<<2][N<<2];
typedef long long ll;
ll sum[N][N];
int a[N][N],maxx,minx,X1,X2,Y1,Y2   ;
inline void buildy(int x,int l,int r,int rtx,int rty)
{
    if (l==r)
    {
        if (x!=-1)
            mx[rtx][rty]=mn[rtx][rty]=a[x][l]; 
        else
        {
            mx[rtx][rty]=max(mx[rtx<<1][rty],mx[rtx<<1|1][rty]);  
            mn[rtx][rty]=min(mn[rtx<<1][rty],mn[rtx<<1|1][rty]);  
        }
        return;
    }
    int mid=(l+r)>>1;
    buildy(x,l,mid,rtx,rty<<1);
    buildy(x,mid+1,r,rtx,rty<<1|1);
    mx[rtx][rty]=max(mx[rtx][rty<<1],mx[rtx][rty<<1|1]);  
    mn[rtx][rty]=min(mn[rtx][rty<<1],mn[rtx][rty<<1|1]);  
}
inline void buildx(int l,int r,int rtx)
{
    if (l==r)
    {
        buildy(l,1,m,rtx,1);
        return;
    }
    int mid=(l+r)>>1;
    buildx(l,mid,rtx<<1);
    buildx(mid+1,r,rtx<<1|1);
    buildy(-1,1,m,rtx,1);
}
inline void queryy(int l,int r,int rtx,int rty)
{
    if (Y1<=l&&r<=Y2)
    {
        maxx=max(maxx,mx[rtx][rty]);
        minx=min(minx,mn[rtx][rty]);
        return;
    }
    int mid=(l+r)>>1;
    if (Y1<=mid)queryy(l,mid,rtx,rty<<1);
    if (Y2>mid)queryy(mid+1,r,rtx,rty<<1|1);
}
inline void queryx(int l,int r,int rtx)
{
    if(X1<=l&&r<=X2)
    {
        queryy(1,m,rtx,1);
        return;
    }
    int mid=(l+r)>>1;
    if (X1<=mid)queryx(l,mid,rtx<<1);
    if (X2>mid)queryx(mid+1,r,rtx<<1|1);
}
int main()
{
    freopen("phalanx.in","r",stdin);
    freopen("phalanx.out","w",stdout);
    scanf("%d%d",&n,&m);
    fo(i,1,n)
    {
        fo(j,1,m)
        {
            scanf("%d",&a[i][j]);
        }
    }
    buildx(1,n,1);
    fo(i,1,n)
    {
        fo(j,1,m)
        sum[i][j]=sum[i-1][j]+sum[i][j-1]-sum[i-1][j-1]+a[i][j];
    }
    scanf("%d",&q);
    fo(i,1,q)
    {
        char ch[5];
        scanf("%s%d%d%d%d",ch,&X1,&Y1,&X2,&Y2);
        X1++,Y1++,X2++,Y2++;
        if (ch[0]=='S')
        {
            printf("%lld\n",sum[X2][Y2]-sum[X2][Y1-1]-sum[X1-1][Y2]+sum[X1-1][Y1-1]);
        }
        else
        {
            int op=0;
            if (ch[1]=='A')op=1;
            if (ch[1]=='I')op=2;
            maxx=-inf,minx=inf;
            queryx(1,n,1);
            if (op==1)
            printf("%d\n",maxx);
            else printf("%d\n",minx);
        }
    }
}