题目大意
一个
n×m
的网格图,其中有的格子有障碍,有的没有。
有
T
个询问,每次询问一个子矩形
1≤n,m≤103,T≤106
题目分析
预处理每隔格子作为左上角能做的最大正方形边长。
一个很显然的想法(和[TJOI&HEOI2016]str/[JZOJ4614]字符串的想法类似):二分边长,转化为判定性问题,查询一个子矩形的最大值是否大于等于二分值。
这里
T
达到了
二维的
ST
表预处理实际上和一维是差不多的。在其中一维指数为
0
时转化为一维问题,否则直接让其中任意一个维一分为二即可。
查询类似一维情况,分成四个子矩形即可。
时间复杂度
代码实现
#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;
}