【POJ】【2420】A Star not a Tree？

模拟退火

　　Orz HZWER

　　这题的题意是在二维平面内找一点，使得这点到给定的n个点的距离和最小……0.0

　　模拟退火算法请戳这里

 //POJ 2420
 #include<ctime>
 #include<cmath>
 #include<cstdio>
 #include<cstring>
 #include<cstdlib>
 #include<iostream>
 #include<algorithm>
 #define rep(i,n) for(int i=0;i<n;++i)
 #define F(i,j,n) for(int i=j;i<=n;++i)
 #define D(i,j,n) for(int i=j;i>=n;--i)
 using namespace std;
 int getint(){
     int v=,sign=; char ch=getchar();
     while(!isdigit(ch)) {if(ch=='-') sign=-; ch=getchar();}
     while(isdigit(ch))  {v=v*+ch-''; ch=getchar();}
     return v*sign;
 }
 /*******************template********************/
 int n;
 double x,y,ans,t;
 struct point{double x,y;} p[];

 double dis(double x,double y,point p){
     return sqrt( (p.x-x)*(p.x-x)+(p.y-y)*(p.y-y) );
 }
 double getsum(double x,double y){
     double temp=;
     F(i,,n) temp+=dis(x,y,p[i]);
     return temp;
 }
 int main(){
     srand(time());
     while(scanf("%d",&n)!=EOF){
         x=y=; ans=1e20; t=1e5;
         F(i,,n){
             scanf("%lf%lf",&p[i].x,&p[i].y);
             x+=p[i].x; y+=p[i].y;
         }
         x/=n; y/=n;
         ans=getsum(x,y);//初始解
         double tmp,tx,ty;
         while(t>0.02){
             tx=ty=;
             F(i,,n){
                 tx+=(p[i].x-x)/dis(x,y,p[i]);
                 ty+=(p[i].y-y)/dis(x,y,p[i]);
             }
             tmp=getsum(x+tx*t,y+ty*t);
             //生成一个新解tmp
             if (tmp<ans){
                 ans=tmp; x+=tx*t; y+=ty*t;
             }//接受较优解
             else if( log((tmp-ans)/t) < (rand()%)/10000.0 ){
                 ans=tmp; x+=tx*t; y+=ty*t;
             }//一定概率接受较差解 （由于t越来越小，所以接受较差解的概率也越来越小……）
             t*=0.9;//“降温”
         }
         printf("%.0f\n",ans);
     }
     return ;
 }

（带注释）