【JZOJ A组】纪念碑

Description

2034年,纪念中学决定修建校庆100周年纪念碑,作为杰出校友的你被找了过来,帮校方确定纪念碑的选址.
纪念中学的土地可以看作是一个长为n,宽为m的矩形.它由n* m个1*1的正方形组成,其中左下角的正方形的坐标为(1,1),右上角的正方形的坐标为(n, m).其中有一些土地已经被用来修建建筑物,每一幢建筑物都可以看做是一个左下角为(x1,y1),右上角为(x2,y2)的矩形.
纪念碑可以看作是一个正方形.校方希望你找出一块最大的正方形区域供他们参考.

Input

每一组数据的第一行包含三个整数n,m和p,分别表示学校的长,宽以及建筑物的数量.
接下来的p行,每行包含四个整数x1,y1,x2,y2,分别表示每一幢建筑物左下角以及右上角的坐标.

Output

输出一个数,表示可能的最大边长.

Sample Input

13 5 8
8 4 10 4
4 3 4 4
10 2 12 2
8 2 8 4
2 4 6 4
10 3 10 4
12 3 12 4
2 2 4 2

Sample Output

3

Data Constraint

对于30%的数据,p<=1000.
对于70%的数据,p<=30000.
对于100%的数据,p<=400000,m,n<=1000000.

思路

现学现卖扫描线！！

维护两个指针l，r表示x轴选取的区间。
用线段树维护y轴的建筑物。
答案是min（建筑物的空隙，r-l+1）

可以发现r向右移答案只会减小不会增加，所以当答案<r-l+1时就移动左端点

代码

#include<iostream>
#include<cstdio>
#include<cstring>
#include<vector>
using namespace std;
const int maxn=1e6+77;
struct T
{
	int l,r,ls,rs,ass;
}tr[maxn*3];
int tag[maxn*3],n,m,p,x[maxn][3],y[maxn][3];
vector<int> ad[maxn],dl[maxn];
void pushup(int u)
{
	int LD=u*2,RD=u*2+1;
    if(tag[u]){tr[u].ls=tr[u].rs=tr[u].ass=0;return;}
    if(tr[u].l==tr[u].r){tr[u].ls=tr[u].rs=tr[u].ass=1;return;}
    tr[u].ls=tr[LD].ls+(tr[LD].ls==tr[LD].r-tr[LD].l+1)*tr[RD].ls;
    tr[u].rs=tr[RD].rs+(tr[RD].rs==tr[RD].r-tr[RD].l+1)*tr[LD].rs;
    tr[u].ass=max(tr[LD].rs+tr[RD].ls,max(tr[LD].ass,tr[RD].ass));
}
void build(int st,int ed,int u)
{
	tr[u].l=st; tr[u].r=ed;
	tr[u].ls=tr[u].rs=tr[u].ass=ed-st+1;
	if(st==ed) return;
	int mid=(st+ed)>>1;
	build(st,mid,u*2); build(mid+1,ed,u*2+1);
}
void ins(int st,int ed,int l,int r,int v,int u)
{
	if(l<=st&&ed<=r)
	{
		tag[u]+=v; pushup(u); return;
	}
    int mid=(st+ed)>>1;
    if(l<=mid) ins(st,mid,l,r,v,u*2);
    if(mid<r) ins(mid+1,ed,l,r,v,u*2+1);
    pushup(u);
}
int main()
{
	scanf("%d%d%d",&n,&m,&p);
	for(int i=1; i<=p; i++)
	{
		scanf("%d%d%d%d",&x[i][0],&y[i][0],&x[i][1],&y[i][1]);
		ad[x[i][0]].push_back(i);
		dl[x[i][1]].push_back(i);
	}
	build(1,m,1);
	int ans=0,l=1,r=1;
	while(r<=n)
	{
		for(int i=0; i<ad[r].size(); i++)
		{
			ins(1,m,y[ad[r][i]][0],y[ad[r][i]][1],1,1);
		}
		ans=max(ans,min(r-l+1,tr[1].ass));
		while(tr[1].ass<r-l+1)
		{
			for(int i=0; i<dl[l].size(); i++)
                ins(1,m,y[dl[l][i]][0],y[dl[l][i]][1],-1,1);
            l++;
		}
		r++;
	}
	printf("%d",ans);
}