ZOJ 3537 Cake 求凸包 区间DP

题意：给出一些点表示多边形顶点的位置（如果多边形是凹多边形就不能切），切多边形时每次只能在顶点和顶点间切，每切一次都有相应的代价。现在已经给出计算代价的公式，问把多边形切成最多个不相交三角形的最小代价是多少。

思路：首先判断多边形是否是凸多边形，之后就是区间dp了。

求出凸包后，按逆时针来看。

设置dp[i][j]为从顶点i到顶点j所围成凸多边形的最优解。

枚举切点k (i < k < j)

dp[i][j] = min(dp[i][k] + dp[k][j] + cost[i][k] + cost[k][j]);

#include<bits/stdc++.h>
using namespace std;
const int MAXN = ;
struct point {
    int x,y;
    friend bool operator < (const point &a,const point &b) {
        if (a.y == b.y) return a.x < b.x;
        return a.y < b.y;
    }
}src[MAXN];
int cost[MAXN][MAXN];
int N,P;
int dp[MAXN][MAXN];
 
point save[MAXN],tmp[MAXN];
int cross(point p0,point p1,point p2) {
    return (p1.x - p0.x) * (p2.y - p0.y) - (p1.y - p0.y) * (p2.x - p0.x);
}
 
int graphm(point * p,int n) {
    sort(p,p + n);
    save[] = p[];
    save[] = p[];
    int top = ;
    for(int i =  ; i < n ; i++){
        while(top && cross(save[top],p[i],save[top-]) >= )top--;
        save[++top] = p[i];
    }
 
    int mid = top;
    for(int i = n -  ; i >=  ; i--){
        while(top > mid && cross(save[top],p[i],save[top-]) >= ) top--;
        save[++top]=p[i];
    }
    return top;
}
 
int getcost(point a,point b) {
    return (abs(a.x + b.x) * abs(a.y+b.y)) % P;
}
 
int main() {
    while (scanf("%d%d",&N,&P) != EOF) {
        for (int i =  ; i < N ; i++) scanf("%d%d",&src[i].x,&src[i].y);
        int num = graphm(src,N);
        if (num < N) {
            printf("I can't cut.\n");
        }
        else {
            memset(cost,,sizeof(cost));
            for (int i =  ; i < N ; i++) {
                for (int j = i +  ; j < N ; j++) cost[i][j] = cost[j][i] = getcost(save[i],save[j]);
            }
            memset(dp,0x3f,sizeof(dp));
            for (int i =  ; i < N ; i++) dp[i][(i + ) % N] = ;
            for (int d =  ; d <= N ; d++) {
                for (int i =  ; i + d -  < N ; i++) {
                    int j = d + i - ;
                    for (int k = i + ; k < j ; k++)
                        dp[i][j] = min(dp[i][j],dp[i][k] + dp[k][j] + cost[i][k] + cost[k][j]);
                }
            }
            printf("%d\n",dp[][N - ]);
        }
    }
    return ;
}