HDU - 1542 - Atlantis （线段树扫描线求矩形并）

题目连接： HDU - 1542 - Atlantis

大佬博客讲解

通过观察可以发现，以y轴为扫描线从下往上扫，若一个区间下边数量大于上边，则扫描的过程中这部分肯定算面积的，因此用线段树维护这样的区间就行了。

#include<bits/stdc++.h>
using namespace std;
typedef double ld;
const int N=207;
int n,m;
struct Line
{
    int x1,x2;
    int down;
    ld y;
    bool operator < (const Line & r) const
    {
        return y<r.y;
    }
};
struct Shape
{
    ld x1,x2,y1,y2;
    Shape(){}
    Shape(ld x1,ld y1,ld x2,ld y2) : x1(x1),y1(y1),x2(x2),y2(y2) { }
};
Shape s[N];
vector<Line> line;
ld xx[N],sum[N<<2];
int c[N<<2];
ld dis(int l,int r) { return xx[r]-xx[l-1]; }
void build(int rt,int l,int r)
{
    c[rt]=0;sum[rt]=0;
    if(l==r) return ;
    int mid=(l+r)>>1;
    build(rt<<1,l,mid);
    build(rt<<1|1,mid+1,r);
}
inline void push_down(int rt,int l,int r)
{
    if(c[rt]!=-1)
    {
        c[rt<<1]=c[rt<<1|1]=c[rt];
        int m=(l+r)>>1;
        sum[rt<<1]=(c[rt]?dis(l,m):0);
        sum[rt<<1|1]=(c[rt]?dis(m+1,r):0);
    }
}
inline void push_up(int rt)
{
    c[rt]=(c[rt<<1]==c[rt<<1|1]&&c[rt<<1]!=-1?c[rt<<1]:-1);
    sum[rt]=sum[rt<<1]+sum[rt<<1|1];
}
void update(int rt,int l,int r,int ql,int qr,int v)
{
    if(ql<=l&&qr>=r)
    {
        if(c[rt]!=-1)
        {
            c[rt]+=v;
            sum[rt]=(c[rt]?dis(l,r):0);
            return;
        }
    }
    push_down(rt,l,r);
    int mid=(l+r)>>1;
    if(ql<=mid) update(rt<<1,l,mid,ql,qr,v);
    if(qr>mid) update(rt<<1|1,mid+1,r,ql,qr,v);
    push_up(rt);
}
int main()
{
    int kase=1;
    while(~scanf("%d",&n))
    {
        if(n==0) break;
        for(int i=0;i<n;++i)
        {
            ld x1,x2,y1,y2;
            scanf("%lf%lf%lf%lf",&x1,&y1,&x2,&y2);
            s[i]=Shape(x1,y1,x2,y2);
            xx[i<<1]=x1;xx[i<<1|1]=x2;
        }
        sort(xx,xx+n*2);
        m=unique(xx,xx+n*2)-xx;
        line.clear();
        for(int i=0;i<n;++i)
        {
            int x1id=lower_bound(xx,xx+m,s[i].x1)-xx+1;
            int x2id=lower_bound(xx,xx+m,s[i].x2)-xx;
            line.push_back({x1id,x2id,1,s[i].y1});
            line.push_back({x1id,x2id,-1,s[i].y2});
        }
        // 线段树上为[l,r]的点 表示[X[l-1],X[r]] 这段区间。
        sort(line.begin(),line.end());
        ld ans=0;
        build(1,1,--m);
        update(1,1,m,line[0].x1,line[0].x2,line[0].down);
        for(int i=1;i<line.size();++i)
        {
            ans+=sum[1]*(line[i].y-line[i-1].y);
            update(1,1,m,line[i].x1,line[i].x2,line[i].down);
        }
        printf("Test case #%d\nTotal explored area: %.2f\n\n",kase++,ans);
    }
    return 0;
}