计算几何 - 二维几何基础 （模板）

最新推荐文章于 2019-08-09 16:23:50 发布
最新推荐文章于 2019-08-09 16:23:50 发布
阅读量1.1k 收藏
点赞数 1
分类专栏： ~~~~~~~~~~计算几何 各种模板 文章标签： ACM 模板 计算几何
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/u014355480/article/details/43605257
版权
这篇博客详述了二维几何的基本概念，包括点积、叉积、旋转、单位法线向量、二直线交点计算、点到直线和线段的距离、点在线段上的投影及线段相交判断等关键知识点，还涵盖了多边形面积的计算。是ACM竞赛和计算几何学习者的重要参考资料。
摘要由CSDN通过智能技术生成



这里总结一下二维几何基础知识！


常用定义：

//定义点的类型 
struct Point {
	double x, y;
	Point(double x = 0, double y = 0) : x(x) , y(y) { }  //构造函数，方便代码编写 
};

typedef Point Vector;  //从程序实现上，Vector只是Point的别名 

//向量 + 向量 = 向量 ，点 + 向量 = 点
Vector operator + (Vector A, Vector B) { return Vector(A.x+B.x, A.y+B.y); }
//点 - 点 = 向量
Vector operator - (Vector A, Vector B) { return Vector(A.x-B.x, A.y-B.y); }
//向量 * 数 = 向量 
Vector operator * (Vector A, double p) { return Vector(A.x*p, A.y*p); }
//向量 / 数 = 向量 
Vector operator / (Vector A, double p) { return Vector(A.x/p, A.y/p); } 

bool operator < (const Point& a, const Point& b) {
	return a.x < b.x || (a.x == b.x && a.y < b.y);
} 

const double eps = 1e-10;
int dcmp(double x) {
	if(fabs(x) < eps) return 0; else return x < 0 ? -1 : 1;
}

bool operator == (const Point& a, const Point& b) {
	return dcmp(a.x - b.x) == 0 && dcmp(a.y - b.y) == 0;
}



点积：

//点积 
double Dot(Vector A, Vector B) { return A.x*B.x + A.y*B.y; } //求点积 
double Length(Vector A) { return sqrt(Dot(A, A)); }			 //求向量长度 
double Angle(Vector A, Vector B) { return acos(Dot(A, B) / Length(A) / Length(B)); }//求向量之间的夹角



叉积：

//叉积 
double Cross(Vector A, Vector B) { return A.x*B.y - A.y*B.x; }//求叉积 
double Area2(Point A, Point B, Point C) { return Cross(B-A, C-A); }//根据叉积求三角形面积的两倍



旋转：

//旋转 
Vector Rotate(Vector A, double rad) {//rad是弧度 
	return Vector(A.x*cos(rad) - A.y*sin(rad), A.x*sin(rad)+A.y*cos(rad) );
}



向量的单位法线：

//向量单位法向量，调用前请确保A不是零向量 
Vector Normal(Vector A) {  
    double L = Length(A);  
    return Vector(-A.y/L, A.x/L);  
}



二直线交点：

//二直线交点（参数式） 
Point GetLineIntersection(Point P, Vector v, Point Q, Vector w) {
	Vector u = P - Q;
	double t = Cross(w, u) / Cross(v, w);
	return P + v * t;
}



点到直线距离 ：

//点到直线距离  
double DistanceToLine(Point P, Point A, Point B) {  
    Vector v1 = B-A, v2 = P - A;  
    return fabs(Cross(v1,v2) / Length(v1));  //如果不取绝对值，得到的是有向距离 
}



点到线段距离  ：

//点到线段距离  
double DistanceToSegment(Point P, Point A, Point B) {  
    if(A==B) return Length(P-A);  
    Vector v1 = B - A, v2 = P - A, v3 = P - B;  
    if(dcmp(Dot(v1, v2)) < 0) return Length(v2);  
    else if(dcmp(Dot(v1, v3)) > 0) return Length(v3);  
    else return fabs(Cross(v1, v2)) / Length(v1);  
}



点在直线上的投影：

//点在直线上的投影
Point GetLineProjection(Point P, Point A, Point B) {
	Vector v = B - A;
	return A + v * ( Dot(v, P-A) / Dot(v, v) ); 
}



线段相交判定：

//线段相交判定
bool SegmentProperIntersection(Point a1, Point a2, Point b1, Point b2) {
	double c1 = Cross(a2 - a1, b1 - a1), c2 = Cross(a2 - a1, b2 - a1),
			c3 = Cross(b2 - b1, a1 - b1), c4 = Cross(b2 - b1, a2 - b1);
	return dcmp(c1) * dcmp(c2) < 0 && dcmp(c3) * dcmp(c4) < 0;
}



判断一个点是否在一条线段上：

//判断一个点是否在一条线段上
bool OnSegment(Point p, Point a1, Point a2) {
	return dcmp(Cross(a1 - p, a2 - p)) == 0 && dcmp(Dot(a1 - p, a2 - p)) < 0;
}



多边形面积：

//多边形面积  
double ConvexPolygonArea(Point* p, int n) {  
    double area = 0;  
    for(int i = 1; i < n-1; i++)  
        area += Cross(p[i] - p[0], p[i + 1] - p[0]);  
    return area / 2;  
}



总结：

//定义点的类型 
struct Point {
	double x, y;
	Point(double x = 0, double y = 0) : x(x) , y(y) { }  //构造函数，方便代码编写 
};

typedef Point Vector;  //从程序实现上，Vector只是Point的别名 

//向量 + 向量 = 向量 ，点 + 向量 = 点
Vector operator + (Vector A, Vector B) { return Vector(A.x+B.x, A.y+B.y); }
//点 - 点 = 向量
Vector operator - (Vector A, Vector B) { return Vector(A.x-B.x, A.y-B.y); }
//向量 * 数 = 向量 
Vector operator * (Vector A, double p) { return Vector(A.x*p, A.y*p); }
//向量 / 数 = 向量 
Vector operator / (Vector A, double p) { return Vector(A.x/p, A.y/p); } 

bool operator < (const Point& a, const Point& b) {
	return a.x < b.x || (a.x == b.x && a.y < b.y);
} 

const double eps = 1e-10;
int dcmp(double x) {
	if(fabs(x) < eps) return 0; else return x < 0 ? -1 : 1;
}

bool operator == (const Point& a, const Point& b) {
	return dcmp(a.x - b.x) == 0 && dcmp(a.y - b.y) == 0;
}

//点积 
double Dot(Vector A, Vector B) { return A.x*B.x + A.y*B.y; } //求点积 
double Length(Vector A) { return sqrt(Dot(A, A)); }			 //求向量长度 
double Angle(Vector A, Vector B) { return acos(Dot(A, B) / Length(A) / Length(B)); }//求向量之间的夹角 

//叉积 
double Cross(Vector A, Vector B) { return A.x*B.y - A.y*B.x; }//求叉积 
double Area2(Point A, Point B, Point C) { return Cross(B-A, C-A); }//根据叉积求三角形面积的两倍 

//旋转 
Vector Rotate(Vector A, double rad) {//rad是弧度 
	return Vector(A.x*cos(rad) - A.y*sin(rad), A.x*sin(rad)+A.y*cos(rad) );
} 

//向量单位法向量，调用前请确保A不是零向量 
Vector Normal(Vector A) {  
    double L = Length(A);  
    return Vector(-A.y/L, A.x/L);  
}

//二直线交点（参数式） 
Point GetLineIntersection(Point P, Vector v, Point Q, Vector w) {
	Vector u = P - Q;
	double t = Cross(w, u) / Cross(v, w);
	return P + v * t;
} 

//点到直线距离  
double DistanceToLine(Point P, Point A, Point B) {  
    Vector v1 = B-A, v2 = P - A;  
    return fabs(Cross(v1,v2) / Length(v1));  //如果不取绝对值，得到的是有向距离 
}  
  
//点到线段距离  
double DistanceToSegment(Point P, Point A, Point B) {  
    if(A==B) return Length(P-A);  
    Vector v1 = B - A, v2 = P - A, v3 = P - B;  
    if(dcmp(Dot(v1, v2)) < 0) return Length(v2);  
    else if(dcmp(Dot(v1, v3)) > 0) return Length(v3);  
    else return fabs(Cross(v1, v2)) / Length(v1);  
}  
  
//点在直线上的投影
Point GetLineProjection(Point P, Point A, Point B) {
	Vector v = B - A;
	return A + v * ( Dot(v, P-A) / Dot(v, v) ); 
}  

//线段相交判定
bool SegmentProperIntersection(Point a1, Point a2, Point b1, Point b2) {
	double c1 = Cross(a2 - a1, b1 - a1), c2 = Cross(a2 - a1, b2 - a1),
			c3 = Cross(b2 - b1, a1 - b1), c4 = Cross(b2 - b1, a2 - b1);
	return dcmp(c1) * dcmp(c2) < 0 && dcmp(c3) * dcmp(c4) < 0;
} 

//判断一个点是否在一条线段上
bool OnSegment(Point p, Point a1, Point a2) {
	return dcmp(Cross(a1 - p, a2 - p)) == 0 && dcmp(Dot(a1 - p, a2 - p)) < 0;
} 

//多边形面积  
double ConvexPolygonArea(Point* p, int n) {  
    double area = 0;  
    for(int i = 1; i < n-1; i++)  
        area += Cross(p[i] - p[0], p[i + 1] - p[0]);  
    return area / 2;  
}




无注释纯净版微笑：

struct Point {
	double x, y;
	Point(double x = 0, double y = 0) : x(x) , y(y) { }  
};

typedef Point Vector;  

Vector operator + (Vector A, Vector B) { return Vector(A.x+B.x, A.y+B.y); }
Vector operator - (Vector A, Vector B) { return Vector(A.x-B.x, A.y-B.y); }
Vector operator * (Vector A, double p) { return Vector(A.x*p, A.y*p); }
Vector operator / (Vector A, double p) { return Vector(A.x/p, A.y/p); } 

bool operator < (const Point& a, const Point& b) {
	return a.x < b.x || (a.x == b.x && a.y < b.y);
} 

const double eps = 1e-10;
int dcmp(double x) {
	if(fabs(x) < eps) return 0; else return x < 0 ? -1 : 1;
}

bool operator == (const Point& a, const Point& b) {
	return dcmp(a.x - b.x) == 0 && dcmp(a.y - b.y) == 0;
}

double Dot(Vector A, Vector B) { return A.x*B.x + A.y*B.y; } 
double Length(Vector A) { return sqrt(Dot(A, A)); }		
double Angle(Vector A, Vector B) { return acos(Dot(A, B) / Length(A) / Length(B)); } 

double Cross(Vector A, Vector B) { return A.x*B.y - A.y*B.x; }
double Area2(Point A, Point B, Point C) { return Cross(B-A, C-A); }

Vector Rotate(Vector A, double rad) {
	return Vector(A.x*cos(rad) - A.y*sin(rad), A.x*sin(rad)+A.y*cos(rad) );
} 

Vector Normal(Vector A) {  
    double L = Length(A);  
    return Vector(-A.y/L, A.x/L);  
}

Point GetLineIntersection(Point P, Vector v, Point Q, Vector w) {
	Vector u = P - Q;
	double t = Cross(w, u) / Cross(v, w);
	return P + v * t;
} 
 
double DistanceToLine(Point P, Point A, Point B) {  
    Vector v1 = B-A, v2 = P - A;  
    return fabs(Cross(v1,v2) / Length(v1)); 
}  

double DistanceToSegment(Point P, Point A, Point B) {  
    if(A==B) return Length(P-A);  
    Vector v1 = B - A, v2 = P - A, v3 = P - B;  
    if(dcmp(Dot(v1, v2)) < 0) return Length(v2);  
    else if(dcmp(Dot(v1, v3)) > 0) return Length(v3);  
    else return fabs(Cross(v1, v2)) / Length(v1);  
}  

Point GetLineProjection(Point P, Point A, Point B) {
	Vector v = B - A;
	return A + v * ( Dot(v, P-A) / Dot(v, v) ); 
}  

bool SegmentProperIntersection(Point a1, Point a2, Point b1, Point b2) {
	double c1 = Cross(a2 - a1, b1 - a1), c2 = Cross(a2 - a1, b2 - a1),
			c3 = Cross(b2 - b1, a1 - b1), c4 = Cross(b2 - b1, a2 - b1);
	return dcmp(c1) * dcmp(c2) < 0 && dcmp(c3) * dcmp(c4) < 0;
} 

bool OnSegment(Point p, Point a1, Point a2) {
	return dcmp(Cross(a1 - p, a2 - p)) == 0 && dcmp(Dot(a1 - p, a2 - p)) < 0;
} 

double ConvexPolygonArea(Point* p, int n) {  
    double area = 0;  
    for(int i = 1; i < n-1; i++)  
        area += Cross(p[i] - p[0], p[i + 1] - p[0]);  
    return area / 2;  
}




















zzuspy
关注
专栏目录
计算几何模板
skyword的专栏
01-07 299
【整理于2016年】 计算几何模板 const double eps = 1e-8; const double pi = acos(-1.0); using namespace std; //采用eps的精度判断大/小于零 int epssgn(double x) { if (fabs(x)&lt;eps) return 0; else return x&lt;0?-1:1; ...
【计算几何-点的定位】
BlessingXRY的博客
09-07 640
1.二维坐标系，输入点的坐标(xi, yi)，判断点是否在三角形内 判断原理：若点p在三角形内，Sabc = Spab + Spac + Spbc； struct point{/*记录点的信息*/ double x; double y; };struct v{/*记录边的信息*/ point star;/*边的起点*/ point en;/*边的终点*/ };str
二维计算几何模板
weixin_33699914的博客
07-29 116
本模板主要来自刘汝佳的《算法竞赛入门经典》，部分缺省的地方我进行了补齐，且对每个函数作出了解释 /***************************************************************** 点和向量的声明，以及加减，数乘比较等基础操作 ****************************************...
ACM计算几何基础模板，杭电老学长精心整理.pdf
03-25
二维点是计算几何的基础，其包含了坐标信息。向量的定义用于表示具有大小和方向的量，不仅用于表示移动方向，也是进行向量运算的基本单位。点和向量的运算包括了向量的内积、外积、长度计算、夹角计算等，这些都是在...
ACM计算几何模板大全 几何 多边形 凸包 三维 圆
05-22
总的来说，ACM计算几何模板大全提供了从基础到高级的几何问题解决方案，帮助参赛者更好地理解和解决实际问题，提升他们在算法竞赛中的表现。通过深入理解和熟练掌握这些知识点，可以在计算几何的复杂世界中游刃有余...
三维几何模板
08-15
根据给定文件的信息，我们可以深入探讨三维几何在计算几何领域中的应用，特别是在ACM竞赛中的三维几何模板。此模板提供了处理三维空间问题的基础结构和算法，对于理解和解决三维空间中的几何问题至关重要。 ### 三...
山东大学ACM计算几何模板：从基础到高级问题
山东大学ACM模板-计算几何是一份针对算法竞赛（ACM）选手准备的模板，特别专注于计算机图形学中的计算几何部分。这份模板涵盖了广泛的问题类型和核心概念，适合在解决几何问题时作为参考。以下是一些关键知识点的...
计算几何模版
01-17
计算几何中大部分问题经典算法模版. 包括弗洛伊德算法，缔结斯特拉算法等等。
计算几何 模板
weixin_30536513的博客
10-13 97
//凸包+最远点对 #include<iostream> #include<cstdio> #include<cstring> #include<cstdlib> #include<cmath> #include<algorithm> #include<vector> using namespace...
计算几何模板（二维）
_rabbit的专栏
02-26 1108
/**  * 二维ACM计算几何模板  * 注意变量类型更改和EPS  * #include  * #include  * By OWenT  */   const double eps = 1e-8; const double pi = std::acos(-1.0); //点 class point { public:     double x, y;
计算几何 - 二维几何(模板)
Sun
08-09 232
几何意义：   向量的点积(a · b )：a · b = |a||b|cos⁡∅{\cos\emptyset }cos∅ ，若cos⁡∅{\cos\emptyset }cos∅ 为正，两向量之间的夹角为锐角；为负，两向量夹角为钝角；为量，两向量夹角为直角。(b 在 a 方向上的投影)。   向量的叉积(a × b )：a · b = |a||b|sin⁡∅{\sin\emptyset }sin∅...
评论 1
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值