- 1、本文档共11页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
漫话二维凸包-Read.doc
PAGE 1
PAGE 11
漫话二维凸包
作者:蓝点^_*
关键字:算法文化 计算机历史 计算几何 凸包
二维凸包的寻找是计算几何学的经典问题之一。
给定平面上的一个点集,找出一个最小点集顺次连结形成一个凸多边形,使得点集中的点皆在此多边形内或此多边形上,这个凸多边形就是给定点集的二维凸包。
例如,下图中有几个点,图中线段围成的凸多边形就是这几个点的二维凸包。
从若干个点中找凸包就是“二维凸包寻找问题”。大批计算机科学家对于这个看似简单的问题进行了研究,取得了丰硕的成果。然而,研究的过程非常曲折,经常有一些算法发布后,很快就被证明出是错误的。*_*
本文中,我准备谈一谈各个时期涌现出的求凸包的算法。
预备篇 点的排序与左转判定
点的排序
找给定点集的凸包,通常需要一些预处理过程,点的排序就是其中之一。下面给出一种将点按照一定规则排序的方法,这个预处理过程在很多凸包寻找算法中都扮演重要角色。
1.找一个必在凸包上的点(这很容易^_^,通常取横坐标或纵坐标最小的点),记为P0,如图:
2.连结P0与其他点,分别计算这些线段与“竖直向下方向”的夹角,按照夹角由小到大的顺序将各线段的另一端(一端是P0)标号为P1、P2、P3……如图:
(为了使图像更清晰,现略去所有线段)
左转判定
这是经典的计算几何学问题,判断向量p1=(x1,y1)到p2=(x2,y2)是否做左转,只需要判断x1*y2-x2*y1的正负,如果结果为正,则从p1到p2做左转。此结论的证明与本文关系不大,有兴趣的读者请参考有关书籍。 ^_*
下面就让我们一起顺着历史的长河探索凸包的奥秘吧。 ^_^
1972 开山鼻祖
下面谈一谈求凸包的鼻祖算法——“三硬币”算法(The Three-Coins Algorithm)。
三硬币算法由斯卡兰斯奇(Sklansky)于1972年提出,我们可以用三个硬币来模拟这个算法。
斯卡兰斯奇三硬币算法:
1.预处理:将各点排序。(参看预备篇)
2.在P0、P1、P2上分别放置一枚硬币;
把这三枚硬币分别命名为“后”、“中”、“前”。
3.反复
如果三枚硬币按“后-中-前”的顺序“做左转”
拿起“后”,放在“前”的前面;
将原先的“后”改名为“前”;
将原先的“前”改名为“中”;
将原先的“中”改名为“后”;
否则
拿起“中”,放在“后”的后面;
移除刚才“中”所在的点;
将原先的“中”改名为“后”;
将原先的“后”改名为“中”;
直到“前”盖在P0上,且三枚硬币“做左转”
4.按照编号大小顺此连结剩下的点(编号最大的点连回P0),
得到的多边形就是给定点集的凸包。
在纸上画几个点,用三枚硬币模拟很容易模拟这个算法的执行过程,而且似乎求得的结果总是正确的,不过,请您不要试图证明它的正确性,因为,事实上它并非正确*_*。
虽然这个算法并非正确,但其中仍有很多地方值得借鉴,后来的许多凸包算法渗透着一些斯卡蓝斯奇算法的内涵。因此,了解这个算法是很有好处的。
1978年,Bykat确定这个算法是错误的,因此有关此算法的反例,有兴趣的读者不妨参考一下Bykat, A.的《Convex hull of a finite set of points in two dimensions》第296-298页。
1975 失传的西莫斯算法
1975年西莫斯(Michael Shamos)提出了他的第一个凸包算法,可悲的是,这个算法在广泛流传之前似乎就被封杀掉了_,毕竟,这也是个错误的算法。因此,原版的算法失传了。
但是,据说日后西莫斯和另一位计算机科学家波珀里达(Franco Preparata)潜心钻研,对原版的西莫斯算法进行了深入探究,严格地分析了误差,并使西莫斯算法以近似算法的面目重新出现了!
算法焕然一新后,代码却臃肿了不少,而且在信息学竞赛中,这种近似算法用处不大,在这里就不具体解释此算法了^_^
真正的开始
1979年,那是一个春天…… ^_^
麦卡仑姆(McCallum)和艾维斯(Avis)在这一年研究出了第一个正确的凸包算法(注:有少数人持不同说法,他们认为Chand和Kapur在1970就提出了后来所说的“贾维斯步进法”)。
这个算法中要用到两个堆栈。堆栈A包含所有目前考虑到的、有可能构成凸包的点,如果算法结束,堆栈A中装的就是构成凸包的点了,各点在对战中按照横坐标由下而上递增的顺序排列。堆栈T 包含目前考虑到的、不在凸包中的点。
将横坐标最小的点和横坐标最大的点放在A中。
将其它点都放在T中,也把横坐标最小的点放在T中,按照横坐标由下而上递增顺序,将T中元素排序。
文档评论(0)