HDU 1255(扫描线)

题目链接：http://acm.hdu.edu.cn/showproblem.php?pid=1255

 

题目大意：求至少两次覆盖的面积大小。

 

题目思路：之前还以为跟那个面积并一样。。直接上来就把cnt改成>=2，样例都没过。。后来看了题解才发现。因为update是到了符合要求的地方就停了，儿子们怎么样它并不知道，所以我们还需要维护一次的，如果这个地方是一次的，那就看它儿子中满足大于等于1次的长度，那那段长度就是已经>=2的长度。太机智了，整体代码都没变，就push_up复制粘贴了一遍改了点。

 

以下是代码：

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
#define inf 0x3f3f3f3f
#define rep(i,a,b) for(int i=a;i<=b;i++)
#define per(i,a,b) for(int i=a;i>=b;i--)
#define ll long long
const int MAXN = 1e5+5;
const int MOD = 1e9+7;
struct node{
    int l,r,cnt;
    double val,val2;
}a[MAXN<<2];
struct seg{
    double x,y1,y2;
    int flag;
}b[MAXN];
double c[MAXN];
bool cmp(seg a,seg b){
    return a.x<b.x;
}
void build(int rt,int l,int r){
    a[rt].l=l,a[rt].r=r;
    a[rt].cnt=a[rt].val=0;
    if(l==r)return;
    int mid=(l+r)>>1;
    build(rt<<1,l,mid);
    build(rt<<1|1,mid+1,r);
}
void push_up(int rt){
    if(a[rt].cnt){
        a[rt].val=c[a[rt].r+1]-c[a[rt].l];
    }
    else if(a[rt].l==a[rt].r)a[rt].val = 0;
    else a[rt].val=a[rt<<1].val+a[rt<<1|1].val;
}
void inpush_up(int rt){
    if(a[rt].cnt>=2){
        a[rt].val2=c[a[rt].r+1]-c[a[rt].l];
    }
    else if(a[rt].l==a[rt].r)a[rt].val2 = 0;
    else if(a[rt].cnt==1)a[rt].val2=a[rt<<1].val+a[rt<<1|1].val;
    else a[rt].val2=a[rt<<1].val2+a[rt<<1|1].val2;
}
void update(int rt,int l,int r,int val){
    if(l<=a[rt].l&&a[rt].r<=r){
        a[rt].cnt+=val;
        push_up(rt);
        inpush_up(rt);
        return;
    }
    int mid=(a[rt].l+a[rt].r)>>1;
    if(l<=mid)update(rt<<1,l,r,val);
    if(r>mid)update(rt<<1|1,l,r,val);
    push_up(rt);
    inpush_up(rt);
}
int main(){
    int n,t;
    scanf("%d",&t);
    while(t--){
        scanf("%d",&n);
        double x1,y1,x2,y2;
        int num=1,p=0;
        rep(i,1,n){
            scanf("%lf%lf%lf%lf",&x1,&y1,&x2,&y2);
            b[num].x=x1,b[num].y1=y1,b[num].y2=y2,b[num++].flag=1;
            b[num].x=x2,b[num].y1=y1,b[num].y2=y2,b[num++].flag=-1;
            c[p++]=x1,c[p++]=y1,c[p++]=x2,c[p++]=y2;
        }
        sort(b+1,b+num,cmp);
        sort(c,c+p);
        p=unique(c,c+p)-c;
        build(1,0,p);
        double ans=0;
        rep(i,1,2*n-1){
            int l=lower_bound(c,c+p,b[i].y1)-c,r=lower_bound(c,c+p,b[i].y2)-c-1;
            update(1,l,r,b[i].flag);
            ans+=(b[i+1].x-b[i].x)*a[1].val2;
        }
        printf("%.2f\n",ans);
    }
    return 0;
}