西安 网站建设外包,微信小程序专业开发公司,网店推广工作内容,wordpress后台无法变中文传送门 文章目录题意#xff1a;思路#xff1a;题意#xff1a; 思路#xff1a;
考虑到最大面积是由四种颜色构成的#xff0c;且四种颜色可以从中心扩展出去#xff0c;所以我们分别维护四种颜色的二维前缀和#xff0c;O(1)O(1)O(1)计算矩阵内颜色的个数。现在我们…传送门
文章目录题意思路题意 思路
考虑到最大面积是由四种颜色构成的且四种颜色可以从中心扩展出去所以我们分别维护四种颜色的二维前缀和O(1)O(1)O(1)计算矩阵内颜色的个数。现在我们只需要找到矩阵左上角即红色的位置让后每次暴力扩展一个长度并且用二维前缀和检查是否合法。记f[i][j][k]f[i][j][k]f[i][j][k]表示以(i,j)(i,j)(i,j)为左上角所要求的矩阵边长一半为kkk的方案是否存在。这个显然我们刚才已经处理出来了。
考虑如何快速查询。
我们可以对fff也求一个前缀和让后从min(n,m)/2min(n,m)/2min(n,m)/2开始往小的跑直到f[x3][y3][k]−f[x1−1][y3][k]−f[x3][y1−1][k]f[x1−1][y1−1][k]0f[x3][y3][k]-f[x1-1][y3][k]-f[x3][y1-1][k]f[x1-1][y1-1][k]0f[x3][y3][k]−f[x1−1][y3][k]−f[x3][y1−1][k]f[x1−1][y1−1][k]0的时候kkk就是最大的了直接输出k∗k∗4k*k*4k∗k∗4即可。
要注意f[i][j][k]f[i][j][k]f[i][j][k]如果每个都开500500500的话会MLEMLEMLE但是其实我们第三维只需要开250250250就好了因为我们存的是边长一半。这个时候下面循环也需要改成min(n,m)/2min(n,m)/2min(n,m)/2不然越界会出现很多奇怪的问题。
复杂度O(qn)O(qn)O(qn)当然可以优化成O(qlogn)O(qlogn)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,tunenative)
//#pragma GCC optimize(2)
#includecstdio
#includeiostream
#includestring
#includecstring
#includemap
#includecmath
#includecctype
#includevector
#includeset
#includequeue
#includealgorithm
#includesstream
#includectime
#includecstdlib
#define X first
#define Y second
#define L (u1)
#define R (u1|1)
#define pb push_back
#define mk make_pair
#define Mid (tr[u].ltr[u].r1)
#define Len(u) (tr[u].r-tr[u].l1)
#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 pairint,int PII;const int N501,mod1e97,INF0x3f3f3f3f;
const double eps1e-6;int n,m,q;
int a[N][N][5],f[N][N][N/220];
char s[N][N];
int x[10],y[10];bool check(int x,int y)
{if(x1||xn||y1||ym) 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 i1;in;i) scanf(%s,s[i]1);for(int i1;in;i)for(int j1;jm;j){for(int k1;k4;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 i1;in;i)for(int j1;jm;j)if(s[i][j]Rs[i][j1]Gs[i1][j]Ys[i1][j1]B){f[i][j][1]1;for(int k2;kmin(n,m)/2;k){int lenk-1;x[1]i-len,y[1]j-len;x[2]i-len,y[2]jlen1;x[3]ilen1,y[3]j-len;x[4]ilen1,y[4]jlen1;if(get(x[1],y[1],i,j,1)!k*k) break;if(get(i-k1,j1,i,j1k-1,2)!k*k||!check(i-k1,j1)||!check(i,j1k-1)) break;if(get(i1,j-k1,i1k-1,j,3)!k*k||!check(i1,j-k1)||!check(i1k-1,j)) break;if(get(i1,j1,x[4],y[4],4)!k*k||!check(i1,j1)) break;f[x[1]][y[1]][k]1;}}//for(int i1;in;i) for(int j1;jm;j) printf(%d %d %d\n,i,j,f[i][j][2]);for(int i1;in;i)for(int j1;jm;j)for(int k1;kmin(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 flag0;for(int kmin(n,m)/2;k1;k--){if(x1kx2||y1ky2) continue;int x3x2-k*21,y3y2-k*21;if(x31||y31||x3x1||y3y1) 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,(kk)*(kk));flag1;break;}}if(!flag) puts(0);}return 0;
}
/**/
