/* 模擬退火法：找到一些隨機點，從這些點出發，隨機的方向坐標向外搜索；最后找到這些隨機點的最大值；坑：//if(xx>-eps&&xx<x+eps&&yy>-eps&&yy<eps+y) {不知道為什么這個判斷方式錯誤？？？？？？ */ #include<iostream> #include<cstdio> #include<math.h> #include<algorithm> #define pi acos(-1.0) #define N 1100 #define inf 1000000000000 #define eps 1e-8 using namespace std; double x,y; int n; struct node { double u,v,dis; }f[N],endd,ff[N]; double diss(double x,double y,int i) {return sqrt((x-f[i].u)*(x-f[i].u)+(y-f[i].v)*(y-f[i].v)); } double Mindis(double x,double y) {double minn=1.0*inf,ss;int i;for(i=1;i<=n;i++) {ss=diss(x,y,i);if(minn>ss)minn=ss;}return minn; } void compute() {int i,j;double xx,yy;for(i=1;i<=20;i++) {//先隨機出來一些點ff[i].u=(rand()%1000+1)/1000.0*x;ff[i].v=(rand()%1000+1)/1000.0*y;ff[i].dis=Mindis(ff[i].u,ff[i].v);}double T=sqrt(1.0*x*x+1.0*y*y);double rate=0.9;while(T>eps) {for(i=1;i<=20;i++)//對于這些點分別向外搜for(j=1;j<=30;j++) {//隨機30個半徑來搜索，更新，得到每個點可以到達的最遠距離double ran=(rand()%1000+1)/1000.0*pi*10;xx=ff[i].u+cos(ran)*T;yy=ff[i].v+sin(ran)*T;//if(xx>-eps&&xx<x+eps&&yy>-eps&&yy<eps+y) {不知道為什么這個判斷方式錯誤？？？？？？if(xx<0.0||xx>x||yy<0.0||yy>y)continue;double di=Mindis(xx,yy);if(ff[i].dis<di) {ff[i].dis=di;ff[i].u=xx;ff[i].v=yy;}}T*=rate;//退火率}for(i=1;i<=20;i++)//找到最大的一個if(endd.dis<ff[i].dis)endd=ff[i];return ; } int main() {// printf("%d\n",(rand()%1000+1)/1000.0*2*pi);int m,i,j,k,t;scanf("%d",&t);while(t--) {scanf("%lf%lf%d",&x,&y,&n);for(i=1;i<=n;i++)scanf("%lf%lf",&f[i].u,&f[i].v);endd.dis=-1;compute();printf("The safest point is (%.1f, %.1f).\n",endd.u,endd.v);} return 0; }

轉載于:https://www.cnblogs.com/thefirstfeeling/p/4410615.html

總結

以上是生活随笔為你收集整理的poj 1379 模拟退火法的全部內容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網站內容還不錯，歡迎將生活随笔推薦給好友。