如何计算可包容一组圆的最小圆及通用算法实现

核心问题定义 #

目标是找到能完全包裹所有给定小圆的最小外接圆，必须满足的硬约束：

• 对每个给定的小圆，包容圆圆心到该小圆圆心的距离 + 该小圆自身半径 ≤ 包容圆半径
• 最终得到的包容圆半径是所有满足上述约束的圆里最小的

你提到的Soddy圆是基于两两相切的小圆生成的相切圆，和本问题的约束场景不同，无法直接套用。

现有代码的问题 #

你当前的实现直接对所有小圆圆心取算术平均作为包容圆圆心、对所有小圆半径取算术平均作为包容圆半径，完全没有匹配上述硬约束，计算结果必然不符合预期：

• 圆心位置没有考虑不同小圆的半径差异
• 半径取值没有结合小圆之间的实际位置距离
• 缺少“小圆是否完全落在包容圆内”的有效性校验

通用计算逻辑 #

这个问题属于带权最小包围圆的经典场景，最小包容圆的参数永远由最多3个边界上的小圆确定，所有完全落在当前圆内部的小圆（比如你提到的可忽略的黄色圆）都不会影响最终结果，和你观察到的特性完全吻合。

不同场景的计算规则 #

• 仅包含1个小圆时，最小包容圆就是这个小圆本身
• 包含2个小圆时，最小包容圆完全由这两个小圆确定：
  1. 半径值为（两小圆圆心的距离 + 两个小圆的半径和）/ 2
  2. 圆心在两个小圆圆心的连线上，位置满足和两个圆心的距离分别适配对应的半径约束
     这对应你左图的等半径小圆场景：当两个小圆半径相等时，圆心为两小圆圆心连线的中点，半径等于两圆心距离的一半加小圆半径，和你手动推导的结果完全一致。
• 包含3个及以上小圆时，逐次校验所有小圆是否在当前已生成的包容圆范围内；遇到落在范围外的小圆时，将其作为新的边界约束重新生成包容圆，直到所有小圆都被完全包裹，收敛得到的就是最终的最小包容圆。

代码修正参考 #

首先实现基础的覆盖关系校验函数：

private static bool IsCovered(Circle outer, Circle inner) {
    float centerDistance = Vector2.Distance(outer.pos, inner.pos);
    // 加入微小浮点容差，避免计算精度问题导致的边界判断错误
    return centerDistance + inner.radius <= outer.radius + 0.0001f;
}

后续迭代过程中，基于边界小圆逐次校验、更新包容圆的位置和半径，替换掉原实现中直接取均值的错误逻辑即可。

说明：你提到的“通过两个已确定的包容圆生成更大包容圆”的逻辑，本质是把每个已生成的包容圆当作一个带半径的点，套用同样的计算流程就能得到覆盖范围更大的最小圆。

内容的提问来源于stack exchange，提问作者Mc Midas