计算几何小技巧

计算几何小技巧

1.类型转换时的精度问题

1.1 加上eps

中间结果中很有可能出现类似0.999999999999这样的数，转换成int之后就会变成0，此时可以将其加上eps再转换即可

2. 三角函数

2.1 值域

$acos：[0,pi] $

$asin： [-pi/2,+pi/2] $

$atan： [-pi/2,+pi/2] $

$atan2(y,x):$返回 y/x 的反正切值，以弧度表示，取值范围为(-π,π]。tan(θ) = y/x，θ = atan2(y, x)。注意，x 和 y 同时为 0 时将导致域错误（domain error），因为此时的角度是不存在的，或者说是没有意义的。

当 (x, y) 在象限中时：

当 (x, y) 在第一象限，0 < θ < π/2

当 (x, y) 在第二象限，π/2 < θ < π

当 (x, y) 在第三象限，-π < θ < -π/2

当 (x, y) 在第四象限，-π/2 < θ < 0

当 (x, y) 在象限的边界（也就是坐标轴）上时：

当 y 是 0，且 x 为非负值，θ = 0

当 y 是 0，且 x 是负值，θ = π

当 y 是正值，且 x 是 0，θ = π/2

当 y 是负值，且 x 是 0，θ = -π/2

由此可知，一般情况下用 atan() 即可，当对所求出角度的象限有特殊要求时，应使用 atan2()。

3.面积公式

3.1 三角形面积公式

$S=(x_1{y}_2-x_1{y}_3+x_2{y}_3-x_2{y}_1+x_3{y}_1-x_2{y}_2)$

3.2 多边形面积公式

$$S=\left|\frac{1}{2}\sum _{i=0}^{n-2}\left(x_i{y}_{i+1}-x_{i+1}{y}_i\right)+\frac{1}{2}\left(x_{n-1}{y}_0-x_0{y}_{n-1}\right)\right|$$

4.体积公式

4.1 棱锥体积公式

$V=\frac{1}{3}Sh$

正n棱锥，边长为$l_n$:

[外链图片转存失败,源站可能有防盗链机制,建议将图片保存下来直接上传(img-o7IF54j8-1625232805646)(https://i.loli.net/2021/04/19/iBS1jayGROFxTPk.png)]

double f(double l, double n){
    return l * l * l * n / (12 * tan(PI / n)) * sqrt(1 - 1 / (4 * sin(PI / n) * sin(PI / n)));
}