BZOJ3680 : 吊打XXX

本题就是找一个受力平衡的点

我们一开始假设这个点是(0,0)

然后求出它受到的力，将合力正交分解后朝着合力的方向走若干步，并不断缩小步长，一步步逼近答案

 

#include<cstdio>

#include<cmath>

#define N 10000

int n,i;double x[N],y[N],w[N],a,b,sumx,sumy,dx,dy,dis,step=1;

double sqr(double x){return x*x;}

int main(){

  for(scanf("%d",&n);i<n;i++)scanf("%lf%lf%lf",x+i,y+i,w+i);

  while(step>1e-8){

    for(sumx=sumy=i=0;i<n;i++){

      if((dis=std::sqrt(sqr(dx=x[i]-a)+sqr(dy=y[i]-b)))<=1e-10)continue;

      sumx+=dx*w[i]/dis,sumy+=dy*w[i]/dis;

    }

    a+=sumx*step,b+=sumy*step,step*=0.97;

  }

  return printf("%.3f %.3f",a,b),0;

}