poj 1151 hdu 1542 atlantis （线段树+扫描线）

盗用一张样例的图，不过x,y轴反了，大家最好照着图先在纸上演一遍。再看代码，代码我注释很详细，希望大家都能看懂。

第二次写这题了，又发现上一次没完全弄懂的东西，所以重新来一遍，写清楚了所有注释。

#include <iostream>
#include <cstring>
#include <cstdio>
#include <algorithm>
#include <cmath>
using namespace std;
#define N 100010
#define lson v<<1
#define rson v<<1|1
struct node{int l,r;double lv,rv,len;int flag; //l,r线段树表示的左右段，lv,rv，所表示左右段对应的y值； 
}t[N<<2];
double ans[N];
struct key{double x,y1,y2;int s;         // 以x排序，y1<y2；存每个x节点对应的段；s就是flag标记它是入边还是出边 
       friend bool operator<(const key &b,const key &c)
       {
             return b.x<c.x;                     // 以x排序； 
       }
}a[N];
const double eps=1e-5;
int f(double a,double b)
{
if(fabs(a-b)<eps)return 0;               //比较函数，返回0说明a==b 
else if(a-b>eps)return 1;                 // a>b
else return -1;                           //a<b ；
                                       // 为以防出现double 比较大小错误故写比较函数，当然这道题精度不大，不写此也可以AC; 
}
void pushup(int v)
{
if(t[v].flag>0)
t[v].len=t[v].rv-t[v].lv;                  //若整个节点的flag都>0，则长度等于其右端-左端； 
else if(t[v].l+1==t[v].r)t[v].len=0;       //因为建树不同，满足左段加1等于右段，则为0；建树如下文注释； 
else
t[v].len=t[lson].len+t[rson].len;      // 不满足以上两个条件，则递归即可； 
}
void build(int v,int l,int r)
{
t[v].l=l,t[v].r=r;t[v].lv=ans[l];t[v].rv=ans[r];t[v].flag=0;t[v].len=0;
if(l==r-1)
{
return ;
}
build(lson,l,(l+r)/2);        //为保证其为满二叉树，及左右两段不能有间隙，故建树有些特别。 
build(rson,(l+r)/2,r);       //左为l to mid, mid to r;   
}
void update(int v,double y1,double y2,int m)
{
if(f(t[v].lv,y1)==0&&f(t[v].rv,y2)==0)
{
  t[v].flag+=m;                        //直接flag加上所扫描段的flag再进行判断，然后pushup执行以上操作 
  pushup(v);      //（上一行）这里得仔细考虑，我的理解是从左扫描到右，必是先入边再出边，故负标记不会对非其对应入边产生影响。 
  return ;
}
if(f(y2,t[rson].lv)!=1)update(lson,y1,y2,m); //如果所扫描段的右端不大于（！=1表示不大于）此节点右儿子的左端，故查询更新左儿子节点 
else if(f(y1,t[lson].rv)!=-1)update(rson,y1,y2,m); //同上； 
else
{
update(lson,y1,t[lson].rv,m);     //如果所扫描线段分部左右儿子两边，则断开它，分别搜索。画画图就可以理解。 
update(rson,t[rson].lv,y2,m);
}
pushup(v);                           //同理； 
}
int main()
{
int n;int q=1;
while(scanf("%d",&n),n!=0)
{
    int p=1;
for(int i=0;i<n;i++)
{
double x1,y1,x2,y2;
   scanf("%lf%lf%lf%lf",&x1,&y1,&x2,&y2); 
   ans[p]=y1,ans[p+1]=y2;                //ans数组是给每个树的节点赋值的，仅保存y值； 
   a[p].x=x1;a[p].y1=y1;a[p].y2=y2;a[p++].s=1; //第一条，入度为1； 
   a[p].x=x2;a[p].y1=y1;a[p].y2=y2;a[p++].s=-1;//第二条，入度为-1； 
    }
    sort(ans+1,ans+p);                      //排序 
    sort(a+1,a+p);
    build(1,1,p-1);                         //p最后多加了1，减掉再建树； 
    update(1,a[1].y1,a[1].y2,a[1].s);       //更新第一条线段；其实最左端的线段，入度必为1； 
    double tot=0;
for(int i=2;i<p;i++)
    {
     tot+=t[1].len*(a[i].x-a[i-1].x);   //开始扫描过去，根节点长度即可，根节点长度及为长，再乘上宽，就是每段扫描需加上的面积； 
     update(1,a[i].y1,a[i].y2,a[i].s);  //更新 
    }
    printf("Test case #%d\n",q++);
    printf("Total explored area: %.2f\n\n",tot);//输出不能用 lf，lf在poj报错。 
}
return 0;
}