看到一道很有意思的題，這個題簡直有毒，是一道物理題。。。好像得用模擬退火。。。但顯然我太弱了不會模擬退火，只能用正交分解暴力。。。

每次沿著力的方向走一定的距離，假如轉頭了，則走的步長就減小一點。

不過這裏有一個坑，就是假如每次二分，就會錯。。。0.7或0.9就可以。有可能走過了回不來了吧。

題面：

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

代碼：

#include<iostream>
#include
 
<cstdio>
#include<algorithm>
#include<cmath>
#include<cstring>
using namespace std;
#define duke(i,a,n) for(int i = a;i <= n;i++)
#define lv(i,a,n) for(int i = a;i >= n;i--)
#define clean(a) memset(a,0,sizeof(a))
const int INF = 1 << 30;
typedef long long ll;
typedef 
 
double db;
template <class T>
void read(T &x)
{
    char c;
    bool op = 0;
    while(c = getchar(), c < ‘0‘ || c > ‘9‘)
        if(c == ‘-‘) op = 1;
    x = c - ‘0‘;
    while(c = getchar(), c >= ‘0‘ && c <= ‘9‘)
        x = x * 10 + c - ‘0‘;
    if(op) x = -x;
}
template <class T>
void write(T x)
{
    if(x < 0) putchar(‘-‘), x = -x;
    if(x >= 10) write(x / 10);
    putchar(‘0‘ + x % 10);
}
struct node
{
    db x,y,power;
}p[50000];
int n;
bool bx = true,by = true;
db x,y;
void solve(db move)
{
    db gx = 0,len = 0,gy = 0;
    duke(i,1,n)
    {
        len = sqrt((x - p[i].x) * (x - p[i].x) + (y - p[i].y) * (y - p[i].y));
        if(len == 0)
        continue;
        gx += p[i].power * (p[i].x - x) / len;
        gy += p[i].power * (p[i].y - y) / len;
    }
    len = sqrt(gx * gx + gy * gy);
    x += move * gx / len;
    y += move * gy / len; 
}
int main()
{
    read(n);
    duke(i,1,n)
    {
        scanf("%lf%lf%lf",&p[i].x,&p[i].y,&p[i].power);
    }
    db move = 5000,tx,ty;
    while(1)
    {
        tx = x;
        ty = y;
        solve(move);
        if(abs(tx - x) < 0.000001 && abs(ty - y) < 0.000001)
            break;
        if(bx != (x > tx) || by != (y > ty))
        {
            bx = !(x > tx);
            by = !(y > ty);
            move = move * 0.9;
        }
    }
    printf("%.3lf %.3lf\n",x,y);
    return 0;
}

B3680 吊打xxx 物理？？？