题意:给定几个点,要求覆盖这些点的最小球半径。
分析:可以采用模拟淬火算法,随机选取一个点作为初始解,然后不断向当前距离最远点靠近,这是一个不断调整的过程,将对应模拟淬火算法中不断向内能最低(半径最小)这一目标函数逼近,温度对应控制变量。
#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;
}