【题解】洛谷P1337[JSOI2004]平衡点/吊打XXX 模拟退火

29 篇文章 0 订阅
2 篇文章 0 订阅

题目链接
在这里插入图片描述
在这里插入图片描述


跟着大佬博客学模拟退火,顺带把这题做()了。

#include<cstdio>
#include<cstdlib>
#include<ctime>
#include<cmath>
const int N=1e3+10;
int n;
const double MAX_TIME=0.777;
double ans=1e18,ansx,ansy;
struct node{
	int x,y,w;
}a[N];
double calc(double x,double y)
{
	double now=0;
	for(int i=1;i<=n;i++)
	{
		double nx=x-a[i].x,ny=y-a[i].y;
		now+=sqrt(nx*nx+ny*ny)*a[i].w;
	}
	return now;
}
void SA()
{
	double t=2000,x=ansx,y=ansy;
	while(t>1e-14)
	{
		double nx=x+(rand()*2-RAND_MAX)*t;
		double ny=y+(rand()*2-RAND_MAX)*t;
		double now=calc(nx,ny);
		if(now<ans)
		    x=nx,y=ny,ansx=nx,ansy=ny,ans=now;
		else if(exp((ans-now)/t)*RAND_MAX>rand())x=nx,y=ny;
		t*=0.996;
	}
}
int main()
{
	//freopen("in.txt","r",stdin);
    srand((unsigned)time(0));
    scanf("%d",&n);
    for(int i=1;i<=n;i++)
        scanf("%d%d%d",&a[i].x,&a[i].y,&a[i].w),ansx+=a[i].x,ansy+=a[i].y;
    ansx/=n;ansy/=n;
    while ((double)clock()/CLOCKS_PER_SEC<MAX_TIME) SA();
    printf("%.3f %.3f\n",ansx,ansy);
    return 0;
}

总结

模拟退火入门

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值