[POJ2420 A Star not a Tree?]

[关键字]:随机化搜索 模拟退火

[题目大意]:给出n个点,找出一个距离所有点之和最小的点,输出距离。

//===============================================================================

[分析]:同样是随机化的方法,只是在把没个点随机移动时有小小的变动。poj1379是使用得随机一个角度x然后利用t*cos(x)和t*sin(x)来处理随机走到的位置,而这道题因为没有了边界的限制,可以用以2t为长以要移动的点位中心的正方形中随机一个移动的位置来求最优值。至于具体退火算法步骤可以看poj1379或http://www.clarkok.com/blog/?p=357

//================================================================================

[代码]:

View Code
#include<iostream>
#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<algorithm>
#include<cmath>
using namespace std;

const int NUM=20;
const int TUM=20;
const int RD=1000;
const int MAXN=150;
const double ZERO=0.1;
const double det=0.7;

struct node
{
double x,y,dat;
}a[MAXN],MIN,MAX,b[MAXN];
int n;
double ans;

double dis(node c)
{
double sum=0;
for (int i=1;i<=n;++i)
sum+=sqrt((a[i].x-c.x)*(a[i].x-c.x)+(a[i].y-c.y)*(a[i].y-c.y));
return sum;
}

double random()
{
double sum=(rand()%(RD+1))*(1.0/RD);
return sum;
}

node GET(node x,node y)
{
node c;
//printf("%lf %lf %lf",x.x,y.x,y.x-x.x);
//system("pause");
c.x=x.x+(y.x-x.x)*random(),c.y=x.y+(y.y-x.y)*random();
//printf("%lf %lf\n",c.x,c.y);
//system("pause");
c.dat=dis(c);
return c;
}

bool Init()
{
if (scanf("%d",&n)==EOF) return 0;
MIN.x=MIN.y=100000;
MAX.x=MAX.y=-100000;
for (int i=1;i<=n;++i)
{
scanf("%lf%lf",&a[i].x,&a[i].y);
if (a[i].x<MIN.x) MIN.x=a[i].x;if (a[i].y<MIN.y) MIN.y=a[i].y;
if (a[i].x>MAX.x) MAX.x=a[i].x;if (a[i].y>MAX.y) MAX.y=a[i].y;
}
for (int i=1;i<=NUM;++i) b[i]=GET(MIN,MAX);//,printf("%lf %lf %lf\n",b[i].x,b[i].y,b[i].dat);
//printf("%lf %lf %lf %lf\n",MIN.x,MIN.y,MAX.x,MAX.y);
return 1;
}

void Solve(double T)
{
node t,x,y;
for (double i=T;i>=ZERO;i*=det)
for (int j=1;j<=NUM;++j)
for (int k=1;k<=TUM;++k)
{
x.x=b[j].x-i,x.y=b[j].y-i;
y.x=b[j].x+i,y.y=b[j].y+i;
t=GET(x,y);
//printf("%lf %lf %lf %lf %lf\n",x.x,x.y,y.x,y.y,i);
if (t.dat<b[j].dat) b[j]=t;
}
ans=1e10;
for (int i=1;i<=NUM;++i)
{
if (ans>b[i].dat) ans=b[i].dat;
//printf("%lf\n",b[i].dat);
}
printf("%.0lf\n",ans);
}

int main()
{
freopen("in.txt","r",stdin);
freopen("out.txt","w",stdout);
while (Init()) Solve(max(MAX.x-MIN.x,MAX.y-MIN.y));
return 0;
}



转载于:https://www.cnblogs.com/procedure2012/archive/2012/03/28/2420816.html

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/news/262827.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

使用file做swap分区

有时候在给系统盘分区时&#xff0c;忘记swap分区&#xff0c;这种情况下&#xff0c;可以用一个文件做swap&#xff0c;效果差不多的。这里就以2Gswap分区为例。步骤&#xff1a;1、创建2G大小的文件dd if/dev/zero of/swap bs1M count20482、格式化swap filemkswap /swap3、把…

react和vue配置本地代理

React 在react中配置开发环境下的本地代理相对比较简单&#xff0c;直接在package.json文件中修改即可。 但是这样做有其局限性&#xff0c;如果开发中代理多个接口的时候将无法满足需求&#xff0c;我们需要的是下面这种的能够代理多个接口的请求方式 proxy(/back, {target: h…

php异步处理,执行系统命令

为什么80%的码农都做不了架构师&#xff1f;>>> 难点1&#xff1a;php执行系统命令&#xff0c;要等到命令执行完成之后才会结束&#xff0c;不然会造成整个php进程挂起&#xff0c;其他php服务僵死&#xff1f; 难点2&#xff1a;php执行的命令完成后&#xff0c;…

axios拦截器的实现

拦截器设计与实现 #需求分析 我们希望能对请求的发送和响应做拦截&#xff0c;也就是在发送请求之前和接收到响应之后做一些额外逻辑。 我们希望设计的拦截器的使用方式如下&#xff1a; // 添加一个请求拦截器 axios.interceptors.request.use(function (config) { // 在发送请…

What can Lines3D Class do?

好几天没来Bolg了&#xff0c;这几天在尝试用3D线条做一个比较优雅的线条效果。文字由3DMAX导出到AS类&#xff0c;其他曲线路径由TweenMax的bezier数组定义。 Preview:niuniuzhu.cn 转载于:https://www.cnblogs.com/niuniuzhu/archive/2008/07/25/1251556.html

10 个最佳的网站分析方法

一个网页设计师在设计网站的各个关键方面时&#xff0c;他们需要了解网站的各种必要细节&#xff0c;可通过市场调研或者是网站本身的统计来获取这些细节信息。 本文介绍的这些工具确实可以帮助网站进行适当和有效的分析。这里面有些是免费的工具&#xff0c;还有的提供非常强大…

C#表达式树浅析

一、前言 在我们日常开发中Lamba 表达式经常会使用&#xff0c;如List.Where(n>Name"abc") 使用起来非常的方便&#xff0c;代码也很简洁&#xff0c;总之一个字就是“爽”。在之前我们总是用硬编码的方式去实现一些底层方法&#xff0c;比如我要查询用户“abc”是…

下载 sdk struts java

<action name"sdkDownload" class"com.curiousby.sdkDownload"><!-- result的Type必须为stream --><result name"success" type"stream"><param name"contentType">application/octet-stream;char…

微信小程序省市区联动,自定义地区字典

最近在做一个项目的时候遇到了这么一个问题&#xff0c;就是省市区的联动呢&#xff0c;我们需要自定义字典来设置&#xff0c;那么微信小程序自带的省市区选择就不能用了&#xff0c;经过三根烟的催化&#xff0c;终于写出来了。下面献上代码示例。 首先是在utils文件夹存入ar…

论文翻译《Object-Level Ranking: Bringing Order to Web Objects》

Object-Level Ranking: Bringing Order to Web Objects Zaiqing Nie Yuanzhi Zhang Jirong Wen Weiying Ma 摘要&#xff1a; 现在的网络搜索方法实际上是做文档级排名和检索&#xff0c;与之相对比&#xff0c;我们在探索一种新的聚合体以实现在对象级的网络检索。我们搜集与某…

前端vscode常用插件

Auto Rename Tag 这是一个html标签的插件&#xff0c;可以让你修改一边标签&#xff0c;另外一边自动改变。 Beautify 格式化代码插件 Braket Pair Colorizer 给js文件中的每一个小括号()花括号{}都配上不同的颜色&#xff0c;方便找到哪一个位置多了少了括号。 Debugger for C…

在线条形码生成器

条形码又称条码、一维码&#xff0c;是将字符按照特定的规则转化成二进制后&#xff0c;描绘成一个宽度不等的多个黑条和空白&#xff0c;按照一定的编码规则排列的图形标识符&#xff0c;条形码现在应用相当广泛&#xff0c;一出门&#xff0c;随便翻一样东西&#xff0c;可能…

[JSOI2008 Prefix火星人]

[关键字]&#xff1a;splay hash 二分 [题目大意]&#xff1a;给出一个字符串&#xff0c;求出给定的两个后缀的的最长公共前缀。在求的过程中会有改变或在某个位置添加字符的操作。 // [分析]&#xff1a;一听最长公共前缀马上想到后缀数组&#xff0c;但因为是动态维护所以后…

ios学习笔记block回调的应用(一个简单的例子)

一、什么是Blocks Block是一个C级别的语法以及运行时的一个特性&#xff0c;和标准C中的函数&#xff08;函数指针&#xff09;类似&#xff0c;但是其运行需要编译器和运行时支持&#xff0c;从ios4.0开始就很好的支持Block。 二、在ios开发中&#xff0c;什么情况下使用…

vue定义global.js,挂载在vue原型上面使用

首先在src目录下创建global目录&#xff0c;在global目录下创建index.js。 export default {install(Vue) {var that this// 1. 添加全局方法或属性// ue.global this// 2. 添加全局资源// 3. 注入组件Vue.mixin({created() {this.global that}})// 大于一的整数验证&#x…

《Windows游戏编程大师技巧》三、Windows高级编程

Windows编程很绝的地方在于&#xff1a;你不用了解太多细节&#xff0c;就可以完成很多工作。使用资源资源就是你的程序代码结合在一起的多块数据&#xff0c;可以被程序本身在运行时加载。资源应当也放在程序的.EXE文件中的原因是&#xff1a;1.同时包含代码和数据的.EXE文件更…

结构型模式--装饰模式

下面先用java&#xff0c;然后用Objective&#xff0d;C行对装饰模式的讲解&#xff1a; 对于java的装饰模式讲解和使用比较详细和难度有点偏高&#xff0c;而对于Objective&#xff0d;C的装饰模式讲解和使用方面比较简单&#xff0c;而且和java的装饰模式略有差异&#xff0c…

ArcGIS.Server.9.2.DotNet自带例子分析(三、一)

目的&#xff1a; 1.arcgis server9.2 ADF的AddGraphics。 准备工作&#xff1a; 1.用ArcGis Server Manager或者ArcCatalog发布一个叫world的Map Service,并且把这个Service启动起来。 2.找到DeveloperKit\SamplesNET\Server\Web_Applications目录下的Common_AddGraphicsCShar…

aspose将datatable导出excel 比自己拼好的多 Bug少-。.net

using System; using System.Collections.Generic; using System.Linq; using System.Web; using System.IO; using System.Data; using Aspose.Cells; /// <summary> ///OutFileDao 的摘要说明 /// </summary> publicclass OutFileDao { public OutFileDa…

整理:Android apk 框架 布局 集锦

2019独角兽企业重金招聘Python工程师标准>>> 看到好的技术教程就想分享一下&#xff0c;不喜勿喷&#xff01;谢谢配合&#xff0c;仅供菜鸟学习研究(^o^)/~ 友情推荐《爱加密》Android apk加密保护视频教程剪辑&#xff1a;http://www.ijiami.cn/Video?v3 Andro…