poj2420 模拟退火

来源：互联网发布：金融数据编辑：程序博客网时间：2024/04/28 23:34

害怕爬山算法挂掉，所以改成了新的模拟退火

模拟退火就是在爬山算法的基础上又新增加了一个以一定概率接受最优解。

就是一个exp 比一个随机在(0,1)的概率小，就接受。

#include<cstdio>#include<iostream>#include<cstring>#include<cstdlib>#include<algorithm>#include<queue>#include<cmath>#include<ctime>using namespace std;const int maxn=105;int n;double xi[maxn],yi[maxn];double dis(double x,double y){double res=0;for(int i=1;i<=n;i++){res+=sqrt((xi[i]-x)*(xi[i]-x)+(yi[i]-y)*(yi[i]-y));}return res;}int main(){srand(time(NULL));while(scanf("%d",&n)!=EOF){for(int i=1;i<=n;i++)scanf("%lf%lf",&xi[i],&yi[i]);double step=100000;double x=xi[1],y=yi[1];double ans=dis(x,y);while(step>1e-10){double rad=rand()%360+1;double nx=x+step*cos(rad);double ny=y+step*sin(rad);double tmp=dis(nx,ny);if(tmp<ans){ans=tmp;x=nx;y=ny;}else{if(exp((tmp-ans)/step)<(rand()%10000)/10000.0){ans=tmp;x=nx;y=ny;}}step*=0.99;}printf("%.0f\n",ans);}return 0;}

0 0