【POJ】2420 A Star not a Tree?（模拟退火）

题目

传送门：QWQ

分析

军训完状态不好QwQ，做不动难题，于是就学了下模拟退火。

之前一直以为是个非常nb的东西，主要原因可能是差不多省选前我试着学一下但是根本看不懂？

骗分利器，但据说由于调参困难，很多情况比不上多点爬山？

总体来说模拟退火的精髓就是：

可以看到T越大，转移的概率越高。 exp是在cmath里面有的，直接用就ok。

本题就是找出n个点的费马点，即找出一个点到这n个点的距离和最小。

代码

 // #include <bits/stdc++.h>
 // POJ 2420
 #include <cstdio>
 #include <cmath>
 #include <cstdlib>
 #include <ctime>
 using namespace std;
 const int maxn=;
 const double pace=0.95;
 struct Point { double x,y; };

 int n;
 Point p[maxn];

 double sqr(double x){ return x*x; }
 double dis(double x,double y,Point a){ return sqrt(sqr(x-a.x)+sqr(y-a.y)); }
 double getnum(double x,double y){
     double ans=;
     for(int i=;i<=n;i++) ans+=dis(x,y,p[i]);
     return ans;
 }
 int main(){
   //  srand(time(NULL));
     srand();
     Point ans;
     scanf("%d",&n);
     for(int i=;i<=n;i++) {
         scanf("%lf%lf",&p[i].x,&p[i].y);
         ans.x+=p[i].x; ans.y+=p[i].y;
     }
     ans.x/=n; ans.y/=n;
     double t=1e5, end=0.02, pre=getnum(ans.x,ans.y);
     while(t>end){

         double xx=, yy=;
         t*=pace;
         for(int i=;i<=n;i++){
             xx+=(p[i].x-ans.x)/dis(ans.x,ans.y,p[i]);
             yy+=(p[i].y-ans.y)/dis(ans.x,ans.y,p[i]);
         }
         double tmp=getnum(ans.x+xx*t,ans.y+yy*t);
         if(tmp<pre){
             pre=tmp; ans.x+=xx*t; ans.y+=yy*t;
         }
         else if(exp((pre-tmp)/t)>(rand()%)/10000.0){
             ans.x+=xx*t; ans.y+=yy*t;
         }
         t*=pace;
     }
     printf("%.0f",pre);
     return ;
 }