Atlantis hdu 1542 线段树扫面线区间合并

最新推荐文章于 2019-08-18 17:01:29 发布

原创最新推荐文章于 2019-08-18 17:01:29 发布 · 604 阅读

0 ·

CC 4.0 BY-SA版权

文章标签：

#ACM #algorithm

本文深入探讨扫描线算法在计算几何中的应用，特别是针对线段树的实现细节进行了详细讲解。文章通过具体实例介绍了如何利用扫描线算法计算平面区域内被线段覆盖的面积，并对关键的数据结构和算法进行了优化。

扫描线第一题,第一次交上去居然PE了,好吧,有傻逼了,少加个了个回车.

注意点:

1. 区间分段的时候要分成 [l,mid] 和[mid,r]不是平常一样的[l,mid]和[mid+1,r](因为我们要算的一段线段的长度,也就是说得连续,比如我们要算这一段[1,10],应该把他分成[1,5][5,10],不能[1,5][6,10]);

2.由于上面的mid和平时的线段树不同,所以判断退出的时候是(l+1==r)才退出,因为我们表示的是一段线段,每个值是一个点,如果l等于r的话,就是一个点了,也没意义了.

#include <iostream>
#include <cstdio>
#include <algorithm>
#include <cstring>
#include <vector>
#include <map>
#define lc idx<<1
#define rc idx<<1|1
#define mid ((l+r)>>1)
#define lson lc,l,mid
#define rson rc,mid,r
#define whole idx,l,r
using namespace std;
const int MAXN=100000+10;
int k;
class Line
{
    public:
    double x,yu,yd;
    int type;
    Line(){}
    Line(double a,double b,double c,int d) :x(a),yu(b),yd(c),type(d) {}
    bool operator <(const Line &a) const
    {
        return x<a.x;
    }
};
vector<Line> line;
vector<double> cor_y;
map<double,int> pj;
class SGtree
{
    public:
    double len[MAXN<<2];
    int cnt[MAXN<<2];
    void build(int idx,int l,int r)
    {
        len[idx]=cnt[idx]=0;
        if(l+1==r) return ;
        build(lson);
        build(rson);
    }
    void maintain(int idx,int l,int r)
    {
        if(cnt[idx]>0) len[idx]=cor_y[r]-cor_y[l];
        else
        {
            if(l+1==r) len[idx]=0;
            else len[idx]=len[lc]+len[rc];
        }
    }
    void update(int idx,int l,int r,int x,int y,int v)
    {
        if(x<=l && r<=y)
            cnt[idx]+=v;
        else
        {
            if(x<mid) update(lson,x,y,v);
            if(y>mid) update(rson,x,y,v);
        }
        maintain(whole);
    }
};
SGtree tree;
int main()
{
//    freopen("in","r",stdin);
    int n,kcase=1;
    while(~scanf("%d",&n) && n)
    {
        line.clear();
        cor_y.clear();
        pj.clear();
        double a,b,c,d;
        for(int i =0; i<n; i++)
        {
            scanf("%lf %lf %lf %lf",&a,&b,&c,&d);
            line.push_back(Line(a,d,b,1));
            line.push_back(Line(c,d,b,-1));
            cor_y.push_back(b);
            cor_y.push_back(d);
        }
        cor_y.push_back(-100000000);
        sort(line.begin(),line.end());
        sort(cor_y.begin(),cor_y.end());
        k=1;
        for(int i=1; i<cor_y.size(); i++)
        {
            if(cor_y[i]!=cor_y[i-1])
            {
                pj[cor_y[i]]=k;
                cor_y[k++]=cor_y[i];
            }
        }
        double ans=0;
        tree.build(1,1,k-1);
        tree.update(1,1,k-1,pj[line[0].yd],pj[line[0].yu],line[0].type);
        for(int i=1; i<line.size(); i++)
        {
            ans+=(line[i].x-line[i-1].x)*(tree.len[1]);
            tree.update(1,1,k-1,pj[line[i].yd],pj[line[i].yu],line[i].type);
        }
        printf("Test case #%d\n",kcase++);
        printf("Total explored area: %.2f\n\n",ans) ;
    }
    return 0;
}