[JSOI2004]平衡点 / 吊打XXX--模拟退火

最新推荐文章于 2022-11-13 13:08:41 发布

小元勋

最新推荐文章于 2022-11-13 13:08:41 发布

阅读量153

点赞数 1

分类专栏：模拟退火文章标签：模拟退火

本文链接：https://blog.csdn.net/weixin_44574520/article/details/89047449

版权

模拟退火专栏收录该内容

1 篇文章 0 订阅

订阅专栏

[JSOI2004]平衡点 / 吊打XXX

在这里插入图片描述

导

Code：

#include <bits/stdc++.h>
using namespace std;
#define maxn 1010
#define down 0.996
int n;
double ansx=0,ansy=0,answ=0;
struct node {
	int x,y,w;
}e[maxn];

inline void init_() {
	freopen("phdiaod.txt","r",stdin);
}

inline int read_() {
	int x=0,f=1;
	char c=getchar();
	while(c<'0'||c>'9') {
		if(c=='-') f=-1;
		c=getchar();
	}
	while(c>='0'&&c<='9') {
		x=x*10+c-'0';
		c=getchar();
	}
	return x*f;
}

inline void readda_() {
	n=read_();
	for(int i=1;i<=n;i++) {
		e[i].x=read_();e[i].y=read_();e[i].w=read_();
		ansx+=e[i].x;ansy+=e[i].y;
	}
}

inline double  energy_(double x,double y) {//根据物理学知识,能量总和越小越稳定 
	double r=0,dx,dy;
	for(int i=1;i<=n;i++) {
		dx=x-e[i].x;
		dy=y-e[i].y;
		r+=sqrt(dx*dx+dy*dy) * e[i].w;
	}
	return r;
}

inline void SA_() {//模拟退火 
	double T=3000;//初始温度
	while(T>1e-15) {//略大于0 
		double ex=ansx+(rand()*2-RAND_MAX)*T;//随机产生新的答案 
		double ey=ansy+(rand()*2-RAND_MAX)*T;
		double ew=energy_(ex,ey);
		double DE=ew-answ;
		if(DE<0) {//如果此答案更优，就接受
			ansx=ex;
			ansy=ey;
			answ=ew;
		}
		else if(exp(-DE/T)*RAND_MAX>rand()) {//以一定概率接受 
			ansx=ex;
			ansy=ey;
		}
		T*=down;
	} 
}

inline void work_() {
	ansx/=n;ansy/=n;//以平均数作为初始答案
	answ=energy_(ansx,ansy);
	SA_();
	SA_();
	SA_();
	SA_();
	printf("%.3lf %.3lf",ansx,ansy);
}

int main() {
	init_();
	readda_();
	work_();
	return 0;
}