[POJ2420]A Star not a Tree?（模拟退火）

2023-07-04 14:22:02

题目链接：http://poj.org/problem?id=2420

求费马点，即到所有其他点总和距离最小的点。

一开始想枚举一个坐标，另一个坐标二分的，但是check的时候还是O(n)的，复杂度相当于O(n^2lgn)，没意义。

学习一种神贪心，模拟退火。感觉和启发式搜索有点像啊，又有点像牛顿迭代。

　　思路就是，固定一个点和一个步长，从这个点开始向四个方向扩展，扩展的步长就是当前步长。如果扩展到的点可以更新答案，那么记住这个点，也就是说这个贪心方向（梯度？）是正确的。就可以拿着这个点继续沿着这个方向走了（剩余的方向可以不考虑了）。

　　如果四个方向都不能走，说明步长过长，这个时候模拟退火，将步长按比率缩小。

 1 #include 
 2 #include 
 3 #include 
 4 #include 
 5 #include 
 6 #include 
 7 #include 
 8 #include 
 9 #include 
10 #include 
11 #include 
12 #include 
13 #include 
14 #include 
15 #include <set>
16 #include 
17 #include 
18 using namespace std;
19 
20 typedef pair<double, double> pdd;
21 #define x first
22 #define y second
23 const double eps = 1e-8;
24 const double delta = 0.98;
25 const int maxn = 10100;
26 const int dx[5] = {1, -1, 0, 0};
27 const int dy[5] = {0, 0, 1, -1};
28 int n;
29 double ret;
30 pdd cur;
31 pdd p[maxn];
32 
33 double dist(pdd a, pdd b) {
34     double p = a.x - b.x;
35     double q = a.y - b.y;
36     return sqrt(p * p + q * q);
37 }
38 
39 int main() {
40     // freopen("in", "r", stdin);
41     while(~scanf("%d", &n)) {
42         for(int i = 0; i < n; i++) {
43             scanf("%lf %lf", &p[i].x, &p[i].y);
44         }
45         int t = 100;
46         ret = 1e100;
47         cur = pdd(.0, .0);
48         while(t > eps) {
49             bool flag = 1;
50             while(flag) {
51                 flag = 0;
52                 for(int i = 0; i < 4; i++) {
53                     pdd pos = pdd(cur.x+dx[i]*t, cur.y+dy[i]*t);
54                     double tmp = .0;
55                     for(int j = 0; j < n; j++) {
56                         tmp += dist(pos, p[j]);
57                     }
58                     if(ret - tmp > eps) {
59                         ret = tmp;
60                         cur = pos;
61                         flag = 1;
62                         break;
63                     }
64                 }
65             }
66             t *= delta;
67         }
68         printf("%.0f\n", ret);
69     }
70     return 0;
71 }

转载于:https://www.cnblogs.com/kirai/p/6229822.html

本文来自互联网用户投稿，文章观点仅代表作者本人，不代表本站立场，不承担相关法律责任。如若转载，请注明出处。 如若内容造成侵权/违法违规/事实不符，请点击【内容举报】进行投诉反馈！

标签：技术

上一篇 > 从统计局抓取2016年最新的全国区县数据！！
下一篇 > window 内核详尽分析

Duilib中list控件支持ctrl和shif多行选中的实现

[ICML2015]Batch Normalization:Accelerating Deep Network Training by Reducing Internal Covariate Shif

win10系统微软输入法于eclipse ctrl+shif+f冲突间接处理办法

Codeforces Round #259 (Div. 2) B. Little Pony and Sort by Shif

读LDD3，内存映射与DMA--PAGE_SHIF…

VMware虚拟机安装XP【要先分区，再设置BOOT 启动CD，shif+上移】

更换iBus五笔的左与右Shif

sublime ctrl+shif+f 没用解决办法

idea 对 ctrl + z 的撤销是 ctrl + shif + z

计算机最早的设计师应用于,计算机应用基础选择题doc.doc

win10自带截图神器：Win+Shift+S

Python基础之文件目录操作

python简述目录_Python基础之文件目录操作(示例代码)

tp5 如何做数据采集

任务2-7(服务器字体+阿里巴巴矢量库)

html标签（1)：h1~h6,p,br,pre,hr

TI 电量计介绍与芯片选型指南

几款TI电源芯片简介

TI DSP芯片C2000系列读取FLASH数据

德州仪器(Ti)平台嵌入式开发基础

TI三相电机智能栅极驱动芯片特点分类

省选模拟（12.08） T3 圈圈圈圈圈圈圈圈

Hadoop生态圈技术栈（上）

大数据开发基础入门与项目实战（三）Hadoop核心及生态圈技术栈之6.Impala交互式查询

小猿圈之Linux下Mysql 操作命令

大数据Hadoop生态圈常用面试题

大数据开发基础入门与项目实战（三）Hadoop核心及生态圈技术栈之4.Hive DDL、DQL和数据操作

备战Noip2018模拟赛11（B组）T3 Monogatari 物语

【智能优化算法-圆圈搜索算法】基于圆圈搜索算法Circle Search Algorithm求解单目标优化问题附matlab代码

NYOJ 78 圈水池

递归问题跑道汽车绕圈问题 Python实现

Hadoop生态圈（三）：MapReduce

[POJ2420]A Star not a Tree?（模拟退火）

相关文章