题目链接

POJ输出不能用%lf!

mmp从4:30改到6:00,把4:30交的一改输出也过了。

于是就有了两份代码。。

//392K	500MS
//用两点构成的矩形更新,就不需要管边界了
#include <cmath>
#include <cstdio>
#include <cctype>
#include <algorithm>
#define gc() getchar()
#define Rate (0.9)
#define Rand() (1.0*((rand()<<1)-RAND_MAX)/RAND_MAX)
const int N=1005; int n,X,Y;
double DIS,ans_dis[20];
struct Node{
double x,y;
Node() {}
Node(double x,double y):x(x),y(y) {}
}p[N],ans[30]; inline int read()
{
int now=0;register char c=gc();
for(;!isdigit(c);c=gc());
for(;isdigit(c);now=now*10+c-'0',c=gc());
return now;
}
inline double Sqr(double x){
return x*x;
}
inline double Dis(double x,double y,int i){
return sqrt(Sqr(x-p[i].x)+Sqr(y-p[i].y));
}
double Calc(double x,double y)
{
double res=DIS;
for(int i=1; i<=n; ++i)
res=std::min(res,Dis(x,y,i));
return res;
}
void Print(){
for(int i=1; i<=15; ++i) printf("%d:%.5lf %.5lf %.5lf %d %d\n",i,ans[i].x,ans[i].y,ans_dis[i],X,Y);
}
void Solve()
{
DIS=sqrt(X*X+Y*Y);
const int tot=15;
for(int i=1; i<=tot; ++i)
ans[i].x=1.0*rand()/RAND_MAX*X,ans[i].y=1.0*rand()/RAND_MAX*Y,ans_dis[i]=Calc(ans[i].x,ans[i].y);
double T=DIS,xx,yy,dis;//this...
while(T>0.01)
{
for(int i=1; i<=tot; ++i)
for(int j=1; j<=30; ++j)
{
xx=ans[i].x+T*Rand(), yy=ans[i].y+T*Rand();
if(xx<0||xx>X||yy<0||yy>Y) continue;//!
dis=Calc(xx,yy);
if(dis>ans_dis[i]) ans[i]=Node(xx,yy),ans_dis[i]=dis;
}
T*=Rate;
}
int res=1;
for(int i=2; i<=tot; ++i)
if(ans_dis[i]>ans_dis[res]) res=i;
printf("The safest point is (%.1f, %.1f).\n",ans[res].x,ans[res].y);
} int main()
{
srand(20180428);
int T=read();
while(T--)
{
X=read(),Y=read(),n=read();
for(int i=1; i<=n; ++i) p[i].x=read(),p[i].y=read();
Solve();
}
return 0;
}

另一种写法:

//392K	1188MS
#include <cmath>
#include <cstdio>
#include <cctype>
#include <algorithm>
#define gc() getchar()
#define Rate (0.97)
//#define Rand() (1.0*((rand()<<1)-RAND_MAX)/RAND_MAX)
#define Rand() (1.0*rand()/RAND_MAX) const int N=1005; int n,X,Y;
double DIS,ans_dis[30];
struct Node{
double x,y;
Node() {}
Node(double x,double y):x(x),y(y) {}
}p[N],ans[30]; inline int read()
{
int now=0;register char c=gc();
for(;!isdigit(c);c=gc());
for(;isdigit(c);now=now*10+c-'0',c=gc());
return now;
}
inline double Sqr(double x){
return x*x;
}
inline double Dis(double x,double y,int i){
return sqrt(Sqr(x-p[i].x)+Sqr(y-p[i].y));
}
double Calc(double x,double y)
{
double res=1e9;//DIS;
for(int i=1; i<=n; ++i)
res=std::min(res,Dis(x,y,i));
return res;
}
void Rand_Point(double &x,double &y,Node a,Node b){
x=a.x+(b.x-a.x)*Rand(), y=a.y+(b.y-a.y)*Rand();
}
void Solve()
{
DIS=sqrt(X*X+Y*Y);
const int tot=15;
ans[1]=Node(0,0), ans[2]=Node(X,Y), ans[3]=Node(X,0), ans[4]=Node(0,Y);//
for(int i=5; i<=tot; ++i)
Rand_Point(ans[i].x,ans[i].y,ans[1],ans[2]);
// ans[i].x=Rand()*X,ans[i].y=Rand()*Y;
for(int i=1; i<=tot; ++i) ans_dis[i]=Calc(ans[i].x,ans[i].y);
double T=std::min(X,Y),xx,yy,dis;//this...
while(T>0.01)
{
for(int i=1; i<=tot; ++i)
for(int j=1; j<=20; ++j)
{
Rand_Point(xx,yy,Node(std::max(0.0,ans[i].x-T),std::max(0.0,ans[i].y-T)),Node(std::min(1.0*X,ans[i].x+T),std::min(1.0*Y,ans[i].y+T)));
dis=Calc(xx,yy);
if(dis>ans_dis[i]) ans[i]=Node(xx,yy),ans_dis[i]=dis;
}
T*=Rate;
}
int res=1;
for(int i=2; i<=tot; ++i)
if(ans_dis[i]>ans_dis[res]) res=i;
printf("The safest point is (%.1f, %.1f).\n",ans[res].x,ans[res].y);
} int main()
{
#ifndef ONLINE_JUDGE
freopen("1.in","r",stdin);
// freopen("1.out","w",stdout);
#endif srand(20180428);
int T=read();
while(T--)
{
X=read(),Y=read(),n=read();
for(int i=1; i<=n; ++i) p[i].x=read(),p[i].y=read();
Solve();
}
return 0;
}

POJ.1379.Run Away(模拟退火)的更多相关文章

  1. poj 1379 Run Away 模拟退火 难度:1

    Run Away Time Limit: 3000MS   Memory Limit: 65536K Total Submissions: 6482   Accepted: 1993 Descript ...

  2. POJ 1379 Run Away 【基础模拟退火】

    题意:找出一点,距离所有所有点的最短距离最大 二维平面内模拟退火即可,同样这题用最小圆覆盖也是可以的. Source Code: //#pragma comment(linker, "/ST ...

  3. PKU 1379 Run Away(模拟退火算法)

    题目大意:原题链接 给出指定的区域,以及平面内的点集,求出一个该区域内一个点的坐标到点集中所有点的最小距离最大. 解题思路:一开始想到用随机化算法解决,但是不知道如何实现.最后看了题解才知道原来是要用 ...

  4. POJ 1379 Run Away

    题意:有n个陷阱,在X,Y范围内要求出一个点使得这个点到陷阱的最小距离最大. 思路:模拟退火,随机撒入40个点,然后模拟退火随机化移动. (这题poj坑爹,加了srand(time(NULL))不能交 ...

  5. POJ 1379 模拟退火

    模拟退火算法,很久之前就写过一篇文章了.双倍经验题(POJ 2420) 题意: 在一个矩形区域内,求一个点的距离到所有点的距离最短的那个,最大. 这个题意,很像二分定义,但是毫无思路,也不能暴力枚举, ...

  6. POJ 1379 (随机算法)模拟退火

    题目大意: 给定一堆点,找到一个点的位置使这个点到所有点中的最小距离最大 这里数据范围很小,精度要求也不高,我们这里可以利用模拟退火的方法,随机找到下一个点,如果下一个点比当前点优秀就更新当前点 参考 ...

  7. poj-1379 Run Away(模拟退火算法)

    题目链接: Run Away Time Limit: 3000MS   Memory Limit: 65536K Total Submissions: 7982   Accepted: 2391 De ...

  8. 【BZOJ1844/2210】Pku1379 Run Away 模拟退火

    [BZOJ1844/2210]Pku1379 Run Away 题意:矩形区域中有一堆点,求矩形中一个位置使得它到所有点的距离的最小值最大. 题解:模拟退火的裸题,再调调调调调参就行了~ #inclu ...

  9. poj 2069 Super Star —— 模拟退火

    题目:http://poj.org/problem?id=2069 仍是随机地模拟退火,然而却WA了: 看看网上的题解,都是另一种做法——向距离最远的点靠近: 于是也改成那样,竟然真的A了...感觉这 ...

随机推荐

  1. JVM调优命令-jps

    JVM Process Status Tool,显示指定系统内所有的HotSpot虚拟机进程. 命令格式 1 jps [options] [hostid] options参数-l : 输出主类全名或j ...

  2. [转]Android ANR 分析解决方法

    一:什么是ANR ANR:Application Not Responding,即应用无响应 二:ANR的类型 ANR一般有三种类型: 1. KeyDispatchTimeout(5 seconds) ...

  3. CSS 编码中超级有用的工具集合

    当你开发网站和 Web 应用时,使用合适的工具可以节省大量的时间.本文我将收集一些非常有用的 CSS 编码工具,希望对你有帮助. Pure CSS Pure 是来自雅虎的 CSS 框架,使用 Norm ...

  4. 利用iis虚拟目录实现文件服务器功能

    要求说明: 通过网站上传文件保存到统一的文件服务器上. 服务器说明: 1.文件服务器以下称为FilesServer,IP地址为:192.168.1.213 2.Web服务器为以下称为WebServer ...

  5. linux的内存文件系统tmpfs

    在centos系统上自带的内存文件系统.这个tmpfs是temporary file system的意思. 一. 使用命令 df -h 查看tmpfs是否正在运行. Filesystem Size U ...

  6. SANS社区帐号邮件激活问题

    注册时,密码需要数字,大写字母,小写字母,符号10位以上才能注册成功    吐槽:谁来爆破一下这种强度的密码,哈哈. 在我的文章中,有 计算机取证 分类,里面的一篇文章 Virtual Worksta ...

  7. MVC常用特性使用

    简介 在以前的文章中,我和大家讨论如何用SingalR和数据库通知来完成一个消息监控应用. 在上一篇文章中,我介绍了如何在MVC中对MongoDB进行CRUD操作. 今天,我将继续介绍一些在开发中非常 ...

  8. 浅介HTML DOM

    什么是DOM? DOM是Document Object Model(文档对象模型)的缩写. DOM是W3C(万维网联盟)的标准. DOM定义了访问HTML和XML文档的标准: “W3C文档对象模型(D ...

  9. testng运行失败,继续执行

    1.重写断言类 public class Verify { public static StringBuffer verificationErrors= new StringBuffer();; pu ...

  10. JavaScript之this学习心得

    this在运行时绑定,并不是在编写时绑定,它的上下文取决于函数调用的各种条件.this既不指向自身,也不指向函数的词法作用域.this是在函数被调用时发生的绑定,指向什么完全取决于函数在哪里被调用. ...