计算几何进阶-扫描线

最新推荐文章于 2024-05-30 22:39:06 发布
最新推荐文章于 2024-05-30 22:39:06 发布
阅读量675 收藏
点赞数
分类专栏： 计算几何 文章标签： 计算几何 扫描线
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/qq_34921856/article/details/80803148
版权

关于扫描线的题目，做法我感觉应该比较容易看出来，但是其中的性质很难找。只要结合题目的背景来找出它其中的数学性质，题目就能迎刃而解了。
A - Coneology
题意:一个笛卡尔坐标系中有许多的圆，有些圆包含其他的圆，一定不存在相交的圆。问有多少个没有被任一个圆包含的圆。
4e4的数据量肯定不允许我们暴力，暴力实际上有许多没必要的计算。因为根据他们的位置关系，有些圆之间不相交能推出另个圆也不和这个圆相交。所以我们得利用他们的位置关系来计算，这里也就是用到了扫描线。
从左往右扫描，左端点进set，右端点出set。并且set是到目前为止右端点还没到的不相交的圆的集合。
碰到左端点的时候，就得将该圆与set里面的圆进行包含判断。而这个判断也并不是与所有set里的圆都得判断一遍，我们只需找到set里面y第一个大于等于它的和y第一个小于它的两个圆，与其判断即可。大家可以画画图想一想。
碰到右端点将set里面的圆删除即可。
所以我们要在set里面存y值及下标。

#include<iostream>
#include<algorithm>
#include<stdio.h>
#include<cstring>
#include<vector>
#include<set>
#include<math.h>
using namespace std;
const int maxn=400005;
const double eps=1e-3;
int sgn(double x)
{
    if(fabs(x)<eps)return 0;
    if(x>0)return 1;
    return -1;
}
double r[maxn],x[maxn],y[maxn];
vector<pair<double,int> >V;
set<pair<double,int> >S;
vector<int>ans; 

bool check(int idx1,int idx2)
{
    double X=x[idx1]-x[idx2];
    double Y=y[idx1]-y[idx2];
    double dis=sqrt(X*X+Y*Y);
    return sgn(dis-r[idx1]-r[idx2])>=0;
}

int main()
{
    int N;
    scanf("%d",&N);
    for(int i=0; i<N; i++)
    {
        scanf("%lf %lf %lf",&r[i],&x[i],&y[i]);
        V.push_back(make_pair(x[i]-r[i],i));
        V.push_back(make_pair(x[i]+r[i],i+N));
    }
    sort(V.begin(),V.end());
    for(int i=0; i<V.size(); i++)
    {
        int idx=V[i].second%N;
        if(V[i].second>=N)
        {
            S.erase(make_pair(y[idx],idx));
        }
        else
        {
            set<pair<double,int> >::iterator tidx=S.lower_bound(make_pair(y[idx],idx));
            if(tidx!=S.end())
                if(!check(tidx->second,idx))
                    continue;
            if(tidx!=S.begin())
                if(!check((--tidx)->second,idx))
                    continue;
            S.insert(make_pair(y[idx],idx));
            ans.push_back(idx);
        }
    }
    sort(ans.begin(),ans.end());
    printf("%d\n",ans.size());
    for(int i=0; i<ans.size(); i++)
    {
        if(i)printf(" ");
        printf("%d",ans[i]+1);
    }
    puts("");
    return 0;
}

B - Moonmist
求最小圆距。
直接讲思路吧，画了半天图总感觉怪怪的。。
求这个圆距，我们可以二分距离num，然后判断是否有圆相交。
判断是否有圆相交可以用两条扫描线，一个枚举左边界i，一个枚举右边界j。
当x[i]-r[i]-num<=x[j]+r[j]+num,则判断i圆并将i圆放入set里，否则将i圆在set里删掉。
这样写可以过题，但是有bug。
大家可以看看这个样例
1
3
-10 20 1
0 0 10
30 40 39
所以我们需要在圆删除之前再判断一次。

#include<bits/stdc++.h>
using namespace std;
const int maxn=50005;
int X[maxn],Y[maxn],R[maxn];
int rL[maxn],rR[maxn],rU[maxn];
bool cmpL(int a,int b)
{
    return X[a]-R[a]<X[b]-R[b];
}
bool cmpR(int a,int b)
{
    return X[a]+R[a]<X[b]+R[b];
}
bool cmpU(int a,int b)
{
    if(Y[a]!=Y[b])return Y[a]<Y[b];
    return X[a]<X[b];
}
int urank[maxn];

double dis(int a,int b)
{
    double x=X[a]-X[b];
    double y=Y[a]-Y[b];
    return sqrt(x*x+y*y);
}
double dis2(int a,int b)
{
    double x=X[a]-X[b];
    double y=Y[a]-Y[b];
    return x*x+y*y;
}
set<int>S;
typedef set<int>::iterator it;
bool check2(int num,double add)
{
    int itrank=urank[num];
    it itor=S.lower_bound(itrank);
    if(itor!=S.end())
    {
        double dist=dis(rU[*itor],num);
        if(dist<=R[rU[*itor]]+R[num]+2*add)return 1;
    }
    if(itor!=S.begin())
    {
        itor--;
        double dist=dis(rU[*itor],num);
        if(dist<=R[rU[*itor]]+R[num]+2*add)return 1;
    }

    return 0;

}
int n;
bool check(double num)
{
    S.clear();
    int i=0,j=0;
    while(i<n||j<n)
    {
        if(j==n||i<n&&X[rL[i]]-R[rL[i]]-num<=X[rR[j]]+R[rR[j]]+num)
        {
            if(check2(rL[i],num))
                return true;
            S.insert(urank[rL[i]]);
            i++;
        }
        else
        {
            S.erase(urank[rR[j]]);
            if(check2(rR[j],num))
                return true;
            j++;
        }
    }
    return false;
}

void solve()
{
    double l=0,r=dis(0,1);
    double mid;
    for(int t=1;t<=40;t++)
    {
        mid=(l+r)/2;
        if(check(mid))
            r=mid;
        else l=mid;
    }
    printf("%.6f\n",mid*2);
}

int main()
{
    int T;
    scanf("%d",&T);
    while(T--)
    {
        scanf("%d",&n);
        for(int i=0;i<n;i++)
        {
            scanf("%d %d %d",&X[i],&Y[i],&R[i]);
            rL[i]=rR[i]=rU[i]=i;
        }
        sort(rL,rL+n,cmpL);
        sort(rR,rR+n,cmpR);
        sort(rU,rU+n,cmpU);
        for(int i=0;i<n;i++)
            urank[rU[i]]=i;
        solve();
    }
}

之后会更详细的讲解该题。。
C - Light and Shadow
这题一开始理解错题意了，该题意为能照到的棍子的个数。
由于题目已经说了，不存在两根棍子相交，所以距离点光源最近的棍子不管点光源怎么射，它都是最前面的。所以我们可以以点光源为原点，以极坐标排序，逆时针扫一下。记录一下每条棍子的起点和终点，当遇到棍子的起点时，我们可以以原点和它的起点作一条射线，计算射线与棍子的交点距原点的距离来重新将set里面的棍子排序。
当然这里有一个特别的地方，就是与射线(-1,0)相交的棍子它们的起点在第二象限，按照我们那种排序方式它会在最后才被加进去，所以我们得预先把它装进来。

#include<bits/stdc++.h>
using namespace std;
const int maxn=10005;
const double eps=1e-8;
int sgn(double x)
{
    if(fabs(x)<eps)return 0;
    if(x>0)return 1;
    return -1;
}
struct Point
{
    double x,y;
    Point(double sx,double sy):x(sx),y(sy){}
    Point(){}
    Point operator-(const Point &b)const
    {
        return Point(x-b.x,y-b.y);
    }
    double operator*(const Point &b)const
    {
        return x*b.x+y*b.y;
    }
    double operator^(const Point &b)const
    {
        return x*b.y-y*b.x;
    }
    bool operator<(const Point &b)const
    {
        return x*b.y<y*b.x;
    }
    void input()
    {
        scanf("%lf %lf",&x,&y);
    }
};
Point cur;
typedef const Point CP;
Point intersection(CP &u1,CP &u2,CP &v1,CP &v2)
{
    Point ret=u1;
    double t=((u1.x-v1.x)*(v1.y-v2.y)-(u1.y-v1.y)*(v1.x-v2.x))
    /((u1.x-u2.x)*(v1.y-v2.y)-(u1.y-u2.y)*(v1.x-v2.x));
    ret.x+=(u2.x-u1.x)*t;
    ret.y+=(u2.y-u1.y)*t;
    return ret;
}
double dis(Point a)
{
    return sqrt(a.x*a.x+a.y*a.y);
}
struct line
{
    Point u,v;
    int idx;
    int flag;
    double ang;
    line(Point _u,Point _v):u(_u),v(_v)
    {
        ang=atan2(u.y,u.x);
    }
    line(){}
    bool operator<(const line &b)const
    {
        Point tmpa=intersection(Point(0,0),cur,u,v);
        Point tmpb=intersection(Point(0,0),cur,b.u,b.v);
        return sgn(dis(tmpa)-dis(tmpb))<0;
    }
}Lines[maxn*2];

bool cmp(line a,line b)
{
    return a.ang<b.ang;
}
set<line>S;
bool ok[maxn];
bool s_intersection(Point a,Point b)
{
    Point tmp=Point(-1,0);
    return sgn((a^tmp)*(a^b))>0&&sgn((a^tmp)*(tmp^b))>0;
}
int n;
void input()
{
    Point O;
    S.clear();
    memset(ok,0,sizeof(ok));
    O.input();
    Point a,b;
    for(int i=0;i<n;i++)
    {
        a.input();
        b.input();
        a=a-O;
        b=b-O;
        if(a<b)swap(a,b);
        Lines[2*i]=line(a,b);
        Lines[2*i].idx=i,Lines[2*i].flag=1;
        Lines[2*i+1]=line(b,a);
        Lines[2*i+1].idx=i,Lines[2*i+1].flag=0;
    }
    sort(Lines,Lines+2*n,cmp);
    for(int i=0;i<2*n;i++)
    {
        if(Lines[i].flag==1&&s_intersection(Lines[i].u,Lines[i].v))
        {
            cur=Lines[i].u;
            S.insert(Lines[i]);
        }
    }
    for(int i=0;i<2*n;i++)
    {
        if(Lines[i].flag==1)
        {
            cur=Lines[i].u;
            S.insert(Lines[i]);
        }
        else
        {
            line tmp=Lines[i];
            tmp.flag=1;
            swap(tmp.u,tmp.v);
            tmp.ang=atan2(tmp.u.y,tmp.u.x);
            S.erase(tmp);
        }
        if(!S.empty())
        {
            set<line>::iterator itor=S.begin();
            ok[(*itor).idx]=1;
        }
    }
    int ans=0;
    for(int i=0;i<n;i++)if(ok[i])ans++;
    printf("%d\n",ans);
}

int main()
{

    while(~scanf("%d",&n))
    {
        input();
    }
    return 0;
}

B题实际上我还没有想清楚，刚刚写的代码还没来得及执行电脑就蓝屏了。。。等想清楚了会更新该博客的。

大学要有梦想
关注
专栏目录
计算几何进阶
supreme567的博客
10-13 186
基础知识 2983. 玩具 #include<bits/stdc++.h> using namespace std; #define int long long typedef pair<int, int> PII; #define x first #define y second //#define d lld int n; PII a[5050],b[5050];//a在下面，b在上面 int ans[5050]; int cross(int x1,int y1,int x2,
计算机视觉基础-图像处理-几何变换
water19111213的博客
04-23 1358
本文来源于Datawhale组队学习的教材手册，供大家参考阅读。 https://github.com/datawhalechina/team-learning/blob/master/%E8%AE%A1%E7%AE%97%E6%9C%BA%E8%A7%86%E8%A7%89%E5%9F%BA%E7%A1%80%EF%BC%9A%E5%9B%BE%E5%83%8F%E5%A4%84%E7%90%8...
计算几何扫描线算法sweep─line算法C++代码实现（CGAL实现的例子）
03-26
基于计算几何算法库实现的Sweep-line算法，能够报告直线的交点（CGAL实现的例子）
计算几何：扫描线
wangqianqianya的博客
08-27 273
【模板】扫描线 theme:给定n个矩阵的左上角和右下角坐标,求覆盖总面积 //给定n个矩阵的左上角和右下角坐标,求覆盖总面积 #include<bits/stdc++.h> using namespace std; #define ls (rt << 1)//左儿子 #define rs (rt << 1 | 1)//右儿子 #define int ...
HDU 4116 Fruit Ninja （ 计算几何 + 扫描线 ）
aegt43853的博客
10-04 142
给你最多1000个圆，问画一条直线最多能与几个圆相交，相切也算。 显然临界条件是这条线是某两圆的公切线，最容易想到的就是每两两圆求出所有公切线，暴力判断一下。 可惜圆有1000个，时间复杂度太高。 网上题解的做法是枚举每个“中心圆”，求出这个圆与其他圆的切线，然后按极角排序，扫一圈。 把每条切线看成扫入线——添加一个圆，或扫出线——删除一个圆。 形象一点就是一条与中心圆相切的直...
HDU 3867 计算几何扫描线
ACM交流帖
08-15 411
#include #include #include #include #include #include #include #define inf 0x7fffffff #define exp 1e-10 #define PI 3.141592654 using namespace std; const int MAXN=31110; struct Point { double x,y;
HDU3867 计算几何扫描线
CandyLien
08-30 1011
转载自AcCry 题意：原子弹爆炸，一些互不相交的线段，求能辐射到的线段（可以将原子弹爆炸点视为泛光源） 以辐射源为中心对周围的点按照极坐标角度进行排序，然后在极坐标上使用扫描线方法。 维护一个集合，集合内的元素是与扫描线相交的线段，排序依据是线段与扫描线的交点到辐
计算几何算法详解：从基础到进阶
"该资源提供计算几何的基础算法，包括点、线段、多边形、圆和矩形的各种运算，采用C语言实现并带有详细注释。" 在计算几何领域，理解和实现这些基础算法至关重要，因为它们是解决更复杂问题的基础。以下是其中的...
ACM计算几何模板
03-31
### ACM计算几何模板知识点概述 本文档是一份详尽的计算几何模板，共计74页，涉及了计算几何领域的多个方面，包括基本的几何公式、直线与线段处理方法、多边形操作、三角形特性分析、圆的相关计算、球面几何以及三...
【免费下载】 FME进阶视频教程： FME常用转换器之几何创建类，讲解如何在FME中创建几何图形
01-21
本视频教程——“FME进阶视频教程：FME常用转换器之几何创建类”专注于介绍如何在FME工作空间中创建和操作几何图形。教程采用的是FME的2020.2中文版，适合对FME有一定基础了解或希望深入学习FME功能的用户。 在FME...
计算几何系列 —— 扫描线算法
Zz_0913的博客
12-21 6259
显然，只要一条线段被覆盖，那么它肯定被图形所截。所以，整个问题就转化为了一个区间查询问题，即：每次将当前扫描线扫到的边对应的信息按照之前赋上的权值更新，然后再查询线段树根节点的信息，最后得到当前扫描线扫过的面积。这就可以用线段树来实现了（线段树：顾名思义，就是树上的结点为一条条线段~），接下来我们来简单看一下模拟过程：
计算几何-扫描线
wyzz_yz的博客
03-27 244
第二步 ：每一个矩形的下边权值为1(代表入边)，上边权值为-1(出边)，结合线段树维护区间，当扫到①时a[1]到a[2]这一段长度更新，扫描线长度a[2]-a[1],当扫到②边a[2]到a[3]更新，a[3]到a[4]，a[2]到a[3]标记累计，扫描线长度a[4]-a[1]。当扫到③边a[1]到a[3]更新，扫描线长度a[4]-a[2]，有一个发现a[2]到a[3]总共被累计两次，所以即便③边出边了一段a[2]到a[3]但是扫描线的长度还是有这段距离的。正解 下面完整代码👇。今日正题基础应用一面积并。
计算几何-扫描线算法
toward_south的博客
05-19 746
计算几何中，扫描线算法（Sweep Line Algorithm）或平面扫描算法（Plane Sweep Algorithm）是一种算法模式，虚拟扫描线或扫描面来解决欧几里德空间中的各种问题，一般被用来解决图形面积，周长等问题，是计算几何中的关键技术之一。这种算法背后的想法是想象一条线（通常是一条垂直线）在平面上扫过或移动，在某些点停止。几何操作仅限于几何对象，无论何时停止，它们都与扫描线相交或紧邻扫描线，并且一旦线穿过所有对象，就可以获得完整的解。
hdu 6127 计算几何模拟题（旋转扫描线）
转角的守望
08-16 719
Hard challenge Time Limit: 4000/2000 MS (Java/Others)    Memory Limit: 524288/524288 K (Java/Others) Total Submission(s): 579    Accepted Submission(s): 235 Problem Description There are n p
【计算几何】扫描线+线段树
gzkeylucky的博客
09-11 390
学习扫描线算法解决矩形面积交、矩形面积并、矩形周长并
【算法】-扫描线
最新发布
knoci的博客
05-30 2699
在计算几何中，扫描线算法（scan line algorithm）一般用来解决几何图形的面积交并，周长交并问题，扫描线算法的核心思想是利用扫描线（通常是水平线或垂直线）在几何空间中“扫描”对象，以确定哪些对象与扫描线相交。下面我们就来通过求矩形的面积并来介绍扫描线算法。先来看看怎么求下面图形的面积并：传统算法是两个矩形面积相加减去重合的面积：但是这样算非常的耗费时间，因为每个矩形都需要两两配对，查看互相之间是否有交集。那么我们接下来就想着把矩形分成三部分：于是现在就变成了。
bzoj 1845: [Cqoi2005] 三角形面积并 （扫描线+计算几何）
clover_hxy的博客
04-25 517
题目描述传送门题目大意：给出n个三角形，求它们并的面积。题解MD被卡了一下午精度，各种不爽！！！ 同样被卡的同学建议使用long double试一试。。。。 我们把三角形的顶点以及所有交点的横坐标离散，那么相邻两列之间的图形要么是梯形要么是三角形，只要统计出两列上的被覆盖区域的长度，就可以用(上底+下底)*高/2来计算。 统计覆盖区域的话，找出所有的区间，左端点排序后扫一遍就可以了。 有一种
计算几何中的扫描（sweeping)算法（转自百度）
CandyLien
09-05 2720
要解决的问题是判断一大堆线段中是否存在一对线段相交。 个人感觉计算几何的扫描算法算是属于解题思想上的东西，类似常规问题的搜索，是遍历处理对象的一种方法。 算法的思想是安照某种方式组织要被处理的对象，比如点、线段。组织方式通常就是从上到下、从左到右等。具体的办法是假想一
【计算几何】【扫描法】The Sky is the Limit
BUAA_Alchemist的博客
02-13 486
 传送门: http://lx.lanqiao.cn/problem.page?gpid=T483 问题描述 　　Banff城雇用了一家广告公司来提升这座城市对潜在的游客的吸引力。其中一个计划中的口号声称延伸在这座城市周围的山脉组成了加拿大最美丽的天际线。但是加拿大消费者保护协会认为“最美丽的天际线”是一种主观的，无法证实的声称，而且可能因此让人误解。 　　然后那个广告公司就想出了一个口号“Ban...
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值