一、问题描述
用多边形顶点的逆时针序列表示凸多边形,即P={v0,v1,…,vn-1}表示具有n条边的凸多边形。 若vi与vj是多边形上不相邻的2个顶点,则线段vivj称为多边形的一条弦。弦将多边形分割成2个多边形{vi,vi+1,…,vj}和{vj,vj+1,…vi}。 多边形的三角剖分是将多边形分割成互不相交的三角形的弦的集合T。 给定凸多边形P,以及定义在由多边形的边和弦组成的三角形上的权函数w。要求确定该凸多边形的三角剖分,使得即该三角剖分中诸三角形上权之和为最小。
二、问题分析
最优子结构性质:
若凸(n+1)边形P={V0,V1……Vn}的最优三角剖分T包含三角形V0VkVn,1<=k<=n,则T的权为三个部分权之和:三角形V0VkVn的权,多边形{V0,V1……Vk}的权和多边形{Vk,Vk+1……Vn}的权之和。由T确定的这两个子多边形的三角剖分也是最优的。
因此,可得递推关系式:
三、代码如下
#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;
}