计算几何基础——operator_

前言

东西好多啊$\cdots$还好都不算难

向量的基本表示与运算

向量的定义：可以理解成只有方向和长度，没有位置的量（亦称矢量）

向量的表示：在平面几何中表示为$\overrightarrow{AB}$，但在程序中一般用平面直角坐标系表示法，即 向量$P(x_1,y_1)$代表从原点指向$P$的向量。用这种定义，显然从$A(x_a,y_a)$指向$B(x_b,y_b)$的向量等价于向量$Q(x_b-x_a,y_b-y_a)$。

$PS.$零向量确实存在，即$\overrightarrow{0}$，大小为$0$，方向为任意方向

向量的加减：向量的加减$⟺$坐标的加减，即$P(x_1,y_1)\pm Q(x_2,y_2)=R(x_1\pm x_2,y_1\pm y_2)$显然向量的加减满足交换律、结合律。

向量的叉乘：$P(x_1,y_1)\times Q(x_2,y_2)=x_1{y}_2-x_2{y}_1$叉乘的核心用途为判断位置关系，即若$P\times Q<0$则$P$在$Q$逆时针方向，反之若$P\times Q》0$则$P$在$Q$顺时针方向，若$P\times Q=0$则$P,Q$共线(但不知道是否同向)。显然向量的叉乘结果为一个数且不满足交换律。

接下来依次讨论两个东西的位置关系（点、线段、矩形、圆、多边形$\cdots$）

点与矩形

点为$P(x_0,y_0)$，矩形左下顶点与右上顶点分别为$A(x_1,y_1),B(x_2,y_2)$

点在矩形内：$x_1\le x_0\le x_2$且$y_1\le y_0\le y_2$

复杂度$O(1)$

点与线段

点为$P(x_0,y_0)$，线段两端点分别为$A(x_1,y_1),B(x_2,y_2)$

点在线段上：$\overrightarrow{AP}\times \overrightarrow{AB}=0$且$x_1\le x_0\le x_2$且$y_1\le y_0\le y_2$

复杂度$O(1)$

线段与线段

两条线段两端点分别为$A(x_1,y_1)$与$B(x_2,y_2)$、$M(x_3,y_3)$与$N(x_4,y_4)$

线段与线段相交：$(\overrightarrow{AM}\times \overrightarrow{AB})\ast (\overrightarrow{AN}\times \overrightarrow{AB})\le 0$且$(\overrightarrow{MA}\times \overrightarrow{MN})\ast (\overrightarrow{MB}\times \overrightarrow{MN})\le 0$

复杂度$O(1)$

注意特判平行和首尾相接

射线、直线均可类比求得（更简单）

点与多边形

点为$P(x_0,y_0)$，多边形有$n$条边

点在多边形内：过点$P$引一条平行于坐标轴的射线，枚举每条边求交点数，交点数之和为奇数则在多边形内，反之不在。

复杂度$O(n)$

注意特判平行和交于顶点

线段与多边形

线段两端点分别为$A(x_1,y_1),B(x_2,y_2)$，多边形有$n$条边

线段在多边形内：$A$在多边形内且$B$在多边形内且$AB$与多边形无内交点（枚举求得）

复杂度$O(n)$

注意特判内交与平行

$PS.$老师有复杂度$O(nlogn)$但不用特判的做法，不过个人感觉没必要。

多边形与多边形

这个反而简单

两个多边形分别有$n$、$m$条边

多边形在多边形内：依次判断每条边是否在另一多边形内即可

复杂度$O(nm)$

点与圆

点为$P(x_1,y_1)$，圆心为$O(x_0,y_0)$，半径为$r$

点在圆内：$PO<r$，若想避免精度误差就用$(x_0-x_1{)}^2+(y_0-y_1{)}^2<r^2$

复杂度$O(1)$

线段与圆

圆是凸的

线段两端点分别为$A(x_1,y_1),B(x_2,y_2)$，圆心为$O(x_0,y_0)$，半径为$r$

线段在圆内：$(x_0-x_1{)}^2+(y_0-y_1{)}^2<r^2$且$(x_0-x_2{)}^2+(y_0-y_2{)}^2<r^2$（即端点都在圆内）

复杂度$O(1)$

多边形与圆

圆是凸的

多边形有$n$条边，圆心为$O(x_0,y_0)$，半径为$r$

多边形在圆内：等价于顶点都在圆内（枚举判断）

复杂度$O(n)$

圆在多边形内：避免不了精度误差了最暴力的求$O$到多边形每条边的端点的距离、用"垂直的一次函数斜率之积为$-1$"求解$O$到该边的垂直距离（如果垂足在线段上的话），求最小值与$r$比较

复杂度$O(n)$

细节很多

圆与圆

圆心分别为$O_1(x_1,y_1)$、$O_2(x_2,y_2)$，半径为$r_1$、$r_2$

圆在圆内：$O_1{O}_2<r_1-r_2$

圆与圆相交：$r_1-r_2<O_1{O}_2<r_1+r_2$

圆与圆相离：$O_1{O}_2>r_1+r_2$

复杂度$O(1)$

特判边界

凸包

定义：形象的将，在木板上有n个钉子，拿根牛皮筋从无限远处开始收缩，问最后成了什么样（可能问周长或面积）

这里只讲单调链算法

对所有点进行$x-y$从小到大排序（以$x$为第一关键字，以$y$为第二关键字）

从最左下点开始维护一个栈，按顺序扫描，当碰到一个新的点时，判断这个点是否在栈顶$2$个点形成的线段的顺时针方向，若是，则弹出栈顶，加入该点，否则不加。

做到最后一个点时，栈中的点就是该凸包的下凸壳，再反向搜一遍求出上凸壳即可。