在二维平面上有 n 个点,第 i 个点的坐标为 (xi,yi)。请你找出一个点,使得该点到这 n个点的距离之和最小。该点可以选择在平面中的任意位置,甚至与这 n个点的位置重合。
输入格式
第一行包含一个整数 n。接下来 n行,每行包含两个整数 xi,yi,表示其中一个点的位置坐标。
输出格式
输出最小距离和,答案四舍五入取整。
数据范围
1≤n≤100,0≤xi,yi≤10000
输入样例:
4 0 0 0 10000 10000 10000 10000 0
输出样例:
28284
详解:代码注释(提前了解模拟退火算法的过程)
AC代码:
#include<iostream> #include<algorithm> #include<cstring> #include<cmath> #include<vector> #include<utility> #include<stdlib.h> #include<stdio.h> using namespace std; typedef long long ll; #define x first //给pair的两个成员起别名 #define y second #define N 110 //最大范围 typedef pair<double, double> pdd; int n; pdd p[N]; //存坐标的变量 double ans = 1e9; //存答案的变量 double dist(pdd a, pdd b) //计算两点距离 { double dx = a.x - b.x; double dy = a.y - b.y; return sqrt(dx*dx+dy*dy); } double s_dist(pdd a) //将随机到的一个点与输入的点计算距离 { double ret = 0; //存随机点的距离 for (int i = 0; i < n; i++)ret += dist(p[i], a); ans = min(ret, ans); //与答案比较 return ret; } double rand(int a, int b) //生成(a,b)的随机数 { return (double)rand() / RAND_MAX * (b - a) + a; } void th() //模拟退火 { pdd a(rand(0, 1e4), rand(0, 1e4)); //随机一个初始点 for (double i = 1e4; i >= 1e-4; i *= 0.99) //退火 for(初始温度;最低温度;衰减因子) { pdd temp(rand(a.x-i,a.x+i), rand(a.y-i,a.y+i)); //新随机点 double sub = s_dist(temp) - s_dist(a); //比较两点 //sub<0时,新点较小,所以必然执行a=temp (exp(-x)>1>rand(0,1)),sub>0时可能执行a=temp if (exp(-sub / i) > rand(0, 1))a = temp; } } int main() { cin >> n; for (int i = 0; i < n; i++)cin >> p[i].x >> p[i].y; for (int i = 0; i < 100; i++)th(); //进行100次模拟退火 printf("%.0f", ans); return 0; }