题解：奶牛吃草

题目大意

给你一个01矩阵，让你维护两个操作
0: 修改一个点
1: 求指定矩形内最大的全1正方形

矩阵大小$size\le 4\times 10^6$，操作数$q\le 2\times 10^3$

闲话

这题部分分对最终方法的求解无太多帮助，故省略

题解

这题看上去挺难 ，实际也挺难。因为要求指定矩形。鉴于数字较小（但是事实上这个数字基本不能再大了），可以带个size再带个log

考虑线段树，求解上下区间中，以$i$为右边界最大的正方形
这样除去上传以外，一切都简单了。
上传中，考虑对于分界线$mid$左右长度的$l$与$r$数组，可以发现，当$j$为左边界，只要满足$min(l[k])+min(r[p])>=i-j+1$且左侧最小即可为最优解（想一想，为什么）
于是题目就解出来了。

题解确实很简单，听着却很容易听不懂。建议与代码一起食用

#include<bits/stdc++.h>
using namespace std;
const int N=2007,M=1.7e7+5;
int n,m,q;
struct node{
	int nm[M];
	int* operator[](int x){
		return nm+x*m;
	}
}mx,sl,sr,a;
int h1,t1,h2,t2,q1[N],q2[N];
void pushup(int rt,int x,int y,int px,int py){
	h1=h2=1;
	t1=t2=0;
	int j=1;
	for(int i=1;i<=m;i++){
		while(h1<=t1&&sr[x][q1[t1]]>sr[x][i])t1--;
		q1[++t1]=i;
		while(h2<=t2&&sl[y][q2[t2]]>sl[y][i])t2--;
		q2[++t2]=i;
		while(j<=i&&sr[x][q1[h1]]+sl[y][q2[h2]]<i-j+1){
			j++;
			if(q1[h1]<j)h1++;
			if(q2[h2]<j)h2++;
		}
		mx[rt][i]=max(i-j+1,max(mx[x][i],mx[y][i]));
	}
	for(int i=1;i<=m;i++){
		sl[rt][i]=(sl[x][i]==px)?px+sl[y][i]:sl[x][i];
		sr[rt][i]=(sr[y][i]==py)?py+sr[x][i]:sr[y][i];
	}
}
void build(int id,int l,int r){
	if(l==r){
		for(int i=1;i<=m;i++)
			sl[id][i]=sr[id][i]=mx[id][i]=a[l][i];
		return;
	}
	int mid=l+r>>1;
	build(id*2,l,mid);
	build(id*2+1,mid+1,r);
	pushup(id,id*2,id*2+1,mid-l+1,r-mid);
}
void change(int id,int l,int r,int x,int y){
	if(l==r){
		a[x][y]^=1;
		sl[id][y]=sr[id][y]=mx[id][y]=a[x][y];
		return;
	}
	int mid=l+r>>1;
	if(x<=mid)
		change(id*2,l,mid,x,y);
	else
		change(id*2+1,mid+1,r,x,y);
	pushup(id,id*2,id*2+1,mid-l+1,r-mid);
}
void find(int id,int l,int r,int x,int y){
	if(x<=l&&r<=y){
		pushup(0,0,id,l-x,r-l+1);
		return;
	}
	int mid=l+r>>1;
	if(x<=mid)
		find(id*2,l,mid,x,y);
	if(y>mid)
		find(id*2+1,mid+1,r,x,y);	
}
main(){
	freopen("car.in","r",stdin);
	freopen("car.out","w",stdout);
	cin>>n>>m>>q;
	for(int i=1;i<=n;i++)
		for(int j=1;j<=m;j++)
			cin>>a[i][j];
	build(1,1,n);
	while(q--){
		int op,x,y,u,v;
		cin>>op>>x>>y;
		if(op==0)
			change(1,1,n,x,y);
		else{
			cin>>u>>v;
			for(int i=1;i<=m;i++)
				mx[0][i]=sl[0][i]=sr[0][i]=0;
			find(1,1,n,x,u);
			int ans=0;
			for(int i=y;i<=v;i++)
				ans=max(ans,min(mx[0][i],i-y+1));
			cout<<ans<<endl;
		}
	}
}//真的很长，能看懂真的很不错，这题难度对标S T4