题目:
题意:
有
n
n
n个重物,分别给出坐标和各自的重量
求一个公共点的坐标能使得所有重物保持一个平衡的状态
分析:
看到是模拟退火经典题就过来了
q
w
q
qwq
qwq
基本都是很常规的操作,对于新解的贡献我们可以通过绳子长度与物体重量来计算稳定程度
代码:
#include<cstdio>
#include<string>
#include<cstring>
#include<iostream>
#include<algorithm>
#include<queue>
#include<cmath>
#include<cstdlib>
#include<ctime>
#define LL long long
#define down 0.9996
using namespace std;
inline LL read() {
LL d=0,f=1;char s=getchar();
while(s<'0'||s>'9'){if(s=='-')f=-1;s=getchar();}
while(s>='0'&&s<='9'){d=d*10+s-'0';s=getchar();}
return d*f;
}
double ansx=0,ansy=0,answ=0;
int n=read();
int x[1005],y[1005],w[1005];
double energy(double xx,double yy)
{
double s=0,dx,dy;
for(int i=1;i<=n;i++)
{
dx=xx-x[i]; dy=yy-y[i];
s+=sqrt(dx*dx+dy*dy)*w[i];
}
return s;
}
void sa()
{
double t=3333;
while(t>1e-15)
{
double ex=ansx+(rand()*2-RAND_MAX)*t,ey=ansy+(rand()*2-RAND_MAX)*t,ew;
ew=energy(ex,ey);
double bc=ew-answ;
if(bc<0) ansx=ex,ansy=ey,answ=ew;
else if(exp(-bc/t)*RAND_MAX>rand()) ansx=ex,ansy=ey;
t*=down;
}
return;
}
int main()
{
srand(1******7);
for(int i=1;i<=n;i++) ansx+=x[i]=read(),ansy+=y[i]=read(),w[i]=read();
ansx/=n;ansy/=n;answ=energy(ansx,ansy);
sa();sa();sa();sa();sa();
printf("%.3lf %.3lf",ansx,ansy);
return 0;
}