看到一道很有意思的题，这个题简直有毒，是一道物理题。。。好像得用模拟退火。。。但显然我太弱了不会模拟退火，只能用正交分解暴力。。。
每次沿着力的方向走一定的距离，假如转头了，则走的步长就减小一点。
不过这里有一个坑，就是假如每次二分，就会错。。。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 ",x,y); return 0; }