1262: 邮局选址问题(模拟退火)

解题思路:

注意事项:

参考代码:

#include<bits/stdc++.h>
using namespace std;
struct poit
{
	int x;
	int y;
}p[10005];
double ansx,ansy,answ;
int n;
double energy(double x,double y)
{
	double ans=0;
	for(int i=0;i<n;i++)
	{
		ans+=abs(x-p[i].x)+abs(y-p[i].y);
	}
	return ans;
}
void sa()//模拟退火就是赌命所以答案不对的话就换换T和delta
{
	double T=3505;//温度
	double delta=0.996;//退火系数
	while(T>1e-15)
	{
		int f[2]={1,-1};
		double newx=ansx+f[rand()%2]*T;//都随机变化
		double newy=ansy+f[rand()%2]*T;
		double neww=energy(newx,newy);//算出随机变化的值
		double flag=neww-answ;
		if(flag<0)/如果新的值更好就用新的值
		{
			ansx=newx;
			ansy=newy;
			answ=neww;
		}
		T*=delta;//退火
	}
}
int main()
{
	cin>>n;
	double x=0,y=0;
	for(int i=0;i<n;i++)
	{
		cin>>p[i].x>>p[i].y;
		x+=p[i].x;
		y+=p[i].y;
	}
	ansx=x/n;//先随便找个值
	ansy=y/n;
	//cout<<ansx<<' '<<ansy;
	answ=energy(ansx,ansy);//将随便找到的x,y的距离和保存
	for(int i=0;i<5;i++)//时间允许可以多退几次，保证答案正确
	{
		sa();//模拟退火
	}
	printf("%.0lf",answ);
	return 0;
}