凸多边形的最优三角剖分

一，题意理解

该题要求给定多边形的边数n，然后给定每一条弦的权重，要求选出一组弦，可以将该多边形分成多个互不相交的三角形，并且使得划分的所有的三角形的总权重最小。

二，解题思路

运用类似于矩阵连乘的思路

矩阵连乘的状态转移方程如下

在矩阵连乘的解决方法中，我们从矩阵数为1开始遍历，逐步求数目为n的矩阵的最优解

在该题中我们同样可以把边数为n的多边形看成n个矩阵相乘，找出最优解，只是需要找出新的状态转移方程。

在此我们同样建立A[i][j]，其表示以i-1，i，i+1,i+2,,,,,,j  为顶点的多边形的最优权重。

W[i][j],记录弦i，j的权重。

由于至少三个点才能构成一个多边形，我们从长度为3（即此时正好为一个三角形）开始遍历逐步求出最优解。

倘若多边形A[i][j]的最优三角剖分T包含三角形vi-1vkvj，i≤k≤j-1，则T的权为3个部分权的和：三角形vi-1vkvj的权，子多边形A[i][k]和A[k+1][j]的权之和。

如此我们便可得状态转移方程

代码如下：

#include<bits/stdc++.h>
using namespace std;
 
const int M =1111;
int n; //顶点数
int s[M][M];//记录最优策略二维数组
double t[M][M];//记录最优值二维数组
double weight[M][M];//记录各顶点之间权值的二维数组
 
void MinWeightTriangulation()
{
    for (int i=1;i<=n;i++) //初始化
    {
        t[i][i]=0;
        s[i][i]=0;
    }
    for(int r=2;r<=n;r++)//r为问题规模，r=2实际上有三个点 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]+weight[i][j]+weight[i-1][j];//将链ij划分为A(i) * ( A[i+1:j] )这里实际上就是k=i
            //策略为k=i的情况
            s[i][j]=i;
            for (int k=i+1;k<j;k++) //枚举划分点i到j所有情况
            {
            	//将链ij划分为( A[i:k] )* (A[k+1:j])
                double u=t[i][k]+t[k+1][j]+weight[i-1][k]+weight[k][j]+weight[i-1][j];
                if(t[i][j]>u)
                {
                    t[i][j]=u;
                    s[i][j]=k;
                }
            }
        }
    }
}
 
void Traceback(int i,int j)  //递归求解所有子问题的弦。
{
    if(i==j)
        return ;
    Traceback(i,s[i][j]);
    Traceback(s[i][j]+1,j);
    cout<<"三角剖分顶点：V"<<i-1<<",V"<<j<<",V"<<s[i][j]<<endl;
}
 
int main()
{
    int i;
    int j;
    cout << "请输入顶点个数n："<< endl;
    cin >> n;
    n--;
    cout << "请依次输入各顶点的连接权值：" << endl;
    for (i=0;i<=n;++i)
    {
        for (j=0;j<=n;++j)
        {
            cin >> weight[i][j];
        }
    }
    MinWeightTriangulation();
    cout << "最优三角剖分的权值和是:" << endl;
    cout << t[1][n]<< endl;
    cout << "三角剖分顶点："<< endl;
    Traceback(1,n);
    return 0;
}