传送门

题意：

思路：

考虑到最大面积是由四种颜色构成的，且四种颜色可以从中心扩展出去，所以我们分别维护四种颜色的二维前缀和，$O(1)$计算矩阵内颜色的个数。现在我们只需要找到矩阵左上角，即红色的位置，让后每次暴力扩展一个长度，并且用二维前缀和检查是否合法。记$f[i][j][k]$表示以$(i,j)$为左上角，所要求的矩阵边长一半为$k$的方案是否存在。这个显然我们刚才已经处理出来了。

考虑如何快速查询。

我们可以对$f$也求一个前缀和，让后从$min(n,m)/2$开始往小的跑，直到$f[x3][y3][k]-f[x1-1][y3][k]-f[x3][y1-1][k]+f[x1-1][y1-1][k]>0$的时候，$k$就是最大的了，直接输出$k\ast k\ast 4$即可。

要注意$f[i][j][k]$如果每个都开$500$的话会$MLE$，但是其实我们第三维只需要开$250$就好了，因为我们存的是边长一半。这个时候下面循环也需要改成$min(n,m)/2$，不然越界会出现很多奇怪的问题。

复杂度$O(qn)$，当然可以优化成$O(qlogn)$，因为其满足单调的性质。

// Problem: E. Nanosoft
// Contest: Codeforces - Codeforces Round #619 (Div. 2)
// URL: https://codeforces.com/contest/1301/problem/E
// Memory Limit: 512 MB
// Time Limit: 2000 ms
// 
// Powered by CP Editor (https://cpeditor.org)//#pragma GCC optimize("Ofast,no-stack-protector,unroll-loops,fast-math")
//#pragma GCC target("sse,sse2,sse3,ssse3,sse4.1,sse4.2,avx,avx2,popcnt,tune=native")
//#pragma GCC optimize(2)
#include<cstdio>
#include<iostream>
#include<string>
#include<cstring>
#include<map>
#include<cmath>
#include<cctype>
#include<vector>
#include<set>
#include<queue>
#include<algorithm>
#include<sstream>
#include<ctime>
#include<cstdlib>
#define X first
#define Y second
#define L (u<<1)
#define R (u<<1|1)
#define pb push_back
#define mk make_pair
#define Mid (tr[u].l+tr[u].r>>1)
#define Len(u) (tr[u].r-tr[u].l+1)
#define random(a,b) ((a)+rand()%((b)-(a)+1))
#define db puts("---")
using namespace std;//void rd_cre() { freopen("d://dp//data.txt","w",stdout); srand(time(NULL)); }
//void rd_ac() { freopen("d://dp//data.txt","r",stdin); freopen("d://dp//AC.txt","w",stdout); }
//void rd_wa() { freopen("d://dp//data.txt","r",stdin); freopen("d://dp//WA.txt","w",stdout); }typedef long long LL;
typedef unsigned long long ULL;
typedef pair<int,int> PII;const int N=501,mod=1e9+7,INF=0x3f3f3f3f;
const double eps=1e-6;int n,m,q;
int a[N][N][5],f[N][N][N/2+20];
char s[N][N];
int x[10],y[10];bool check(int x,int y)
{if(x<1||x>n||y<1||y>m) return false;return true;
}int get(int x1,int y1,int x2,int y2,int k)
{return a[x2][y2][k]-a[x1-1][y2][k]-a[x2][y1-1][k]+a[x1-1][y1-1][k];
}int main()
{
//	ios::sync_with_stdio(false);
//	cin.tie(0);scanf("%d%d%d",&n,&m,&q);for(int i=1;i<=n;i++) scanf("%s",s[i]+1);for(int i=1;i<=n;i++)for(int j=1;j<=m;j++){for(int k=1;k<=4;k++) a[i][j][k]=a[i-1][j][k]+a[i][j-1][k]-a[i-1][j-1][k];a[i][j][1]+=s[i][j]=='R';a[i][j][2]+=s[i][j]=='G';a[i][j][3]+=s[i][j]=='Y';a[i][j][4]+=s[i][j]=='B';}for(int i=1;i<=n;i++)for(int j=1;j<=m;j++)if(s[i][j]=='R'&&s[i][j+1]=='G'&&s[i+1][j]=='Y'&&s[i+1][j+1]=='B'){f[i][j][1]=1;for(int k=2;k<=min(n,m)/2;k++){int len=k-1;x[1]=i-len,y[1]=j-len;x[2]=i-len,y[2]=j+len+1;x[3]=i+len+1,y[3]=j-len;x[4]=i+len+1,y[4]=j+len+1;if(get(x[1],y[1],i,j,1)!=k*k) break;if(get(i-k+1,j+1,i,j+1+k-1,2)!=k*k||!check(i-k+1,j+1)||!check(i,j+1+k-1)) break;if(get(i+1,j-k+1,i+1+k-1,j,3)!=k*k||!check(i+1,j-k+1)||!check(i+1+k-1,j)) break;if(get(i+1,j+1,x[4],y[4],4)!=k*k||!check(i+1,j+1)) break;f[x[1]][y[1]][k]=1;}}//for(int i=1;i<=n;i++) for(int j=1;j<=m;j++) printf("%d %d %d\n",i,j,f[i][j][2]);for(int i=1;i<=n;i++)for(int j=1;j<=m;j++)for(int k=1;k<=min(n,m)/2;k++)	f[i][j][k]+=f[i-1][j][k]+f[i][j-1][k]-f[i-1][j-1][k];while(q--){int x1,y1,x2,y2;scanf("%d%d%d%d",&x1,&y1,&x2,&y2);int flag=0;for(int k=min(n,m)/2;k>=1;k--){if(x1+k>x2||y1+k>y2) continue;int x3=x2-k*2+1,y3=y2-k*2+1;if(x3<1||y3<1||x3<x1||y3<y1) continue;if(f[x3][y3][k]-f[x1-1][y3][k]-f[x3][y1-1][k]+f[x1-1][y1-1][k]>0){printf("%d\n",(k+k)*(k+k));flag=1;break;}}if(!flag) puts("0");}return 0;
}
/**/