poj 2069 Super Star —— 模拟退火

题目：http://poj.org/problem?id=2069

仍是随机地模拟退火，然而却WA了；

看看网上的题解，都是另一种做法——向距离最远的点靠近；

于是也改成那样，竟然真的A了...感觉这个做法的随机因素好像很小？不过似乎也是一种套路。

代码如下：

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
#include<cstdlib>
#include<ctime>
#define eps 1e-17
#define dt 0.997
using namespace std;
typedef double db;
int const xn=;
int n;
db xx[xn],yy[xn],zz[xn],ans,ansx,ansy,ansz;
int rd()
{
  int ret=,f=; char ch=getchar();
  while(ch<''||ch>''){if(ch=='-')f=; ch=getchar();}
  while(ch>=''&&ch<='')ret=(ret<<)+(ret<<)+ch-'',ch=getchar();
  return f?ret:-ret;
}
db dist(db x,db y,db z,db a,db b,db c){return sqrt((x-a)*(x-a)+(y-b)*(y-b)+(z-c)*(z-c));}
db cal(db x,db y,db z)
{
  db ret=;
  for(int i=;i<=n;i++)ret=max(ret,dist(x,y,z,xx[i],yy[i],zz[i]));
  return ret;
}
void SA()
{
  db T=,x=ansx,y=ansy,z=ansz,lst=ans,tx,ty,tz,tmp;
  while(T>eps)
    {
      int id; db ret=,tt;
      for(int i=;i<=n;i++)
    if((tt=dist(x,y,z,xx[i],yy[i],zz[i]))>ret)ret=tt,id=i;
      ans=min(ans,ret);
      x+=(xx[id]-x)/ret*T;
      y+=(yy[id]-y)/ret*T;
      z+=(zz[id]-z)/ret*T;
      T*=dt;
    }
}
int main()
{
  srand(time());
  while()
    {
      n=rd(); if(!n)return ; db sx=,sy=,sz=;
      for(int i=;i<=n;i++)
    scanf("%lf%lf%lf",&xx[i],&yy[i],&zz[i]),sx+=xx[i],sy+=yy[i],sz+=zz[i];
      ansx=sx/n; ansy=sy/n; ansz=sz/n;
      ans=cal(ansx,ansy,ansz);
      SA();
      printf("%.5lf\n",ans);
    }
  return ;
}