草庐IT

「分治法」最近点对距离问题

安和橋北 2023-06-17 原文

一、问题描述

给定平面S上n个点,找其中的一对点,使得在n个点组成的 所有点对中,该点对间的距离最小。

假设所讨论的点是以标准笛卡儿坐标形式(x, y)给出的。因 此,在两个点Pi =(xi , yi )和Pj =(xj , yj )之间的距离是标准的欧几里德 距离: d = 根号下( x i − x j ) 2 + ( y i − y j )2

二、问题分析

直接用暴力解法很简单,用结构体把每个点的x、y坐标都存下来之后,用一个二重循环即可获得距离的最大值。但是这样做的时间复杂度也达到了O(n的平方)

那么有没有什么办法,可以把时间复杂度降下来,达到nlogn呢?

思路

首先,利用分治法,将集合S分成两个相等的子集S1、S2。然后在这两个子集中递归的求出距离最近的点对。

这样做,如果最近点对都是分别在两边子集中,那还好说。如果不是分别在两边子集,而是这两个点一个在左边,一个在右边怎么办?这就不是单纯的递归了,子问题并不互斥,存在了overlap的问题。

所以我们把问题划分成两个子集S1、S2之后,出现三种情况

  1. 最近点对的两个点都在S1中
  2. 最近点对的两个点都在S2中
  3. 最近点对的两个点一个在S1中,一个在S2中

情况1和2可以递归求解,但是3不行。那为了解决情况3,我们应该怎么做呢?

精髓

首先我们假设点集S已经按照x坐标升序排列,并且点集S中x坐标的中值就是中间点x的坐标。

我们还是在左边点集S1中和右边点集S2中递归求解,假设S1和S2的最短点对距离分别是d1和d2,然后令d = min { d1, d2 }。

如果说,最后求解出来的S的最短点对距离是小于d的,也就是答案还不在左右两个子集中产生,那么说明这两个点肯定一个在左边点集S1中,一个在右边点集S2中。

这个可以用反证法来证明:

如果不是左右各一个,而是都在一边的话,那么由于这个距离是小于d的,而d又是min { d1, d2 },所以直接矛盾了,左边点集或者右边点集的最小点对距离肯定不是d1或者d2了。因此得证。

现在我们可以将求解限制在以点集S中线x = m为中心,宽度为2d的垂直地带P1和P2中,垂直带之外的任何点对的距离肯定都大于d。

假设点p(x,y)是集合P1和P2中y坐标最小的点,那么实际上点p可能在P1也可能在P2,但这并不重要。我们现在的目的就是找出在这个垂直地带中,有没有点对的距离比d还小的了。如果存在,那么这个点对就是答案;如果不存在,那么答案就是d。

那么以点p(x,y)为基准,那么我们要找的这个点,肯定是在垂直域[ y , y + d ]这个范围里面。而且在这个范围里面,点的个数不会超过8个!

这又是为什么呢?课堂上老师只是简单的说这可以用鸽笼原理来证明,让同学们自己下去看,接下来我简略解释一下:

看这个正方形,假设这个是中线x = m右边的这半矩形。为了保证右半边内的点对的距离不小于d(因为如果超过了,那S2点集的最短点对距离就不是d了!),只能把四个点放在四个角。

这时但凡再出现第五个点,那么这个第五个点的距离和另外一个点的距离,肯定是小于d的。也就证明一边最多只有四个点。所以总的来说,在这个矩形内,最多只有8个点。

所以,到目前为止,我们就限定了情况3的可能答案的范围。接下来就可以进行代码书写了。

算法时间复杂度分析

那么在这种情况下,时间复杂度是多少呢?

  1. 当点数n为1,2,3的时候,可以不用划分,直接求解即可,这时是O(1)
  2. 当点数大于3的时候,就需要划分子问题,进行递归求解。
    递归方程为T(n) = 2T(n/2) + O(n)
    为什么呢?因为划分和合并虽然用的是一个二重循环,但是实际上内层循环的最大值也就是8,所以相当于是O(an),也就是O(n)。

那么求解递归方程之后即可得出时间复杂度是nlogn。

三、代码

由于本题是书上内容,没有题目和测试用例,所以此处直接附上书上的代码。书上的代码笔者认为写的还是很清晰明了的。

有关「分治法」最近点对距离问题的更多相关文章

  1. ruby - 在 64 位 Snow Leopard 上使用 rvm、postgres 9.0、ruby 1.9.2-p136 安装 pg gem 时出现问题 - 2

    我想为Heroku构建一个Rails3应用程序。他们使用Postgres作为他们的数据库,所以我通过MacPorts安装了postgres9.0。现在我需要一个postgresgem并且共识是出于性能原因你想要pggem。但是我对我得到的错误感到非常困惑当我尝试在rvm下通过geminstall安装pg时。我已经非常明确地指定了所有postgres目录的位置可以找到但仍然无法完成安装:$envARCHFLAGS='-archx86_64'geminstallpg--\--with-pg-config=/opt/local/var/db/postgresql90/defaultdb/po

  2. ruby - 通过 rvm 升级 ruby​​gems 的问题 - 2

    尝试通过RVM将RubyGems升级到版本1.8.10并出现此错误:$rvmrubygemslatestRemovingoldRubygemsfiles...Installingrubygems-1.8.10forruby-1.9.2-p180...ERROR:Errorrunning'GEM_PATH="/Users/foo/.rvm/gems/ruby-1.9.2-p180:/Users/foo/.rvm/gems/ruby-1.9.2-p180@global:/Users/foo/.rvm/gems/ruby-1.9.2-p180:/Users/foo/.rvm/gems/rub

  3. ruby - 通过 RVM (OSX Mountain Lion) 安装 Ruby 2.0.0-p247 时遇到问题 - 2

    我的最终目标是安装当前版本的RubyonRails。我在OSXMountainLion上运行。到目前为止,这是我的过程:已安装的RVM$\curl-Lhttps://get.rvm.io|bash-sstable检查已知(我假设已批准)安装$rvmlistknown我看到当前的稳定版本可用[ruby-]2.0.0[-p247]输入命令安装$rvminstall2.0.0-p247注意:我也试过这些安装命令$rvminstallruby-2.0.0-p247$rvminstallruby=2.0.0-p247我很快就无处可去了。结果:$rvminstall2.0.0-p247Search

  4. ruby - Fast-stemmer 安装问题 - 2

    由于fast-stemmer的问题,我很难安装我想要的任何ruby​​gem。我把我得到的错误放在下面。Buildingnativeextensions.Thiscouldtakeawhile...ERROR:Errorinstallingfast-stemmer:ERROR:Failedtobuildgemnativeextension./System/Library/Frameworks/Ruby.framework/Versions/2.0/usr/bin/rubyextconf.rbcreatingMakefilemake"DESTDIR="cleanmake"DESTDIR=

  5. ruby - 安装 Ruby 时遇到问题(无法下载资源 "readline--patch") - 2

    当我尝试安装Ruby时遇到此错误。我试过查看this和this但无济于事➜~brewinstallrubyWarning:YouareusingOSX10.12.Wedonotprovidesupportforthispre-releaseversion.Youmayencounterbuildfailuresorotherbreakages.Pleasecreatepull-requestsinsteadoffilingissues.==>Installingdependenciesforruby:readline,libyaml,makedepend==>Installingrub

  6. java - 从 JRuby 调用 Java 类的问题 - 2

    我正在尝试使用boilerpipe来自JRuby。我看过guide从JRuby调用Java,并成功地将它与另一个Java包一起使用,但无法弄清楚为什么同样的东西不能用于boilerpipe。我正在尝试基本上从JRuby中执行与此Java等效的操作:URLurl=newURL("http://www.example.com/some-location/index.html");Stringtext=ArticleExtractor.INSTANCE.getText(url);在JRuby中试过这个:require'java'url=java.net.URL.new("http://www

  7. ruby-on-rails - 简单的 Ruby on Rails 问题——如何将评论附加到用户和文章? - 2

    我意识到这可能是一个非常基本的问题,但我现在已经花了几天时间回过头来解决这个问题,但出于某种原因,Google就是没有帮助我。(我认为部分问题在于我是一个初学者,我不知道该问什么......)我也看过O'Reilly的RubyCookbook和RailsAPI,但我仍然停留在这个问题上.我找到了一些关于多态关系的信息,但它似乎不是我需要的(尽管如果我错了请告诉我)。我正在尝试调整MichaelHartl'stutorial创建一个包含用户、文章和评论的博客应用程序(不使用脚手架)。我希望评论既属于用户又属于文章。我的主要问题是:我不知道如何将当前文章的ID放入评论Controller。

  8. 【高数】用拉格朗日中值定理解决极限问题 - 2

    首先回顾一下拉格朗日定理的内容:函数f(x)是在闭区间[a,b]上连续、开区间(a,b)上可导的函数,那么至少存在一个,使得:通过这个表达式我们可以知道,f(x)是函数的主体,a和b可以看作是主体函数f(x)中所取的两个值。那么可以有,  也就意味着我们可以用来替换 这种替换可以用在求某些多项式差的极限中。方法: 外层函数f(x)是一致的,并且h(x)和g(x)是等价无穷小。此时,利用拉格朗日定理,将原式替换为 ,再进行求解,往往会省去复合函数求极限的很多麻烦。使用要注意:1.要先找到主体函数f(x),即外层函数必须相同。2.f(x)找到后,复合部分是等价无穷小。3.要满足作差的形式。如果是加

  9. SPI接收数据异常问题总结 - 2

    SPI接收数据左移一位问题目录SPI接收数据左移一位问题一、问题描述二、问题分析三、探究原理四、经验总结最近在工作在学习调试SPI的过程中遇到一个问题——接收数据整体向左移了一位(1bit)。SPI数据收发是数据交换,因此接收数据时从第二个字节开始才是有效数据,也就是数据整体向右移一个字节(1byte)。请教前辈之后也没有得到解决,通过在网上查阅前人经验终于解决问题,所以写一个避坑经验总结。实际背景:MCU与一款芯片使用spi通信,MCU作为主机,芯片作为从机。这款芯片采用的是它规定的六线SPI,多了两根线:RDY和INT,这样从机就可以主动请求主机给主机发送数据了。一、问题描述根据从机芯片手

  10. git使用常见问题(提交代码,合并冲突) - 2

    文章目录git常用命令(简介,详细参数往下看)Git提交代码步骤gitpullgitstatusgitaddgitcommitgitpushgit代码冲突合并问题方法一:放弃本地代码方法二:合并代码常用命令以及详细参数gitadd将文件添加到仓库:gitdiff比较文件异同gitlog查看历史记录gitreset代码回滚版本库相关操作远程仓库相关操作分支相关操作创建分支查看分支:gitbranch合并分支:gitmerge删除分支:gitbranch-ddev查看分支合并图:gitlog–graph–pretty=oneline–abbrev-commit撤消某次提交git用户名密码相关配置g

随机推荐