题目大意：有个人（gty）被吊打，他机智的使用了分身，但是每个分身有他的重力，把这些gty的分身绑起来，经过一个公共的绳结，求这个绳结最后在哪里。
思路：其实这个题就转化成了：定义一个点到一个分身的距离是两点间的距离 * 分身的重力。求平面内到这些点的距离的和的最小值。和poj2420差不多，这个题只需要在统计的时候吧权值乘上每个分身的重量就可以了。
值得一提的是这个题要求精度到1e-3，写的时候还是有点怕的，但是其实精度把握的合适还是能切这个题的。
一开始MAX开了10000，怎么改参数都是挂。。血的教训啊。。。
CODE：
#include
#include
#include
#include
#include
#define MAX 10010 #define INF 1e17 #define EPS 1e-3 #define PI acos(-1.0) using namespace std; struct Point{ double x,y; double weight; double total_weight; Point(double _x,double _y):x(_x),y(_y) {} Point() {} void Read() { scanf("%lf%lf%lf",&x,&y,&weight); } }point[MAX],now,ans; int points; inline double Calc(Point p1,Point p2); inline double Statistic(Point p); inline double Rand(); void SinulatedAnnealing(); int main() { srand(19980531); cin >> points; for(int i = 1;i <= points; ++i) { point[i].Read(); now.x += point[i].x; now.y += point[i].y; } now.x /= points,now.y /= points; ans.total_weight = INF; SinulatedAnnealing(); printf("%.3lf %.3lf\n",ans.x,ans.y); return 0; } inline double Calc(Point p1,Point p2) { return sqrt((p1.x - p2.x) * (p1.x - p2.x) + (p1.y - p2.y) * (p1.y - p2.y)); } inline double Statistic(Point p) { double re = 0.0; for(int i = 1;i <= points; ++i) re += Calc(p,point[i]) * point[i].weight; if(re < ans.total_weight) ans = p,ans.total_weight = re; return re; } void SinulatedAnnealing() { double T = 100000.0; while(T > EPS) { double alpha = 2.0 * PI * Rand(); Point temp(now.x + T * cos(alpha),now.y + T * sin(alpha)); double dE = Statistic(now) - Statistic(temp); if(dE >= 0 || exp(dE / T) >= Rand()) now = temp; T *= .99; } T = .001; for(int i = 1;i <= 1000; ++i) { double alpha = 2.0 * PI * Rand(); Point temp(ans.x + T * cos(alpha) * Rand(),ans.y + T * sin(alpha) * Rand()); Statistic(temp); } } inline double Rand() { return (rand() % 1000 + 1) / 1000.0; }