POJ 2420 A Star not a Tree? (爬山法||模拟退火)

最新推荐文章于 2022-11-18 19:36:55 发布

原创最新推荐文章于 2022-11-18 19:36:55 发布 · 596 阅读

1 ·

CC 4.0 BY-SA版权

文章标签：

#模拟退火

模板同时被 2 个专栏收录

65 篇文章

订阅专栏

爬山法||模拟退火

1 篇文章

订阅专栏

本文介绍了一种结合爬山法与模拟退火算法的优化方法，并提供了一个具体的实现示例。该方法用于寻找一组给定点到某特定位置的距离之和最小化的解决方案。通过迭代调整候选解并接受一定概率的较差解来跳出局部最优。

爬山法+模拟退火:

#include<cstdio>
#include<cstring>
#include<iostream>
#include<queue>
#include<cmath>
#include<math.h>
#include<ctime>
#include<cstdlib>
using namespace std ;
struct node
{
    double x,y;
}q[120];
double xx,yy;
int n;
double dis(double x,double y,node p)
{
    return sqrt( (x-p.x)*(x-p.x)+(y-p.y)*(y-p.y) );
}
double getsum(double a,double b)
{
    double tmp=0;
    for(int i=0;i<n;i++)
        tmp+=dis(a,b,q[i]);
    return tmp;
}
int main()
{
    int m,i,j,k;
    while(~scanf("%d",&n))
    {
        double ans;
        xx=0,yy=0;double t=10000;
        for(i=0;i<n;i++)
        {
            scanf("%lf%lf",&q[i].x,&q[i].y);
            xx+=q[i].x;yy+=q[i].y;
        }
        xx/=n;yy/=n;
        ans=getsum(xx,yy);
        double tmp,x,y;
        while(t>0.02)
        {
            x=y=0;
            for(i=0;i<n;i++)
            {
                x+=(q[i].x-xx)/dis(xx,yy,q[i]);
                y+=(q[i].y-yy)/dis(xx,yy,q[i]);
            }
            tmp=getsum(xx+x*t,yy+y*t);
            if(tmp<ans)
            {
                ans=tmp;
                xx+=x*t;yy+=y*t;
            }
            /*else if(log( (tmp-ans)/t)<(rand()%10000/10000.0) )//注意使用的cstdlib/ctime头文件
            {
                ans=tmp;
                xx+=x*t;yy+=y*t;
            }*////模拟退火和爬山就是一个else if(log(tmp-ans )/t<(rand()%10000)/10000.0 )的区别
            t*=0.9;
        }
        printf("%.0lf\n",ans);
    }
    return 0;
}