码迷,mamicode.com
首页 > 其他好文 > 详细

hdu1542线段树(扫描线+离散化)

时间:2018-02-19 19:19:57      阅读:181      评论:0      收藏:0      [点我收藏+]

标签:xpl   根据   扫描   log   min   复杂   tar   数组   ons   

题目链接

要求矩形的面积并

代码不复杂,主要要理解扫描线的思想以及一些细节的处理。

首先需要将接收到的x坐标离散化,方法就是排序去重。接下来的线段树建立在这个

关于x坐标的数组上,这很关键。线段树的节点代表一段区间,这个区间是由x坐标数组的下标

来构成的。更新的时候就根据水平线段的左右x坐标获得区间,然后更新区间。

看了这么多,离散化之后关于r+1,r-1的问题还是很模糊。

我尽力理解一下,有一个[0,4]的区间,叶节点分别为[0,0],[1,1],[2,2],[3,3],[4,4],

那么[0,0]节点表示的是区间[0,1];[1,1]节点表示的区间[1,2];[2,2]节点表示的区间[2,3];[3,3]节点表示区间[3,4];

当向上传递时[0,1]节点表示区间[0,2],以此类推

#include<cstdio>
#include<algorithm>
#include<vector>
using namespace std;
const int maxn=110;//最多矩形个数
struct seg{
    double lx,rx,y;
    short f;
    seg(){}
    seg(double lx_,double rx_,double y_,short f_):lx(lx_),rx(rx_),y(y_),f(f_){}
    bool operator <(const seg &a)const 
    {
            return y<a.y;
    }
};
int sgsNum;
seg sgs[maxn*2];
int x2nNum;
//x2n
double x2n[maxn*2];
double len[maxn*8];
int num[maxn*8];
void build(int root,int l,int r)
{
    len[root]=num[root]=0;
    if(l==r)return;
    int mid=(l+r)/2;
    build(root*2,l,mid);
    build(root*2+1,mid+1,r);
}
void pushUp(int root,int l,int r)
{
    if(num[root]){
        len[root]=x2n[r+1]-x2n[l];
    }else if(l==r){
        len[root]=0;
    }
    else {
        len[root]=len[root*2]+len[root*2+1];
    }
}
void update(int root,int L,int R,int f,int l,int r)
{
    if(L<=l&&r<=R){
        num[root]+=f;
        pushUp(root,l,r);
        return;
    }
    int mid=l+(r-l)/2;
    if(L<=mid)update(root*2,L,R,f,l,mid);
    if(mid<R)update(root*2+1,L,R,f,mid+1,r);
    pushUp(root,l,r);
}
int bin(double k)
{
    int m,l=0,r=x2nNum-1;
    while(l<=r){
        m=(l+r)/2;
        if(x2n[m]==k)return m;
        else if(x2n[m]>k)r=m-1;
        else l=m+1;
    }
    return -1;
}
int main()
{
    //freopen("in.txt","r",stdin);
    int n,p=0;
    while(scanf("%d",&n)!=EOF&&n!=0){
        sgsNum=0,x2nNum=0;
        p++;
        for(int i=0;i<n;i++){
            double lx,ly,rx,ry;
            scanf("%lf%lf%lf%lf",&lx,&ly,&rx,&ry);
            sgs[sgsNum++]=seg(lx,rx,ly,1);
            sgs[sgsNum++]=seg(lx,rx,ry,-1);
            x2n[x2nNum++]=lx;
            x2n[x2nNum++]=rx;
        }
        double area=0;
        sort(x2n,x2n+x2nNum);
        sort(sgs,sgs+sgsNum);
        x2nNum=unique(x2n,x2n+x2nNum)-x2n;
        build(1,0,x2nNum-1);
        for(int i=0;i<sgsNum-1;i++){
            int l=bin(sgs[i].lx);
            int r=bin(sgs[i].rx)-1;
            update(1,l,r,sgs[i].f,0,x2nNum-1);
            area+=len[1]*(sgs[i+1].y-sgs[i].y);
        }
        printf("Test case #%d\n",p);
        printf("Total explored area: %.2lf\n\n",area);
    }
    //while(1);
}

参考资料

http://www.cnblogs.com/zhangmingcheng/p/3907072.html

http://blog.csdn.net/qq_37497322/article/details/75126251

hdu1542线段树(扫描线+离散化)

标签:xpl   根据   扫描   log   min   复杂   tar   数组   ons   

原文地址:https://www.cnblogs.com/MalcolmMeng/p/8454230.html

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!