标签：const Point Meow Country int 模拟退火 double start

题意

三维空间有\(n\)个点，找到另外一个点，离所有点的最大距离最小。求这个距离。

题解

\(1\)、最小球覆盖，要找的点为球心。

\(2\)、模拟退火。

还是补一下模拟退火的介绍吧。

模拟退火有一个初始温度，温度越高，接受较差的解的可能性就越大。每次走完后，都会降低温度，使得接受较差解的可能性变小。在走的过程中，更新最优解的值。

对这个题目来说，从空间中某一个点出发，如果每次都找离当前点最远的点，往那个点的方向走，大概率可以使结果变得更优。

随便设了个温度下降速率为\(0.97\)，一遍就AC了。参数都不用调。

#include <bits/stdc++.h>

#define FOPI freopen("in.txt", "r", stdin)
#define FOPO freopen("out.txt", "w", stdout)

using namespace std;
typedef long long LL;
const int maxn = 100 + 5;
const double eps = 1e-8;
const int inf = 0x3f3f3f3f;
const double start_T = 1000;

struct Point
{
    double x, y, z;
    Point() {}
    Point(int _x, int _y, int _z):x(_x), y(_y), z(_z) {}
}a[maxn];

double dist(Point a, Point b)
{
    return sqrt((a.x-b.x)*(a.x-b.x) + (a.y-b.y)*(a.y-b.y) + (a.z-b.z)*(a.z-b.z));
}

int n;
double SA()
{
    Point p = Point(0,0,0);
    int s = 0;
    double ans = inf, rate = 0.97, T = start_T;
    while(T > eps)
    {
        for (int i = 1; i <= n; i++)
            if(dist(p, a[i]) > dist(p, a[s])) s = i;
        double d = dist(p, a[s]);
        ans = min(ans, d);
        p.x += (a[s].x - p.x) * T / start_T;
        p.y += (a[s].y - p.y) * T / start_T;
        p.z += (a[s].z - p.z) * T / start_T;
        T *= rate;
    }
    return ans;
}

int main()
{
//    FOPI;
    scanf("%d", &n);
    for (int i = 1; i <= n; i++)
        scanf("%lf%lf%lf", &a[i].x, &a[i].y, &a[i].z);
    double ans = SA();
    printf("%.8f\n", ans);
}

标签：const,Point,Meow,Country,int,模拟退火,double,start
来源： https://www.cnblogs.com/ruthank/p/10910513.html

本站声明： 1. iCode9 技术分享网（下文简称本站）提供的所有内容，仅供技术学习、探讨和分享；
2. 关于本站的所有留言、评论、转载及引用，纯属内容发起人的个人观点，与本站观点和立场无关；
3. 关于本站的所有言论和文字，纯属内容发起人的个人观点，与本站观点和立场无关；
4. 本站文章均是网友提供，不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属；如您发现该文章侵犯了您的权益，可联系我们第一时间进行删除；
5. 本站为非盈利性的个人网站，所有内容不会用来进行牟利，也不会利用任何形式的广告来间接获益，纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。