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

最新推荐文章于 2022-04-05 12:14:37 发布

滑稽大佬

最新推荐文章于 2022-04-05 12:14:37 发布

阅读量123

点赞数 1

分类专栏：模拟

本文链接：https://blog.csdn.net/qq_35786326/article/details/109154423

版权

模拟专栏收录该内容

23 篇文章 0 订阅

订阅专栏

玄学

题目：
题意：
分析：
代码：

题目：

传送门

题意：

有 $n$ 个重物，分别给出坐标和各自的重量
求一个公共点的坐标能使得所有重物保持一个平衡的状态

分析：

看到是模拟退火经典题就过来了 $q w q$
基本都是很常规的操作，对于新解的贡献我们可以通过绳子长度与物体重量来计算稳定程度

代码：

#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;
}