洛谷 P1337 [JSOI2004]平衡点 / 吊打XXX

最新推荐文章于 2022-11-23 19:25:41 发布

verdin黄大锤

最新推荐文章于 2022-11-23 19:25:41 发布

阅读量264

点赞数

CC 4.0 BY-SA版权

分类专栏： # 模拟退火

本文链接：https://blog.youkuaiyun.com/rabbit_ZAR/article/details/86673439

模拟退火专栏收录该内容

1 篇文章

订阅专栏

本文探讨了在给定点集上寻找平衡点的问题，通过模拟退火算法进行优化，实现对复杂数据集的高效处理。文章详细介绍了模拟退火算法的原理及其在平衡点求解中的应用。

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

题目：平衡点 / 吊打XXX

思路：
退个火？
关于退火，安利。

代码：

#include<bits/stdc++.h>
using namespace std;

#define maxn 1000
#define read(x) scanf("%d",&x)
#define db double

#define eps 1e-15

int n;
int X[maxn+5],Y[maxn+5],w[maxn+5];

db ansx=0,ansy=0;

void readin() {
	read(n);
	for(int i=1;i<=n;i++) {
		read(X[i]),read(Y[i]),read(w[i]);
		ansx+=X[i],ansy+=Y[i];
	}
	ansx/=n,ansy/=n;
}

db slv(db x,db y) {
	db s=0;
	for(int i=1;i<=n;i++) {
		db delx=x-X[i],dely=y-Y[i];
		s+=sqrt(delx*delx+dely*dely)*w[i];
	}
	return s;
}

void mnth() {
	db T=200;
	while(T>eps) {
		db nowx=ansx+(rand()*2-RAND_MAX)*T;
		db nowy=ansy+(rand()*2-RAND_MAX)*T;
		db delta=slv(nowx,nowy)-slv(ansx,ansy);
		if(delta<0) ansx=nowx,ansy=nowy;
		else if(exp(-delta/T)*RAND_MAX>rand()) ansx=nowx,ansy=nowy;
		T*=0.99;
	}
}

int main() {
	srand((int)time(NULL));
	
	readin();
	mnth();
	printf("%.3lf %.3lf",ansx,ansy);
	
	return 0;
}