凸多边形的最优三角剖分问题

最新推荐文章于 2022-04-16 18:49:15 发布

hechenghai

最新推荐文章于 2022-04-16 18:49:15 发布

阅读量2.5k

点赞数 1

分类专栏：学习笔记文章标签：算法 matrix qt n2

学习笔记专栏收录该内容

77 篇文章 0 订阅

订阅专栏

凸多边形的最优三角剖分问题

问题描述

多边形是平面上一条分段线性的闭曲线。也就是说，多边形是由一系列首尾相接的直线段组成的。组成多边形的各直线段称为该多边形的边。多边形相接两条边的连接点称为多边形的顶点。若多边形的边之间除了连接顶点外没有别的公共点，则称该多边形为简单多边形。一个简单多边形将平面分为3个部分：被包围在多边形内的所有点构成了多边形的内部；多边形本身构成多边形的边界；而平面上其余的点构成了多边形的外部。当一个简单多边形及其内部构成一个闭凸集时，称该简单多边形为凸多边形。也就是说凸多边形边界上或内部的任意两点所连成的直线段上所有的点均在该凸多边形的内部或边界上。

通常，用多边形顶点的逆时针序列来表示一个凸多边形，即P=<v₀,v₁,… ,v_n-1>表示具有n条边v₀v₁，v₁v₂，… ,v_n-1v_n的一个凸多边形，其中，约定v₀=v_{n 。}

若v_i与v_j是多边形上不相邻的两个顶点，则线段v_iv_j称为多边形的一条弦。弦将多边形分割成凸的两个子多边形<v_i,v_i+1,… ,v_j>和<v_j,v_j+1,… ,v_i>。多边形的三角剖分是一个将多边形分割成互不重迭的三角形的弦的集合T。图1是一个凸多边形的两个不同的三角剖分。


(a)	(b)

图1 一个凸多边形的2个不同的三角剖分

在凸多边形P的一个三角剖分T中，各弦互不相交且弦数已达到最大，即P的任一不在T中的弦必与T中某一弦相交。在一个有n个顶点的凸多边形的三角刮分中，恰好有n-3条弦和n-2个三角形。

凸多边形最优三角剖分的问题是：给定一个凸多边形P=<v₀,v₁,… ,v_n-1>以及定义在由多边形的边和弦组成的三角形上的权函数ω。要求确定该凸多边形的一个三角剖分，使得该三角剖分对应的权即剖分中诸三角形上的权之和为最小。

注意：解决此问题的算法必须适用于任意的权函数。

参考解答

用动态规划算法也能有效地求解凸多边形的最优三角剖分问题。尽管这是一个计算几何学问题，但在本质上，它与矩阵连乘积的最优计算次序问题极为相似。

1.三角剖分的结构及其相关问题

凸多边形的三角剖分与表达式的完全加括号方式之间具有十分紧密的联系。正如所看到过的，矩阵连乘积的最优计算次序问题等价于矩阵链的完全加括号方式。这些问题之间的相关性可从它们所对应的完全二叉树的同构性看出。这里的所谓完全二叉树是指叶结点以外的所有结点的度数都为2的二叉树(注意与满二叉树和近似满二叉树的区别)。

一个表达式的完全加括号方式对应于一棵完全二叉树，人们称这棵二叉树为表达式的语法树。例如，与完全加括号的矩阵连乘积((A₁(A₂A₃))(A₄(A₅A₆)))相对应的语法树如图3(a)所示。

图3 表达式语法树与三角剖分的对应

语法树中每一个叶子表示表达式中一个原子。在语法树中，若一结点有一个表示表达式E₁的左子树，以及一个表示表达式E_r的右子树，则以该结点为根的子树表示表达式(E₁E_r)。因此，有n个原子的完全加括号表达式对应于唯一的一棵有n个叶结点的语法树，反之亦然。

凸多边形<v₀,v₁,… ,v_n-1>的三角剖分也可以用语法树来表示。例如，图3(a)中凸多边形的三角剖分可用图3(b)所示的语法树来表示。该语法树的根结点为边v₀v₆，三角剖分中的弦组成其余的内部结点。多边形中除v₀v₆边外的每一条边是语法树的一个叶结点。树根v₀v₆是三角形v₀v₃v₆的一条边，该三角形将原多边形分为3个部分：三角形v₀v₃v₆，凸多边形<v₀,v₁,… ,v₃>和凸多边形<v₃,v₄,… ,v₆>。三角形v₀v₃v₆的另外两条边，即弦v₃v₆和v₀v₃为根的两个儿子。以它们为根的子树分别表示凸多边形<v₀,v₁,… ,v₃>和凸多边形<v₃,v₄,… ,v₆>的三角剖分。

在一般情况下，一个凸n边形的三角剖分对应于一棵有n-1个叶子的语法树。反之，也可根据一棵有n-1个叶子的语法树产生相应的一个凸n边形的三角剖分。也就是说，凸n边形的三角剖分与n-1个叶子的语法树之间存在一一对应关系。由于n个矩阵的完全加括号乘积与n个叶子的语法树之间存在一一对应关系，因此n个矩阵的完全加括号乘积也与凸(n+1)边形的三角剖分之间存在一一对应关系。图3的(a)和(b)表示出了这种对应关系，这时n=6。矩阵连乘积A₁A₂..A₆中的每个矩阵A_i对应于凸(n+1)边形中的一条边v_i-1v_i。三角剖分中的一条弦v_iv_j，i<j-1，对应于矩阵连乘积A_i+1..j 。

事实上，矩阵连乘积的最优计算次序问题是凸多边形最优三角剖分问题的一个特殊情形。对于给定的矩阵链A₁A₂..A_n，定义一个与之相应的凸(n+1)边形P=<v₀,v₁,… ,v_n>，使得矩阵A_i与凸多边形的边v_i-1v_i一一对应。若矩阵A_i的维数为p_i-1×p_i,i=1,2,…,n，则定义三角形v_iv_jv_k上的权函数值为： ω(△v_iv_jv_k)=p_ip_jp_k。依此权函数的定义，凸多边形P的最优三角剖分所对应的语法树给出矩阵链A₁A₂..A_n的最优完全加括号方式。

2.最优子结构性质

凸多边形的最优三角剖分问题有最优子结构性质。事实上，若凸(n+1)边形P=<v₀,v₁,… ,v_n>的一个最优三角剖分T包含三角形v₀v_kv_{n ，}1≤k≤n-1，则T的权为3个部分权的和，即三角形v₀v_kv_n的权，子多边形<v₀,v₁,… ,v_k>的权和<v_k,v_k+1,… ,v_n>的权之和。可以断言由T所确定的这两个子多边形的三角剖分也是最优的，因为若有<v₀,v₁,… ,v_k>或<v_k,v_k+1,… ,v_n>的更小权的三角剖分，将会导致T不是最优三角剖分的矛盾。

3.最优三角剖分对应的权的递归结构

首先，定义t[i,j]，1≤i<j≤n，为凸子多边形<v_i-1,v_i,… ,v_j>的最优三角剖分所对应的权值，即最优值。为方便起见，设退化的多边形<V_i-1 ,v_i>具有权值0。据此定义，要计算的凸(n+1)边多边形P对应的权的最优值为t[1,n]。

t[i,j]的值可以利用最优子结构性质递归地计算。由于退化的2顶点多边形的权值为0，所以t[i,i]=0，i=1,2,…,n 。当j一i≥1时，子多边形<v_i-1,v_i,… ,v_j>至少有3个顶点。由最优于结构性质,t[i,j]的值应为t[i,k]的值加上t[k+1,j]的值，再加上△v_i-1v_kv_j的权值，并在i≤k≤j-1的范围内取最小。由此，t[i,j]可递归地定义为：

4.计算最优值

将(2.3)式与矩阵连乘积的最优计算次序问题中计算m[i,j]的(2.l)式进行比较容易看出，除了权函数的定义外，两个递归式是完全一样的。因此只要对计算m[i,j]的算法MATRIX_CHAIN_ORDER做很小的修改就完全适用于计算t[i,j]。

下面描述的计算凸(n+1)边形P=<v₀,v₁,… ,v_n>的三角剖分最优权值的动态规划算法MINIMUM_WEIGHT_TRIANGULATION，输入是凸多边形P=<v₀,v₁,… ,v_n>的权函数ω，输出是最优值t[i,j]和使得t[i,k]+t[k+1,j]+ω(△v_i-1v_kv_j)达到最优的位置(k=)s[i,j]，1≤i≤j≤n 。

Procedure MINIMUM_WEIGHT_TRIANGULATION(P,w)；
begin
  n:=length[p]-1;
  for i:=1 to n do t[i,i]:=0;
    for ll:=2 to n do
      for i:=1 to n-ll+1 do
        begin
          j:=i+ll-1;
          t[i,j]:=∞;
          for k:=i to j-1 do
            begin
              q:=t[i,k]+t[k+1,j]+ω(△vi-1vkvj);
              if q<t[i,j] then
                begin
                  t[i,j]:=q;
                  s[i,j]:=k;
                end;
            end;
        end;
  return(t,s);
end;

与MATRIX_CHAIN_ORDER一样，算法MINIMUM_WEIGHT_TRIANGULATION占用θ(n²)空间，耗时θ(n³)。

5.构造最优三角剖分

如我们所看到的，对于任意的1≤i≤j≤n ，算法MINIMUM_WEIGHT_TRIANGULATION在计算每一个子多边形<v_i-1,v_i,… ,v_j>的最优三角剖分所对应的权值t[i,j]的同时，还在s[i,j]中记录了此最优三角剖分中与边(或弦)v_i-1v_j构成的三角形的第三个顶点的位置。因此，利用最优子结构性质并借助于s[i,j]，1≤i≤j≤n ，凸(n+l)边形P=<v₀,v₁,… ,v_n>的最优三角剖分可容易地在Ο(n)时间内构造出来。