BZOJ3356 : [Usaco2004 Jan]禁闭围栏

首先将坐标离散化，考虑从左往右扫描线

碰到插入操作则插入

碰到删除操作的：

当前包含i的矩形数=y1在[1,y2[i]]之间的矩形数-y2在[1,y1[i]-1]之间的矩形数

用两棵树状数组维护即可，时间复杂度$O(n\log n)$。

 

#include<cstdio>

#include<algorithm>

#define N 500010

int n,m,i,x1,y1,x2,y2,b[N],bl[N],br[N],now,ans=-1,cnt;

struct P{int x,l,r,t;}a[N];

inline bool cmp(P a,P b){return a.x<b.x;}

inline void addl(int x,int y){for(;x<=m;x+=x&-x)bl[x]+=y;}

inline void addr(int x,int y){for(;x<=m;x+=x&-x)br[x]+=y;}

inline int askl(int x){int t=0;for(;x;x-=x&-x)t+=bl[x];return t;}

inline int askr(int x){int t=0;for(;x;x-=x&-x)t+=br[x];return t;}

inline int lower(int x){

  int l=1,r=m,mid,t;

  while(l<=r)if(b[mid=(l+r)>>1]<=x)l=(t=mid)+1;else r=mid-1;

  return t;

}

inline void read(int&a){char c;while(!(((c=getchar())>='0')&&(c<='9')));a=c-'0';while(((c=getchar())>='0')&&(c<='9'))(a*=10)+=c-'0';}

int main(){

  read(n);

  while(n--){

    read(x1),read(y1),read(x2),read(y2);

    a[++m].x=x1,a[m].l=b[m]=y1,a[m].r=y2,a[m].t=1;

    a[++m].x=x2,a[m].l=y1,a[m].r=b[m]=y2;

  }

  for(std::sort(b+1,b+m+1),std::sort(a+1,a+m+1,cmp),i=1;i<=m;i++){

    a[i].l=lower(a[i].l),a[i].r=lower(a[i].r);

    if(a[i].t)addl(a[i].l,1),addr(a[i].r,1);else{

      now=askl(a[i].r)-askr(a[i].l-1);

      if(now>ans)ans=now,cnt=1;else if(now==ans)cnt++;

      addl(a[i].l,-1),addr(a[i].r,-1);

    }

  }

  return printf("%d %d",ans,cnt),0;

}