动态规划-凸多边形最优三角形剖分

一.题目概述

①通常，用多边形顶点的逆时针序列表示凸多边形，即P={V0, VI, … Vn 1}表示具有n条V0V1、 V1V2、 … Vn-1Vn的凸多边形。其中，约定V0=Vn。
②多边形的三角剖分是指将多边形分割成互不相交的三角形的弦的集合T。

凸多边形最优三角剖分问题:给定凸多边形P={V0,VI, … Vn},以及定义在由凸多边形的边和弦组成的三角形上的权函数w,要求确定该凸多边形的三角剖分，使得该三角剖分所对应的权，即三角剖分中诸三角形上权之和为最小。

二.问题分析

1.分析：这个问题类似于矩阵连乘的最小乘积次数的问题。
1.矩阵连乘：
①对于矩阵 1~n 数组元素a[i] 表示第i+1个矩阵的行，其中a[n] 表示第n个矩阵的列。
②每每将一个连乘的式子划分成两个部分a[i:k] a[k+1:j]，产生的“权值”（乘法计算次数）是 a[i-1] a[k] a[j]
2.凸多边形剖分：
①对于每个多边形 v[i:j]，i表示第一个顶点，j表示最后一个顶点。
②每每划一条弦ViVk，将多边形分成两个子多边形v[i:k] v[k:j] 产生的权值是w(i,k,j)
可以用下面的树结构表示这个过程：

2.最优子结构：
假设 v0 vn vk 是包含在最优解中， 那么划分出来的两个子多边形的最优解也是包含在最优解中的。
在这个问题中k的选取范围：[1:n+1]

3.状态转移方程：
首先，定义t[i][i] (1≤i<j≤n)为凸子多边形{Vi-1; Vi, ……Vj}的最优三角剖分对应的权函数值，即其最优值。为方便起见，设退化的多边形{V 1, v}具有权值0。据此定义，要计算的凸n+1边形P的最优权值为 t [ 1 ] [ n ]。

三、代码

template<class Type>
void MinWeightTriangulation(int n, Type **t, int **s) {
	for (int i=1; i <= n; i++)//vivi=0
		t[i][i] = 0;
	for (int r=2; r <= n; r++) {//从2到n规模的多边形
		for (int 1=1; i <= n-r+1; i++) {//每个规模为r的多边形
			int j = i+r-1;
			t[i][j] = t[i+1][j] + W(i-1, i, j);
			s[i][j] = i;
			for (int k=i+1; k < i+r-1; k++) {//对于某个规模为r的多边形，求得它的最小权值,k属于【l+1:r-1】
				int u = t[i][k] + t[k+1][j] + W(i-1, k, j);
				if (u<t[i][j]) {
					t[i][j] = u;
					s[i][j] = k;
				}
			}
		}
	}
}