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

一、问题描述

        给定一个凸多边形P={v0,v1,…,vn-1}，以及定义在由多边形的边和弦组成的三角形上的权函数w。要求确定该凸多边形的三角剖分，使得即该三角剖分中诸三角形上权之和为最小。

二、问题分析

        最优子结构：若凸(n+1)边形p={V0,V1,…,VN}的最优三角剖分T包含三角形V0VkVn, 1<=k<n。则T的权为3个部分权的和：三角形V0VkVn的权，子多边形{V0,V1,…,Vk}和{Vk,Vk+1,…,Vn}的权之和。可断言，有T确定的这2个子多女性的三角剖析也是最优。因为若有更小的三角剖析将导致T不是最优的三角剖析的矛盾。

        递归结构：定义t[i][j],1<=i<j<=n为凸子多边形{Vi-1，Vi,…,Vj}的最优三角剖分所对应的权函数值，即其最优值。最优剖分包含三角形Vi-1VkVj的权，子多女性{Vi-1,Vi,…,Vk}和{Vk,Vk+1,…,Vj}的权之和。凸(n+1)边形P的最优权值为t[1][n]。T[1][n]=min{t[1][k]+t[k+1][n]+w(V0VkVn)|1<=k<n}设退化的多边形{Vi-1,Vi}具有的权值为0，因为此刻够不成多边形。

        多边形的点数需要大于3，等于3，可以构成一个三角形，并返回它们的权值和。对于点数大于3，循环遍历所有的划分方式，递归计算左右两侧的最小权值和，并计算当前划分下的权值和。同时记录使得权值和最小的划分方式，并将其对应的点索引保存下来。  最后，输出权值和最小的划分方式，并返回最小权值和。

代码：

#include <iostream> 
using namespace std;
const int N = 7;//凸多边形边数+1
int weight[][N] = { {0,2,2,3,1,4},{2,0,1,5,2,3},{2,1,0,2,1,4},{3,5,2,0,6,2},{1,2,1,6,0,1},{4,3,4,2,1,0} };//凸多边形的权

int MinWeightTriangulation(int n, int** t, int** s);
void Traceback(int i, int j, int** s);//构造最优解
int Weight(int a, int b, int c);//权函数

int main()
{
    int** s = new int* [N];
    int** t = new int* [N];
    for (int i = 0; i < N; i++)
    {
        s[i] = new int[N];
        t[i] = new int[N];
    }

    cout << "此多边形的最优三角剖分值为：" << MinWeightTriangulation(N - 1, t, s) << endl;
    cout << "最优三角剖分结构为：" << endl;
    Traceback(1, 5, s); //s[i][j]记录了Vi-1和Vj构成三角形的第3个顶点的位置
    system("pause");
    return 0;
}

int MinWeightTriangulation(int n, int** t, int** s)
{
    for (int i = 1; i <= n; i++)
    {
        t[i][i] = 0;
    }
    for (int r = 2; r <= n; r++) //r为当前计算的链长（子问题规模）  
    {
        for (int i = 1; i <= n - r + 1; i++)//n-r+1为最后一个r链的前边界  
        {
            int j = i + r - 1;//计算前边界为r，链长为r的链的后边界  
            t[i][j] = t[i + 1][j] + Weight(i - 1, i, j);
            s[i][j] = i;
            for (int k = i + 1; k < j; k++)
            {
                //将链ij划分为( A[i:k] )* (A[k+1:j])   
                int u = t[i][k] + t[k + 1][j] + Weight(i - 1, k, j);
                if (u < t[i][j])
                {
                    t[i][j] = u;
                    s[i][j] = k;
                }
            }
        }
    }
    return t[1][N - 2];
}

void Traceback(int i, int j, int** s)
{
    if (i == j) return;
    Traceback(i, s[i][j], s);
    Traceback(s[i][j] + 1, j, s);
    cout << "三角剖分顶点：V" << i - 1 << ",V" << j << ",V" << s[i][j] << endl;
}

int Weight(int a, int b, int c)
{
    return weight[a][b] + weight[b][c] + weight[a][c];
}