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
0 0 1
0 2 1
1 1 1
Sample Output
0.577 1.000
HINT
Source
By wangxz
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
模拟退火~
推荐一篇讲得相当完美的博客：http://www.cnblogs.com/heaad/archive/2010/12/20/1911614.html~
注意：
1.最大值要取得足够大，否则……
2.所有函数返回值都是double！
#include<cstdio>
#include<cstdlib>
#include<iostream>
#include<cmath>
#include<ctime>
using namespace std;int n;
double minn=99999999999999999ll;struct node{double x,y,p;
}a[10001],ans;double dis(node u,node v)
{return sqrt((u.x-v.x)*(u.x-v.x)+(u.y-v.y)*(u.y-v.y));
}double ran()
{return (rand()%1000)/1000.0;
}double jud(node now)
{double re=0;for(int i=1;i<=n;i++) re+=dis(a[i],now)*a[i].p;if(re<minn) minn=re,ans=now;return re;
}void cal()
{node now=ans,ne;double T=1e6,dE;while(T>1e-3){ne.x=now.x+T*(ran()*2-1);ne.y=now.y+T*(ran()*2-1);dE=jud(now)-jud(ne);if(dE>0 || exp(dE/T)>ran()) now=ne;T*=0.993;}for(int i=1;i<=1000;i++){ne.x=ans.x+T*(ran()*2-1);ne.y=ans.y+T*(ran()*2-1);jud(ne);}
}int main()
{srand(20002333);scanf("%d",&n);for(int i=1;i<=n;i++){scanf("%lf%lf%lf",&a[i].x,&a[i].y,&a[i].p);ans.x+=a[i].x;ans.y+=a[i].y;}ans.x/=n;ans.y/=n;cal();printf("%.3lf %.3lf\n",ans.x,ans.y);return 0;
}