POJ2420 A Star not a Tree?【模拟退火】

题目描述：

在平面直角坐标系中求一个点使得它到$n$个已知点的距离之和最小。
$n\le 100$

题目分析：

模拟退火模板，求费马点。
模拟退火的学习文章
这题网上大多数的解法是往八个方向尝试然后缩短步长。
也可以采用向量的方式，相互抵消最后哪个剩的多就往哪边走。

Code1：

#include<cstdio>
#include<cmath>
#include<algorithm>
#define maxn 105
using namespace std;
int n;
struct Point{double x,y;}p[maxn];
double x,y,ans,Ans;
inline double calc(double x,double y){
	double s=0;
	for(int i=1;i<=n;i++) s+=sqrt((x-p[i].x)*(x-p[i].x)+(y-p[i].y)*(y-p[i].y));
	return s;
}
int main()
{
	scanf("%d",&n);
	for(int i=1;i<=n;i++) scanf("%lf%lf",&p[i].x,&p[i].y),x+=p[i].x,y+=p[i].y;
	x/=n,y/=n,Ans=ans=calc(x,y);
	for(double T=10000;T>1e-4;T*=0.95)
		for(int i=-1;i<=1;i++)
			for(int j=-1;j<=1;j++){
				double xx=x+i*T,yy=y+j*T,tmp=calc(xx,yy);
				if(tmp<Ans) Ans=tmp;
				if(tmp<ans||exp((ans-tmp)/T)>1.*rand()/RAND_MAX) ans=tmp,x=xx,y=yy;
				//结合此题实际情况，exp中乘10更优秀。
			}
	//printf("%.0f\n",ans);
	printf("%.0f\n",Ans);
}

Code2：

#include<cstdio>
#include<cmath>
#include<algorithm>
#define maxn 105
using namespace std;
int n;
struct Point{double x,y;}p[maxn];
double x,y,ans,Ans;
inline double dis(double x,double y,int i){return sqrt((x-p[i].x)*(x-p[i].x)+(y-p[i].y)*(y-p[i].y));}
inline double calc(double x,double y){double s=0;for(int i=1;i<=n;i++) s+=dis(x,y,i);return s;}
int main()
{
	scanf("%d",&n);
	for(int i=1;i<=n;i++) scanf("%lf%lf",&p[i].x,&p[i].y),x+=p[i].x,y+=p[i].y;
	x/=n,y/=n,Ans=ans=calc(x,y);
	for(double T=10000;T>1e-2;T*=0.97){
		double xx=x,yy=y;
		for(int i=1;i<=n;i++) xx+=(p[i].x-x)/dis(x,y,i),yy+=(p[i].y-y)/dis(x,y,i);
		double tmp=calc(xx,yy);
		if(tmp<Ans) Ans=tmp;
		if(tmp<ans||exp((ans-tmp)/T)>1.0*rand()/RAND_MAX) ans=tmp,x=xx,y=yy;
	}
	//printf("%.4f %.4f\n%.0f\n",x,y,ans);
	printf("%.0f\n",Ans);
}