poj1379 模擬退火

【題意】

地圖中有N個陷阱，給出他們的座標，求一個點，使得這個點到所有陷阱的最小距離最大。

【題解】

乍一看是二分，但沒有好方法。

看到題目要求的精度爲0.1，則自然想到了模擬退火。

模擬退火的流程大致如下

1、隨即生成若干解

2、定下步長d，和降溫速度

3、開始模擬退火過程

【代碼】

#include <iostream>
#include <cmath>
#include <cstdlib>
#include <cstdio>
#include <ctime>
#include <cstring>
using namespace std;
const double pi=acos(-1.0);
const double INF=1e20;
const double eps=1e-3;

double px[31],py[31],d[31],x[1005],y[1005];

double dist(double x1,double y1,double x2,double y2)
{
    return sqrt((x1-x2)*(x1-x2)+(y1-y2)*(y1-y2));
}

int main()
{
    freopen("in.txt","r",stdin);
    int cc,X,Y,n,i,j,k;
    double theta,delta,tx,ty,td,dx,dy,ans,qx,qy;
    cin >> cc;
    srand((unsigned)time(NULL));
    while (cc--)
    {
        scanf("%d%d%d",&X,&Y,&n);
        for (i=1;i<=n;i++)
            scanf("%lf%lf",&x[i],&y[i]);
        for (i=1;i<=30;i++)
        {
            px[i]=double(rand()%1000+1)/1000.000*X;
            py[i]=double(rand()%1000+1)/1000.000*Y;
            d[i]=INF;
            for (j=1;j<=n;j++)
                d[i]=min(d[i],dist(px[i],py[i],x[j],y[j]));
        }
        delta=double(max(X,Y))/(sqrt(1.0*n));
        while (delta>eps)
        {
            for (i=1;i<=30;i++)
            {
                qx=px[i];qy=py[i];
                for (j=1;j<=30;j++)
                {
                    theta=double(rand()%1000+1)/1000.000*10*pi;
                    dx=delta*cos(theta);
                    dy=delta*sin(theta);
                    tx=qx+dx;ty=qy+dy;
                    if (tx<0 || tx>X || ty<0 || ty>Y)
                        continue;
                    td=INF;
                    for (k=1;k<=n;k++)
                        td=min(td,dist(tx,ty,x[k],y[k]));
                    if (td>d[i])
                    {
                        d[i]=td;px[i]=tx;py[i]=ty;
                    }
                }
            }
            delta*=0.8;
        }
        ans=0;
        for (i=1;i<=30;i++)
            if (d[i]>ans)
            {
                k=i;
                ans=d[i];
            }
        printf("The safest point is (%.1lf, %.1lf).\n",px[k],py[k]);
    }
    return 0;
}