A Star not a Tree? 模拟退火

136 篇文章 0 订阅
8 篇文章 0 订阅
/*模拟退火的思路非常巧妙。从大处定位,再慢慢细微。相当于先粗条再微调。*/
#include <stdio.h>
#include <cmath>
struct point
{
    double x,y;
}p[101];
int n;
double dis(point p1,point p2)
{
    return sqrt((p1.x - p2.x) * (p1.x - p2.x) + (p1.y - p2.y) * (p1.y - p2.y));
}
point getpoint(double a,double b)
{
    point t;
    t.x=a;
    t.y=b;
    return t;
}
double alldis(point a)
{
   double sum=0;
   for(int i=0;i<n;i++)
   sum+=dis(a,p[i]);
   return sum;
}
int main()
{
    while(scanf("%d",&n)==1)
    {
        for(int i=0;i<n;i++)
        scanf("%lf%lf",&p[i].x,&p[i].y);
        double step=100.0;
        point tt,qq,pp;
        pp=p[0];
        double ans=alldis(pp);
        while(step>0.2)
        {
            bool ret=true;
            while(ret)
            {
                ret=false;
                double tmp=0;
                qq=getpoint(pp.x,pp.y+step);//tt为原始点,pp为以tt为中心更新后的点
                tt=pp;
                tmp=alldis(qq);
                if(tmp<ans)
                {
                    ans=tmp;
                    ret=true;
                    tt=qq;
                }
                qq=getpoint(pp.x,pp.y-step);
                tmp=alldis(qq);
                if(tmp<ans)
                {
                    ans=tmp;
                    ret=true;
                    tt=qq;
                }
                qq=getpoint(pp.x+step,pp.y);
                tmp=alldis(qq);
                if(tmp<ans)
                {
                    ans=tmp;
                    ret=true;
                    tt=qq;
                }
                qq=getpoint(pp.x-step,pp.y);
                tmp=alldis(qq);
                if(tmp<ans)
                {
                    ans=tmp;
                    ret=true;
                    tt=qq;
                }
                pp=tt;
            }
            step=step/2.0;
        }
        printf("%d\n",int((ans+0.5)*100/100));
    }
    return 0;
}


                
  • 1
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值