18 南京 D

裸的最小球覆盖。

坐标范围大一些所以我们把初始的温度也设置的大一些。

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
typedef double db;
const db INF = 1e100;
const db delta = 0.98;
const db eps = 1e-8;
struct p3{
    db x,y,z;
    p3 operator +(const  p3&k1)const { return p3{x+k1.x,y+k1.y,z+k1.z};}
    p3 operator -(const p3&k1)const { return p3{x-k1.x,y-k1.y,z-k1.z};}
    p3 operator *(const db k1)const { return p3{x*k1,y*k1,z*k1};}
    p3 operator /(const db k1)const { return p3{x/k1,y/k1,z/k1};}
    db abs(){ return sqrt(x*x+y*y+z*z);}
    db dis(p3 k1){ return (*this-k1).abs();}
};
int n;
p3 p[105];
int main(){
    scanf("%d",&n);
    for(int i=1;i<=n;i++){
        scanf("%lf%lf%lf",&p[i].x,&p[i].y,&p[i].z);
    }
    db T = 100000;
    db ans = INF;
    p3 x = {0,0,0};
    while (T>eps){
        db mx =0;int id=-1;
        for(int i=1;i<=n;i++){
            if((p[i]-x).abs()>mx){
                mx=(p[i]-x).abs();
                id=i;
            }
        }
        ans = min(ans,mx);
        x=x+(p[id]-x)/mx*T;
        T*=delta;
    }
    printf("%.11f
",ans);
}