计算几何基础

zxy巨巨
神仙
ak了winter camp回来给我们讲知识点

开始吧
计算及处理几何问题尽量减小误差
涉及到double，计算少用除法、开根（牛顿迭代法三次除法）

精度：

妙啊
向量表示线段（有向线段）
还可以用美丽的文玉的方法
$对于\forall C\in AB,\exists p\in [0,1],C=pA+(1-p)B$

一个方法

直线方程$Ax+By+C=0$可以用

向量

struct vec{
	double x,y;
};
vec operator+(const vec a,const vec b){
	return (vec){a.x+b.x,a.y+b.y};
}
vec a,b,c;
	a.x=in,a.y=in,b.x=in,b.y=in;
	c=a+b;
	printf("%lf %lf",c.x,c.y);
	return 0;

重载运算符是个好东西
但是你要用就用你重载的东西啊！（定义了+不要有-）
重载后不影响原来对于其他数据类型的用法
函数也可以解决这个问题

点积

详见数学书，必修四欢迎您

叉乘

叉乘结果的大小是两向量构成的平行四边形的有向面积，用行列式算（详见线形代数）（证明显然）
方向利用右手定则判断（你发现你比good的时候就是正的，比bad的时候就是负的）
向量的叉乘有有序性，不满足交换律
$\vec{a}\times \vec{b}=-\vec{b}\times \vec{a}$

（期末考前玩着玩着搞出来了向量叉乘坐标运算的一个严格的证明，我是不是有点不务正业T_T）
向量叉乘坐标运算：
仿照数学书上点乘的坐标运算，我们需要知道向量的叉乘是否遵循分配率
即$(\vec{a}+\vec{b})\times \vec{c}=?=\vec{a}\times \vec{c}+\vec{b}\times \vec{c}$
做出几何图形如下，易证

$(\vec{a}+\vec{b})\times \vec{c}=S_{AC{C}^{\prime }D}$
$\vec{a}\times \vec{c}=S_{AB{B}^{\prime }D}$
$\vec{b}\times \vec{c}=S_{BC{C}^{\prime }{B}^{\prime }}$
由初中几何知识易知：$(\vec{a}+\vec{b})\times \vec{c}=\vec{a}\times \vec{c}+\vec{b}\times \vec{c}$

接着，我们就可以放心大胆地使用分配率了（下面有点不严谨：$i和j$上面应该打向量的箭头，但是咕了）
$\vec{a}=(x_1,y_1)=x_{1}i+y_{1}j,\vec{b}=(x_2,y_2)=x_{2}i+y_{2}j$
$\vec{a}\times \vec{b}=x_1{x}_2(i\times i)+y_1{y}_2(j\times j)+x_1{y}_2(i\times j)+y_1{x}_2(j\times i)$
$=x_1{y}_2-x_2{y}_1$（ixi=0，jxj=0，ixj=1，jxi=-1）

妙不妙？

向量的旋转

$\vec{a}=(x,y)=x\vec{i}+y\vec{j}$
旋转角$\theta$
各位应该很明显的看到直接把$x\vec{i}$和$y\vec{j}$转了不就行了？

叉乘

我是不是又打了一个一样的标题？
那是因为我懒得想这个标题究竟应该写什么

· 点乘和叉乘用

void dottimes(vec a,vec b,int &c){
	c=a.x*b.x+a.y*b.y;
}
void crosstimes(vec a,vec b,int &c){
	c=a.x*b.y-b.x*a.y;
}

给一点叉乘的小结论
· 求有向面积：废话，用叉乘$S_{\Delta OAB}=\overrightarrow{OA}\times \overrightarrow{OB}$

· 又看数学书：三点共线的时候是不是围成三角形的面积为0？是的，可以根据叉乘$x_1{y}_2-x_2{y}_1=0$

· 还可以依据三角形的有向面积判断角的大小（优劣角）——
称$\alpha$为以$\vec{a}$为始边、以$\vec{b}$为终边的角
可以发现：
$\vec{a}\times \vec{b}>0\Rightarrow \alpha \in (0,\pi )$
$\vec{a}\times \vec{b}<0\Rightarrow \alpha \in (\pi ,2\pi )$
一会儿会用到它来做凸包

· 另：点乘的正负由$cos<\vec{a},\vec{b}>$提供，因此可以判断夹角为钝角或锐角

· 求距离：三角形面积与高，他不容易吗？$d=\frac{S_{\Delta PAB}}{|\overrightarrow{AB}|}=\frac{\overrightarrow{OA}\times \overrightarrow{OB}}{2|\overrightarrow{AB}|}$
提醒！叉乘算的是平行四边形的面积！

· 判断是否相交：
端点分布在两侧，利用面积，叉乘的结果同号（注意边的顺序）
还可以用来求交点——知道面积知道面积比知道边上比例

· 点与多边形的关系：在内或在外
做一条直线，记穿边的个数，奇内偶外
如果交了顶点，算远不算近（不要问我，对于顶点肯定不能说交了1个点对吧，那么交0个或者交2个对奇偶性是没有影响的，算不算无所谓）

· 求多边形面积就很妙了

叉乘自动给你改变正负
这就是向量的魅力，所以我第一眼看到向量的时候就爱上他了

· 还有一种判断点是否在多边形的方法
把多边形放在极坐标系里，以一个点作为原点，易知每条对角线的极角是单增的，知道点的极角，二分它所在的区域，再用叉乘判断点在不在里面就行了

【模板】二维凸包 / [USACO5.1]圈奶牛Fencing the Cows

求凸包：
step1 取极点，选左下角的那个
step2 排序，根据极角单增与距离单增（叉乘实现）
step3 挨个连，用栈维护，弹优角
现在栈里面就是凸包的边界了

码的时候遇到了一点问题：
运算符之类的符号是不能对结构体运行
在结构体里面要么重载运算符，要么return true（不改变顺序）| return false（改变顺序）