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;
}


评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值