考虑静态怎么做：枚举右边界，然后枚举上边界，对应的下边界一定单调不降，单调栈维护每一列从当前枚举的右边界向左最长空位的长度，这样是O(nm)的
注意到n>=m，所以m<=2000，可以枚举右边界，然后考虑怎么快速知道当前枚举的右边界向左最长空位的长度
用线段树维护行，每个节点都维护一段连续的列，p[ro][i]表示当i列从ro的区间最上面开始有多少行是全空的，q[ro][i]表示从下，v[ro][i]表示i列向左扩展最大的最大全空正方形的边长，大概是下面这种感觉：
然后p和q合并的时候类似HOTEL那题，看看左右儿子是否全空来决定直接继承还是加上另一段（因为不涉及其他列所以比较好写）
然后合并l行到r行的v的时候用单调递增的单调队列，分别维护左儿子的q和右儿子的p，从左到右扫，维护当前[l,r]里的最大正方形
询问的时候是按顺序把查询区间里的值都合并到一起，注意是按顺序！
#include<iostream>
#include<cstdio>
using namespace std;
const int N=4000005;
int n,m,Q,len[N<<2],ql[N],qr[N],ll,lr,rl,rr;
struct qwe
{
int f[N<<2];
int* operator [](int x)
{
return f+x*m;
}
}a,p,q,v;
int read()
{
int r=0,f=1;
char p=getchar();
while(p>'9'||p<'0')
{
if(p=='-')
f=-1;
p=getchar();
}
while(p>='0'&&p<='9')
{
r=r*10+p-48;
p=getchar();
}
return r*f;
}
void ud(int ro,int ls,int rs)
{
ll=rl=1,lr=rr=0;
for(int i=1,j=1;i<=m;i++)
{
while(rl<=rr&&q[ls][qr[rr]]>q[ls][i])
rr--;
qr[++rr]=i;
while(ll<=lr&&p[rs][ql[lr]]>p[rs][i])
lr--;
ql[++lr]=i;
while(ll<=lr&&rl<=rr&&q[ls][qr[rl]]+p[rs][ql[ll]]<i-j+1)
{
if(qr[rl]<=j)
rl++;
if(ql[ll]<=j)
ll++;
j++;
}
v[ro][i]=max(i-j+1,max(v[ls][i],v[rs][i]));
}
for(int i=1;i<=m;i++)
p[ro][i]=(p[ls][i]==len[ls])?len[ls]+p[rs][i]:p[ls][i];
for(int i=1;i<=m;i++)
q[ro][i]=(q[rs][i]==len[rs])?len[rs]+q[ls][i]:q[rs][i];
}
void build(int ro,int l,int r)
{
len[ro]=r-l+1;
if(l==r)
{
for(int i=1;i<=m;i++)
p[ro][i]=q[ro][i]=v[ro][i]=a[l][i];
return;
}
int mid=(l+r)>>1;
build(ro<<1,l,mid);
build(ro<<1|1,mid+1,r);
ud(ro,ro<<1,ro<<1|1);
}
void update(int ro,int l,int r,int x,int y)
{
if(l==r)
{
p[ro][y]=q[ro][y]=v[ro][y]=a[x][y];
return;
}
int mid=(l+r)>>1;
if(x<=mid)
update(ro<<1,l,mid,x,y);
else
update(ro<<1|1,mid+1,r,x,y);
ud(ro,ro<<1,ro<<1|1);
}
int hb(int la,int ro,int l,int r)
{
ll=rl=1,lr=rr=0;
int nw=0;
for(int i=l,j=l;i<=r;i++)
{
while(rl<=rr&&q[la][qr[rr]]>q[la][i])
rr--;
qr[++rr]=i;
while(ll<=lr&&p[ro][ql[lr]]>p[ro][i])
lr--;
ql[++lr]=i;
while(ll<=lr&&rl<=rr&&q[la][qr[rl]]+p[ro][ql[ll]]<i-j+1)
{
if(qr[rl]<=j)
rl++;
if(ql[ll]<=j)
ll++;
j++;
}
nw=max(nw,i-j+1);
}
for(int i=l;i<=r;i++)
p[la][i]=(p[la][i]==len[la])?len[la]+p[ro][i]:p[la][i];
for(int i=l;i<=r;i++)
q[la][i]=(q[ro][i]==len[ro])?len[ro]+q[la][i]:q[ro][i];
len[la]+=len[ro];
return nw;
}
int ques(int ro,int l,int r,int x,int xx,int y,int yy)
{
if(l==x&&r==xx)
{
int nw=hb(0,ro,y,yy);
for(int i=y;i<=yy;i++)
nw=max(nw,min(i-y+1,v[ro][i]));
return nw;
}
int mid=(l+r)>>1;
if(xx<=mid)
return ques(ro<<1,l,mid,x,xx,y,yy);
else if(x>mid)
return ques(ro<<1|1,mid+1,r,x,xx,y,yy);
else
{
int nw=ques(ro<<1,l,mid,x,mid,y,yy);
return max(nw,ques(ro<<1|1,mid+1,r,mid+1,xx,y,yy));//因为要求按从前往后的顺序合并所以不能直接max两个函数，max的时候处理顺序可能会反
}
}
int main()
{
n=read(),m=read(),Q=read();
for(int i=1;i<=n;i++)
for(int j=1;j<=m;j++)
a[i][j]=read();
build(1,1,n);
while(Q--)
{
int o=read();
if(o==0)
{
int x=read(),y=read();
a[x][y]^=1;
update(1,1,n,x,y);
}
else
{
int x=read(),y=read(),xx=read(),yy=read();
for(int i=1;i<=m;i++)
p[0][i]=q[0][i]=v[0][i]=0;
len[0]=0;
printf("%d
",ques(1,1,n,x,xx,y,yy));
}
}
return 0;
}