【题目大意】:找出一个点,使得到多边形每个点的距离和最小。输出最小距离和。
【解题思路】:求的是一个多边形的费马点。三边以上没有公式,我们可以使用随机化变步长贪心法来解这道题。
【随机化变步长贪心法】:随机选取一个点,再取一个步长,朝这个方向走,如果新位置到各点距离比原来小,则走过去。直到走不动为止,再缩小步长。直到步长小于题目精度。
【代码】:
#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
#include <vector>
#include <queue>
#include <cmath>
#include <string>
#include <cctype>
#include <map>
#include <iomanip>
using namespace std;
#define eps 1e-8
#define pi acos(-1.0)
#define inf 1<<30
#define linf 1LL<<60
#define pb push_back
#define lc(x) (x << 1)
#define rc(x) (x << 1 | 1)
#define lowbit(x) (x & (-x))
#define ll long long
struct Point{
double x,y;
Point(){}
Point(double a,double b){
x=a,y=b;
}
}point[101];
double pt_distance(const Point &p1,const Point &p2){
return sqrt((p1.x-p2.x)*(p1.x-p2.x)+(p1.y-p2.y)*(p1.y-p2.y));
}
double get_all_dis(const Point &p,int n){
double ans=0.0;
for (int i=0; i<n; i++) ans+=pt_distance(point[i],p);
return ans;
}
int main() {
int n;
while (~scanf("%d",&n)){
for (int i=0; i<n; i++)
scanf("%lf%lf",&point[i].x,&point[i].y);
Point st=point[0];
double step=100,mind=get_all_dis(st,n);
while (step>0.2){
int ok=1;
while (ok){
Point tmp,nt;
double t;
ok=0,nt=st;
tmp=Point(st.x,st.y+step);
t=get_all_dis(tmp,n);
if (t<mind) mind=t,ok=1,nt=tmp;
tmp=Point(st.x,st.y-step);
t=get_all_dis(tmp,n);
if (t<mind) mind=t,ok=1,nt=tmp;
tmp=Point(st.x+step,st.y);
t=get_all_dis(tmp,n);
if (t<mind) mind=t,ok=1,nt=tmp;
tmp=Point(st.x-step,st.y);
t=get_all_dis(tmp,n);
if (t<mind) mind=t,ok=1,nt=tmp;
st=nt;
}
step=step/2.0;
}
int ans=(int)(mind+0.5)*100/100;
printf("%d\n",ans);
}
return 0;
}