[JZOJ4817]square

题目大意

一个$n\times m$的网格图，其中有的格子有障碍，有的没有。
有$T$个询问，每次询问一个子矩形$(x_1,y_1)-(x_2,y_2)$。查询该子矩形内边长最大的格点正方形边长为多少。一个正方形不能包括有障碍的格子。

$1\le n,m\le10^3,T\le10^6$

---

题目分析

预处理每隔格子作为左上角能做的最大正方形边长。
一个很显然的想法（和[TJOI&HEOI2016]str/[JZOJ4614]字符串的想法类似）：二分边长，转化为判定性问题，查询一个子矩形的最大值是否大于等于二分值。
这里$T$达到了$10^6$规模，显然我们需要一个查询$\mathrm O(1)$的二维$\mathrm{RMQ}$算法。二维的$\mathrm{ST}$表可以胜任这一工作。
二维的$\mathrm{ST}$表预处理实际上和一维是差不多的。在其中一维指数为$0$时转化为一维问题，否则直接让其中任意一个维一分为二即可。
查询类似一维情况，分成四个子矩形即可。
时间复杂度$\mathrm O(nm\log_2n\log_2m+T\log_2n)$，卡卡常就过了。

---

代码实现

#include <iostream>
#include <cstdio>
#include <cctype>
#include <cmath>

using namespace std;

typedef pair<int,int> P;
#define ft first
#define sd second
#define mkp(a,b) make_pair(a,b)

inline int read()
{
    int x=0,f=1;
    char ch=getchar();
    while (!isdigit(ch)) f=ch=='-'?-1:f,ch=getchar();
    while (isdigit(ch)) x=x*10+ch-'0',ch=getchar();
    return x*f;
}

int buf[30];

inline void write(int x)
{
    if (x<0) putchar('-'),x=-x;
    buf[0]=0;
    while (x) buf[++buf[0]]=x%10,x/=10;
    if (!buf[0]) buf[0]=1,buf[1]=0;
    while (buf[0]) putchar('0'+buf[buf[0]--]);
}

const int N=1050;
const int LGN=10;

int down[N][N],Right[N][N],dwl[N][N],rgl[N][N];
int n,m,T,lgn,head,tail;
int f[LGN][LGN][N][N];
bool a[N][N];
int pw[LGN];
int LOG[N];
P que[N];

void add(P x)
{
    while (head<=tail&&que[tail].ft>=x.ft) tail--;
    que[++tail]=x;
}

int getmin(int x)
{
    while (head<=tail&&que[head].sd>x) head++;
    return que[head].ft;
}

void pre()
{
    for (int i=1;i<=n;i++)
        for (int j=m;j>=1;j--)
            Right[i][j]=a[i][j]?Right[i][j+1]+1:0;
    for (int j=1;j<=m;j++)
        for (int i=n;i>=1;i--)
            down[i][j]=a[i][j]?down[i+1][j]+1:0;
    for (int i=1;i<=n;i++)
    {
        head=1,tail=0;
        for (int cur=m,j=m;j>=1;j--)
        {
            add(mkp(down[i][j],j));
            while (cur>=j&&getmin(cur)<cur-j+1) cur--;
            rgl[i][j]=cur;
        }
    }
    for (int j=1;j<=m;j++)
    {
        head=1,tail=0;
        for (int cur=n,i=n;i>=1;i--)
        {
            add(mkp(Right[i][j],i));
            while (cur>=i&&getmin(cur)<cur-i+1) cur--;
            dwl[i][j]=cur;
        }
    }
    lgn=trunc(log(n)/log(2));
    pw[0]=1;
    for (int i=1;i<=lgn;i++) pw[i]=pw[i-1]<<1;
    for (int i=1;i<=n;i++) LOG[i]=trunc(log(i)/log(2));
    for (int i=1;i<=n;i++)
        for (int j=1;j<=m;j++)
            f[0][0][i][j]=min(rgl[i][j]-j+1,dwl[i][j]-i+1);
    for (int k=0;k<=lgn;k++)
        for (int l=0;l<=lgn;l++)
            if (k||l)
                for (int i=1;i+pw[k]-1<=n;i++)
                    for (int j=1;j+pw[l]-1<=m;j++)
                        if (!k) f[k][l][i][j]=max(f[k][l-1][i][j],f[k][l-1][i][j+pw[l-1]]);
                        else if (!l) f[k][l][i][j]=max(f[k-1][l][i][j],f[k-1][l][i+pw[k-1]][j]);
                        else f[k][l][i][j]=max(f[k][l-1][i][j],f[k][l-1][i][j+pw[l-1]]);
}

int query(int x1,int y1,int x2,int y2)
{
    int lg1=LOG[x2-x1+1],lg2=LOG[y2-y1+1];
    return max(max(f[lg1][lg2][x1][y1],f[lg1][lg2][x2-pw[lg1]+1][y2-pw[lg2]+1]),max(f[lg1][lg2][x2-pw[lg1]+1][y1],f[lg1][lg2][x1][y2-pw[lg2]+1]));
}

int solve(int x1,int y1,int x2,int y2)
{
    int l=0,r=min(x2-x1+1,y2-y1+1),mid,ret=0;
    while (l<=r)
    {
        mid=l+r>>1;
        if (query(x1,y1,x2-mid+1,y2-mid+1)>=mid) l=(ret=mid)+1;
        else r=mid-1;
    }
    return ret;
}

int main()
{
    freopen("square.in","r",stdin),freopen("square.out","w",stdout);
    n=read(),m=read();
    for (int i=1;i<=n;i++)
        for (int j=1;j<=m;j++)
            a[i][j]=read();
    pre();
    for (T=read();T--;)
    {
        int x1=read(),y1=read(),x2=read(),y2=read();
        write(solve(x1,y1,x2,y2)),putchar('\n');
    }
    fclose(stdin),fclose(stdout);
    return 0;
}