BZOJ_P3680 吊打XXX(模拟退火)

BZOJ传送门
Time Limit: 10 Sec Memory Limit: 128 MBSec Special Judge
Submit: 1505 Solved: 504
[Submit][Status][Discuss]
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

HINT

Source
By wangxz

模拟退火找重心，%gty学长

#include<cstdio>
#include<cmath>
#include<iostream>
using namespace std;
#define N 10005
inline int in(int x=0,char ch=getchar(),int v=1){
    while(ch!='-'&&(ch>'9'||ch<'0')) ch=getchar();if(ch=='-') v=-1,ch=getchar();
    while(ch>='0'&&ch<='9') x=(x<<3)+(x<<1)+ch-'0',ch=getchar();return x*v;}
int x[N],y[N],w[N];int n;double ansx,ansy,tmpx,tmpy,tmp,T;
inline double calc(double x1,double x2,double y1,double y2){return sqrt((x1-x2)*(x1-x2)+(y1-y2)*(y1-y2));}
int main(){
    n=in();
    for(int i=1;i<=n;i++) x[i]=in(),y[i]=in(),w[i]=in(),ansx+=x[i]*w[i],ansy+=y[i]*w[i];
    ansx/=n,ansy/=n;T=1;
    while(T>1e-8){
        tmpx=tmpy=0;
        for(int i=1;i<=n;i++){
            if((tmp=calc(ansx,x[i],ansy,y[i]))<=1e-10) continue;
            tmpx+=(x[i]-ansx)*w[i]/tmp,tmpy+=(y[i]-ansy)*w[i]/tmp;
        }

        ansx+=tmpx*T,ansy+=tmpy*T,T*=0.97;
    }
    printf("%.3lf %.3lf\n",ansx,ansy);
    return 0;
}