最小外接球是包含所有点且半径最小的球，球心未必是点集均值，由1~4个边界点唯一确定；需用Welzl算法（随机打乱+递归支撑集）求解，避免暴力枚举和浮点退化问题。
最小外接球是什么，为什么不能直接用平均点当球心
最小外接球（Minimum Bounding Sphere）不是点集的几何中心，也不是坐标均值——std::accumulate算出的平均点很可能离最远点很远。它必须满足：球心到所有点的欧氏距离 ≤ 半径，且半径尽可能小。这个球由最多 4 个点“支撑”（即在球面上），其中 1~4 个点决定球的唯一解。暴力枚举所有四元组太慢（O(n⁴)），实际要用迭代算法。
推荐用 Welzl’s 算法，C++ 实现要点
Welzl’s 是目前最实用的随机增量算法，期望时间复杂度 O(n)，常数小，代码可控。核心是递归函数 welzl，接受点集和“边界点集（support set）”，后者大小不超过 4。关键细节：
- 输入点必须用
std::vector<:array>></:array>或自定义struct Point3D { double x,y,z; };，避免std::vector<:vector>></:vector>带来额外拷贝开销 - 务必先对点做随机打乱：
std::shuffle(points.begin(), points.end(), std::mt19937{std::random_device{}()});，否则最坏情况退化为 O(n²) - 递归终止条件：support set 大小为 0 → 返回空球；大小为 1 → 球心=该点，半径=0；大小为 2 → 球心=中点，半径=距离一半；大小为 3 或 4 → 调用
circumsphere_3d或circumsphere_4points
三点确定球面？不，三维里三点只能确定一个圆，需要第四点或约束
三维空间中，任意不共线三点确定唯一圆（位于某平面内），但无法唯一确定球——有无穷多个球包含该圆。真正能唯一确定最小外接球的，是 2~4 个点，且它们必须处于球表面。实操中：
C++ Code Review Master
下载
组合式C++代码评审方案，融合静态分析、AI推理、多轮迭代评审和C++专项检查，适用于PR审查、增量代码审查、全项目评审和代码质量评分，触发词包括review cpp、cpp代码评审、C++review、代码审查。
- 2 点 → 直径端点，球心 =
(p0 + p1) / 2.0，半径 =distance(p0, p1) / 2.0 - 3 点 → 先求三点所在平面的外接圆圆心（二维投影），再沿法向偏移得到球心；更稳的做法是解方程组：
|c - p0|² = |c - p1|² = |c - p2|²，转为线性系统求解 - 4 点 → 解
|c - p0|² = |c - p1|² = |c - p2|² = |c - p3|²，得唯一球心（前提是四点不共面）；可用 Cramer 法则或 QR 分解，但注意数值稳定性，建议用Eigen::Vector3d配合Eigen::FullPivLU
容易崩的坑：浮点精度、退化点集、重复点
真实数据常含重复点、近似共面/共线点，导致矩阵奇异或除零。应对方式：
立即学习“C++免费学习笔记（深入）”；
- 预处理去重：用
std::set配合自定义比较（加 eps 比较，如fabs(a.x-b.x) ） - 计算距离平方时别开根号，所有比较用
sq_distance，避免 sqrt 误差累积 - 解线性系统前检查行列式绝对值是否
，若是则降维处理（如三点近似共线，退化为两点方案） - 最终验证：遍历所有原始点，确认
sq_distance(center, p) ，否则重新运行或换初始随机种子
Welzl 实现不到 150 行，但边界 case 处理占一半工作量；别省略验证步，尤其当输入含激光扫描噪声点时。
