poj 1151 hdu 1542 Atlantis 线段树扫描线，详细讲解，(*^__^*) 嘻嘻……

我的扫面线第一题，一开始看网上讲的都好抽象，最后还是研究别人代码整明白的，所以我要写一个直观的，哈哈哈！！希望大家都能看懂

如图虚线将整个图型分成三个矩形，我们现将每个点的x进行排序，也就是（10，15，20，25.5）

这样就知道矩形的长了，还差宽，然后我们就要开始扫描线啦

先每个点的y值离散话，因为我们要求的是映射在x轴的线段，然后建树，共有t哥点，样立t-1等于4

下面的线段树我每个离散花后的点对应的y值我都标记上了

将一号边加入线段，先加入一号边，二号边三号边四号边具体过程建议自己模拟一下

我每个线段树的点用s标记这条边被覆盖还是没被都盖，每条边的flag表示这天边是入度还是出度

注意啊！！我这个图x，y轴画反啦

#include<iostream>#include<cstdio>#include<algorithm>using namespace std;#define N 300struct Node{    double x;double y1;double y2;    int flag;}node[N];bool cmp(Node a,Node b){    return a.x-b.x<0.0000001;}double y[N];struct node{    int l;int r;double ml;double mr;int s;double len;}a[N*3];void build(int i,int left,int right){    a[i].l=left;    a[i].r=right;    a[i].ml=y[left];    a[i].mr=y[right];    a[i].s=0;    a[i].len=0;    if(a[i].l+1==a[i].r){        return ;    }    int mid=(left+right)>>1;    build(i*2,left,mid);    build(i*2+1,mid,right);//建树时注意这里不是mid+1，因为做相减的时候如果mid+1这么建回到值左孩子的右边与有孩子的左边无法进行运算}void callen(int i){    if(a[i].s>0){//注意这里不是所有边都是左孩子的长度加上右孩子的长度，他存在一个覆盖问题        a[i].len=a[i].mr-a[i].ml;    }else if(a[i].r-a[i].l==1){        a[i].len=0;    }else{        a[i].len=a[i*2].len+a[i*2+1].len;    }    return ;}void updata(int i,Node b){    if(a[i].ml==b.y1&&a[i].mr==b.y2){        a[i].s+=b.flag;        callen(i);        return ;    }    if(b.y2<=a[i*2].mr) updata(i*2,b);    else if(b.y1>=a[i*2+1].ml) updata(i*2+1,b);    else{        Node temp=b;        temp.y2=a[i*2].mr;        updata(i*2,temp);        temp=b;        temp.y1=a[i*2+1].ml;        updata(i*2+1,temp);    }    callen(i);    return ;}int main(){//    freopen("in.txt","r",stdin);    int n,t,p=1,te;    double x1,x2,y1,y2;    while(scanf("%d",&n),n){            t=1;            for(int i=0;i<n;i++){            scanf("%lf%lf%lf%lf",&x1,&y1,&x2,&y2);            node[t].x=x1;            node[t].y1=y1;            node[t].y2=y2;            node[t].flag=1;//入边            y[t++]=y1;            node[t].x=x2;            node[t].y1=y1;            node[t].y2=y2;            node[t].flag=-1;//出边            y[t++]=y2;        }        sort(node+1,node+t,cmp);        sort(y+1,y+t);        build(1,1,t-1);        updata(1,node[1]);        double sum=0;        for(int i=2;i<t;i++){            sum+=a[1].len*(node[i].x-node[i-1].x);            updata(1,node[i]);        }        printf("Test case #%d\n",p++);        printf("Total explored area: %.2lf\n\n",sum);    }    return 0;}