Weird Flecks, But OK（最小圆覆盖问题）

题目链接：点击这里

题目大意：
给出三维坐标中的 $n$ 个点，求一个圆柱的最小直径，该圆柱垂直于坐标平面且能覆盖住所有点

题目分析：
可以将这 $n$ 个点投影到三个坐标平面上，在三个坐标面上分别做期望时间复杂度为 $O(n)$ 的最小圆覆盖算法，最后取最小值即可
下面详细讲一下此最小圆覆盖算法（随机化增量算法）求半径，其实现步骤如下：
我们依次遍历所有点，当遍历到第 $i$ 个点时，我们判断一下当前点是否在前 $i-1$ 个点的最小圆之中，若不在我们将最小圆重新定义为以点 $i$ 为圆心，$r=0$ 为半径的圆，再遍历前 $i-1$ 个点，若找到一个点 $j$ 不在当前新的最小圆中我们就以点 $i,j$ 为最小圆的直径更新一下圆，然后再在前 $j-1$ 个点中找不在新的最小圆中的点 $k$ ，若找到了这样的点 $k$ 则用 $i,j,k$ 三点构成的三角形的外接圆的半径更新答案即可

具体细节见代码：

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
#include<vector>
#include<set>
#include<map>
#define ll long long
#define inf 0x3f3f3f3f
#define Inf 0x3f3f3f3f3f3f3f3f
#define int ll
using namespace std;
int read()
{
	int res = 0,flag = 1;
	char ch = getchar();
	while(ch<'0' || ch>'9')
	{
		if(ch == '-') flag = -1;
		ch = getchar();
	}
	while(ch>='0' && ch<='9')
	{
		res = (res<<3)+(res<<1)+(ch^48);//res*10+ch-'0';
		ch = getchar();
	}
	return res*flag;
}
const int maxn = 2e5+5;
const int maxm = 3e4+5;
const int mod = 998244353;
const double pi = acos(-1);
const double eps = 1e-8;
int n;
struct Point{
	double x,y;
}a[maxn],b[maxn],c[maxn];
double get_dis(Point a,Point b)
{
	return hypot(a.x-b.x,a.y-b.y);
}
void updat(Point p1,Point p2,Point p3,Point &o,double &r) //用确定的三点更新外接圆圆心和半径 
{
	double xx1 = 2*(p2.x-p1.x),yy1 = 2*(p2.y-p1.y);
    double c1 = p2.x*p2.x+p2.y*p2.y-p1.x*p1.x-p1.y*p1.y;
    double xx2 = 2*(p3.x-p1.x),yy2 = 2*(p3.y-p1.y);
    double c2 = p3.x*p3.x+p3.y*p3.y-p1.x*p1.x-p1.y*p1.y;
    o.x = (yy1*c2-yy2*c1)/(yy1*xx2-yy2*xx1);
    o.y = (xx2*c1-xx1*c2)/(yy1*xx2-yy2*xx1);
    r=get_dis(o,p1);
}
double get_circle(Point a[]) //求能覆盖a中所有点的最小外接圆直径 
{
	random_shuffle(a+1,a+n+1);
	Point o = a[1];double r = 0;
	for(int i = 2;i <= n;i++)
	  	if(get_dis(o,a[i]) > r+eps)
	  	{
	  		o = a[i],r = 0;
	  		for(int j = 1;j < i;j++)
	  			if(get_dis(o,a[j]) > r+eps)
	  			{
	  				o.x = (a[i].x+a[j].x)/2;
	  				o.y = (a[i].y+a[j].y)/2;
	  				r = get_dis(a[i],a[j])/2;
	  				for(int k = 1;k < j;k++)
	  					if(get_dis(o,a[k]) > r+eps) 
	  						updat(a[i],a[j],a[k],o,r);
				}
	 	}
	return 2*r;
}
signed main()
{
	n = read();
	for(int i = 1;i <= n;i++)
	{
		double x,y,z;
		scanf("%lf%lf%lf",&x,&y,&z);
		a[i] = {x,y};b[i] = {y,z};c[i] = {x,z};
	}
	double ans;
	ans = min({get_circle(a),get_circle(b),get_circle(c)});
	printf("%.10f\n",ans);
	return 0;
}