[IOI1998] Pictures

用线段树维护区间最小值和最小值个数来求一段区间里0的个数，把横的和竖的边分别拿出来，排序，然后每次查一下重复部分的长度即可

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#define ll long long
#define MAXN 5000
#define MAXL 10000
#define mid (T[k].l+T[k].r)/2
using namespace std;
inline int read()
{
    int x=,f=;char ch=getchar();
    while(ch<''||ch>''){if(ch=='-') f=-;ch=getchar();}
    while(ch>=''&&ch<=''){x=x*+ch-''; ch=getchar();}
    return x*f;
} 

int n,cnt=,ans=;
struct sq{
    int x1,y1,x2,y2;
}s[MAXN+];
struct line{
    int k,l,r,x;
}l[MAXL+];
struct TREE{
    int l,r,x,num,val;
}T[MAXL*];

bool cmp(line x,line y){return x.x<y.x||(x.x==y.x&&x.k<y.k);}

void pushdown(int k)
{
    int ls=k<<,rs=(k<<)+;
    T[ls].x+=T[k].val;T[rs].x+=T[k].val;
    T[ls].val+=T[k].val;T[rs].val+=T[k].val;
    T[k].val=;
}

void combine(int k)
{
    int ls=k<<,rs=(k<<)+;
    if(T[ls].x==T[rs].x) T[k].x=T[ls].x,T[k].num=T[ls].num+T[rs].num;
    else if(T[ls].x<T[rs].x) T[k].x=T[ls].x,T[k].num=T[ls].num;
    else T[k].x=T[rs].x,T[k].num=T[rs].num;
}

void renew(int k,int l,int r,int ad)
{
    if(T[k].l==l&&T[k].r==r)
    {T[k].x+=ad;T[k].val+=ad;return;}
    if(T[k].val) pushdown(k);
    if(r<=mid) renew(k<<,l,r,ad);
    else if(l>mid) renew((k<<)+,l,r,ad);
    else {renew(k<<,l,mid,ad);renew((k<<)+,mid+,r,ad);}
    combine(k);
}

int query(int k,int l,int r)
{
//  cout<<"q"<<k<<" "<<l<<" "<<r<<endl;
    if(T[k].l==l&&T[k].r==r)
    {
        if(T[k].x>) return T[k].r-T[k].l+;
        else return T[k].r-T[k].l+-T[k].num;
    }
    if(T[k].val) pushdown(k);
    if(r<=mid) return query(k<<,l,r);
    else if(l>mid) return query((k<<)+,l,r);
    else return query((k<<)+,mid+,r)+query(k<<,l,mid);
}

void init(int k,int l,int r)
{
    //cout<<"init"<<k<<" "<<r<<endl;
    T[k].l=l;T[k].r=r;T[k].val=;
    if(l==r){T[k].x=;T[k].num=;return;}
    init(k<<,l,mid);init((k<<)+,mid+,r);
    combine(k);
}

int main()
{
    n=read();
    for(int i=;i<=n;++i)
    {
        s[i].x1=read()+;s[i].y1=read()+;
        s[i].x2=read()+;s[i].y2=read()+;
        ans+=(s[i].x2+s[i].y2-s[i].y1-s[i].x1)*;
        l[++cnt].k=;l[cnt].x=s[i].x1;l[cnt].l=s[i].y1+;l[cnt].r=s[i].y2;
        l[++cnt].k=;l[cnt].x=s[i].x2;l[cnt].l=s[i].y1+;l[cnt].r=s[i].y2;
    }
    sort(l+,l+cnt+,cmp);init(,,);
//  for(int i=1;i<=cnt;i++)
//  {
//      cout<<l[i].k<<" "<<l[i].x<<" "<<l[i].l<<" "<<l[i].r<<endl;
//  }
    for(register int i=;i<=cnt;i++)
    {
        if(!l[i].k){ans-=query(,l[i].l,l[i].r);renew(,l[i].l,l[i].r,);}
        else       {renew(,l[i].l,l[i].r,-);ans-=query(,l[i].l,l[i].r);}

    }
    cnt=;
    for(register int i=;i<=n;++i)
    {
        l[++cnt].k=;l[cnt].x=s[i].y1;l[cnt].l=s[i].x1+;l[cnt].r=s[i].x2;
        l[++cnt].k=;l[cnt].x=s[i].y2;l[cnt].l=s[i].x1+;l[cnt].r=s[i].x2;
    }
    sort(l+,l+cnt+,cmp);init(,,);
//  for(int i=1;i<=cnt;i++)
//  {
//      cout<<l[i].k<<" "<<l[i].x<<" "<<l[i].l<<" "<<l[i].r<<endl;
//  }
    for(register int i=;i<=cnt;i++)
    {
        if(!l[i].k){ans-=query(,l[i].l,l[i].r);renew(,l[i].l,l[i].r,);}
        else       {renew(,l[i].l,l[i].r,-);ans-=query(,l[i].l,l[i].r);}
    }
    printf("%d",ans);
    return ;
}