HDU 1255 覆盖的面积 线段树+扫描线求面积并

基本和上一题类似...只是push_up函数中多了一部分东西，len1代表被至少覆盖一次的长度，len2代表至少被覆盖两次的长度。一段区间内被覆盖两次的长度有三种情况，一种是整个区间都被覆盖了两次，另外一个是整个区间被覆盖一次，然后左右孩子也有部分被覆盖一次的经历...还有就是左右孩子有部分被覆盖两次...然后和上一题基本不变..

#include<iostream>
#include<algorithm>
#include<cstdio>
#include<cstring>
#include<cmath>
#include<cstdlib>
#include<set>
#include<vector>
#include<iterator>
#include<map>
using namespace std;
#define LL long long

vector<double> v;
vector<double>::iterator it1;

struct node
{
    int l,r,cover;
    double len1,len2;
}data[20000];

struct lin
{
    double yu,yd,x;
    int cover;
}line[20000];

int cmp(lin a,lin b)
{
    return a.x<b.x;
}

void build(int l,int r,int k)
{
    data[k].l=l;
    data[k].r=r;
    data[k].len1=data[k].len2=0;
    data[k].cover=0;
    if(l+1==r)
        return ;
    int mid=(l+r)/2;
    build(l,mid,k*2);
    build(mid,r,k*2+1);
}

void push_up(int k)
{
    if(data[k].cover)
        data[k].len1=v[data[k].r]-v[data[k].l];
    else if(data[k].l+1==data[k].r)
        data[k].len1=0;
    else
        data[k].len1=data[k*2].len1+data[k*2+1].len1;

    if(data[k].cover>1)
        data[k].len2=v[data[k].r]-v[data[k].l];
    else if(data[k].l+1==data[k].r)
        data[k].len2=0;
    else if(data[k].cover==1)
        data[k].len2=data[k*2].len1+data[k*2+1].len1;
    else
        data[k].len2=data[k*2].len2+data[k*2+1].len2;
}

void updata(int l,int r,int k,int x)
{
    if(data[k].l>r || data[k].r<l)
        return;
    if(data[k].l>=l&&data[k].r<=r)
    {
        data[k].cover+=x;
        push_up(k);
        return ;
    }
    updata(l,r,k*2,x);
    updata(l,r,k*2+1,x);
    push_up(k);
}

int main()
{
    int n;
    int t;
    scanf("%d",&t);
    while(t--)
    {
        v.clear();
        v.push_back(-1);
        scanf("%d",&n);
        int m=0;
        double x1,x2,y1,y2;
        for(int i=0;i<n;i++)
        {
            scanf("%lf %lf %lf %lf",&x1,&y1,&x2,&y2);
            v.push_back(y1);
            line[m].cover=1;
            line[m].yd=y1;
            line[m].yu=y2;
            line[m].x=x1;
            m++;
            v.push_back(y2);
            line[m].cover=-1;
            line[m].yd=y1;
            line[m].yu=y2;
            line[m].x=x2;
            m++;
        }
        sort(v.begin(),v.end());
        it1=unique(v.begin(),v.end());
        v.erase(it1,v.end());
        int k=v.size();
        build(1,k-1,1);
        sort(line,line+m,cmp);
        double ans=0;
        for(int i=0;i<m-1;i++)
        {
            int x=lower_bound(v.begin(),v.end(),line[i].yd)-v.begin();
            int y=lower_bound(v.begin(),v.end(),line[i].yu)-v.begin();
            updata(x,y,1,line[i].cover);
            ans+=data[1].len2*(line[i+1].x-line[i].x);
        }
        printf("%0.2f\n",ans);
    }
    return 0;
}