HDU-1255

覆盖的面积

链接 HDU-1255
亚特兰蒂斯升级版，就是求的是覆盖两次的面积而不是并面积了，就在原来体面上用线段树多维护一个len2信息就好了

代码

---

#include<iostream>
#include<algorithm>
#include<cstring>
#include<vector>
#include<cstdio>
using namespace std;
const int N=1e4+10;
struct Segment{
    double y,x1,x2;
    int k;
    bool operator<(const Segment &s) const
    {
        return y<s.y;
    }
}seg[N*2];

//n个输入，每次输入四个坐标即两条我们需要存信息的边，也就是两个线段
//我们线段树存的区间，这样下来一共2n-1个区间，本来N*4是对于n来说的，现在2n-1所以就是n*4*2=8*n;
struct node{
    int l,r;
    int cnt;
    double len,len2;
}tr[N<<3];

int n;

//用来存放离散化后的存的y轴的信息（我们按照上面的图模型来说，这里我们建立的x,y与题目给出的x,y轴不一样，但是不影响后续的做题，和输入对应起来就行）;
vector<double> vs;
void build(int u,int l,int r)
{
    tr[u]={l,r,0,0,0};
    //printf("t %d %d %d\n",u,l,r);
    if(l==r) return ;
    else
    {
        int mid=l+r>>1;
        build(u<<1,l,mid);
        build(u<<1|1,mid+1,r);
    }
    
}

int find(double x)//和它相等的最左边的下标
{
    int l=0,r=vs.size();
    while(l<r)
    {
        int mid=l+r>>1;
        if(vs[mid]>=x)
        {
            r=mid;//缩小右边，不断找到最左边的它
        }
        else
        {
            l=mid+1;
        }
    }
    return l;
}

void pushup(int u)
{
    if(tr[u].cnt)
    {
    tr[u].len=vs[tr[u].r+1]-vs[tr[u].l];	
	}
	else if(tr[u].l!=tr[u].r)
	{
		tr[u].len=tr[u<<1].len+tr[u<<1|1].len;
	}
	else
	{
		tr[u].len=0;
	}
	if(tr[u].cnt>=2)
	{
		tr[u].len2=vs[tr[u].r+1]-vs[tr[u].l];
	}
	else if(tr[u].l!=tr[u].r&&tr[u].cnt==1)
	{
		tr[u].len2=tr[u<<1].len+tr[u<<1|1].len;
	}
	else if(tr[u].l!=tr[u].r)
	{
		tr[u].len2=tr[u<<1].len2+tr[u<<1|1].len2;
	}
	else tr[u].len2=0;

}

void modify(int u,int l,int r,int k)
{
    if(tr[u].l>=l&&tr[u].r<=r)
    {
        tr[u].cnt+=k;
        pushup(u);
    }
    else
    {
        int mid=tr[u].l+tr[u].r>>1;
        if(l<=mid) modify(u<<1,l,r,k);
        if(r>mid) modify(u<<1|1,l,r,k);
        pushup(u);
    }
    
}

int main()
{
    int t;
    cin>>t;
    while(t--)
    {
     cin>>n;
     int idx=0;
     vs.clear();
     for(int i=0;i<n;i++)
     {
     double x1,y1,x2,y2;
     scanf("%lf%lf%lf%lf",&x1,&y1,&x2,&y2);
     seg[idx++]={y1,x1,x2,1};
     seg[idx++]={y2,x1,x2,-1};
     vs.push_back(x1);
     vs.push_back(x2);
    } 
    sort(seg,seg+2*n);
    // for(int i=0;i<2*n;i++)
    // {
    //     cout<<seg[i].x<<" "<<seg[i].y1<<" "<<seg[i].y2<<" "<<seg[i].k<<endl;
    // }
    sort(vs.begin(),vs.end());
    //unique是给定排好序的开头结尾，返回去重后的结尾
    //erase是给定开头结尾然后删去，就是从重复的地方开始到最后
    vs.erase(unique(vs.begin(),vs.end()),vs.end());
    // for(int i=0;i<vs.size();i++)
    // {
    //     cout<<vs[i]<<" ";
    // }
    //cout<<endl;
    build(1,0,vs.size()-2);//一共2n-1个区间，维护的区间从vs[0]开始,维护(vs[0]--vs[1]),(vs[1]--vs[2])的一个个区间
    double res=0;//存一下答案
    for(int i=0;i<2*n;i++)
    {
    if(i)//第一个线段无法求面积，从第二个开始
    {
        //cout<<"len: "<<tr[1].len<<endl;
        res+=tr[1].len2*(seg[i].y-seg[i-1].y);
    }
    
    //因为线段树维护的是vs[tr[i].l],vs[tr[i].r+1],所以我们传入l,r-1就可以了
    modify(1,find(seg[i].x1),find(seg[i].x2)-1,seg[i].k);
    }
    printf("%.2f\n",res);
    }
    return 0;
}

---