BZOJ 3680 吊打XXX 模拟退火

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

PoPoQQQ

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

阅读量4.2k

点赞数 1

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

本文链接：https://blog.csdn.net/PoPoQQQ/article/details/39322569

版权

博客介绍了BZOJ 3680题目的解决过程，该题是一个模拟退火算法的应用问题，要求找到n个质点的重心，使得所有质点到重心距离乘以其权重g[i]之和最小。博主在实现过程中遇到了INF值设置不当导致的错误，最终通过调整INF值至极大整数附近成功解决问题，但程序运行时间达到了10188MS，成为BZOJ上的倒数第一。

摘要由CSDN通过智能技术生成

首先这题应该改名叫吊打出题人

题目大意：给定n个质点，求重心

这n个质点的重心满足Σ(重心到点i的距离)*g[i]最小

模拟退火的裸题

尼玛交了两篇死活过不去各种改参数最后发现是我的INF不够大尼玛！

这题INF开0x3f妥妥过不去。。。起码要max_of _long_long附近才可以

最后写了10188MS，BZOJ倒数第一……这也是种艺术啊0.0

#include<cmath>
#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
#define M 10100
using namespace std;
struct point{
	double x,y,g;
}points[M],ans;
int n;
double minans=23333333333333333ll;
double dis(const point &x,const point &y)
{
	return sqrt( (x.x-y.x)*(x.x-y.x) + (x.y-y.y)*(x.y-y.y) );
}
double Judge(const point &p)
{
	int i;
	double re=0;
	for(i=1;i<=n;i++)
		re+=points[i].g*dis(p,points[i]);
	if(re<minans)
		ans=p,minans=re;
	return re;
}
double Rand()
{
	return rand()%1000/1000.0;
}
void SA(double T)
{
	int i;
	point Now=ans;
	while(T>0.001)
	{
		point Ne