HDU-1542：Atlantis（扫描线模板题）

题目链接：点击打开链接

题目大意：

题目给你一些矩形的左下角坐标和右上角坐标，然后计算形成的新图形的面积。中间重叠的部分算作一部分。

题意解析：

要用到线段树扫描线的知识，从坐标系的下面往上面一次扫，每次维护当前高度的有效线段长。然后就可以通过高度差还有效线段长计算有效面积。

代码：

#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
#include<queue>
#include<vector>
using namespace std;

struct node     //记录线段左端点右端点和当前线段高度
{
    double l,r,h;
    int flag;       //记录当前线段为上位边和下位边
    bool operator<(const node &b)      //通过线段高度对线段排序
    {
        return h<b.h;
    }
}nodes[500000];
struct tree
{
    int l,r,mid;
    int cnt;    //下位边减上位边的数
    double sum;  //本节点控制区域内cnt不为0的总长度
}t[500000];
int n;
double x[500000];
void pushdown(int rt)       //从上到下维护
{
    if(t[rt].cnt!=-1)
    {
        t[rt<<1].cnt=t[rt<<1|1].cnt=t[rt].cnt;
        if(t[rt].cnt)
        {
            t[rt<<1].sum=x[t[rt].mid+1]-x[t[rt].l];
            t[rt<<1|1].sum=x[t[rt].r+1]-x[t[rt].mid+1];
        }
        else
        {
            t[rt<<1].sum=0;
            t[rt<<1|1].sum=0;
        }
    }
}
void pushup(int rt)     //自上而下维护，cnt=-1代表当前区间左右儿子的上下位边数不同
{
    if(t[rt<<1].cnt==-1||t[rt<<1|1].cnt==-1)
        t[rt].cnt=-1;
    else if(t[rt<<1].cnt!=t[rt<<1|1].cnt)
        t[rt].cnt=-1;
    else
        t[rt].cnt=t[rt<<1].cnt;
    t[rt].sum=t[rt<<1].sum+t[rt<<1|1].sum;
    //printf("----%lf %d %lf %d\n",t[rt<<1].sum,t[rt<<1].l,t[rt<<1|1].sum,t[rt<<1|1].l);
}
void build(int l,int r,int rt)
{
    if(l==r)
    {
        t[rt].l=t[rt].r=l;
        t[rt].mid=(t[rt].l+t[rt].r)/2;
        t[rt].cnt=0;
        t[rt].sum=0.0;
        return ;
    }
    int mid=(l+r)/2;
    t[rt].l=l;
    t[rt].r=r;
    t[rt].mid=mid;
    t[rt].cnt=0;
    t[rt].sum=0.0;
    build(l,mid,rt<<1);
    build(mid+1,r,rt<<1|1);
    pushup(rt);
}
void update(int l,int r,int flag,int rt)    //更新结果
{
    if(t[rt].l>=l&&t[rt].r<=r)
    {
        if(t[rt].cnt!=-1)
        {
            t[rt].cnt=t[rt].cnt+flag;
            if(t[rt].cnt)
                t[rt].sum=x[t[rt].r+1]-x[t[rt].l];
            else
                t[rt].sum=0;
            return ;
        }
    }
    pushdown(rt);
    if(l<=t[rt].mid)
        update(l,r,flag,rt<<1);
    if(r>t[rt].mid)
        update(l,r,flag,rt<<1|1);
    pushup(rt);
}
int bin(double key,int n,double d[])
{
    int l=1,r=n;
    while(l<=r)
    {
        int m=(l+r)/2;
        if(d[m]==key)
            return m;
        else if(d[m]>key)
            r=m-1;
        else
            l=m+1;
    }
    return -1;
}
int main()
{
    int kase=0;
    while(scanf("%d",&n)!=EOF)
    {
        int gg=0;
        if(n==0)
            break;
        for(int i=0;i<n;i++)
        {
            double x1,y1,x2,y2;
            scanf("%lf%lf%lf%lf",&x1,&y1,&x2,&y2);
            x[++gg]=x1;
            nodes[gg].l=x1;
            nodes[gg].r=x2;
            nodes[gg].h=y1;
            nodes[gg].flag=1;
            x[++gg]=x2;
            nodes[gg].l=x1;
            nodes[gg].r=x2;
            nodes[gg].h=y2;
            nodes[gg].flag=-1;
        }
        sort(x+1,x+gg+1);
        sort(nodes+1,nodes+gg+1);
        int k=1;
        for(int i=2;i<=gg;i++)
        {
            if(x[i-1]!=x[i])
                x[++k]=x[i];
        }
        build(1,k-1,1);     //用每个节点代表一段区间
        double ans=0.0;
        for(int i=1;i<gg;i++)
        {
            int l=bin(nodes[i].l,k,x);
            int r=bin(nodes[i].r,k,x)-1;
            if(l<=r)
                update(l,r,nodes[i].flag,1);
            //printf("%lf\n",t[1].sum);
            ans=ans+t[1].sum*(nodes[i+1].h-nodes[i].h);
        }
        printf("Test case #%d\nTotal explored area: %.2lf\n\n",++kase,ans);
    }
}