传送门
解题思路
特别狗的一道题
根据重力势能最小时平衡（？？），可以进行模拟退火。
其他的套板子行了。
然后我就写挂了（交了24次）
AC代码
#include<cstdio>
#include<iostream>
#include<cstring>
#include<cmath>
#include<algorithm>
#include<ctime>
#include<iomanip>
using namespace std;
const int maxn=1005;
const double delta=0.996;
int n;
double ans,ansx,ansy,sumx,sumy;
struct node{
double x,y,w;
}a[maxn];
double cal(double x,double y){
double ans=0;
for(int i=1;i<=n;i++){
double d=(x-a[i].x)*(x-a[i].x)+(y-a[i].y)*(y-a[i].y);
d=sqrt(d);
ans+=d*a[i].w;
}
return ans;
}
void SA(){
double t=3000,x=ansx,y=ansy;
while(t>1e-15){
double newx=x+(rand()*2-RAND_MAX)*t,newy=y+(rand()*2-RAND_MAX)*t;
double res=cal(newx,newy);
double cha=res-ans;
if(cha<0) ans=res,ansx=x=newx,ansy=y=newy;
else if(exp(-cha/t)*RAND_MAX>rand()) x=newx,y=newy;
t*=delta;
}
}
int main(){
ios::sync_with_stdio(false);
srand(time(NULL));
cin>>n;
for(int i=1;i<=n;i++){
cin>>a[i].x>>a[i].y>>a[i].w;
sumx+=a[i].x;
sumy+=a[i].y;
}
ansx=sumx/n;
ansy=sumy/n;
ans=1e18;
while((double)clock()/CLOCKS_PER_SEC<=0.85) SA();
cout<<fixed<<setprecision(3)<<ansx<<" "<<ansy;
return 0;
}