zoj3537 cake（凸包+区间dp）

首先利用凸包判断给出的点是否正好组成一个凸多边形，如果可以的话，定义状态dp[i][j]为从i点到j点构成的凸多边形的最优解。因为凸多边形是一个环，所以我们将节点后面再复制一遍节点，这样方便计算。然后就可以按照区间dp去做了。

转移方程： dp[i][j]=min(dp[i][j],dp[i][k]+dp[k][j]+temp);k为枚举的切割点，temp为割i到j点的代价。当i和j相邻时，显然temp为0，因为它是一条边，不需要我们割。

代码：

#include<stdio.h>
#include<iostream>
#include<algorithm>
#include<string.h>
#include<math.h>
#include<stack>
using namespace std;
#define inf 0x3f3f3f3f
#define pi acos(-1)
#define maxn 650
int xx,yy;
int dp[maxn][maxn];
struct Node
{
    int x;
    int y;
} node[maxn];
Node st[maxn];
bool cmp1(Node a,Node b)
{
    if(a.y==b.y)return a.x<b.x;
    return a.y<b.y;
}
bool cmp2(Node a,Node b)
{
    if(atan2(a.y-yy,a.x-xx)==atan2(b.y-yy,b.x-xx))
        return a.x<b.x;
    return atan2(a.y-yy,a.x-xx)<atan2(b.y-yy,b.x-xx);
}
int chaji(Node a,Node b,Node c)
{
    return (b.x-a.x)*(c.y-a.y)-(b.y-a.y)*(c.x-a.x);
}
int main()
{
    int n,p;
    while(~scanf("%d%d",&n,&p))
    {
        for(int i=1; i<=n; i++)
        {
            scanf("%d%d",&node[i].x,&node[i].y);

        }
        sort(node+1,node+n+1,cmp1);
        st[1]=node[1];
        xx=node[1].x;
        yy=node[1].y;
        sort(node+1,node+n+1,cmp2);
        st[2]=node[2];
        int top=2;
        for(int i=3; i<=n; i++)
        {
            while(chaji(st[top-1],st[top],node[i])<=0)//注意此处还要判断等于
            {
                top--;
            }
            st[++top]=node[i];
        }
//        for(int i=1;i<=top;i++)
//        {
//            cout<<st[i].x<<" "<<st[i].y<<endl;
//        }
        if(top!=n)//全部节点不能构成一个凸包
        {
            printf("I can't cut.\n");
            continue;
        }
        for(int i=1; i<=n; i++)//因为是一个环，我们在复制一遍节点，然后则这个数组上dp
        {
            node[i+n]=node[i];
        }
        memset(dp,inf,sizeof(dp));
        int temp;
        for(int len=3; len<=n; len++)
        {
            for(int i=1; i<=2*n-len+1; i++)
            {
                int j=i+len-1;
                for(int k=i+1; k<j; k++)
                {
                    if(k-i+1<3)dp[i][k]=0;//这部分的代价为0，因为不能构成一个三角形
                    if(j-k+1<3)dp[k][j]=0;//同理
                    if(abs(i-j)==1||abs(i-j)==n-1)temp=0;//i和j相邻，不需要割
                    else
                    {
                        temp=(abs(node[i].x+node[j].x)*abs(node[i].y+node[j].y))%p;
                    }
                    dp[i][j]=min(dp[i][j],dp[i][k]+dp[k][j]+temp);//维护最小值
                }
            }
        }
        int ans=inf;
        for(int i=1; i<=n; i++)
        {
            ans=min(ans,dp[i][i+n-1]);//找出最优解
        }
        printf("%d\n",ans);


    }
}
/*
4 6
6 1
7 9
4 12
1 9
0
*/