【学习笔记】计算几何
只写我会的。
封装好再用，注意精度。
这里定义角度 \(\pi=180^{\circ}\)，角度的值域为 \([-\pi,\pi]\)。规定从初始边顺时针旋转的角为负角，即度数为负数。
记一条直线 \(PQ\) 所夹矩形为以 \(PQ\) 为对角线的矩形。
其实就是点，线段，直线，多边形，圆的组合。
因为本身自带大常数，一般来说理论复杂度能过就不会超时。
- 正弦定理：令 \(R\) 为三角形外接圆半径，则 \(\frac{a}{\sin A}=2R\)，\(a\) 为一条边，\(A\) 为边 \(a\) 的对角。
- 余弦定理：令 \(a,b,c\) 为三角形三边，\(A\) 为 \(a\) 的对角，则 \(a^2=b^2+c^2-2bc\cos A\)。这给了我们角度与长度的换算。\(\cos A=\frac{b^2+c^2-a^2}{2bc}\)。
极坐标系
类似平面直角坐标系，同样有原点 \(O\) 和单位长度 \(1\)，其中定义一条自极点引出的射线极轴 \(OX\)，就建立了一个极坐标系。
设 \(A\) 为极坐标系上一点，极径 \(|OA|\) 定义为 \(A\) 到原点的距离，记为 \(\rho\)。以极轴为始边 \(OA\) 为终边的角 \(\angle XOA\) 即为极角，记为 \(\varphi\)。
如果限制 \(0\le\varphi\le 2\pi\)，那么除原点外极坐标系中每一个点都可表示为 \((\rho,\varphi)\) 的形式。
与平面直角坐标系的转换
点 \(A(\rho,\varphi)\) 在平面直角坐标系可表示为 \(A(\rho\cos\varphi,\rho\sin\varphi)\)。
点 \(A(x,y)\) 在极坐标系可表示为 \(A(\sqrt{x^2+y^2},\text{atan2}(y,x))\)。其中 atan2(y,x) 是 STL 中自带的函数。
向量
既有大小又有方向的量称为向量。或者称作一个有向线段，即带有起点，长度，方向的线段。一般表示为 \(\vec a\) 或者 \(\vec{AB}\)。也记作向量 \((x,y)\) 表示从 \((0,0)\) 到 \((x,y)\) 的向量。向量的大小也称向量的模，定义为有向线段 \(\vec{AB}\) 的长度。令 \(A(x,y)\)，则 \(|A|=\sqrt{x^2+y^2}\)。
特殊的向量：
- 零向量：模长为 \(0\)，记作 \(\vec 0\) 。方向任意。
- 单位向量：大小为 \(1\) 的向量为该方向的单位向量。
平行向量：方向相同或者相反的两个非零向量。记作 \(a\mid \mid b\)。也叫共性向量。
相等向量：模相等且方向相同的向量。相对地，相反向量指模相等但方向相反的向量。
向量计算：
加减：定义向量 \(a(x_a,y_a),b(x_b,y_b)\)，\(a+b=(x_a+x_b,y_a+y_b),a-b=(x_a-x_b,y_a-y_b)\)。得到的结果仍是向量。
数乘：定义向量 \(a(x_a,y_a)\) 和常数 \(k\) ，\(ka=(kx_a,ky_a)\)。得到的结果仍是向量。
点乘：定义向量 \(a(x_a,y_a),b(x_b,y_b)\)，\(a·b=x_a\times x_b+y_a\times y_b=|a|\times|b|\cos\theta(\theta=\langle a,b\rangle)\)。得到的结果是数值，表示 \(a\) 在 \(b\) 上的投影长度与 \(b\) 的长度乘积。所以计算两向量的夹角 \(\cos\theta=\frac{a\cdot b}{|a|\times|b|}\)。用 acos 得到 \(\theta\)。
叉乘：定义向量 \(a(x_a,y_a),b(x_b,y_b)\)，\(a\times b=x_a\times y_b-x_b\times y_a=|a|\times|b|\sin\theta(\theta=\langle a,b\rangle)\)。得到的结果是数值，表示 \(a\) 和 \(b\) 所夹的平行四边形的有向面积（即可以为负）
移动向量：朝一个方向向量 \(U\) 移动 \(k\) 个单位长度。算出单位向量的角度（与 \(x\) 正半轴即 \((0,1)\)）\(\theta\)，\(x\leftarrow x+k\cos\theta,y\leftarrow y+k\sin\theta\)。
有向量 \(A,B\)，则向量 \(\vec{AB}=B-A\)。
叉乘没有交换律，但是 \(a\times b=-(b\times a)\)。
角度计算中，\(a\) 和 \(b\) 的夹角为负，当且仅当 \(a\) 可以通过顺时针旋转 \(\le 180^{\circ}\) 得到 \(b\)。
旋转向量
对于向量 \((x,y)\)，其绕原点顺时 针旋转 \(\theta\) 后得到的向量为 \((x\cos\theta+y\sin\theta,y\cos\theta-x\sin\theta)\)。会掉非常大的精度，慎用。
存储（表示）
存储一个点（向量）：直接存储坐标。
存储一条直线：以直线上任意一点 \(P\) 和一个方向向量 \(U\) 表示。直线上任意一个点可以表示为 \(P+kU\)，其中 \(k\) 为常数。
判断问题：
两向量垂直：
\(a\cdot b=0\) 或 \(\cos\theta=0\)。\(O(1)\)。
两向量平行或相同：
\(A\times B=0\)。
折线段拐向：
对于点 \(A(x_a,y_a),B(x_b,y_b),C(x_c,y_c)\)，判断 \(A\to B\to C\) 的拐弯方向。
记 \(x=(B-A)\times(C-A)\)，若 \(x=0\) 则三点共线，若 \(x>0\) 则往左拐，若 \(x<0\) 则往右拐。\(O(1)\)。
点是否在线段上：
对于点 \(Q\) 和线段 \(P_1(x_1,y_1)P_2(x_2,y_2)\)。先判断 \(Q,P_1,P_2\) 是否三点共线，然后判断 \(Q\) 是否在 \(P_1P_2\) 的范围内，即 \(x_1\le x_q\le x_2,y_1\le y_q\le y_2\)。\(O(1)\)。
两线段是否相交：
对于线段 \(P_1(x_1,y_1)P_2(x_2,y_2)\) 和 \(Q_1(x_3,y_3)Q_2(x_4,y_4)\)。先判断 \(P_1P_2\) 所夹矩形是否和 \(Q_1Q_2\) 所夹矩形相交，若不相交则线段必定不相交。否则进行跨立实验。对于跨立实验，若 \([(P_1-P_2)\times(Q_1-P_2)][(P_1-P_2)\times(Q_2-P_2)]\le0\) 则称 \(P_1P_2\) 跨立 \(Q_1Q_2\)。判断 \(P_1P_2,Q_1Q_2\) 是否互相跨立。可以自行画图，原理就是看另一条线段的两个端点是否在这条线段的左右两侧或者位于这个线段上。
跨立实验可以解决很多线段直线之类的相交问题，这里不多写。\(O(1)\)。
判断点是否在圆内
判断点到圆心的距离是否小于等于半径即可。
判断线段、折线、多边形是否在凸多边形内
判断所有端点是否都在该多边形内。\(O(nm)\)，\(n\) 为待判断图形点数，\(m\) 为凸多边形边数。特殊地，对于圆 \(m=1\)。
判断圆是否在多边形内
两个条件，圆心在多边形内，圆心到多边形的最短距离大于等于半径，如何计算在后面有。\(O(n)\)，\(n\) 是边数。
判断点是否在多边形内
先判断该点是否为端点或者在边上面。如果可以就能结束了。令 \(P(p_x,p_y)\)。
引一条向左的射线出来，实现可以用一个保证不在多边形内的点 \(Q(-Inf,p_y)\) 变为线段，然后对每条边做跨立，记录与多边形边的相交次数 \(t\)。如果 \(t\) 是奇数代表在多边形内。然后需要剔除一些特殊情况，如不考虑多边形的水平边（判断端点纵坐标是否相等）以及如果 \(PQ\) 与一条边的交点为端点（因为 \(PQ\) 水平所以很好判）那么仅当这个端点为这条边中纵坐标更大的点才算做统计。\(O(n)\)，\(n\) 为边数。
判断线段是否在多边形内
有 \(O(n)\) 做法，但是不会。红温了等会再补。
判断折线、多边形是否在多边形内
对于每条边判断即可。
判断圆是否在圆内
对于圆 \(C_1(P_1,r_1),C_2(P_2,r_2)\)，判断 \(C_2\) 是否在 \(C_1\) 内。如果 \(r_2>r_1\) 显然不可能。否则若 \(dis(P_1,P_2)\le r_1-r_2\) 则满足。注意如果取等代表两圆相切。
计算问题
三角形重心公式
\(\Delta ABC\) 的重心为 \((\frac{x_A+x_B+x_C}{3},\frac{y_A+y_B+y_C}{3})\)。挺好记的。
求点到直线的距离
在直线上任取两点 \(A,B\)，计算向量 \(X=P-A,Y=B-A\)，则距离 \(\Delta=\frac{abs(X\times Y)}{|Y|}\)。相当于平行四边形的面积除以底得到高。
求点到线段的距离
特判一下取端点问题，可以用一个类似跨立的东西判断是否取端点最优，如果不是，那么问题转化为求点到直线的距离。
求点到直线的垂足
求距离，然后随便推一推得到与向量 \((x,y)\) 垂直的向量是 \((y,-x)\)，朝这个方向移动距离个单位得到垂足。
求两直线交点
先判断两直线是否平行或相等。如果平行，没有交点；在平行的基础上过同一点即为重合，有无数个交点。
令两条直线分别为 \((P_1,U_1),(P_2,U_2)\)。设交点坐标为 \(P_1+kU_1\)，有 \((P_1+kU_1-P_2)\times U_2=0\)，所以交点 \(O=P_1+\frac{U_2\times(P_2-P_1)}{U_2\times U_1}U_1\)。
求直线与圆的交点
先判断不交，若圆心到直线距离大于半径则不交。算出垂足坐标，若圆心到直线距离等于半径则只有一个交点即为垂足，否则根据勾股定理，以半径为斜边，垂径为一条直角边算出偏移量，在直线上移动偏移量得到两个交点。
求三角形面积
叉积是平行四边形的面积，我们只要求任意两个边的向量的叉积除以二就是三角形面积。
求多边形面积
\(S=\sum_{i=1}^n P_i\times P_{(i\bmod n)+1]}\)。可以画图辅助理解。
pick 定理
对于一个顶点在格点上的多边形，令面积为 \(S\)，\(a\) 为多边形内部的点，\(b\) 为多边形边上的点，有 \(S=a+\frac{b}{2}-1\)。一般是知道面积回推用。
距离
- 曼哈顿距离：\(dis(A,B)=|x_A-x_B|+|y_A-y_B|\)。
- 欧几里得距离：\(dis(A,B)=\sqrt{(x_A-x_B)^2+(y_A-y_B)^2}\)。
- 切比雪夫距离：\(dis(A,B)=\max(|x_A-x_B|,|y_A-y_B|)\)。
曼哈顿距离和切比雪夫距离可以互相转化。偷懒不写。
凸包
一个包含给定点集的最小凸多边形。首先显然所有端点都被包含在点集之中。一个凸包由上下两凸壳构成。
凸壳就是把凸包切开，然后一个凸壳内的斜率单调不增（上凸壳）或单调不减（下凸壳）。
静态求凸包
不算排序的话是 \(O(n)\) 的。
以 \(x\) 为第一关键字，\(y\) 为第二关键字排序，用单调栈从前往后扫一遍维护出下凸壳，然后从后往前扫一遍维护出上凸壳。注意维护上凸壳的时候不要破坏了下凸壳，具体地可以把栈顶设置成下凸壳的末尾，这样就不会访问到下凸壳部分。注意在栈内元素个数小于两个时新加入的点一定能够加入凸壳。
关于维护凸壳，我们发现凸壳的斜率是单调的，所以我们维护斜率或者说这个向量和上个向量的拐向是否符合现在维护的凸壳性质，如果不是就弹出。
带增点凸包
set 维护凸包内的点，新增点的时候判断一下在不在凸包里面，然后看能不能更新凸包，往前往后扫一下踢掉不满足凸性的点。每个点只会被加一次删一次所以还是 \(O(n\log n)\)。还是注意上下凸包的影响要开两个 set 分开维护。
带增点删点离线凸包
学长说在线不是很可做。
cdq 分治，但是不会。咕咕咕。
平面最近点对
人类的智慧是无敌的
P7883 平面最近点对（加强加强版）里最高赞题解写道：
我们充分发扬人类智慧：
将所有点全部绕原点旋转同一个角度，然后按 \(x\times y\) 排序
根据数学直觉，在随机旋转后，答案中的两个点在数组中肯定不会离得太远
所以我们只取每个点向前的 \(50\) 个点来计算答案
这样速度快得飞起，在 \(n=400000\) 时都可以在 124ms 内卡过
当然，这是对的。
分治解法
复杂度 \(O(n\log n)\)。常数大。
对于一张图，按照 \(x\) 坐标排序，我们分治求解。对于一次 \(solve(l,r)\)，我们解决范围在 \([l,r]\) 区间内的最近点对距离。
容易发现 \(l=r\) 时没有答案，\(l=r-1\) 时答案可以直接出来。首先答案应该现在左右最小值，令这个值为 \(d\)。考虑 \([l,mid]\) 和 \([mid+1,r]\) 的合并。考虑直线 \(x=X_{mid}\)，首先 \(x\) 坐标在 \([x-d,x+d]\) 外的可以直接排除了，因为它们与另一边的距离一定大于 \(d\)。然后考虑对于剩下的点，每个点考虑它纵坐标往上 \(d\) 的范围围成的矩形，显然如果纵坐标相差大于 \(d\) 也可以直接排除。对于这个范围内的所有其他点暴力更新 \(d\)。我们考虑按照 \(y\) 坐标归并排序，这样每次考虑纵坐标往上的 \(d\) 时可以通过双指针做到 \(O(n)\)。
正确性显然。考虑复杂度分析。对于每个点可行范围暴力更新的复杂度为 \(O(m)\)，\(m\) 为范围内可行点数量，那么答案就是 \(O(nm\log n)\)。可以证明 \(m\le5\)，是常数级别的。不证哦，可以自己画图感性理解。
旋转卡壳
求凸包的直径。凸包的直径也可以理解成平面最远点对的距离。可以证明，对于凸包的每一条边，其与凸包点的关系是单峰的，也就是说我们可以用双指针对于每条边求出与其距离最远的点然后以这个点向两端点更新答案。
注意此时最好不要把共线的点放到凸包里面，并且特判凸包内只有两个点的情况。
半平面交
咕咕咕。
闵可夫斯基和
咕咕咕。
