模拟退火

好文推荐

（我就懒得再讲了）

1. 模拟退火总结（模拟退火）
2. 最优化算法之模拟退火算法（SA）

例题

吊打XXX luogu1337

题目传送门：吊打XXX/平衡点

分析

到达平衡，也就是总势能最小时

此时系统的总能量是等于各个物体的重力势能，在质量一定时，即要求物体离地最近，离桌子最远。

那么，也就是绳子在桌子上的距离尽量的小。及每个洞到平衡点的距离与相应重物的重量的乘积之和最小

代码

/*
User:Mandy
Language:c++
Problem:1337
*/
//emmmm,注意double,然后参数可以调大点 ，注意时间 
#include<bits/stdc++.h>
#define RD T*(rand()*2-RAND_MAX)

using namespace std;

const int maxn=1005;
const long double EPS=1e-16;
const long double D=0.98;

int n;
long double ansx,ansy,minans,curans;
long double x[maxn],y[maxn],w[maxn];

void readdata()
{
	scanf("%d",&n);
	for(int i=1;i<=n;++i)
	{
		scanf("%Lf%Lf%Lf",&x[i],&y[i],&w[i]);
		ansx+=x[i];
		ansy+=y[i];
	}
}

long double get(long double x0,long double y0)
{
	long double sum=0;
	for(int i=1;i<=n;++i)
	{
		long double nx=x0-x[i];
		long double ny=y0-y[i];
		sum+=(sqrt(nx*nx+ny*ny))*w[i];
	}
	return sum;
}

void work()
{
	srand(20030110);
	ansx/=n;ansy/=n;
	minans=get(ansx,ansy);
	int Time=2;
	long double curx,cury;
	while(Time--)
	{
		
		curans=minans;curx=ansx;cury=ansy;
		for(long double T=100000;T>EPS;T*=D)
		{
			
			long double xx=curx+RD,yy=cury+RD;
			long double sum=get(xx,yy);
			//printf("%.5Lf %.5Lf %.5Lf\n",xx,yy,sum);
			if(minans>sum)
			{
				minans=sum;ansx=xx;ansy=yy;
			}
			if(curans>sum||exp((curans-sum)/T)>(long double)rand()/RAND_MAX)
			{
				curans=sum;curx=xx;cury=yy;
			}
		}
	}
	printf("%.3Lf %.3Lf",ansx,ansy);
}

int main()
{
	readdata();
	work();
	return 0;
}