zoukankan      html  css  js  c++  java
  • Country Meow(模拟退火入门)

    题目链接:http://codeforces.com/gym/101981/attachments

    题意:给你n个三维空间点(n<=100)要求找到一个离所有点最远的距离最小的点,并输出这个距离

    思路:模拟退火,首先取一个值,代表答案点在的位置,然后暴力一遍所有的点,会得到里当前点最远的点是哪个

    然后更新当前ans ,并且往这个最远的点移动,贪心,最后就是答案了

    看代码:

    #include<iostream>
    #include<cstdio>
    #include<cstring>
    #include<algorithm>
    #include<cmath>
    #include<vector>
    #include<map>
    #include<cstdlib>
    using namespace std;
    typedef long long LL;
    const LL mod=1e9+7;
    const LL INF=1e9+7;
    const int maxn=1e2+50;
    const double StartT=50000;//初始温度 1000也可以过这道题
    const double rate=0.99;//退火系数 一般0.97-0.99
    const double eps=1e-6;//精度 可以在题目要求的小数点后面位数多加3左右
    const double FINF=1e18;
    /**
    没有算退火的复杂度
    rate=0.98
    T=1000 1026次
    T=20000 1175次
    T=50000 1220次
    
    rate=0.99
    T=50000 2452
    **/
    int N;
    struct Node
    {
        double x,y,z;
        Node(){}
        Node(double x1,double y1,double z1)
        {
            x=x1;y=y1;z=z1;
        }
    }a[maxn];
    double dist(Node n1,Node n2)
    {
        return sqrt((n1.x-n2.x)*(n1.x-n2.x)+(n1.y-n2.y)*(n1.y-n2.y)+(n1.z-n2.z)*(n1.z-n2.z));
    }
    void solve()
    {
    //    int sum=0;
        double ans=FINF;
    //    cout<<ans<<endl;
        double T=StartT;
        Node now=Node(0,0,0);//随机选取一个点
        Node max_p=a[1];//离当前最远的点
        while(T>eps)
        {
    //        sum++;
            for(int i=1;i<=N;i++)//寻找里当前最远的点
            {
                if(dist(a[i],now)>=dist(max_p,now))
                {
                    max_p=a[i];
                }
            }
            ans=min(ans,dist(now,max_p));
            now.x=(now.x+(max_p.x-now.x)*(T/StartT));
            now.y=(now.y+(max_p.y-now.y)*(T/StartT));
            now.z=(now.z+(max_p.z-now.z)*(T/StartT));
            T=T*rate;
        }
    //    cout<<"sum:"<<sum<<endl;
        printf("%.15lf
    ",ans);
    }
    int main()
    {
        scanf("%d",&N);
        for(int i=1;i<=N;i++)
        {
            scanf("%lf%lf%lf",&a[i].x,&a[i].y,&a[i].z);
        }
        solve();
        return 0;
    }
    当初的梦想实现了吗,事到如今只好放弃吗~
  • 相关阅读:
    使用air16sdk打包ipa报错
    笔试题目
    网络编程Socket之UDP
    网络编程Socket之TCP
    Spring的HelloWorld
    一起学makefile
    使用Eclipse创建Hibernate工程
    log4j日志输出使用教程
    监听tomcat服务器启动/关闭并从配置文件中读取参数进行初始化
    java加载配置文件
  • 原文地址:https://www.cnblogs.com/caijiaming/p/11674182.html
Copyright © 2011-2022 走看看