利用分治来求平面最近点对

只需要查后面6个点就好了

原因在于https://blog.csdn.net/liufeng_king/article/details/8484284

两个集合的话就把不同集合的设为inf就好了

 1 const db inf = 1e100;
 2
 3 ll sqr(ll x) { return x*x; }
 4
 5 struct point {
 6     int x,y;
 7     int f;
 8     void in() {
 9         scanf("%d%d",&x,&y);
10     }
11     point operator-(const point &b) const {
12         point ret;
13         ret.x=x-b.x;
14         ret.y=y-b.y;
15         return ret;
16     }
17     db mean() {
18         ll ret=0;
19         ret=sqr(x)+sqr(y);
20         return sqrt(ret);
21     }
22 };
23
24 point a[N];
25
26 db dis(point a,point b) {
27     db ret=0;
28     if(a.f==b.f) return inf;
29     ret=(a-b).mean();
30     return ret;
31 }
32
33 bool cmp_x(point a,point b) {
34     return a.x<b.x;
35 }
36
37 int t[N];
38 bool cmp_t(int x,int y) {
39     return a[t[x]].y<a[t[y]].y;
40 }
41
42 db F(int l,int r) {
43     if(l==r) return inf;
44     if(l+1==r) return dis(a[r],a[l]);
45     int mid=(l+r)>>1;
46     db ans1=F(l,mid);
47     db ans2=F(mid+1,r);
48     db ans=min(ans1,ans2);
49     if(ans==0) return ans;
50     int nm=0;
51     for(int i=l;i<=r;i++) if(a[mid].x-ans<=a[i].x && a[i].x<=a[mid].x+ans) {
52         t[++nm]=i;
53     }
54     sort(t+1,t+nm+1,cmp_t);
55     for(int i=1;i<=nm;i++) {
56         for(int j=i+1;j<=min(nm,i+6);j++) {
57             db anss=dis(a[t[i]],a[t[j]]);
58             ans=min(ans,anss);
59         }
60     }
61     return ans;
62 }
63
64 int main() {
65     int T;
66     scanf("%d",&T);
67     while(T--) {
68         int n;
69         scanf("%d",&n);
70         for(int i=1;i<=n;i++) {
71             a[i].in();
72             a[i].f=1;
73         }
74         for(int i=n+1;i<=2*n;i++) {
75             a[i].in();
76             a[i].f=2;
77         }
78         sort(a+1,a+2*n+1,cmp_x);
79         db ans=F(1,2*n);
80         printf("%.3f\n",ans);
81     }
82     return 0;
83 }

View Code

转载于:https://www.cnblogs.com/MyGirlfriends/p/9317655.html

POJ3714 Raid 平面最近点对相关推荐

  1. POJ - 3714 Raid(平面最近点对模板题,几何)

    题目链接:点击查看 题目大意:给出两个含有n个点的集合,在两个集合中分别任选一点,使得这两个点之间的距离最小 题目分析:因为n给到了1e5,所以n*n的暴力肯定是不行了,直接从网上copy了个分治优化 ...

  2. CF429D Tricky Function(求解公式、经分析转为求平面最近点对、思维)

    整理的算法模板合集: ACM模板 目录 CF429D Tricky Function 题意实际上就是给定长度为 nnn 的一串序列a1,a2,...,ana_1, a_2,...,a_na1​,a2​ ...

  3. HDU-4631 Sad Love Story 平面最近点对

    题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=4631 数据是随机的,没有极端数据,所以可以分段考虑,最小值是一个单调不增的函数,然后每次分治算平面最近 ...

  4. POJ3714 Raid 分治/K-D Tree

    VJ传送门 简要题意:给出两个大小均为\(N\)的点集\(A,B\),试在\(A\)中选择一个点,在\(B\)中选择一个点,使得它们在所有可能的选择方案中欧几里得距离最小,求出这个距离 下面给出的两种 ...

  5. 平面最近点对问题(分治)

    题目描述 在与联盟的战斗中屡战屡败后,帝国撤退到了最后一个据点. 依靠其强大的防御系统,帝国击退了联盟的六波猛烈进攻. 经过几天的苦思冥想,联盟将军亚瑟终于注意到帝国防御系统唯一的弱点就是能源供应. ...

  6. 平面最近点对问题求解—基于Java语言

    平面最近点对问题求解-基于Java语言 1. 问题描述: 本问题来自<编程之美2.11-寻找最近点对>,文中给出了两种解法:暴力解法,分治解法.其中,暴力解法很简单,求出所有点之间的距离并 ...

  7. 分治法求平面最近点对入门

    一.平面最近点对问题. 平面最近点对:在一个平面上有 n n n个点,求出距离最近的两个点. 平面最近点对是计算几何中一个十分经典且基础的问题,通常采用分治法来解决. 二.直线最近点对的分治法. 在用 ...

  8. HDU 5721 Palace(平面最近点对(分治))

    http://acm.hdu.edu.cn/showproblem.php?pid=5721 n个点,去掉一个点的情况下,最近距离平方之和. 平面最近点对模版题,先求出最近距离,然后找到是哪两个点,然 ...

  9. 【算法分析与设计】平面最近点对(含最近距离、最近点对、第一次分割点集合的输出)

    [问题描述] 给定二维平面上n个点,找其中的一对点,使得在n个点组成的所有点对中,该点对间的距离最小.使用递归与分治策略求解二维平面上的最接近点对问题.假设所有点的集合为S,m为S中所有点的x坐标的中 ...

最新文章

  1. 焦李成教授谈深度神经网络发展历程
  2. C# 多线程读取数据并保存在DataGridView中,如果操作?代码是什么?谢谢
  3. 搜推广遇上用户画像:Lookalike相似人群拓展算法
  4. PostgreSQL 10.1 手册_部分 III. 服务器管理_第 16 章 从源代码安装_16.5. 安装后设置...
  5. squirrelsql安装
  6. 算法基础:递归算法知识笔记
  7. Kafka是什么、主要应用在什么场景?
  8. git 应用 rebase
  9. 的概念产生于计算机芯片,是指采用光导原位合成或微量点样等方法,,微阵列芯片...
  10. 机器学习——共享单车数据集预测
  11. SpringCloud微服务面试题
  12. 51单片机智能远程遥控温控PWM电风扇系统红外遥控温度速度定时关机
  13. WAS下Sanp、heapdump、javacore
  14. [工具]更新音乐下载软件工具音乐下载网站,MP3音乐无损音乐下载器
  15. 计算机英语课外知识竞赛,高2012级英语课外阅读知识竞赛.docx
  16. IDEATerminate vs Disconnect
  17. android10的coloros,OPPOReno10倍变焦版ColorOSV7刷机包(官方最新固件升级包安卓10ColorOS V7)...
  18. 您访问的网页出错了!
  19. 视频质量分析工具介绍
  20. linux 系统黑洞简单使用案例

热门文章

  1. vissim跟驰模型_VISSIM是什么工具?你对VISSIM了解多少?
  2. python写网站 linux_linux下利用nginx部署python网站
  3. java 方法特点是什么_Java String 的特点是什么?它有哪些重要的方法?
  4. 关于长江的题目_长江流域综合治理练习题
  5. 与word2vec_NLP--Word2Vec详解
  6. java并发编程详解,Java架构师成长路线
  7. 三面蚂蚁金服(交叉面)定级阿里P6
  8. mysql1045错误如何解决,Java面试真题精选
  9. 别再说自己不会了!最新高频Java笔试题分享
  10. 解决mac找不到隐藏的local文件夹