【动态规划】UVa 1331 最大面积最小三角形剖分

题目

请点击

题目大意

将一个多边形用它不相交的对角线将它分成若干个三角形，使得最大的三角形面积最小，求最大三角形的面积。如图是一个六边形的几种剖分：

思路

记由点u,u+1,…,v-1,v(u< v)组成的多边形为$F(i,j)$

首先，总的多边形为$F(1,N)$，
对于$F(1,N)$的一个子多边形$F(i,j)(1<=i<j<=N)$，可以在点i+1-点j-1中任意找到一个点k，将其分割为：$\Delta ijk,F(i,k),F(j,k)$，如图:

设d[i][j]表示$F(i,j)$中的最优解，很显然，先要找到$max \{ S\Delta ijk,d[i][k],d[k][j]\}$（$S\Delta ijk$为$\Delta ijk$的面积），也就是找到以k为分割点，$F(i,j)$中的最大的三角形的面积，然后枚举k点，使这个最大的三角形面积最小即可。

状态转移方程：$d[i][j]=min\{ max \{ S\Delta ijk,d[i][k],d[k][j]\}\}(i<k<j)$

这样就完了吗？不是的，还有一个问题：$\Delta ijk$合法吗？例如这样的“分割”：

当$k$点在如图位置时，$\Delta ijk$显然不合法，因为线段$ik$不是$F(1,7)$的对角线！

所以怎么判断$\Delta ijk$是否合法呢？
换个思路，判断$\Delta ijk$是否合法实际上就是看$\Delta ijk$内有没有其他的点。然后枚举除ijk外每一个点呗，再看这点是不是在$\Delta ijk$内。

怎么判断点是否在$\Delta ijk$内呢？
“面积法”是很好的办法：设现在枚举到的点为$x$，计算$T=S\Delta ijx,S\Delta ikx,S\Delta kjx$，再把它们加起来，如果$x$不在$\Delta ijk$，显然$T==S\Delta ijk$，所以判断$T$与$S\Delta ijk$的大小关系即可。不明白的可以画画图，我懒得再画了……

至于计算面积，我用的海伦公式。

代码

#include<cstdio>
#include<cmath>
#include<cstring>
#include<algorithm>
using namespace std;

#define MAXN 50
#define INF 0x7fffffff
#define eps 0.0001
struct point
{
    double x,y;
}p[MAXN+5];//点
int N;
double d[MAXN+5][MAXN+5];//动归

double dis(int x,int y){return sqrt((p[x].x-p[y].x)*(p[x].x-p[y].x)+(p[x].y-p[y].y)*(p[x].y-p[y].y));}//计算两点距离
double area(int x,int y,int z)//计算由点x,y,z构成的三角形的面积（海伦公式）
{
    double a=dis(x,y),b=dis(y,z),c=dis(x,z);
    double p=(a+b+c)/2;
    return sqrt(p*(p-a)*(p-b)*(p-c));
}
bool check(int x,int y,int z)//判断由点x,y,z构成的三角形中有没有点
{
    double tarea=area(x,y,z);
    for(int i=1;i<=N;i++)
    {
        if(i==x||i==y||i==z) continue;
        double a=area(i,x,y),b=area(i,y,z),c=area(i,x,z);
        if(fabs(a+b+c-tarea)<eps)//double计算有精度误差，不要用==
            return 0;
    }
    return 1;
}

int main()
{
    int T;
    scanf("%d",&T);
    while(T--)
    {
        scanf("%d",&N);
        for(int i=1;i<=N;i++)
            scanf("%lf%lf",&p[i].x,&p[i].y);
        for(int i=N-2;i>=1;i--)
            for(int j=i+2;j<=N;j++)
            {
                d[i][j]=INF;
                for(int k=i+1;k<j;k++)//选择分割点
                    if(check(i,j,k))
                        d[i][j]=min(d[i][j],max(area(i,j,k),max(d[i][k],d[k][j])));
            }
        printf("%.1lf\n",d[1][N]);
        memset(d,0,sizeof(d));
    }
}