欢迎关注更多精彩
关注我,学习常用算法与数据结构,一题多解,降维打击。
算法提出
在论文中提到一个结论:
对于任意点集P(n个点), 里面的点不同时分部在一条直线上,则至少有一个三角剖分解。
证明:可以使用构造法来证明。
先对所有点进行字典序排序(先按照X坐标排序,X坐标相同则按照Y坐标排序)。
取出最少的前m个点,使得m个点至少有一个与其他m-1个点不共线。
接下来依次加入下面的n-m个点。每加入一个点pi,就去连接之前已经构建好的凸包的边界顶点(可以从pi观察可见如上图(b)所示)。加完所有剩下的点,三角化就完成了。
正确性证明
对于上述步骤需要说明3点:
- 剩下n-m个点,每次加入的点都在之前凸包的外面。
红点为已经构成凸包的点,pi-1是上次加入的点,黑点为新加入的点。
可以分成2种情况x轴与pi-1相同,x轴与pi-1不同。
不管哪种情况,总是可以找到一条分割线,使得黑点与红点不在同一边。
- 新加入的点必定在新构建的凸包的边界上。
结论1可知,新加入的点不在凸包内,那么想要组成新的凸包,新的点必须在凸包边界上。
- 新形成的三角面与原来的三角形刚好形成新的凸包
先连接最远的两个点,形成凸包。里同再依次按照最小三角形面积部分即可。
判断“可以看到”
上述算法步骤中有一步是连接“可以看到”的点,普通做法是把每个点都连接一遍判断是否有与其他线相交。这样复杂度是O(n^2)。
我们可以发现,每次连接的点都是在凸包的边界上的。只要拿边界的点判断就可以,可且还到两端最远端是存在一个临界点的。
优化:可以维护一个边界表,可以快速判断某条边是否为当前凸包的边界。
可以从上一次加入的点开始连边,依次向两边扩展,直到某个边界点被某条边界挡住就停止。
整体复杂度为排序nlog(n) + 边数量|E|, E为最终边的集合。|E| = 3n-3-h(h是凸包边的数量)
上图列举了几种(不代表全部情况)新点与原凸包对就关系时的搜索路径。
本人码农,希望通过自己的分享,让大家更容易学懂计算机知识。创作不易,帮忙点击公众号的链接。