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 300
struct 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;
}