平衡点(balance)题解

题目描述

现有一张无限大的桌子（二维平面），桌面上有n个洞，每个洞有一根绳子穿过挂着一个重物，质量为wi，洞的坐标为（xi，yi）。这n根绳子有个公共的绳结，求系统平衡后绳结的坐标。

输入

输入第一行为一个正整数n(1<=n<=10000)，表示洞的数目。
接下来n行,每行三个整数xi，yi，wi，表示第i个重物的横坐标，纵坐标和重力。

输出

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

样例输入

3
0 0 1
0 2 1
1 1 1

样例输出

0.577 1.000

提示

对于20%的数据，桌面上的洞排列成一条直线。

对于另外40%的数据，1<=n<=1000，-10000<=xi,yi<=10000。

对于100%的数据，1<=n<=10000,-100000<=xi,yi<=100000,wi<=1000

想法

• 考场上只想出20%的写法 一条直线。。。
• 分析以后不难发现该题要求一个广义的费马点(完全没看出来。。)
• 想不出有什么靠谱的算法(还是太弱。。)

算法

• 实际上该题完全可以用一个类似随机化算法的思想来做
• 大题思路有点像模拟退火(而不是爬山！)
• 先求出合力的方向 再按照合力的方向迭代
• 据某神犇(isdkfj)验证迭代次数大概只有几百。。

代码

#include <iostream>
#include <cstdio>
#include <cstdlib>
#include <algorithm>
#include <cstring>
#include <string>
#include <cmath>
#define MAXN 10005
#define eps 1e-7
using namespace std;
int n,cnt;
double d=10000.0;
struct Tpoint
{
    double x,y,w;
}point[MAXN],pos,now;
inline double cal(Tpoint a,Tpoint b)
{
    return sqrt((a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y));
}
int main()
{
    scanf("%d",&n);
    for (int i=1;i<=n;i++)
    {
        scanf("%lf%lf%lf",&point[i].x,&point[i].y,&point[i].w);
        pos.x+=point[i].x,pos.y+=point[i].y,cnt+=point[i].w;
    }
    pos.x/=cnt,pos.y/=cnt;
    while(d-eps>=0)
    {
        now.x=now.y=0;
        for (int i=1;i<=n;i++)
        {
            double dis=cal(pos,point[i]);
            //printf("%lf\n",dis);
            now.x+=(point[i].x-pos.x)/dis*point[i].w;
            now.y+=(point[i].y-pos.y)/dis*point[i].w;
        }
        double mc=sqrt(now.x*now.x+now.y*now.y);
        if(mc<eps)break;
        pos.x+=now.x/mc*d;
        pos.y+=now.y/mc*d;
        d*=0.65;
    }
    printf("%.3lf %.3lf\n",pos.x,pos.y);
    return 0;
}