【算法设计与分析 李春葆】计算几何（一）

于 2022-09-11 10:42:05 发布
阅读量262 收藏 1
点赞数 1
分类专栏： # ACM进阶训练（进阶） 文章标签： 算法 c++
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/xjsc01/article/details/126803148
版权

这里讨论的是二维平面中的算法

基础常识

这里一改之前草率的作风，使用类来进行构造。

在这里默认把一个点看做还是一个以原点为起点的向量。

最开始的类

#include <bits/stdc++.h>
using namespace std;
class Point{
public:
    double x, y;
    Point(){
        x = 0, y = 0;
    }//默认构造函数
    Point(double x1, double y1){//重载构造函数
        x = x1;
        y = y1;
    }
    void disp()
    {
        printf("(%g, %g)", x, y);
    }
};
int main()
{
    Point p1(1, 2);
    p1.disp();   
    return 0;
}

实现向量的相关简单运算

//加运算
Point operator +(const Point &p1, const Point &p2)
{
    return Point(p1.x + p2.x, p1.y + p2.y);
}

//减运算
Point operator -(const Point &p1, const Point &p2)
{
    return Point(p1.x - p2.x, p1.y - p2.y);
}

//相等
bool operator ==(const Point &p1, const Point &p2)
{
    if(p1.x == p2.x && p1.y == p2.y)
        return true;
    else 
        return false;
}

//点积运算
double Dot(const Point &p1,const Point &p2)
{
    return p1.x*p2.x + p1.y * p2.y;
}
/*
通过点积可以判断两个向量的夹角
*/

//求向量的模
double Length(const Point &p)
{
    return sqrt(Dot(p, p));
}

//对于具有相同起点的两个向量的夹角的计算
int Angle(const Point &p0, const Point &p1, const Point &p2)
{
    double t = Dot((p1-p0), (p2-p0));
    if(t > 1e-8) return 1;//锐角
    else if(t < -(1e-8)) return -1;//钝角
    else return 0;//直角
}

//pay attention to：
//如果是一个二维向量，那么叉积为一个标量
double Det(const Point &p1, const Point p2)
{
    return p1.x * p2.y - p1.y * p2.x;
}
//两点之间的距离（欧式距离）
double Distance(const Point &p1, const Point &p2)
{
    return Length(p2-p1);
}

叉积详解

叉积在上面已经提到过。叉积在运算中具有很重要的地位。

叉积的算法：p1.x * p2.y - p1.y * p2.x;

叉积的意义：

  1. 绝对值大小：表示(0, 0), p1, p2, p1+p2这四个点所组成的平行四边形的面积。
  2. 符号：（前提：Det(p1, p2)）
    正：p1逆时针到p2
    负：p1逆时针到p2
    0：共线，但是可能同向，也可能反向

如果对于p0->p1, p0->p2

p0, p1, p2在右手螺旋方向上Det( (p1-p0), (p2-p0) ) > 0

p0, p1, p2在左手螺旋方向上Det( (p1-p0), (p2-p0) ) < 0

当仅仅使用叉积的大小时，别忘记加绝对值

点到线段的距离

有以下三种情况

image

double DistPtoSegment(const Point &p0, const Point &p1, const Point &p2)
{
    if(p1 == p2)
    {
        return Distance(p1, p0);
    }
    if(Dot((p0 - p1), (p2-p1)) < 0)
    {
        return Distance(p0, p1);
    }
    else if(Dot(p1-p2, p0-p2) < 0)
    {
        return Distance(p0, p2);
    }
    return (fabs(Det(p2-p1, p0-p1))/Distance(p1, p2));
}

判断一个点是不是在一个矩形里面

假设p1以及p2是矩形的左上角以及右下角

bool InRectangle(const Point &p0, const Point &p1, const Point &p2)
{
   return Dot(p2-p0, p1-p0) <= 0;
}

判断一个点是不是在一条线段上

注意：限制一个点的范围可以使用矩形框起来！！！

bool OnSegment(const Point &p0, const Point &p1, const Point &p2)
{  
    return Det(p1-p0, p2-p0) == 0 && Dot(p1-p0, p2-p0) <= 0;
}

实际情况下，等于号应该是取一个精度

判断两条线段是否平行

bool Parallel(const Point &p1, const Point &p2, const Point &p3, const Point &p4)
{
    return Det(p2-p1, p4-p3) == 0;
}

判断两条线段是否相交

这里这一个思路比较巧妙，见代码

  1. p1,p2在线段p3p4的两端，并且p3,p4在线段p1p2的两端
  2. 点在线上

image

bool SegIntersect(const Point &p1, const Point &p2, const Point &p3, const Point &p4)
{
    if(Det(p3-p1, p2-p1)*Det(p4-p1, p2-p1) < 0 && Det(p1-p4, p3-p4)*Det(p2-p4, p3-p4)<0)
        return true;
    else if(Det(p3-p1, p2-p1) == 0 && OnSegment(p3, p1, p2))
        return true;
    else if(Det(p4-p1, p2-p1) == 0 && OnSegment(p4, p1, p2))
        return true;
    else if(Det(p1-p4, p3-p4) == 0 && OnSegment(p1, p2, p4))
        return true;
    else if(Det(p2-p4, p3-p4) == 0 && OnSegment(p2, p3, p4))
        return true;
    else
        return false;
}

判断一个点是否在多边形内

简单多边形：所有的边不相交。（这里讨论的情况就是简单多边形）

要求：判断p0是否在a[0...n]中（a[0] == a[n]）[包含边界]

我们需要在这一个点做一条向右的射线，如果与多边形的边有奇数个交点，那么这一个点就在多边形的内部，否则就在外部。

教课书是一般不会错的，每一个异样处都有原因！

double x = (p0.y - p1.y) * (p2.x - p1.x) / (p2.y - p1.y) + p1.x;

这是算出p0向右发出的射线与多边形的这一条边的交点

bool PointInPolygon(Point p0, vector<Point> a)
{
    int cnt = 0;
    for(int i = 0; i < a.size()-1; i++)//我改了一下，感觉书上有越界现象
    {
        Point p1 = a[i], p2 = a[i+1];
        if(OnSegment(p0, p1, p2)) return true;//如果在多边形的边上
        if(p1.y == p2.y) continue;//跳过水平线
        if(p0.y < p1.y && p0.y < p2.y) continue;
        if(p0.y >= p1.y && p0.y >= p2.y) continue;//这里有玄机，等于号不要删除
        double x = (p0.y - p1.y) * (p2.x - p1.x) / (p2.y - p1.y) + p1.x;
        if(x > p0.x) cnt++;
    }
    return (cnt & 1);
}

对于“玄机”的解释：

先上一段代码：

bool PointInPolygon(Point p0, vector<Point> a)
{
    int cnt = 0;
    for(int i = 0; i < a.size()-1; i++)
    {
        Point p1 = a[i], p2 = a[i+1];
        if(OnSegment(p0, p1, p2)) return true;
        if(p1.y == p2.y) continue;
        if(p0.y < p1.y && p0.y < p2.y) continue;
        if(p0.y > p1.y && p0.y > p2.y) continue;//等于号已经删除！！！！！
        double x = (p0.y - p1.y) * (p2.x - p1.x) / (p2.y - p1.y) + p1.x;
        if(x > p0.x) cnt++;
    }
    return (cnt & 1);
}


int main()
{
    vector<Point> a;
    a.push_back(Point(-2, 0));
    a.push_back(Point(0, 2));
    a.push_back(Point(2, 0));
    a.push_back(Point(1, 0));
    a.push_back(Point(0, -1));//更改后面的坐标为-1和1，观察结果
    a.push_back(Point(-1, 0));
    a.push_back(Point(-2, 0));
    bool ret = PointInPolygon(Point(0, 0), a);
    cout << ret;
    return 0;
}

正是由于如果一个image

三个点构成的三角形的面积

利用向量的叉积处理

double TriangleArea(const Point &p0, const Point &p1, const Point &p2)
{
    return fabs(Det(p1-p0, p2-p0))/2;
}

求一个多边形的面积

这里把这一个多边形分成多块三角形，然后分别计算面积。

注意：采用这一种方法可以有效地处理凹形的问题。

double PolyArea(vector<Point> p)
{
    double ans = 0.0;
    for(int i = 1; i < p.size()-1; i++)
    {
        ans += Det(p[i] - p[0], p[i+1] - p[0]);
        //不可以使用 TriangleArea(p[0], p[i], p[i+1]);，因为三角形的面积是带有绝对值的
    }
    return fabs(ans) / 2;
}

最终代码

#include <bits/stdc++.h>
using namespace std;
#include <bits/stdc++.h>
using namespace std;
class Point{
public:
    double x, y;
    Point(){
        x = 0, y = 0;
    }//默认构造函数
    Point(double x1, double y1){//重载构造函数
        x = x1;
        y = y1;
    }
    void disp()
    {
        printf("(%g, %g)", x, y);
    }
};
//加运算
Point operator +(const Point &p1, const Point &p2)
{
    return Point(p1.x + p2.x, p1.y + p2.y);
}

//减运算
Point operator -(const Point &p1, const Point &p2)
{
    return Point(p1.x - p2.x, p1.y - p2.y);
}

//相等
bool operator ==(const Point &p1, const Point &p2)
{
    if(p1.x == p2.x && p1.y == p2.y)
        return true;
    else 
        return false;
}

//点积运算
double Dot(const Point &p1,const Point &p2)
{
    return p1.x*p2.x + p1.y * p2.y;
}
/*
通过点积可以判断两个向量的夹角
*/

//求向量的模
double Length(const Point &p)
{
    return sqrt(Dot(p, p));
}

//对于具有相同起点的两个向量的夹角的计算
int Angle(const Point &p0, const Point &p1, const Point &p2)
{
    double t = Dot((p1-p0), (p2-p0));
    if(t > 1e-8) return 1;//锐角
    else if(t < -(1e-8)) return -1;//钝角
    else return 0;//直角
}

//pay attention to：
//如果是一个二维向量，那么叉积为一个标量
double Det(const Point &p1, const Point p2)
{
    return p1.x * p2.y - p1.y * p2.x;
}

//两点之间的距离（欧式距离）
double Distance(const Point &p1, const Point &p2)
{
    return Length(p2-p1);
}

//点到线段的距离
double DistPtoSegment(const Point &p0, const Point &p1, const Point &p2)
{
    if(p1 == p2)
    {
        return Distance(p1, p0);
    }
    if(Dot((p0 - p1), (p2-p1)) < 0)
    {
        return Distance(p0, p1);
    }
    else if(Dot(p1-p2, p0-p2) < 0)
    {
        return Distance(p0, p2);
    }
    return (fabs(Det(p2-p1, p0-p1))/Distance(p1, p2));
}

//判断一个点是不是在一个矩形里面,假设`p1`以及`p2`是矩形的左上角以及右下角
bool InRectangle(const Point &p0, const Point &p1, const Point &p2)
{
   return Dot(p2-p0, p1-p0) <= 0;
}

//判断一个点是不是在一条线段上
bool OnSegment(const Point &p0, const Point &p1, const Point &p2)
{  
    return Det(p1-p0, p2-p0) == 0 && Dot(p1-p0, p2-p0) <= 0;
}

//判断两条线段是否平行
bool Parallel(const Point &p1, const Point &p2, const Point &p3, const Point &p4)
{
    return Det(p2-p1, p4-p3) == 0;
}

//判断两条线段是否相交
bool SegIntersect(const Point &p1, const Point &p2, const Point &p3, const Point &p4)
{
    if(Det(p3-p1, p2-p1)*Det(p4-p1, p2-p1) < 0 && Det(p1-p4, p3-p4)*Det(p2-p4, p3-p4)<0)
        return true;
    else if(Det(p3-p1, p2-p1) == 0 && OnSegment(p3, p1, p2))
        return true;
    else if(Det(p4-p1, p2-p1) == 0 && OnSegment(p4, p1, p2))
        return true;
    else if(Det(p1-p4, p3-p4) == 0 && OnSegment(p1, p2, p4))
        return true;
    else if(Det(p2-p4, p3-p4) == 0 && OnSegment(p2, p3, p4))
        return true;
    else
        return false;
}

//判断一个点是否在多边形内
bool PointInPolygon(Point p0, vector<Point> a)
{
    int cnt = 0;
    for(int i = 0; i < a.size()-1; i++)//我改了一下，感觉书上有越界现象
    {
        Point p1 = a[i], p2 = a[i+1];
        if(OnSegment(p0, p1, p2)) return true;//如果在多边形的边上
        if(p1.y == p2.y) continue;//跳过水平线
        if(p0.y < p1.y && p0.y < p2.y) continue;
        if(p0.y >= p1.y && p0.y >= p2.y) continue;//这里有玄机，等于号不要删除
        double x = (p0.y - p1.y) * (p2.x - p1.x) / (p2.y - p1.y) + p1.x;
        if(x > p0.x) cnt++;
    }
    return (cnt & 1);
}

//三个点构成的三角形的面积
double TriangleArea(const Point &p0, const Point &p1, const Point &p2)
{
    return fabs(Det(p1-p0, p2-p0))/2;
}


//求一个多边形的面积
double PolyArea(vector<Point> p)
{
    double ans = 0.0;
    for(int i = 1; i < p.size()-1; i++)
    {
        ans += Det(p[i] - p[0], p[i+1] - p[0]);
        //不可以使用 TriangleArea(p[0], p[i], p[i+1]);，因为三角形的面积是带有绝对值的
    }
    return fabs(ans) / 2;
}
int main()
{
    vector<Point> a;
    a.push_back(Point(-2, 0));
    a.push_back(Point(0, 2));
    a.push_back(Point(2, 0));
    a.push_back(Point(1, 0));
    a.push_back(Point(0, 1));//
    a.push_back(Point(-1, 0));
    a.push_back(Point(-2, 0));
    double ret = TriangleArea(a[0], a[1], a[2]);
    cout << ret;
    return 0;
}
xjsc01
关注
  • 1
    点赞
  • 踩
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
专栏目录
计算几何——算法设计与分析(第3版)
09-02
老周的《计算几何——算法设计与分析》第三版，有算法思路的分析，学习计算几何相当应手。
计算几何学习总结（使用教材算法设计与分析（第二版））
qq_44934958的博客
07-05 752
** 计算几何总结（使用教材算法设计与分析（第二版）李春葆 清华大学出版社） ** （代码大部分为书中原有代码，如有雷同，实属正常。） #include<bits/stdc++.h> using namespace std;//万能头 #include<fundament.h>//向量计算的函数库 本次学习主要主要有以下几个方面：1、向量计算  2...
算法设计与分析：图与计算几何算法
qq_42339734的博客
12-18 343
图和计算几何 有n个点、m条无向边的图，每条边都有长度d和花费p,再给出一个七点s和一个重点t，要求输出起点到终点的最短距离及其花费，如果最短距离有多条路线，则输出花费最少的。 给定三角形的3条边，a、b、c，判断该三角形的类型。 直接利用狄克斯特拉算法求从顶点s到顶点t花费最小的最短路径。在狄克斯特拉算法中要做两点修改，一是增加记录路径最小花费的数字cost，cost[j]表示从顶点s到顶点j...
计算几何-算法设计与分析
01-14
《计算几何：算法设计与分析》系统地介绍了计算几何中的基本概念、求解诸多问题的算法及复杂性分析，概括了求解几何问题所特有的许多思想方法、几何结构与数据结构。全书共分10章，包括：预备知识，几何查找（检索），多边形，凸壳及其应用，Voronoi图、三角剖分及其应用，交与并及其应用，多边形的获取及相关问题，几何体的划分与等分，路径与回路，几何拓扑网络设计等。 《计算几何：算法设计与分析》可作为高等院校计算机、自动化等专业研究生或本科高年级学生的教材或教学参考书，也可供软件开发人员、相关专业科技工作者参考。
计算几何——算法分析与设计周培德.pdf
05-13
计算几何——算法分析与设计周培德.pdf 计算几何——算法分析与设计周培德.pdf 计算几何——算法分析与设计周培德.pdf
李春葆.zip 李春葆 算法设计与分析2nd习题答案代码课件
12-20
在《李春葆 算法设计与分析2nd》中，读者将学习到如何设计高效、准确的算法，并学习分析算法的时间复杂度和空间复杂度，这对于优化程序性能和理解算法效率至关重要。 课件部分可能包含PPT或PDF文档，涵盖了课程中的...
算法设计与分析习题答案 .pdf
06-02
总结的关于中科大研究生课程算法设计与分析习题答案，包括分治法、动态规划、贪心算法、回溯、分支限界等章节内容
算法设计与分析课后习题答案（李春保）
05-07
算法设计与分析（李春葆）版的课后习题答案，我在网上找了好久都没有找到的
计算机算法设计和分析习题及答案解析.pdf
05-24
计算机算法设计和分析是计算机科学中一个重要的领域，它涉及到算法的设计、分析和实现。本文档提供了计算机算法设计和分析习题及答案解析，涵盖了算法设计和分析的基础知识和重要概念。 选择题： 1. 二分搜索算法...
算法设计与分析(第二版）李春保课后习题答案
12-29
算法设计与分析是计算机科学中的核心课程之一，它主要研究如何高效地解决问题并实现这些解决方案。本章概论部分的习题涵盖了算法的基本概念、效率评估以及算法设计的基础技巧。 1. 算法的基本特性包括：有穷性...
算法分析与设计pdf
04-19
算法分析与设计内容包含：递归与分治，动态规划，贪心算法
算法设计与分析（详细解析（含源代码））
05-21
常用算法设计方法详细解析（含源代码） 算法是问题求解过程的精确描述，一个算法由有限条可完全机械地执行的、有确定结果的指令组成。指令正确地描述了要完成的任务和它们被执行的顺序。计算机按算法指令所描述的顺序执行算法的指令能在有限的步骤内终止，或终止于给出问题的解，或终止于指出问题对此输入数据无解。 通常求解一个问题可能会有多种算法可供选择，选择的主要标准是算法的正确性和可靠性，简单性和易理解性。其次是算法所需要的存储空间少和执行更快等。 算法设计是一件非常困难的工作，经常采用的算法设计技术主要有迭代法、穷举搜索法、递推法、贪婪法、回溯法、分治法、动态规划法等等。另外，为了更简洁的形式设计和藐视算法，在算法设计时又常常采用递归技术，用递归描述算法。 一、迭代法 二、穷举搜索法 三、递推法 四、递归 五、回溯法 六、贪婪法 七、分治法 八、动态规划法
算法设计与分析 清华大学出版社
11-24
内容包括ACM程序中设计的很多算法 递归 排序 贪婪 随机 动态规划等基本算法
计算机算法设计与分析(第2版) 习题解答
05-28
计算机算法设计与分析(第2版) 习题解答
算法设计与分析之计算几何
qq_40831898的博客
01-05 414
向量的加法运算 对于两个点表示的向量p1和p2（起点均为原点(0，0)），向量加法定义为：p1+p2=(p1.x+p2.x，p1.y+p2.y)，其结果仍为一个向量。 向量加法一般可用平行四边形法则，如下图所示，两个向量为p1(2，-1)，p2(3，3)，则p3=p1+p2=(5，2)。 算法思路： Point operator +(const Point &p1，const Point...
计算机图形学常用算法
03-20 3190
（引自：http://blog.csdn.net/orbit/article/details/7082678  作者：orbit） 我的专业是计算机辅助设计（CAD），算是一半机械一半软件，《计算机图形学》是必修课，也是我最喜欢的课程。热衷于用代码摆平一切的我几乎将这本教科书上的每种算法都实现了一遍，这种重复劳动虽然意义不大，但是收获很多，特别是丢弃了多年的数学又重新回到了脑袋中，算是最
算法设计与分析李春葆pdf
最新发布
07-16
《算法设计与分析李春葆pdf》是由李春葆所编著的一本关于算法设计与分析的参考书籍，该书以系统全面的方式介绍了算法设计与分析的相关概念、技巧和方法。 该书第一部分主要介绍了算法与问题求解的基础知识，包括...
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值