AcWing算法进阶课 计算几何 基础知识

 

AcWing算法进阶课 计算几何 基础知识

1.前置知识点

（1） pi = acos(-1);

（2）余弦定理：c ^ 2 = a ^ 2+ b ^ 2 - 2abcos(theta)

2.浮点数的比较

const double eps = 1e-8;
int sign(double x)//符号函数，返回一个数的正负或0
{
    if (fabs(x) < eps) return 0;
    if (x < 0) return -1;
    return 1;
}
int cmp(double x, double y) //比较函数，比较两个数的大小
{
    if (fabs(x - y) < eps) return 0;
    if (x < y) return -1;
    return 1;
}

3.向量

3.1向量的加减法和数乘运算

3.2向量（点积）A * B = |A||B|cos(C)

(1).几何意义：向量A在向量B上的投影与B的长度的乘积。

(2).代码实现

double dot(Point a, Point b)
{
    return a.x * b.x + a.y * b.y;
}

3.3外积（叉积）AxB = |A||B|sin(C)

(1).几何意义：向量A与B张成的平行四边形的有向面积。B在A的逆时针方向为正。

(2).代码实现：

double cross(Point a, Point b)
{
    return a.x * b.y - b.x * a.y;
}

3.4常用函数

（1).取模

double get_lenth(Point a)
{
    return sqrt(dot(a, a));
}

(2).计算向量夹角

double get_angle(Point a, Point b)
{
    return acos(dot(a,b) / get_length(a) / get_length(b));
}

(3).计算两个向量构成的平行四边形有向面积

double area(Point a, Point b, Point c)
{
    return cross(b - a, c - a);
}

(4).向量A顺时针旋转C的角度

Point rotate(Point a, double angle)
{
    return Point(a.x * cos(angle) + a.y * sin(angle), -a.x * sin(angle) * cos(angle));
}

4.点与线

4.1 直线定理

（1）一般式： ax + by + c = 0;

（2）点向式：p0 + vt

（3）斜截式：y = kx + b

4.2 常用操作

（1）判断点在直线上 A x B = 0

（2）两直线相交

//cross(v, w) == 0 则两直线平行或重合

//就两条直线的交点
//v和w都是向量 
Point get_line_intersection(Point p, Vector v, Point q, vector w)
{
    vector u = p - q;
    double t = cross(w, u) / cross(v, w);
    return p + v * t;
}

（3）点到直线的距离

//平行四边形面积除以底边长度
double distance_to_line(Point p, Point a, Point b)
{
    vector v1 = b - a, v2 = p - a;
    return fabs(cross(v1, v2)) / get_length(v1);
}

（4）点到线段的距离

double distance_to_segment(Point p, Point a, Point b)
{
    if (a == b) return get_length(p - a);
    Vector v1 = b - a, v2 = p - a, v3 = p - b;
    //点乘小于0，说明夹角大于90度
    if (sign(dot(v1, v2)) < 0) return get_length(v2);
    if (sign(dot(v1, v3)) > 0) return get_length(v3);
    return distance_to_line(p, a, b);
}

（5）点在直线上的投影

double get_line_projection(Point p, Point a, Point b)
{
    Vector v = b - a;
    return a + v * (dot(v, p - a) / dot(v, v));
}

（6）点是否在线段上

bool on_segment(Point p, Point a, Point b)
{
    //共线， 点在线内，点积之积为0 
    return sign(cross(p - a, p - b) == 0 && sign(dot(p - a, p - b))) <= 0;
}

（7）判断两线段是否相交

跨立实验

//判断a1和a2是不是在b1b2线段的两侧，再判断b1和b2是不是在a1a2的两侧
bool segment_intersection(Point a1, Point a2, Point b1, Point b2)
{
    double c1 = cross(a2 - a1, b1 - a1), c2 = cross(a2 - a1, b2 - a1);
    double c3 = cross(b2 - b1, a2 - b1), c4 = cross(b2 - b1, a1 - b1);
    return sign(c1) * sign(c2) <= 0 && sign(c3) * sign(c4) <= 0;
}

5.多边形

5.1 三角形

（1）面积

（1）叉积

（2）海伦公式

p = (a + b + c) / 2;

s = sqrt(p * (p - a) * (p - b) * (p - c))

（2）三角形四心

a.外心，外接圆圆心

三边中垂线交点，到三角形三个顶点的距离相等

b.内心，内切圆圆心

角平分线交点，到三边距离相等

c.垂心

三条垂线交点

d.重心

三条中线交点（到三角形三顶点距离平方和最小的点，三角形内到是三边距离之积最大的点。

5.2 普通多边形

通常按逆时针存储所有点

（1）定义

（1）【简单】多边形

由在同一平面且不在同一直线上的多条线段首尾顺次连接且不相交组成的图形叫多边形

（2）简单多边形

简单多边形是除相邻边外其他边不相交的多边形

（3）凸多边形

过多边形的任意一边多一条直线，如果其他各个顶点都在这条直线的同侧，则把这个多边形叫做凸多边形

任意凸多边形外角和均为360°

任意凸多边形内角和为(n - 2)180°

（2）常用函数

（1）求多边形面积（不一定是凸多边形）

我们可以从第一个顶点出发把凸多边形分成n-2个三角形，然后把面积加起来。

当然更一般的情况是，任意选择一个点，然后将n边形顺次选择每条边加上我们任选的这个点，组成了n-1个三角形，算一下叉积之和

double ploygon_area(Point p[], int n)
{
    double s = 0;
    for (int i = 1; i + 1 < n; i ++ )
        s += cross(p[i] - p[0], p[i + 1] - p[i]);
    //每个叉积算的是平行四边形，是原来面积的两倍
    return s / 2;
}

（2）判断点是否在多边形内（不一定是凸多边形）

a.射线法，从盖蒂埃纳任意做一条和所有边都平行的射线。交点个数为偶数，则在多边形外，为奇数，则在多边形内部

b.转角法

（3）判断点是否在凸多边形内

判断点是否在所有边的左边（逆时针存储多边形）

5.3 皮克定理

皮克定理是指一个计算点阵中定点在格点上的多边形面积公式，该公式可以表示为：

S = a + b / 2 + 1

其中a表示多边形内部的点数，b表示多边形边界上的点数，s表示多边形的面积

6.圆

（1）圆于直线交点

（2）两圆交点

（3）点到圆的切线

（4）两圆公切线

（5）两圆相交面积