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 2010struct 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;}