【JZOJ 4238】纪念碑

Description

给你一个大矩形空地，其中有很多小矩形障碍在里面（有重叠），让你在里面找一个最大的正方形的边长（里面不包含任何小矩形）
对于100%的数据,p<=400000,m,n<=1000000.

Analysis

经典的扫描线题目
一般一个坐标系里有些矩形，就考虑扫描线了
一条扫描线从左到右扫过去，中途会遇到矩形的左右边界
x轴上维护两个指针l,r，表示正方形的左右边界
对于正方形，其x轴上的长度就是r-l+1
y轴上的长度怎么算？显然是当前y轴上没被覆盖的最大长度
对于一个矩形，设其左右边界为x1,x2，上下边界为y1,y2
当且仅当r>=x1且l<=x2时，y轴上y1~y2会被覆盖
那么用线段树维护覆盖情况，扫描的时候就区间修改
由于加入与删除的区间一一对应，所以可以标记不下传

Code

#include <ctime>
#include <cstdio>
#include <cstring>
#include <algorithm>
#define fo(i,a,b) for(int i=a;i<=b;i++)
#define fd(i,a,b) for(int i=a;i>=b;i--)
#define efo(i,v) for(int i=last[v];i;i=next[i])
using namespace std;
typedef long long ll;
const int P=800005,N=1000010;
int n,m,num,p,tr[N*4],mx[N*4],le[N*4],ri[N*4];
struct line
{
    int x,y1,y2,z;
}a[P];
void read(int &t)
{
    t=0;char ch;
    for(ch=getchar();ch<'0' || ch>'9';ch=getchar());
    for(;'0'<=ch && ch<='9';ch=getchar()) t=t*10+ch-'0';
}
bool cmp(line a,line b)
{
    return a.x<b.x;
}
void build(int v,int l,int r)
{
    mx[v]=le[v]=ri[v]=r-l+1;
    if(l==r) return;
    int mid=(l+r)>>1;
    build(v+v,l,mid);build(v+v+1,mid+1,r);
}
void up(int v,int l,int r)
{
    int mid=(l+r)>>1;
    le[v]=le[v+v];
    if(le[v]==mid-l+1) le[v]+=le[v+v+1];
    ri[v]=ri[v+v+1];
    if(ri[v]==r-mid) ri[v]+=ri[v+v];
    mx[v]=max(mx[v+v],mx[v+v+1]);
    mx[v]=max(mx[v],ri[v+v]+le[v+v+1]);
}
void change(int v,int l,int r,int x,int y,int z)
{
    if(l==x && r==y)
    {
        tr[v]+=z;
        if(tr[v]) mx[v]=le[v]=ri[v]=0;
        else
        {
            if(l==r) mx[v]=le[v]=ri[v]=1;
            else up(v,l,r);
        }
        return;
    }
    int mid=(l+r)>>1;
    if(y<=mid) change(v+v,l,mid,x,y,z);
    else
    if(x>mid) change(v+v+1,mid+1,r,x,y,z);
    else change(v+v,l,mid,x,mid,z),change(v+v+1,mid+1,r,mid+1,y,z);
    if(tr[v]) mx[v]=le[v]=ri[v]=0;
    else up(v,l,r);
}
int main()
{
    int x1,y1,x2,y2;
    read(n),read(m),read(p);
    fo(i,1,p)
    {
        read(x1),read(y1),read(x2),read(y2);
        a[++num].x=x1,a[num].y1=y1,a[num].y2=y2,a[num].z=1;
        a[++num].x=x2+1,a[num].y1=y1,a[num].y2=y2,a[num].z=-1;
    }
    sort(a+1,a+num+1,cmp);
    build(1,1,m);
    int l=1,r=1,pl=1,pr=1,ans=0;
    for(;a[pl].x==l;pl++)
        if(a[pl].z==1) change(1,1,m,a[pl].y1,a[pl].y2,a[pl].z);
    for(;a[pr].x==r;pr++)
        if(a[pr].z==-1) change(1,1,m,a[pr].y1,a[pr].y2,a[pr].z);
    while(r<=n)
    {
        ans=max(ans,min(r-l+1,mx[1]));
        if(mx[1]>=r-l+1)
        {
            ++r;
            for(;a[pr].x==r;pr++)
                if(a[pr].z==1) change(1,1,m,a[pr].y1,a[pr].y2,a[pr].z);
        }
        else
        {
            ++l;
            for(;a[pl].x==l;pl++)
                if(a[pl].z==-1) change(1,1,m,a[pl].y1,a[pl].y2,a[pl].z);
        }
    }
    printf("%d\n",ans);
    return 0;
}