凸多边形最优三角剖分

凸多边形最优三角剖分一、问题描述多边形是平面上一条分段线性的闭曲线。

也就是说，多边形是由一系列首尾相接的直线段组成的。

组成多边形的各直线段称为该多边形的边。

多边形相接两条边的连接点称为多边形的顶点。

若多边形的边之间除了连接顶点外没有别的公共点，则称该多边形为简单多边形。

一个简单多边形将平面分为3个部分：被包围在多边形内的所有点构成了多边形的内部；多边形本身构成多边形的边界；而平面上其余的点构成了多边形的外部。

当一个简单多边形及其内部构成一个闭凸集时，称该简单多边形为凸多边形。

也就是说凸多边形边界上或内部的任意两点所连成的直线段上所有的点均在该凸多边形的内部或边界上。

通常，用多边形顶点的逆时针序列来表示一个凸多边形，即P={v0 ,v1 ,… ,v n-1}表示具有n条边v0v1，v1v2，… ,v n-1v n的一个凸多边形，其中，约定v0=v n。

若v i与v j是多边形上不相邻的两个顶点，则线段v i v j称为多边形的一条弦。

弦将多边形分割成凸的两个子多边形{v i ,v i+1 ,… ,v j}和{v j ,v j+1 ,… ,v i}。

多边形的三角剖分是一个将多边形分割成互不相交的三角形的弦的集合T。

图1是一个凸多边形的两个不同的三角剖分。

图1 一个凸多边形的2个不同的三角剖分在凸多边形P的一个三角剖分T中，各弦互不相交，且弦数已达到最大，即P的任一不在T中的弦必与T 中某一弦相交。

在一个有n个顶点的凸多边形的三角剖分中，恰好有n-3条弦和n-2个三角形。

凸多边形最优三角剖分的问题是：给定一个凸多边形P={v0 ,v1 ,… ,v n-1}以及定义在由多边形的边和弦组成的三角形上的权函数ω。

要求确定该凸多边形的一个三角剖分，使得该三角剖分对应的权即剖分中诸三角形上的权之和为最小。

可以定义三角形上各种各样的权函数ω。

例如：定义ω(v i v j v k)=|v i v j|+|v i v k|+|v k v j|，其中，|v i v j|是点v i到v j 的欧氏距离。

相应于此权函数的最优三角剖分即为最小弦长三角剖分。

二、算法思路凸多边形的三角剖分与表达式的完全加括号方式之间具有十分紧密的联系。

正如所看到过的，矩阵连乘积的最优计算次序问题等价于矩阵链的完全加括号方式。

这些问题之间的相关性可从它们所对应的完全二叉树的同构性看出。

一个表达式的完全加括号方式对应于一棵完全二叉树，人们称这棵二叉树为表达式的语法树。

例如，与完全加括号的矩阵连乘积((A1(A2A3))(A4(A5A6)))相对应的语法树如图2(a)所示。

图2 表达式语法树与三角剖分的对应语法树中每一个叶子表示表达式中一个原子。

在语法树中，若一结点有一个表示表达式E1的左子树，以及一个表示表达式E r的右子树，则以该结点为根的子树表示表达式(E1E r)。

因此，有n个原子的完全加括号表达式对应于惟一的一棵有n个叶结点的语法树，反之亦然。

凸多边形{v0 ,v1 ,… ,v n-1}的三角剖分也可以用语法树来表示。

例如，图1(a)中凸多边形的三角剖分可用图2(b)所示的语法树来表示。

该语法树的根结点为边v0v6，三角剖分中的弦组成其余的内部结点。

多边形中除v0v6边外的每一条边是语法树的一个叶结点。

树根v0v6是三角形v0v3v6的一条边，该三角形将原多边形分为3个部分：三角形v0v3v6，凸多边形{v0 ,v1 ,… ,v3}和凸多边形{v3 ,v4 ,… ,v6}。

三角形v0v3v6的另外两条边，即弦v3v6和v0v3为根的两个儿子。

以它们为根的子树分别表示凸多边形{v0 ,v1 ,… ,v3}和凸多边形{v3 ,v4 ,… ,v6}的三角剖分。

在一般情况下，一个凸n边形的三角剖分对应于一棵有n-1个叶子的语法树。

反之，也可根据一棵有n-1个叶子的语法树产生相应的一个凸n边形的三角剖分。

也就是说，凸n边形的三角剖分与n-1个叶子的语法树之间存在一一对应关系。

由于n个矩阵的完全加括号乘积与n个叶子的语法树之间存在一一对应关系，因此n个矩阵的完全加括号乘积也与凸(n+1)边形的三角剖分之间存在一一对应关系。

图2的(a)和(b)表示出了这种对应关系，这时n=6。

矩阵连乘积A1A2..A6中的每个矩阵A i对应于凸(n+1)边形中的一条边v i-1v i。

三角剖分中的一条弦v i v j，i<j，对应于矩阵连乘积A[i+1:j ]。

事实上，矩阵连乘积的最优计算次序问题是凸多边形最优三角剖分问题的一个特殊情形。

对于给定的矩阵链A1A2..A n，定义一个与之相应的凸(n+1)边形P={v0 ,v1 ,… ,v n}，使得矩阵A i与凸多边形的边v i-1v i一一对应。

若矩阵A i的维数为p i-1×p i,i=1,2,…,n，则定义三角形v i v j v k上的权函数值为：ω(v i v j v k)=p i p j p k。

依此权函数的定义，凸多边形P的最优三角剖分所对应的语法树给出矩阵链A1A2..A n的最优完全加括号方式。

三、实验源程序新建一个类CTriangle，Class type选择Generic Class。

Triangle.h代码typedef struct{int x;int y;}point;class CTriangle{public:bool Run();CTriangle();virtual ~CTriangle();private://用递归的方法输出剖分后的各个三角形void Traceback(int i,int j,int **s);//计算最优值算法bool minWeightTriangulation();//处理键盘输入，同时判断能否构成一个凸多边形 bool Input();//计算三角形权值的函数int w(point X,point Y,point Z);//计算平面上任意两点间距离的函数int distance(point X,point Y); //记录最优三角剖分中所有三角形信息 int **s;//记录最优三角剖分所对应的权函数值 int **t;//记录凸多边形各顶点坐标point *v;//记录坐标在直线方程中的值int *total;int M;};Triangle.cpp代码#define N 50#include <iostream.h>#include <math.h>#include <stdlib.h>#include "Triangle.h"CTriangle::CTriangle(){M = 0;t = new int *[N];s = new int *[N];for(int i=0 ; i<N ; i++){t[i] = new int[N];s[i] = new int[N];}v = new point[N];total = new int[N];}CTriangle::~CTriangle(){for(int i=0 ; i<N ; i++){delete []t[i];delete []s[i];}delete []t;delete []s;delete []v;delete []total;}int CTriangle::distance(point X, point Y) {int dis = (Y.x-X.x)*(Y.x-X.x) +(Y.y-X.y)*(Y.y-X.y);return (int)sqrt(dis);}int CTriangle::w(point X, point Y, point Z) {return distance(X,Y) + distance(Y,Z) + distance(Z,X); }bool CTriangle::Input(){int m;int a,b,c;cout<<"请输入凸多边形顶点个数：";cin>>m;M = m-1;for(int i=0 ; i<m ; i++){cout<<"输入顶点v"<<i<<"的坐标："; cin>>v[i].x>>v[i].y;}//根据顶点坐标判断是否能构成一个凸多边形 for(int j=0 ; j<m ; j++){int p = 0;int q = 0;if(m-1 == j){a = v[m-1].y - v[0].y;b = v[m-1].x - v[0].y;c = b * v[m-1].y - a * v[m-1].x;}else{a = v[j].y - v[j+1].y;b = v[j].x - v[j+1].x;c = b * v[j].y - a * v[j].x;}for(int k=0 ; k<m ; k++){total[k] = a * v[k].x - b * v[k].y + c;if(total[k] > 0){p = p+1;}else if(total[k] < 0){q = q+1;}}if((p>0 && q>0) || (p==0 && q==0)){cout<<"无法构成凸多边形！"<<endl; exit(1);}}if(NULL != v) return true;elsereturn false;}bool CTriangle::minWeightTriangulation() {if(NULL == v)return false;for(int i=1 ; i<=M ; i++)t[i][i] = 0;for(int r=2 ; r<=M ; r++)for(int i=1 ; i<=M-r+1 ; i++){int j = i+r-1;t[i][j] = t[i+1][j] + w(v[i-1],v[i],v[j]); s[i][j] = i;for(int k=i+1 ; k<i+r-1 ; k++){int u = t[i][k] + t[k+1][j] +w(v[i-1],v[k],v[j]);if(u < t[i][j]){t[i][j] = u;s[i][j] = k;}}}return true;}void CTriangle::Traceback(int i, int j, int **s) {if(i == j)return;/*{cout<<"分成的三角形依次为："<<endl; cout<<"v"<<i-1<<"v"<<i<<"v"<<j<<endl; }else*/Traceback(i,s[i][j],s);Traceback(s[i][j]+1,j,s);cout<<"三角形：v"<<i-1<<"v"<<s[i][j]<<"v"<<j<<endl;}bool CTriangle::Run(){if(Input()){ if(CTriangle::minWeightTriangulation()) {CTriangle::Traceback(1,M,s);cout<<endl;cout<<"最优权值之和为："<<t[1][M]<<endl;return true;}elsereturn false;}elsereturn false;}main.cpp代码#include "Triangle.h"void main(){CTriangle triangle;triangle.Run();}四、测试结果给定的七个顶点坐标分别为(8,26)、(0,20)、(0,10)、(10,0)、(22,12)、(27,21)、(15,26)，测试结果如下：该程序有一定的灵活性，用户可以自己选择顶点的个数以及顶点坐标，并对顶点坐标进行分析，判断这些顶点能否构成一个凸多边形，例如输入四个点坐标(0,0)、(0,10)、(4,4)、(10,0)，这四个点构成的多边形显然是一个凹多边形，测试结果为：五、总结（1）凸多边形的最优三角剖分本来是一个几何问题，但通过分析，它本质上于矩阵连乘积的最优计算次序问题极为相似，从而可以将问题进行转化，用动态规划算法有效的解决这个问题。