直接用距离公式等于半径判定圆上点几乎总会失败，因浮点舍入误差使理论在圆上的点计算距离不等于r；应改用容差比较，如abs(dist - r)

浮点数比较导致的“看似在圆上却判定失败”

直接用 sqrt((x - cx)*(x - cx) + (y - cy)*(y - cy)) == r 判定几乎总会失败——浮点运算存在舍入误差，哪怕点理论上在圆上，计算出的距离也大概率不等于 r。哪怕用 abs(dist - r) ，仍可能因尺度失配误判（比如半径是 <code>1e6 时，1e-9 远小于机器精度相对误差）。

正确做法是避免开方，改用平方距离与半径平方比较，并引入相对容差：

• 计算平方距离：(x - cx)*(x - cx) + (y - cy)*(y - cy)
• 计算 r * r（注意：若 r 是浮点数，确保不重复计算）
• 判断：abs(sq_dist - r_sq) ，其中 <code>epsilon 取 1e-10 或更小，取决于坐标量级

C++中推荐的健壮判定函数写法

封装成函数时，需处理边界情况（如 r 、无穷坐标），并统一使用 <code>double 避免单精度累积误差：

bool onCircle(double x, double y, double cx, double cy, double r) {
    if (r < 0) return false; // 负半径无意义
    if (!std::isfinite(x) || !std::isfinite(y) || !std::isfinite(cx) || 
        !std::isfinite(cy) || !std::isfinite(r)) return false;
<pre class="brush:php;toolbar:false;">double dx = x - cx;
double dy = y - cy;
double sq_dist = dx*dx + dy*dy;
double r_sq = r * r;

double diff = std::abs(sq_dist - r_sq);
double scale = std::max({sq_dist, r_sq, 1.0});
return diff <= 1e-12 * scale;

}

立即学习“C++免费学习笔记（深入）”；

关键点：

C++14

C++14 对 C++11 的修正与增强版本，适合旧系统维护和较老工具链兼容。

下载

• 不用 sqrt，避免额外误差和性能损耗
• 容差随尺度自适应（1e-12 * scale），比固定阈值更可靠
• 显式检查 std::isfinite，防止 NaN 或 inf 传播

当需要更高精度或整数坐标时的特殊处理

如果所有输入都是整数（如像素坐标、网格点），且半径平方也是整数，可完全避免浮点误差：

• 用 long long 存储平方距离和 r*r，防止溢出（注意：坐标范围大时需检查溢出）
• 直接判断 sq_dist == r_sq，无需容差
• 示例：onCircleInt(3, 4, 0, 0, 5) → (3-0)^2 + (4-0)^2 == 25 → true

但一旦涉及除法、旋转、缩放等操作，坐标大概率转为浮点，此时必须回归带容差的平方比较。

常见误判场景与调试建议

实际调试时，常发现“明明画出来就在圆上，代码却返回 false”。典型原因包括：

• 传入的 r 是直径而非半径（检查变量命名和单位一致性）
• 坐标系 y 轴方向反了（如 OpenGL 与图像坐标系差异），导致 y - cy 符号错
• 使用了 float 类型，中间计算溢出或精度丢失（尤其大坐标时）
• 调试打印用了 std::cout 默认只显示 6 位小数，掩盖了真实误差值

建议打印完整精度：std::cout

真正难的不是公式本身，而是意识到“精确位于周界线”在计算机里从来不是数学意义上的相等，而是一个尺度敏感的近似区间。