最小球覆盖

题意：给n个点的坐标(xi,yi,zi)(xi,yi,zi)，求覆盖这n个点的最小球的半径rr。
分析：最小球覆盖模板题，模拟退火或者三分套三分再套三分，模拟退火法就是纯贪心的修改版。在搜索到某个局部最优解之后，普通的贪心算法就结束搜索了，但是模拟退火法会以一定的概率进行下一次搜索，这个概率随着时间的变化而减小。选择任意一个点作为初始解，然后将所有点按到当前圆/球心的距离排序，向距离最远的点靠近。每次靠近的距离都会减小（相当于模拟退火法的概率逐渐减小），最后圆/球心原来越靠近最优解，半径趋近稳定。

附上AC的模拟退火解法，要适当设置rate,太大则超时，太小则不精确

#include<cstdio>
#include<cmath>
#include<algorithm>
using namespace std;
const double eps=1e-3;
int n;
struct Point{
    double x;
    double y;
    double z;
    Point(double _x=0.0,double _y=0.0,double _z=0.0){x=_x;y=_y;z=_z;}
};
Point points[105];
inline double dis(Point a,Point b)
{
    double norm=(a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y)+(a.z-b.z)*(a.z-b.z);
    return sqrt(norm);
}
double solve()
{
    double T=1000.0; //初始温度
    double rate=0.99999; //温度下降系数
    Point ans_p;    //初始圆心
    int cur;
    Point max_p; //距离圆心最远的点
    double ans=1e99;
    while(T>eps) //模拟降温
    {
        double max_dis=0.0; //选定圆心到最远点距离
        for(int i=1;i<=n;++i)
        { //最远的点即points[cur]
            if(dis(ans_p,points[i])>max_dis){
                max_dis=dis(ans_p,points[i]);
                cur=i;
            }
        }
        ans=min(ans,max_dis);
        ans_p.x+=(points[cur].x-ans_p.x)*(T/1000.0);
        ans_p.y+=(points[cur].y-ans_p.y)*(T/1000.0);
        ans_p.z+=(points[cur].z-ans_p.z)*(T/1000.0);
        T*=rate;
    }
    return ans;
}
int main()
{
    while(scanf("%d",&n)!=EOF)
    {
        for(int i=1;i<=n;++i)
        {
            scanf("%lf%lf%lf",&points[i].x,&points[i].y,&points[i].z);
        }
        printf("%.15f\n",solve());
    }
    return 0;
}