poj 1379 Run Away

题意就是寻找最安全的点，也就是找到一点使得他在题目给的区域之内与所有点最短距离最大，很容易想到这题可以用不断趋近的模拟退火来做，先随机一部分点集，对每个点分别点的扩展，选取较优点来替代原来的点，然后不断缩小步长，点的扩展时试了几种方法，发现还是使用正余弦比较靠谱。

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<algorithm>
#include<cmath>
using namespace std;
struct point
{
double x,y,dis;
}a[2000],b[100],ans;
int n,X,Y;

double dis(point c)			//计算最短距离
{
double sum=1e8;
for (int i=1;i<=n;++i)
if (sum>sqrt((c.x-a[i].x)*(c.x-a[i].x)+(c.y-a[i].y)*(c.y-a[i].y)))
sum=sqrt((c.x-a[i].x)*(c.x-a[i].x)+(c.y-a[i].y)*(c.y-a[i].y));
return sum;
}

void Solve()
{
for (int i=1;i<=30;++i) 		//随机点集
{
b[i].x=rand()%X+1,b[i].y=rand()%Y+1;
b[i].dis=dis(b[i]);
}
point t;
for (double i=max(X,Y);i>=1e-2;i*=0.9)		//步长缩短
for (int j=1;j<=30;j++)					//固定步长，对每个点进行扩展
for (int k=1;k<=30;k++)				//扩展的尝试
{
double ran=rand();
t.x=b[j].x+cos(ran)*i;
t.y=b[j].y+sin(ran)*i;
if (0>t.x || t.x>X || 0>t.y || t.y>Y)
continue;
t.dis=dis(t);
if (t.dis>b[j].dis) b[j]=t;
}
ans.dis=-1;
for (int i=1;i<=30;++i)
if (ans.dis<b[i].dis) ans=b[i];			//取距离最远的点
}

int main()
{
int test;
scanf("%d",&test);
while (test--)
{
scanf("%d %d %d",&X,&Y,&n);
for (int i=1;i<=n;i++)
scanf("%lf %lf",&a[i].x,&a[i].y);
Solve();
printf("The safest point is (%.1f, %.1f).\n",ans.x,ans.y);
}
return 0;
}