hdu1255 覆盖的面积(线段树+扫描线+离散化，求矩形面积并)

给定平面上若干矩形,求出被这些矩形覆盖过至少两次的区域的面积。和周长并套路差不多。。甚至简单了些我觉得。重点是维护重复了两次以上的线段长度。分类讨论即可。
以coverlen[2]为例：
如果cover==0 则coverlen[2]利用子树的coverlen[2]来更新
如果cover==1 则coverlen[2]利用子树的coverlen[1]来更新
如果cover==2 则coverlen[2]利用子树的coverlen[0]来更新
其中coverlen[0]就是每个节点的总长度。
注意一下叶子节点特殊处理，不能再往下找儿子了。
tips:MLE N次，以为是写的太丑内存爆了？结果是数组开小了？？？无奈.jpg

#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std;
#define N 2010
struct Scanline{
    double x,y1,y2;bool f;
}line[N];
inline bool cmp(Scanline x,Scanline y){
    return x.x==y.x?x.f>y.f:x.x<y.x;
}
struct node{
    int cover;double coverlen[3];//cover[i]--覆盖了i次及以上的长度 
}tree[N<<2];
int tst,n;
double a[N];
inline void build(int p,int l,int r){
    memset(tree[p].coverlen,0,sizeof(tree[p].coverlen));
    tree[p].cover=0;
    if(l==r){tree[p].coverlen[0]=a[l+1]-a[l];return;}
    int mid=l+r>>1;
    build(p<<1,l,mid);build(p<<1|1,mid+1,r);
    tree[p].coverlen[0]=tree[p<<1].coverlen[0]+tree[p<<1|1].coverlen[0];
}
inline void pushup(int p){
    if(tree[p].cover>=2){
        tree[p].coverlen[2]=tree[p].coverlen[0];
        tree[p].coverlen[1]=tree[p].coverlen[0];return;
    }
    else if(tree[p].cover==1){
        tree[p].coverlen[1]=tree[p].coverlen[0];
        tree[p].coverlen[2]=tree[p<<1].coverlen[1]+tree[p<<1|1].coverlen[1];
        return;
    }
    else{
        tree[p].coverlen[2]=tree[p<<1].coverlen[2]+tree[p<<1|1].coverlen[2];
        tree[p].coverlen[1]=tree[p<<1].coverlen[1]+tree[p<<1|1].coverlen[1];
    }
}
inline void update(int p,int l,int r,int x,int y,int val){
    if(x<=l&&r<=y){
        tree[p].cover+=val;
        if(l==r){
            tree[p].coverlen[1]= tree[p].cover>=1?tree[p].coverlen[0]:0;
            tree[p].coverlen[2]= tree[p].cover>=2?tree[p].coverlen[0]:0;
        }
        else pushup(p);return;
    }
    int mid=l+r>>1;
    if(x<=mid) update(p<<1,l,mid,x,y,val);
    if(y>mid) update(p<<1|1,mid+1,r,x,y,val);
    pushup(p);
}
int main(){
//  freopen("a.in","r",stdin);
    scanf("%d",&tst);
    while(tst--){
        int num=0,m=0;double ans=0;
        scanf("%d",&n);
        while(n--){
            double x1,y1,x2,y2;
            scanf("%lf%lf%lf%lf",&x1,&y1,&x2,&y2);
            line[++num].x=x1;line[num].y1=y1;line[num].y2=y2;line[num].f=1;
            line[++num].x=x2;line[num].y1=y1;line[num].y2=y2;line[num].f=0;
            a[++m]=y1;a[++m]=y2;
        }
        sort(a+1,a+m+1);m=unique(a+1,a+m+1)-a-1;
        build(1,1,m-1);
        sort(line+1,line+num+1,cmp);
        line[0]=line[1];
        for(int i=1;i<=num;++i){
            ans+=tree[1].coverlen[2]*(line[i].x-line[i-1].x);
            update(1,1,m-1,lower_bound(a+1,a+m+1,line[i].y1)-a,lower_bound(a+1,a+m+1,line[i].y2)-a-1,line[i].f?1:-1);
        }
        printf("%.2lf\n",ans);
    }
    return 0;
}