本文主要是介绍POJ 2069 Super Star（最小球覆盖），希望对大家解决编程问题提供一定的参考价值，需要的开发者们随着小编来一起学习吧！
题目链接：http://poj.org/problem?id=2069
wori，随机退火根本满足不了精度。。。
要每次向半径最大的那个点的方向前进，感觉像是梯度下降啊？？？
然后要过程中记录最小值，如果最后计算的话好像精度还是有问题？？？
再毒瘤一点
代码：
#include<cstdio>
#include<cmath>
#include<algorithm>
#include<time.h>
using namespace std;
namespace SA
{const double EPS=1e-8;const double PI=acos(-1.0);int n;double Delta,ans;struct Point{double x,y,z;}P[55];struct Solution{double x,y,z;}S;double Dis(double x1,double y1,double z1,double x2,double y2,double z2){return sqrt((x1-x2)*(x1-x2)+(y1-y2)*(y1-y2)+(z1-z2)*(z1-z2));}void Trans(){int id=1;double mx=0;for(int k=1;k<=n;k++){double dis=Dis(S.x,S.y,S.z,P[k].x,P[k].y,P[k].z);if(dis>mx) mx=dis,id=k;}S.x=S.x+(P[id].x-S.x)/mx*Delta;S.y=S.y+(P[id].y-S.y)/mx*Delta;S.z=S.z+(P[id].z-S.z)/mx*Delta;ans=min(ans,mx);}void solve(){for(int i=1;i<=n;i++){scanf("%lf%lf%lf",&P[i].x,&P[i].y,&P[i].z);}S.x=S.y=S.z=0;Delta=100;ans=1e10;while(Delta>EPS){Trans();Delta=Delta*0.98;//0.95 0.9}printf("%.5lf\n",ans);}
}
int main()
{while(scanf("%d",&SA::n)==1&&SA::n){SA::solve();}return 0;
}这篇关于POJ 2069 Super Star（最小球覆盖）的文章就介绍到这儿，希望我们推荐的文章对编程师们有所帮助！
