BZOJ3680 吊打XXX

Description

gty又虐了一场比赛，被虐的蒟蒻们决定吊打gty。gty见大势不好机智的分出了n个分身，但还是被人多势众的蒟蒻抓住了。蒟蒻们将
n个gty吊在n根绳子上，每根绳子穿过天台的一个洞。这n根绳子有一个公共的绳结x。吊好gty后蒟蒻们发现由于每个gty重力不同，绳
结x在移动。蒟蒻wangxz脑洞大开的决定计算出x最后停留处的坐标，由于他太弱了决定向你求助。
不计摩擦，不计能量损失，由于gty足够矮所以不会掉到地上。

Input

输入第一行为一个正整数n(1<=n<=10000)，表示gty的数目。
接下来n行,每行三个整数xi，yi，wi，表示第i个gty的横坐标，纵坐标和重力。
对于20%的数据，gty排列成一条直线。
对于50%的数据，1<=n<=1000。
对于100%的数据，1<=n<=10000,-100000<=xi,yi<=100000

Output

输出1行两个浮点数（保留到小数点后3位），表示最终x的横、纵坐标。
Sample Input

3

0 0 1

0 2 1

1 1 1

Sample Output

0.577 1.000

PoPoQQQ老师使用模拟退火解决的经典问题，不过我没能理解为什么重心到各点距离乘以质量之和最小。。。希望知道的人能够赐教
通过向量加法，可以不断沿合力方向移动来逼近最终解。开始时每一次移动很大距离，之后越来越小，直到停止。
注意要把合向量化为那个方向的单位向量然后再乘上本次要移动的距离。
自己的代码

#include<cstdio>
#include<cmath>
using namespace std;
struct vec
{
    double x,y;
    vec(double x,double y):x(x),y(y){}
    vec operator +(const vec &b)const
    {return vec(x+b.x,y+b.y);}
    vec operator *(const double &b)const
    {return vec(x*b,y*b);}
    vec operator /(const double &b)const
    {return vec(x/b,y/b);}
    double len()
    {return sqrt(x*x+y*y);}
    vec one(double k)
    {return *this/len()*k;}
}v(0,0);
struct pnt
{
    double x,y;
    pnt():x(),y(){}
    pnt(double x,double y):x(x),y(y){}
    vec operator -(const pnt &b)const
    {return vec(x-b.x,y-b.y);}
    pnt operator +(const vec &b)const
    {return pnt(x+b.x,y+b.y);}
}t[10000],m;
int n,g[10001];
#define getv\
    v.x=v.y=0;\
    for(int zy=0;zy<n;zy++)\
    v=v+(t[zy]-m).one(g[zy])
int main()
{
    scanf("%d",&n);
    for(int i=0;i<n;i++)
    {
        scanf("%lf%lf%d",&t[i].x,&t[i].y,&g[i]);
        m.x+=t[i].x;
        m.y+=t[i].y;
    }
    m.x/=n;
    m.y/=n;
    getv;
    for(double i=1000;i>=0.0001;i*=0.5)
    {
        m=m+(v.one(1)*i);
        getv;
    }
    printf("%.3lf %.3lf",m.x,m.y);
    return 0;
}