计算几何小技巧
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] $
a
t
a
n
2
(
y
,
x
)
:
atan2(y,x):
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 ) S=(x_1y_2-x_1y_3+x_2y_3-x_2y_1+x_3y_1-x_2y_2) S=(x1y2−x1y3+x2y3−x2y1+x3y1−x2y2)
3.2 多边形面积公式
S = ∣ 1 2 ∑ i = 0 n − 2 ( x i y i + 1 − x i + 1 y i ) + 1 2 ( x n − 1 y 0 − x 0 y n − 1 ) ∣ 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| S=∣∣∣∣∣21i=0∑n−2(xiyi+1−xi+1yi)+21(xn−1y0−x0yn−1)∣∣∣∣∣
4.体积公式
4.1 棱锥体积公式
V = 1 3 S h V=\frac{1}{3}Sh V=31Sh
正n棱锥,边长为 l n l_n ln:
[外链图片转存失败,源站可能有防盗链机制,建议将图片保存下来直接上传(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)));
}