乱搞-模拟退火

bzoj3680

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

这道题我有一点点思路，就是大概是力矩最小的时候比较靠谱吧，物理题

然后就不会了

貌似是模拟退火啊

每次高出答案后算一下剩余点的力的矢量合成

然后答案向这个方向偏移——感性的说，就是这个终点被拉过去了，然后拉不动的时候就停啦

顺便讨一个模拟退火的模板

hzwer的模板

#include<iostream>
#include<cstdio>
#include<cmath>
using namespace std;
int n;
double ansx,ansy;
struct data{double x,y;int w;}p[10005];
inline double sqr(double x){return x*x;}
inline double dis(double x,double y,data p)
{
    return sqrt(sqr(x-p.x)+sqr(y-p.y));
}
void hillclimb()
{
     double t=1000,x,y;
     for(int i=1;i<=n;i++)
         ansx+=p[i].x*p[i].w,ansy+=p[i].y*p[i].w;
     ansx/=n;ansy/=n;
     while(t>0.00000001)
     {
         x=y=0;
         for(int i=1;i<=n;i++)
         {
             x+=(p[i].x-ansx)*p[i].w/dis(ansx,ansy,p[i]);
             y+=(p[i].y-ansy)*p[i].w/dis(ansx,ansy,p[i]);
         }
         ansx+=x*t;ansy+=y*t;
         if(t>0.5)t*=0.5;
         else t*=0.97;
     }
     printf("%.3lf %.3lf",ansx,ansy);
}
int main()
{
     scanf("%d",&n);
     for(int i=1;i<=n;i++)
         scanf("%lf%lf%d",&p[i].x,&p[i].y,&p[i].w);
     hillclimb();
     return 0;
}