发现可以在 x 方向找到一个点 mid ，使得 [1,mid],[mid+1,r] 的点分别用一个圆来覆盖
然后可以二分 + 最小圆覆盖来解决
但这样显然是错的，因为划分可能不沿 x 轴，于是我们可以将坐标轴旋转
考虑 n≤1000，加上二分，数据组数，大概可以旋转个 100 次，每种旋转求一次即可
顺便一提最小圆覆盖：
考虑增量法：
假设当前加入第 i 个点，前 i−1 个点求出的圆是 Ci−1
如果 i 在圆中，那么可以跳过，如果 i 不在圆中，那么 i 一定在圆上
于是我们固定 i 点，枚举 j，假设当前加入 j 点，前 j−1 个点求出的最小圆是 Ci,j−1
如果 j 在圆中那么可以跳过，否则 j 一定在圆上，这个时候我们以 i,j 为直径作圆
这个圆 Ci,j 不一定包涵 [1,j−1] 的点，于是我们继续枚举 k∈[1,j−1]，如果 k 不在圆中那么 k 在圆上，三点定圆，这样求出的圆一定包涵 {i}∪[1,j] 的所有点
考虑复杂度，j 这一步是要枚举满的，但是进入 k 的循环当且仅当 j 在圆上，概率为 3/j
第 i 个点不在圆中的概率是 3/i，这个时候才会进入 j，于是复杂度为 O(n)
#include<bits/stdc++.h>
#define cs const
using namespace std;
cs int N = 1e3 + 5;
cs double eps = 1e-8, PI = acos(-1.0);
cs double angle = PI/100, si=sin(angle), co=cos(angle);
cs double INF = 1e9;
int sgn(double a){ if(fabs(a)<eps) return 0; return a<0?-1:1; }
struct P{
double x,y; P(double _x=0, double _y=0){ x=_x; y=_y; }
friend P operator + (cs P &a, cs P &b){ return P(a.x+b.x,a.y+b.y); }
friend P operator - (cs P &a, cs P &b){ return P(a.x-b.x,a.y-b.y); }
friend double operator * (cs P &a, cs P &b){ return a.x*b.y-a.y*b.x; }
P operator * (cs double &k){ return P(x*k, y*k); }
double norm(){ return x*x+y*y; }
P rot(){ return P(x*co-y*si, x*si+y*co); }
bool operator < (cs P &a)cs{ return x<a.x; }
}p[N],q[N];
double dis(P a,P b){ return sqrt((a-b).norm()); }
struct Cir{
P o; double r; Cir(){}
Cir(cs P &_o, cs double &_r):o(_o),r(_r){}
bool in(P a){ return r+eps>dis(a,o); }
}C;
int n;
Cir getCir(P a,P b,P c){
if(!sgn((a-b)*(b-c))){
if(a.x==b.x){
if((a.y<b.y)==(b.y<c.y)) return Cir((a+c)*0.5,dis(a,c)*0.5);
if((a.y<c.y)==(c.y<b.y)) return Cir((a+b)*0.5,dis(a,b)*0.5);
return Cir((b+c)*0.5,dis(b,c)*0.5);
}
else{
if((a.x<b.x)==(b.x<c.x)) return Cir((a+c)*0.5,dis(a,c)*0.5);
if((a.x<c.x)==(c.x<b.x)) return Cir((a+b)*0.5,dis(a,b)*0.5);
return Cir((b+c)*0.5,dis(b,c)*0.5);
}
}
double A=b.x-a.x,B=b.y-a.y,C=c.x-b.x,D=c.y-b.y,E=b.norm()-a.norm(),F=c.norm()-b.norm();
P o(0.5*(B*F-D*E)/(B*C-A*D),0.5*(A*F-C*E)/(A*D-B*C));
return Cir(o,dis(a,o));
}
double Solve(int l,int r){
if(l>r) return 0; int n=0;
for(int i=l;i<=r;i++) q[++n]=p[i];
random_shuffle(q+1,q+n+1);
C=Cir(q[1],0);
for(int i=1;i<=n;i++) if(!C.in(q[i])){
C.o=q[i]; C.r=0;
for(int j=1; j<i; j++) if(!C.in(q[j])){
C.o=(q[i]+q[j])*0.5;
C.r=dis(C.o,q[i]);
for(int k=1;k<j;k++) if(!C.in(q[k])) C=getCir(q[i],q[j],q[k]);
}
} return C.r;
}
int main(){
srand(time(0));
while(scanf("%d",&n)&&n){
double ans=INF;
for(int i=1; i<=n; i++) scanf("%lf%lf",&p[i].x,&p[i].y);
for(int T=0; T<100; T++){
for(int i=1; i<=n; i++) p[i]=p[i].rot();
sort(p+1, p+n+1);
int l=1, r=n;
while(l<=r){
int mid=(l+r)>>1; double r1=Solve(1,mid),r2=Solve(mid+1,n);
if(min(r1,r2)>=ans) break;
ans = min(ans,max(r1,r2));
if(r1>r2) r=mid-1; else l=mid+1;
}
} printf("%.2lf\n", ans);
} return 0;
}