模拟退火算法

http://poj.org/problem?id=2420
题目大意：
给出平面上n个点，求到所有点距离和最短的点，输出该距离和。

#include <stdio.h>
#include <stdlib.h>
#include <math.h>
#include <string.h>
#include <time.h>
struct pos{double x,y;}p[101];
int n;
///////// 计算点x到各点距离之和 /////////////
double dis(pos x) {
    double d=0;
    for(int i=0;i<n;i++)
        d+=hypot(p[i].x-x.x,p[i].y-x.y);
    return d;
}
///////// 模拟退火 /////////////
double SA() {
    double t,mi,tmp;// 温度，最小距离，临时距离
    pos p,tp;       // 最优点，临时点
    int k=0;        // 记录连续未接受点数
    p.x=p.y=5000;   // 初始化最优点为中心点
    mi=dis(p);
    srand(1);
    ///// 初始化温度，温度以一定速率减小，温度小于临界（1） 或 连续超过一定数量（20）个临时点未被接受时 终止 /////
    for(t=5000;t>1&&k++<20;t*=0.93){
        ///// 在p周围一定 范围 内生成临时点， 该范围随温度t降低逐渐减小 /////
        tp.x=p.x+(rand()&1?1:-1)*rand()*t/RAND_MAX;
        tp.y=p.y+(rand()&1?1:-1)*rand()*t/RAND_MAX;
        tmp=dis(tp);
        ///// 若生成的临时点距离小于当前最小，更新最优点 /////
        ///// 否则，以一定 概率 接受该临时点，该概率随温度t降低逐渐减小 /////
        if(tmp<mi||rand()*10000.0/RAND_MAX<t){p=tp;mi=tmp;k=0;}
    }return mi;
}
int main(){
    scanf("%d",&n);
    for(int i=0;i<n;i++)
        scanf("%lf%lf",&p[i].x,&p[i].y);
    printf("%.0f\n",SA());
    return 0;
}

题意：在三维空间内给出一些行星的坐标，求出把所有点都包含在内的球的最小半径。

#include<iostream>  
#include<cstdio>  
#include<algorithm>  
#include<cstring>  
#include<cmath>  
using namespace std;  
const double eps=1e-7;  
struct node {double x,y,z;}A[35],p;  
int n;  
double dis(node a,node b)  {  
    return sqrt((a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y)+(a.z-b.z)*(a.z-b.z));  
}  
double solve()  {  
    double step=100,ans=1e30,mt;  
    p=(node){50,50,50};
    int s=0,k=0;  
    for(double t=50;t>eps;t*=0.98){  
        for(int i=0;i<n;i++)  
            if(dis(p,A[s])<dis(p,A[i]))s=i;  
        mt=dis(p,A[s]);  
        ans=min(ans,mt);
        //以一定的比例向它移动 
        p.x+=(A[s].x-p.x)/mt*t;  
        p.y+=(A[s].y-p.y)/mt*t;  
        p.z+=(A[s].z-p.z)/mt*t;  
    }printf("%.5f\n",ans);
}  
int main(){  
    while(~scanf("%d",&n),n){  
        for(int i=0;i<n;i++)  
            scanf("%lf%lf%lf",&A[i].x,&A[i].y,&A[i].z);  
        solve();
    }return 0;  
}