poj 2069 Super Star (模拟退火)

poj 2069 题目链接

题意:给出n个点,找出包含这n个点的最小半径的外接球。求球的半径。

思路:模拟退火。不过在走的时候,不是随机上下左右前后6个方向走,而是每次往距离当前球心最远的点的方向走。这样才能通过(随机6个方向的写法样例也是可以通过的)

所以模拟退火的精髓大概是“概率减小” 而不是“随机”?

以及,我看到的资料中,对于模拟退火的介绍,都有一部分是“允许一定概率向着不优的方向移动,但是这个概率会越来越小(因此叫退火)”

但是做了几道所谓的“模拟退火”的题目,发现并没有这部分。。。?

那样的话应该是爬山法?

以及看到说模拟退火适合并行计算。。。

是不是体现在我可以放k个初始解,然后每次最优的答案是在k个答案中取最优?

  1/* ***********************************************
  2Author :111qqz
  3Created Time :2016年08月30日 星期二 14时31分15秒
  4File Name :code/poj/2069.cpp
  5************************************************ */
  6#include <cstdio>
  7#include <cstring>
  8#include <iostream>
  9#include <algorithm>
 10#include <vector>
 11#include <queue>
 12#include <set>
 13#include <map>
 14#include <string>
 15#include <cmath>
 16#include <cstdlib>
 17#include <ctime>
 18#define fst first
 19#define sec second
 20#define lson l,m,rt<<1
 21#define rson m+1,r,rt<<1|1
 22#define ms(a,x) memset(a,x,sizeof(a))
 23typedef long long LL;
 24#define pi pair < int ,int >
 25#define MP make_pair
 26#define INF 1E20
 27#define MAX 1E6
 28using namespace std;
 29const double eps = 1E-7;
 30const int dx6[6]={1,0,0,-1,0,0};
 31const int dy6[6]={0,1,0,0,-1,0};
 32const int dz6[6]={0,0,1,0,0,-1};
 33const int inf = 0x3f3f3f3f;
 34int n;
 35int dblcmp(double d)
 36{
 37    return d<-eps?-1:d>eps;
 38}
 39struct point
 40{
 41    double x,y,z;
 42    void input()
 43    {
 44	scanf("%lf%lf%lf",&x,&y,&z);
 45    }
 46    point(){}
 47    point (double _x,double _y,double _z):x(_x),y(_y),z(_z){}
 48    void look()
 49    {
 50//	cout<<"x:"<<x<<" y:"<<y<<" z:"<<z<<" ";
 51	printf("x:%.5f y:%.5f z:%.5f ",x,y,z);
 52    }
 53    double dis(point b)
 54    {
 55	double res = 0 ;
 56	res = (x-b.x)*(x-b.x)+(y-b.y)*(y-b.y)+(z-b.z)*(z-b.z);
 57//	if (res<0.000001){
 58//	    look();
 59//	    cout<<"res:"<<res<<endl;
 60//	}
 61	return sqrt(res);
 62    }
 63    bool ok ()
 64    {
 65	if (x<0||y<0||z<0||x>100||y>100||z>100) return false;
 66	return true;
 67    }
 68}p[35];
 69double ax,ay,az;
 70double SA(const point *p,int n,point &ret)
 71{
 72    double t,cur,ans = INF;
 73    const double delta = 0.98;
 74    int mp = 0; //记录距离当前圆心最远的点,每次不是随机移动而是朝着最远的点的方向移动。
 75    for (ret = point(ax,ay,az) ,t = 100 ; t>eps ; t*=delta){
 76	for ( int i = 0 ; i < n ; i++)
 77	    if (ret.dis(p[mp])<ret.dis(p[i]))  mp = i ;
 78	cur = ret.dis(p[mp]);
 79	if (dblcmp(cur-ans)<0) ans = cur;
 80	ret.x+=(p[mp].x-ret.x)/cur*t;
 81	ret.y+=(p[mp].y-ret.y)/cur*t;
 82	ret.z+=(p[mp].z-ret.z)/cur*t;
 83    }
 84    return ans;
 85}
 86int main()
 87{
 88	#ifndef  ONLINE_JUDGE 
 89	freopen("code/in.txt","r",stdin);
 90  #endif
 91      while (~scanf("%d",&n))
 92      {
 93	  if (n==0) continue;
 94	  ax = ay = az = 0 ;
 95	  for ( int i = 0 ; i < n ; i++)
 96	    {
 97		p[i].input();
 98		ax+=p[i].x;
 99		ay+=p[i].y;
100		az+=p[i].z;
101	    }
102	  ax = ax/n;
103	  ay = ay/n;
104	  az = az/n;
105	  point ret;
106	  double ans = SA(p,n,ret);
107	  printf("%.5f\n",ans);
108      }
109  #ifndef ONLINE_JUDGE  
110  fclose(stdin);
111  #endif
112    return 0;
113}