POJ, 2069 Super Star(模拟退火算法)

题意：给定几个点，要求覆盖这些点的最小球半径。

分析：可以采用模拟淬火算法，随机选取一个点作为初始解，然后不断向当前距离最远点靠近，这是一个不断调整的过程，将对应模拟淬火算法中不断向内能最低（半径最小）这一目标函数逼近，温度对应控制变量。

#include <iostream>
#include<cstdio>
#include<ctime>
#include<cmath>
#include<algorithm>
#include<cstdlib>
using namespace std;
#define N 35
#define INF 999999
#define EPS 1e-8
double a[N],b[N],c[N];
int main()
{
    int m;
    while(~scanf("%d",&m)&m)
    {
        for(int i=0;i<m;i++)
        {
            scanf("%lf%lf%lf",&a[i],&b[i],&c[i]);
        }
        double T=100;
        double  ansx=0;
        double ansy=0;
        double  ansz=0;
            int lm=0;
           double ans=INF;
        while(T>EPS)
        {
              for(int i=0;i<m;++i){
                if((ansx-a[i])*(ansx-a[i])+(ansy-b[i])*(ansy-b[i])+(ansz-c[i])*(ansz-c[i])>((ansx-a[lm])*(ansx-a[lm])+(ansy-b[lm])*(ansy-b[lm])+(ansz-c[lm])*(ansz-c[lm])))
                    lm=i;
            }
            double d=sqrt((ansx-a[lm])*(ansx-a[lm])+(ansy-b[lm])*(ansy-b[lm])+(ansz-c[lm])*(ansz-c[lm]));
            ans=min(ans,d);
            ansx+=(a[lm]-ansx)/d*T;
            ansy+=(b[lm]-ansy)/d*T;
            ansz+=(c[lm]-ansz)/d*T;
            T=T*0.98;
        }
        printf("%.5f\n",ans);

    }
    return 0;
}