[POJ][2069][Super Star][模拟退火]

  
    

   
     /*
题目：Super Star
题目来源：POJ 2069
题目难度：中等偏难
题目内容或思路：
 模拟退火
 一开始一直徘徊在WA和TLE之间，提交了有四十多次，步长和点数改多了就TLE,
 改少了就WA。无奈之下只得网上看别人的代码。发现别人代码中的点不是随机
 移动的，而是朝最远的点移动，这时有些明白了，朝最远的点移动应该会更好
 一些，但是还是WA。再仔细看了看，发现我是先判移动后的地方是不是更优，
 再移动，如果不是更优就不会移动，这样最远的点不变，那个点又不动，这就
 造成了很多次无效的循环，这样反而不好。于是去掉了判断更优的条件。
 
 看来模拟退火也不是那么容易的，其中的技巧还是要好好参透一翻啊。
做题日期：2011.3.8
*/
#include <cstdio>
#include <cstdlib>
#include <climits>
#include <iostream>
#include <algorithm>
#include <cstring>
#include <string>
#include <queue>
#include <map>
#include <vector>
#include <bitset>
#include <cmath>
#include <set>
#include <utility>
#include <ctime>
#define sqr(x) ((x)*(x)) 
using namespace std;

const int N = 33, nt = 1, L = 30;
const double inf = 1e100;
const double eps = 1e-8;
const double pi = acos(-1.0);
int n;
double delta, maxv, minv;

struct vpoint{
 double x, y, z, d;
 int id;
}vp[N], test[N];

double dis(vpoint a, vpoint b) {
 return sqrt(sqr(a.x - b.x) + sqr(a.y - b.y) + sqr(a.z - b.z));
}

void calc(vpoint &p) {
 p.d = 0;
 for (int i = 0; i < n; ++i) {
 double t = dis(p, vp[i]);
 if (p.d < t) {
 p.d = t;
 p.id = i;
 }
 }
}

void updata(int id) {
 double d = test[id].d;
 test[id].x += (vp[test[id].id].x - test[id].x) / d * delta;
 test[id].y += (vp[test[id].id].y - test[id].y) / d * delta;
 test[id].z += (vp[test[id].id].z - test[id].z) / d * delta;
 calc(test[id]);
}

void solve() {
 maxv = -inf; minv = inf;
 for (int i = 0; i < n; ++i) {
 scanf("%lf%lf%lf", &vp[i].x, &vp[i].y, &vp[i].z);
 maxv = max(maxv, vp[i].x);
 maxv = max(maxv, vp[i].y);
 minv = min(minv, vp[i].x);
 minv = min(minv, vp[i].y);
 }
 for (int i = 0; i < nt; ++i) {
 test[i].x = minv + (rand() % (int(maxv - minv + 1) * 1000)) / 1000.0;
 test[i].y = minv + (rand() % (int(maxv - minv + 1) * 1000)) / 1000.0;
 test[i].z = minv + (rand() % (int(maxv - minv + 1) * 1000)) / 1000.0;
 calc(test[i]);
 }
 double r = 0.98;
 for (delta = 100; delta > eps; delta *= r) {
 for (int i = 0; i < nt; ++i)
 for (int j = 0; j < L; ++j)
 updata(i);
 }
 double res = inf;
 for (int i = 0; i < nt; ++i)
 res = min(res, test[i].d);
 printf("%lf\n", res);
}

int main() {
#ifndef ONLINE_JUDGE
 freopen("D:\\in.txt", "r", stdin);
#endif
 while (scanf("%d", &n), n) {
 solve();
 }
 return 0;
}