点P在圆内(含边界)当且仅当distance(P,O)≤r;用std::hypot计算欧氏距离可避免溢出,或直接比较距离平方以提升性能。
std::hypot 或距离平方比较判断点在圆内核心就一条:点 P 在圆心为 O、半径为 r 的圆内(含边界),当且仅当 distance(P, O) 。但直接开方有性能和精度风险,实际应优先用距离平方比较。
(px - ox) * (px - ox) + (py - oy) * (py - oy) —— 无浮点误差、无函数调用开销、支持整数坐标
std::hypot(px - ox, py - oy) ,但注意 std::hypot 在 C++11 起才保证不溢出,且比平方比较慢
r)是否算“内”,取决于业务逻辑,通常包含(即用 )
当 px, oy, r 是 float 或 double 时,直接比较 distance 可能因舍入误差失败,尤其在临界点附近。
std::sqrt((px-ox)*(px-ox) + (py-oy)*(py-oy)) —— 开方引入额外误差,且效率低
std::hypot(px-ox, py-oy)
1e-9:对 double 坐标且 r ~ 1e6 时,1e-9 远小于机器精度,应改用相对容差,例如 1e-12 * std::max({std::abs(r), std::abs(px), std::abs(ox), ...})
避免重复写距离平方逻辑,建议封装为带类型推导的函数,兼容 int、float、double 甚至自定义向量类。
templatebool is_point_in_circle(T px, T py, T cx, T cy, T r) { static_assert(std::is_arithmetic_v , "T must be arithmetic"); auto dx = px - cx; auto dy = py - cy; return dx * dx + dy * dy <= r * r; }
bool,语义清晰;使用 包含圆周
,无运行时开销,编译期全内联struct Point { double x, y; };),可重载为 is_point_in_circle(Point p
, Point c, double r)
算法本身简单,但工程中出错多因数据来源混乱。
py - cy 符号错,结果完全相反几何计算真正难的从来不是公式,而是确认所有坐标在同一参考系、同一单位、同一手性约定。写完 is_point_in_circle 后,务必用已知的边界点(圆心、圆上点、圆外点)手动验算三组输入。