【JZOJ 5417】【NOIP2017提高A组集训10.24】方阵

Description

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

Solution

用分块+RMQ，
块的大小设为16，这样可以用倍增来预处理，
复杂度：$O(n^2\log(n)*4+16q)$

Code

#include <cstdio>
#include <algorithm>
#include <cmath>
#define fo(i,a,b) for(int i=a;i<=b;i++)
#define fod(i,a,b) for(int i=a;i>=b;i--)
#define min(q,w) ((q)>(w)?(w):(q))
#define max(q,w) ((q)<(w)?(w):(q))
using namespace std;
const int N=1000;
int read(int &n)
{
    char ch=' ';int q=0,w=1;
    for(;(ch!='-')&&((ch<'0')||(ch>'9'));ch=getchar());
    if(ch=='-')w=-1,ch=getchar();
    for(;ch>='0' && ch<='9';ch=getchar())q=q*10+ch-48;n=q*w;return n;
}
int m,n,ans,m1;
int a[N][N],sum[N][N];
int lg[N],er[11];
int MI[10][N][N],MX[10][N][N];
int MI1[10][N][N],MX1[10][N][N];
int main()
{
    freopen("phalanx.in","r",stdin);
    freopen("phalanx.out","w",stdout);
    int q,w,x,y,_;
    read(n),read(m);
    m1=16;
    fo(i,1,n)fo(j,1,m)
    {
        sum[i][j]=read(a[i][j])+sum[i-1][j]-sum[i-1][j-1]+sum[i][j-1];
        MX1[0][i][j]=MI1[0][i][j]=MI[0][i][j]=MX[0][i][j]=a[i][j];
    }
    fo(i,1,n)lg[i]=log2(i);
    er[1]=1;fo(i,2,10)er[i]=er[i-1]<<1;
    fo(k,1,9)fo(i,1,n-er[k])fo(j,1,m)
    {
        MI1[k][i][j]=MI[k][i][j]=min(MI[k-1][i][j],MI[k-1][i+er[k]][j]);
        MX1[k][i][j]=MX[k][i][j]=max(MX[k-1][i][j],MX[k-1][i+er[k]][j]);
    }
    fo(k,0,9)fo(l,1,4)fo(i,1,n-er[k])fo(j,1,m-er[l])
    {
        MI1[k][i][j]=min(MI1[k][i][j],MI1[k][i][j+er[l]]);
        MX1[k][i][j]=max(MX1[k][i][j],MX1[k][i][j+er[l]]);
    }
    read(_);
    fo(I,1,_)
    {
        char ch=' ';
        for(;ch<'A'||ch>'Z';ch=getchar());
        ch=getchar();
        read(q),read(w),read(x),read(y);
        q++,w++,x++,y++;
        if(ch=='U')
        {
            ans=sum[x][y]-sum[x][w-1]+sum[q-1][w-1]-sum[q-1][y];
        }else if(ch=='I')
        {
            ans=1e9;
            int e=lg[x-q+1];
            int t=x+1-er[e+1],i;
            for(i=w;i+m1<=y;i+=m1)ans=min(ans,min(MI1[e][q][i],MI1[e][t][i]));
            for(;i<=y;i++)ans=min(ans,min(MI[e][q][i],MI[e][t][i]));
        }else 
        {
            ans=-1;
            int e=lg[x-q+1];
            int t=x+1-er[e+1],i;
            for(i=w;i+m1<=y;i+=m1)ans=max(ans,max(MX1[e][q][i],MX1[e][t][i]));
            for(;i<=y;i++)ans=max(ans,max(MX[e][q][i],MX[e][t][i]));
        }
        printf("%d\n",ans);
    }
    return 0;
}