凸多边形最优三角形剖分 区间DP

题目描述

众所周知，sciorz会画画。某天，sciorz画了一个凸多边形，这个多边形的每个顶点都有一个权值a[i]。sciorz觉得这个凸多边形不够美丽，于是他决定在n个点之间连线，最终用n-3条不相交的线将这个凸n边形分割成n-2个三角形。sciorz认为，一个三角形的美丽值是三个顶点权值的乘积，凸多边形的美丽值是其内部三角形的美丽值的和。sciorz想找到一种分割方案，使得这个凸多边形的美丽值最大。sciorz忙着刷难题，所以他随手就把这个签到题扔给你，希望你帮sciorz算出最大的美丽值。

输入
第一行一个t,表示有t组样例。
每组样例的第一行是一个n,表示多边形的边数。
第二行n个数,第i个数表示多边形第i个顶点的权值a[i]，按逆时针顺序给出。
输出
对于每组样例，输出一行。格式为"Case #x: y"，x为样例编号，y为答案。

#include<cstdio>
#include<iostream>
#include<algorithm>
#include<cstring>
#include<cmath>
#include<vector>
#include<set>
#include<map>
#include<queue>
using namespace std;
#define maxn 100005
#define ll long long
#define pi acos(-1.0)
#define inf 0x3f3f3f3f
int t,n;
int a[105];
int b[105];
ll dp[105][105];
int main()
{
    scanf("%d",&t);
    for(int w=1;w<=t;w++)
    {
        scanf("%d",&n);
        for(int i=1;i<=n;i++)
        {
            scanf("%d",&a[i]);
        }
        for(int i=0;i<=100;i++)
        {
            for(int j=0;j<=100;j++)
            {
                dp[i][j]=0;
            }
        }
        for(int i=n;i>=1;i--)
        {
            for(int j=i+1;j<=n;j++)
            {
                if(j-i==1)
                    dp[i][j]=0;
                else if(j-i==2)
                {
                    dp[i][j]=a[i]*a[i+1]*a[i+2];
                }
                else
                {
                    for(int k=i+1;k<=j-1;k++)
                    {
                        dp[i][j] = max(dp[i][j],(dp[i][k]+dp[k][j]+a[i]*a[j]*a[k]));
                    }
                }
            }
        }
        printf("Case #%d: %d\n",w,dp[1][n]);
    }
}

dp[i][j] 表示从第i到第j个点最优剖分
的答案。 当 j=i 或 j=i+1 时，无法构成三角形， dp[i][j]=0 。 当 j>=i+2 时，考虑一个分割点k, ijk
三点构成一个三角形，此时划分为i到k部分，k到j部分，以及三角形ijk部分这三个部分。i到k部分和k到j
部分是子问题，可以利用之前求出的结果，三角形ijk的值是三点权值之积。 状态转移方程为DP(i,j)=max(DP(i,j),DP(i,k)+DP(k,j)+ant[i]*ant[k]*ant[j])
当 j=i 或 j=i+1 时，无法构成三角形， dp[i][j]=0 。
当 j=i+2 时，三个相邻的点构成了一个三角形。