BZOJ3680 吊打XXX（模拟退火算法）

最新推荐文章于 2021-10-21 14:49:42 发布

原创最新推荐文章于 2021-10-21 14:49:42 发布 · 945 阅读

0 ·

CC 4.0 BY-SA版权

文章标签：

#模拟退火

近似算法专栏收录该内容

1 篇文章

订阅专栏

本文探讨了一种使用力学原理进行路径优化的算法，通过模拟物理运动来解决实际问题，实现路径的有效优化。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

其实我并不确定这是什么算法。

代码参考自hzwer巨神

#include<stdio.h>
#include<stdlib.h>
#include<math.h>
#define eps 1e-6
double x[10005],y[10005],w[10005];
double dis(double x1,double y1,double x2,double y2)
{
	double d=sqrt( (x1-x2)*(x1-x2) + (y1-y2)*(y1-y2) );
	if(d>-eps&&d<eps) return eps;
	return d;
}
int main()
{
	double t,xx=0.0,yy=0.0,fx,fy,F;//设初始位置为(0,0)
	int n,i;
	scanf("%d",&n);
	for(i=1;i<=n;i++)
		scanf("%lf%lf%lf",&x[i],&y[i],&w[i]);
	for(t=10000;t>0.00001;t*=0.99)//t:步长 
	{
		fx=fy=0;
		for(i=1;i<=n;i++)
		{
			fx+=w[i]*(x[i]-xx)/dis(x[i],y[i],xx,yy);//fx:水平合外力 
			fy+=w[i]*(y[i]-yy)/dis(x[i],y[i],xx,yy);//fy:竖直合外力 
		}
		F=dis(fx,fy,0,0);//F:合外力，(fx/F,fy/F)为合外力的方向向量 
		if(F<0.0001) break;
		xx+=t*fx/F;
		yy+=t*fy/F;
	}
	printf("%.3lf %.3lf",xx,yy);
	return 0;
}