UVA 1331 【区间DP+三角剖分】

题目文档

 题目大意：按照顺时针或者逆时针的顺序给出多边的点，要将这个多边形分解成n-2个三角形，要求使得这些三角行中面积最大的三角形面积尽量小，求最小值。

分析：设dp[i][j]为i,i+1,i+2,....j-1这个多边形的最大三角形面积的最小值，则很容易得出dp[i][j]=min(dp[i][j],max(S(i,j,k),dp[i][k],dp[k][j]));(i<k<j);（假设多边形为凸多边形）。但是实际多边形不一定是凸多变形，这个时候就对每个三角形加以判断，对这个三角形，如果没有其它点在它内部，就是多边形内部的三角形。

祥见代码

#include<cstring>
#include<cstdio>
#include<iostream>
#include<algorithm>
#include<cmath>
using namespace std;
const double INF=1e30;
const double eps=1e-6;
struct edge
{
    double x,y;
}a[55];
int n;
double dp[55][55];
double S(edge a,edge b,edge c)
{
    return fabs((b.x - a.x) * (c.y - a.y) - (c.x - a.x) * (b.y - a.y)) / 2.0 ;
}
bool judge(int aa,int b,int c)
{
    for(int i=1;i<=n;i++)
    {
        if(aa==i||b==i||c==i)continue;
        double s=S(a[i],a[aa],a[b])+S(a[i],a[b],a[c])+S(a[i],a[aa],a[c]);
        if((s-S(a[aa],a[b],a[c]))<eps)return true;
    }
    return false;
}
int main()
{
    int t;
    scanf("%d",&t);
    while(t--)
    {
        scanf("%d",&n);
        for(int i=1;i<=n;i++)
        {
            scanf("%lf%lf",&a[i].x,&a[i].y);
        }
        memset(dp,0,sizeof(dp));
        for(int l=2;l<=n;l++)
        {
            for(int i=1;i<=n-l;i++)
            {
                int j=i+l;
                dp[i][j]=INF;
                for(int k=i+1;k<j;k++)
                {    if(judge(i,k,j))continue;
                    dp[i][j]=min(dp[i][j],max(S(a[i],a[k],a[j]),max(dp[i][k],dp[k][j])));
                }
               // printf("%lf   ",S(i,j,i+1));
                //printf("%d---%d %lf\n",i,j,dp[i][j]);
            }
        }
        printf("%.1f\n",dp[1][n]);
    }
    return 0;
}