[LUOGU] P2716 和谐的雪花

https://www.luogu.org/problemnew/show/P2716

给出一个n*m的矩形，求里面边长最小的正方形，使得该正方形内最大值与最小值的差大于等于给定的K。

第一反应是二分答案+二维RMQ，时间复杂度是能过的，O(nlognlognlogn)，空间O(n*nlogn)。

（待补充）

考虑一个看起来不怎么优秀的O(n^3)做法，实际上是可以通过本题的。

设左上角在(i,j)的边长为k的正方形的最大值为f[i][j][k]，最小值g[i][j][k]，转移方程显然有f[i][j][k]=max{f[i][j][k-1],f[i+1][j][k-1],f[i][j+1][k-1],f[i+1][j+1][k-1]}，g同理。

如果在顺推时找到了一个f[i][j][k]-g[i][j][k]>=K，那么这个边长k一定是最小的了，

但是这样空间受限是吃不消的，注意到转移中第三维k是单调转移的，可以用适当的转移方式（也就是正常顺推）缩掉这一维。

空间复杂度O(n^2)，时间复杂度O(n^3)，卡常可不开O2通过。

#include<iostream>
#include<cstdio>
using namespace std;

inline int rd(){
  int ret=,f=;char c;
  while(c=getchar(),!isdigit(c))f=c=='-'?-:;
  while(isdigit(c))ret=ret*+c-'',c=getchar();
  return ret*f;
}

const int MAXN=;

int n,m,q;
int a[MAXN*MAXN];
int f[MAXN*MAXN],g[MAXN*MAXN];
int main(){
  n=rd();m=rd();q=rd();
  for(register int i=;i<=n;i++)
    for(register int j=;j<=m;j++)
      f[(i<<)+j]=g[(i<<)+j]=a[(i<<)+j]=rd();
  int up=min(n,m)-;
  for(register int k=;k<=up;k++)
    for(register int i=;i<=n-k;i++)
      for(register int j=;j<=m-k;j++){
        f[(i<<)+j]=max(f[(i<<)+j],max(f[(i+<<)+j],max(f[(i<<)+j+],f[(i+<<)+j+])));
        g[(i<<)+j]=min(g[(i<<)+j],min(g[(i+<<)+j],min(g[(i<<)+j+],g[(i+<<)+j+])));
        if(f[(i<<)+j]-g[(i<<)+j]>=q) return printf("%d\n",k+),;
      }
  puts("-1");
  return ;
}