【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=2000;
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=0;
    for(int i=1;i<=n;i++) ans+=dis(x,y,p[i]);
    return ans;
}
int main(){
  //  srand(time(NULL));
    srand(0);
    Point ans;
    scanf("%d",&n);
    for(int i=1;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=0, yy=0;
        t*=pace;
        for(int i=1;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)/10000.0){
            ans.x+=xx*t; ans.y+=yy*t;
        }
        t*=pace;
    }
    printf("%.0f",pre);
    return 0;
}