漫话二维凸包-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中元素排序。