HDU

题目连接： 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;}