版权声明：本文为博主原创文章，未经博主允许不得转载。 https://blog.csdn.net/sxy201658506207/article/details/83590140
题意：平面上给你n个点，让你求一个点，到这n点的距离和最小
板子得模拟退火，系数也是很好控制，对数器和ACdream得代码跑了1884组样例，结果只有一组不一样,可是我得解更优啊, 交上去RE,
不想再改了,
真心感觉或者模拟退火算法是个假得算法,a不a多数靠运气啊,
脱坑
6
6878 4617
2040 3626
2067 9225
9729 4725
8635 7922
5480 3415
20010//我的跑出来的结果
20040//ACdream跑出来的结果
哇,A了A了，第二天，想着用VS2017debug一番,果真有不通过的地方 原来是他妈~scanf("%d",&n); 把他改成EOF就过了 操
re到哭，长记性了
#include <iostream>
#include <cmath>
#include <cstdio>
#include <stdlib.h>
#include <ctime>
#define mp make_pair
#define sz(x) int((x).size())
#define fin freopen("in.txt","r",stdin)
#define fout freopen("out.txt","w",stdout)
#define io ios::sync_with_stdio(0),cin.tie(0)
using namespace std;
typedef long long ll;
typedef pair<int, int> PII;
const int inf = 0x3f3f3f3f;
//const int mod = 1e9 + 7;
//const int maxn = 1e6 + 5;
int n, x, y;
struct point { double x, y, w; } p[1000006], now, nex, ansp;
double dis(point a, point b) { return sqrt((a.x - b.x)*(a.x - b.x) + (a.y - b.y)*(a.y - b.y)); }
double f(point x)
{ //评估函数
double res = 0;
for (int i = 1;i <= n;i++) res += dis(x, p[i]);
return res;
}
double ans = 1e20;//最开始的能量值,初始很大就可以,不用修改
void sa()
{
ans = 1e20;
double T = 100; //初始温度, (可以适当修改,最好和给的数据最大范围相同,或者缩小其原来0.1)
double d = 0.99; //降温系数 (可以适当修改,影响结果的精度和循环的次数,)
double eps = 1e-8; //最终温度 (要是因为精度问题，可以适当减小最终温度)
double TT = 1.0; //采纳新解的初始概率
double dd = 0.98; //(可以适当修改,采纳新解变更的概率)(这个概率下面新解更新的时候,最好和未采纳的新解更新的次数是一半一半)
double res = f(now); //传入的初始默认解(now)下得到的评估能量值
if (res < ans) ans = res, ansp = now;//ansp终解
while (T > eps)
{
for (int i = -1;i <= 1;++i)
for (int j = -1;j <= 1;++j)
if ((i || j) && (now.x + T * i <= 10000) && (now.x + T * i >= 0) && (now.y + T * j <= 10000) && (now.y + T * j >= 0))
{
nex.x = now.x + T * i, nex.y = now.y + T * j;//新解
double tmp = f(nex);//新解下的评估能量值
if (tmp < ans) ans = tmp, ansp = nex;//降温成功,更新当前最优解
if (tmp < res) res = tmp, now = nex;// 降温成功,采纳新解
else if (TT > rand() % 10000 / 10000.0) res = tmp, now = nex;//,cout<<"======"<<endl;//没有 降温成功，但是以一定的概率采纳新解
//else cout<<"="<<endl;//用于测试，设定的采纳新解的概率,是否为一半一半,可以适当修改降温参数dd
}
T *= d; TT *= dd;
}
}
int main()
{
srand(time(0));
while (scanf("%d", &n!=EOF)
{
now.x = now.y = 0;
for (int i = 1;i <= n;++i)
{
cin >> p[i].x >> p[i].y;//scanf("%lf%lf%lf",&p[i].x,&p[i].y),
now.x += p[i].x, now.y += p[i].y;
}
now.x /= n, now.y /= n; sa();
printf("%.0f\n", ans);//cout<<ans<<endl;
}
return 0;
}