Codeforces713D（二维RMQ）

区间最大子正方形问题。
预处理后二分答案验证

AC Code：

#include<bits/stdc++.h>
#define maxn 1005
#define lim 10
using namespace std;

int n,m,dp[maxn][maxn],s[maxn][maxn],f[10][10][maxn][maxn],lg[1<<10];
inline int Query(int a,int b,int c,int d)
{
	if(a>c || b>d) return 0;
	int l1 = lg[c - a + 1]  ,l2=lg[d - b + 1];
	return max(
		max(f[l1][l2][a][b],f[l1][l2][c-(1<<l1)+1][b]),
		max(f[l1][l2][a][d-(1<<l2)+1],f[l1][l2][c-(1<<l1)+1][d-(1<<l2)+1])
	);
}

int main()
{
	scanf("%d%d",&n,&m);
	for(int i=1;i<=n;i++)
		for(int j=1;j<=m;j++)
			scanf("%d",&s[i][j]);
	for(int i=1;i<=n;i++)
		for(int j=1;j<=m;j++)
			if(s[i][j])
				if(i==1 || j==1) dp[i][j] = 1;
				else dp[i][j] = min(dp[i-1][j],min(dp[i][j-1],dp[i-1][j-1]))+1;
	for(int a=0;a<lim;a++)
		for(int b=0;b<lim;b++)
			for(int i=1;i<=n;i++)
				for(int j=1;j<=m;j++)
				{
					if(a==0 && b==0) f[a][b][i][j] = dp[i][j];
					else if(a==0)
					{
						int v = j + (1<<b-1);
						f[a][b][i][j] = max(f[a][b-1][i][j],f[a][b-1][i][v]);
					}
					else if(b==0)
					{
						int u = i + (1<<a-1);
						f[a][b][i][j] = max(f[a-1][b][i][j],f[a-1][b][u][j]);
					}
					else 
					{
						int u = i + (1<<a-1) , v = j + (1<<b-1);
						f[a][b][i][j] = max(
							max(f[a-1][b-1][i][j],f[a-1][b-1][u][j]) , 
							max(f[a-1][b-1][i][v],f[a-1][b-1][u][v])
						);
					}
				}
	for(int i=2;i<maxn;i++) lg[i] = lg[i>>1] + 1;
		
	int q,a,b,c,d;
	for(scanf("%d",&q);q--;)
	{
		scanf("%d%d%d%d",&a,&b,&c,&d);
		int L = 0 , R = min(c-a+1,d-b+1) , mid;
		for(;L<R;)
		{	mid = (L+R+1) >> 1;
			if(Query(a+mid-1,b+mid-1,c,d) >= mid) L = mid;
			else R = mid - 1;
		}
		printf("%d\n",L);
	}
}