[论文解读] Algorithms for Sampling 3-Orientations of Planar Triangulations

Sarah Miracle, Dana Randall+2 more|arXiv (Cornell University)|Feb 22, 2012

Computational Geometry and Mesh Generation参考文献 27被引用 4

一句话总结

本文分析了两种马尔可夫链——三角形反转和边翻转——用于对平面三角剖分的3-定向进行采样。研究证明，对于所有具有n个内部顶点的平面三角剖分，边翻转链都是快速混合的；而三角形反转链仅在最大内部度数不超过六时是快速混合的，但在高阶三角剖分中可能变得缓慢。

ABSTRACT

Given a planar triangulation, a 3-orientation is an orientation of the internal edges so all internal vertices have out-degree three. Each 3-orientation gives rise to a unique edge coloring known as a Schnyder wood that has proven powerful for various computing and combinatorics applications. We consider natural Markov chains for sampling uniformly from the set of 3-orientations. First, we study a "triangle-reversing" chain on the space of 3-orientations of a fixed triangulation that reverses the orientation of the edges around a triangle in each move. It was shown previously that this chain connects the state space and we show that (i) when restricted to planar triangulations of maximum degree six, the Markov chain is rapidly mixing, and (ii) there exists a triangulation with high degree on which this Markov chain mixes slowly. Next, we consider an "edge-flipping" chain on the larger state space consisting of 3-orientations of all planar triangulations on a fixed number of vertices. It was also shown previously that this chain connects the state space and we prove that the chain is always rapidly mixing. The triangle-reversing and edge-flipping Markov chains both arise in the context of sampling other combinatorial structures, such as Eulerian orientations and triangulations of planar point sets, so our results here may shed light on the mixing rate of these related chains as well.

研究动机与目标

• 分析用于对平面三角剖分中3-定向进行采样的马尔可夫链的混合时间。
• 解决Felsner和Zickfeld提出的关于3-定向高效采样的开放问题。
• 确定自然局部操作（三角形反转与边翻转）在何种条件下能产生快速混合的链。
• 为使用马尔可夫链蒙特卡洛方法对3-定向进行采样建立收敛速度的理论界限。
• 扩展对相关组合结构（如欧拉定向和三角剖分）中混合速率的理解。

提出的方法

• 提出一种三角形反转马尔可夫链，通过反转有向三角形周围的边方向，保持顶点出度不变。
• 引入一种边翻转马尔可夫链，适用于所有具有n个顶点的平面三角剖分，允许在不同三角剖分之间转移。
• 使用规范路径分析来界定边翻转链的混合时间，建立多项式时间运行保证。
• 应用三角形反转链与环路反转链之间的比较技术，推导混合时间界限。
• 将转移路径分解为多个阶段（如红色/绿色边交换）以界定每条边所用规范路径数的上界。
• 利用3-定向与二分图中完美匹配之间已知的双射关系，指导采样复杂度的分析。

实验结果

研究问题

• RQ1 在何种条件下，三角形反转马尔可夫链对固定平面三角剖分的3-定向采样是快速混合的？
• RQ2 当在所有具有n个内部顶点的平面三角剖分中采样时，边翻转马尔可夫链是否快速混合？
• RQ3 能否为最大内部度数Δ ≤ 6的三角剖分，对三角形反转链的混合时间进行界定？
• RQ4 当最大度数无界时，三角形反转链的混合时间是否会变为指数级？
• RQ5 规范路径方法能否有效应用于边翻转链在更大状态空间中的混合时间界定？

主要发现

• 边翻转马尔可夫链在所有具有n个内部顶点的平面三角剖分的并集中是快速混合的，混合时间为O(n^10 log(n/ε))。
• 对于最大内部度数不超过六的三角剖分，三角形反转链是快速混合的，混合时间为O(n^5 log(ε⁻¹))。
• 对于最大内部度数超过六的三角剖分，三角形反转链可能混合缓慢，需要指数时间。
• 规范路径分析表明，任意给定边最多被O(n^4)条规范路径使用。
• 边翻转链的混合时间界限是通过路径分解、路径长度和平稳分布界限的组合推导得出的。
• 结果表明，边翻转链为固定大小的所有平面三角剖分提供了高效的3-定向采样算法。