（磨人的小妖精）hdu 1542 Atlantis：线段树扫描线，离散化

http://acm.hdu.edu.cn/showproblem.php?pid=1542

啊~~~~~做完这道题感觉智商被掏空！！！！！
1，首先由于横纵坐标不一定是整数，需要根据实际的数据进行离散化；
2，由于我的扫描是从下往上的，所以将横坐标进行离散化；
3，从下到上进行扫描的时候，当遇到一条底边则将覆盖数目加一，否则减一；
4，每次枚举每一条边，累加将前后两条边的高度差乘以区间覆盖长度；

#include <cstdio>
#include <iostream>
#include <algorithm>
#include <cstring>
using namespace std;
const int maxn=1009;
double pos[maxn*2];//保留离散化的结果
struct line
{
    double l,r,h;
    int f;//标记是上面（1）的边还是下面的边（-1）
} seg[maxn<<1];
struct node
{
    int l,r,cnt;
    double len;
} ac[maxn<<3];
int cmp(line a,line b)
{
    return a.h<b.h;//按照线段按高度排序
}
void build(int l,int r,int rt)//build维护的是一个闭区间
{
    ac[rt].l=l,ac[rt].r=r,ac[rt].len=ac[rt].cnt=0;
    if(l==r)
        return;
    int m=ac[rt].l+ac[rt].r>>1;
    build(l,m,rt<<1);
    build(m+1,r,rt<<1|1);
}
void pushup(int rt)//为什么要加一呢，因为build 维护的是闭区间！而更新的时候我们更新的是一个左闭右开的区间
{
    if(ac[rt].cnt) ac[rt].len=pos[ac[rt].r+1]-pos[ac[rt].l];//cnt不为零则全部被覆盖
    else if(ac[rt].l==ac[rt].r) ac[rt].len=0;//此时已不再是线段
    else ac[rt].len=ac[rt<<1].len+ac[rt<<1|1].len;//如果区间覆盖不连续，则从子区间找值
}
int fin(double key,int low,int high)//找出离散化后的位置
{
    while(low<=high)
    {
        int mid=high+low>>1;
        if(pos[mid]==key)
            return mid;
        else if(pos[mid]<key) low=mid+1;
        else high=mid-1;
    }
    return -1;
}
void update(int l,int r,int val,int rt)
{
    if(l<=ac[rt].l&&r>=ac[rt].r)
    {
        ac[rt].cnt+=val;
        pushup(rt);
        return;
    }
    int m=ac[rt].l+ac[rt].r>>1;
    if(m>=r)
        update(l,r,val,rt<<1);
    else if(m<l) update(l,r,val,rt<<1|1);
    else
    {
        update(l,r,val,rt<<1);
        update(l,r,val,rt<<1|1);
    }
    pushup(rt);
}
int main()
{
    int ase=1;
    int n;
    while(scanf("%d",&n)&&n)
    {
        double x1,x2,y1,y2;
        int m=0;
        for(int i=0; i<n; i++)
        {
            scanf("%lf%lf%lf%lf",&x1,&y1,&x2,&y2);
            seg[m].l=x1,seg[m].r=x2,seg[m].h=y1,seg[m].f=1;//遇到下面的边则加一
            seg[m+1].l=x1,seg[m+1].r=x2,seg[m+1].h=y2,seg[m+1].f=-1;//遇到上面的边则减一
            pos[m]=x1,pos[m+1]=x2;
            m+=2;
        }
        sort(seg,seg+m,cmp);
        sort(pos,pos+m);
        int cnt=1;
        for(int i=1; i<m; i++)//离散化后去重
            if(pos[i]!=pos[i-1])
                pos[cnt++]=pos[i];
        build(0,cnt-1,1);
        double ans=0;
        for(int i=0; i<m; i++)
        {
            int left=fin(seg[i].l,0,cnt-1);
            int right=fin(seg[i].r,0,cnt-1)-1;//要减一
            update(left,right,seg[i].f,1);
            ans+=(seg[i+1].h-seg[i].h)*ac[1].len;
        }
        printf("Test case #%d\n",ase++);
        printf("Total explored area: %.2f\n\n",ans);
    }
    return 0;
}