HDU 3932 Groundhog Build Home （模拟退火算法）

来源：互联网发布：红色高棉知乎编辑：程序博客网时间：2024/06/05 02:31

题意：找出一个点使得这个点到n个点的最长距离最短。

分析：模拟退火算法。

#include <iostream>#include<cstdio>#include<ctime>#include<cmath>#include<algorithm>#include<cstdlib>using namespace std;#define N 1005#define NUM 20#define D 50#define EPS 1e-3double a[N],b[N];int mo[8][2]={1,0,0,1,-1,0,0,-1,1,-1,-1,1,-1,-1,1,1};int main(){    int x,y,m;    srand((unsigned int)time(NULL));    while(~scanf("%d%d%d",&x,&y,&m))    {        for(int i=0;i<m;i++)        {            scanf("%lf%lf",&a[i],&b[i]);        }        double T=sqrt((double)x*x+y*y)+10;         double  ansx=(double)(rand()%1000+1)/1000.000*x;        double ansy=(double)(rand()%1000+1)/1000.000*y;        double ansdis=0;        for(int i=0;i<m;i++)            {             ansdis=max(ansdis,sqrt((ansx-a[i])*(ansx-a[i])+(ansy-b[i])*(ansy-b[i])));            }        while(T>EPS)        {            for(int k=0;k<D;k++){            for(int j=0;j<8;j++){            double tempx=ansx+(double)(rand()%1000+1)/1000.000*T*mo[j][0];            double tempy=ansy+(double)(rand()%1000+1)/1000.000*T*mo[j][1];            if(tempx<0||tempx>x||tempy<0||tempy>y)                continue;            double tempdis=0;            for(int i=0;i<m;i++)            {               tempdis=max(tempdis,sqrt((tempx-a[i])*(tempx-a[i])+(tempy-b[i])*(tempy-b[i])));            }            if(tempdis<ansdis)            {                ansx=tempx;                ansy=tempy;                ansdis=tempdis;            }            }            }            T=T*0.9;        }        printf("(%.1f,%.1f).\n%.1f\n",ansx,ansy,ansdis);    }    return 0;}

0 0