覆盖的面积（线段树-扫描线）

覆盖的面积

题目要求的时统计覆盖过两次以上部分的面积，query判断修改下就行了，除此之外和Atlantis完全一样

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<map>
using namespace std;
#define MAXN 1005
#define ll p<<1
#define rr p<<1|1
struct node
{
    int l,r;
    int cover;
}t[MAXN*8];
struct edge
{
    double x,yd,yu;
    int lr;
    bool operator < (const edge a) const {return x<a.x;}
}e[MAXN*2];
double yy[MAXN*2];
void construct(int l,int r,int p)
{
    t[p].l=l,t[p].r=r;
    t[p].cover=0;
    if(r-l<=1) return ;
    int m=(l+r)>>1;
    construct(l,m,ll);
    construct(m,r,rr);
}
void modify(int l,int r,int val,int p)
{
    if(t[p].l==l&&t[p].r==r)
    {
        t[p].cover+=val;
        return ;
    }
    int m=(t[p].l+t[p].r)>>1;
    if(r<=m) modify(l,r,val,ll);
    else if(l>=m) modify(l,r,val,rr);
    else modify(l,m,val,ll),modify(m,r,val,rr);
}
void query(double &ans,int p)
{
    if(t[p].cover>1)
    {
        ans+=yy[t[p].r]-yy[t[p].l];
        return ;
    }
    else if(t[p].r-t[p].l>1)
    {
        t[ll].cover+=t[p].cover,t[rr].cover+=t[p].cover;
        query(ans,ll);
        query(ans,rr);
        t[ll].cover-=t[p].cover,t[rr].cover-=t[p].cover;
    }
}
int main()
{
    map<double,int> hash;
    int cas,n,i;
    double x1,x2,yd,yu;
    scanf("%d",&cas);
    while(cas--)
    {
        hash.clear();
        scanf("%d",&n);
        for(i=0;i<n*2;i+=2)
        {
            scanf("%lf%lf%lf%lf",&x1,&yd,&x2,&yu);
            e[i].x=x1,e[i].yd=yd,e[i].yu=yu,e[i].lr=1;
            e[i+1].x=x2,e[i+1].yd=yd,e[i+1].yu=yu,e[i+1].lr=-1;
            yy[i]=yd,yy[i+1]=yu;
        }
        sort(e,e+i);
        sort(yy,yy+i);
        int up=int (unique(yy,yy+i)-yy);
        for(i=0;i<up;i++)
            hash[yy[i]]=i;
        construct(0,up-1,1);
        double ans=0,dy;
        for(i=0;i<n*2-1;i++)
        {
            modify(hash[e[i].yd],hash[e[i].yu],e[i].lr,1);
            dy=0;
            query(dy,1);
            ans+=dy*(e[i+1].x-e[i].x);
        }
        printf("%.2lf\n",ans);
    }
    return 0;
}