题目描述:
在平面直角坐标系中求一个点使得它到
n
n
n个已知点的距离之和最小。
n
≤
100
n\le100
n≤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);
}