最小圆覆盖

易证，最后 最小的圆一定有三个点在圆的边界上（只有两个点的时候会有两个点）

所以我们把点rand_shuffle随机排一下

那么现在这些点就都不确定了

然后按照最小圆覆盖的算法，下一个点在前i个点中是作为边界的概率只有3/i

而一旦概率高了过后，那么能够一次覆盖的点又会变多

所以期望的复杂度其实只是O(nk)的

#include<bits/stdc++.h>
using namespace std;
const double eps=1e-12;
struct point
{
	double x,y;
}a[500005],o;
int n;
double r;
double dist(point x,point y)
{
	return sqrt((x.x-y.x)*(x.x-y.x)+(x.y-y.y)*(x.y-y.y));
}
void circle(point x,point y,point z)
{
    double a,b,c,d,e,f;
    a=x.x-y.x;
    b=x.y-y.y;
    c=x.x-z.x;
    d=x.y-z.y;
    e=(x.x*x.x-y.x*y.x+x.y*x.y-y.y*y.y)/2;
    f=(x.x*x.x-z.x*z.x+x.y*x.y-z.y*z.y)/2;
    o.x=(b*f-d*e)/(b*c-a*d);
    o.y=(a*f-c*e)/(a*d-b*c);
    r=dist(o,x);
}
int main()
{
	cin>>n;
	for(int i=1;i<=n;i++)
	{
		scanf("%lf%lf",&a[i].x,&a[i].y);
	}
	random_shuffle(a+1,a+n+1);
	o=a[1],r=0;
	for(int i=2;i<=n;i++)
	{
		if(dist(a[i],o)>r+eps)
		{
			o=a[i],r=0;
			for(int j=1;j<i;j++)
			{
				if(dist(o,a[j])>r+eps)
				{
					o.x=(a[i].x+a[j].x)/2;
					o.y=(a[i].y+a[j].y)/2;
					r=dist(o,a[i]);
					for(int k=1;k<j;k++)
					{
						if(dist(a[k],o)>r+eps)
						{
							circle(a[i],a[j],a[k]);
						}
					}
				}
			}
		}
	}
	printf("%.10lf\n%.10lf %.10lf",r,o.x,o.y);
	return 0;
}