# zoj 3537 cake 切蛋糕 区间DP+凸包+递归 最优三角形剖分

/*
*  Author:      illuz <iilluzen[at]gmail.com>
*  Blog:        http://blog.csdn.net/hcbbt
*  File:        zoj3537.cpp
*  Create Date: 2013-12-04 16:35:27
*  Descripton:  convex hull + intervel dp
*/

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

#define sqr(a) ((a) * (a))
#define dis(a, b) sqrt(sqr(a.x - b.x) + sqr(a.y - b.y))

const int MAXN = 305;
const int INF = 0x3c3c3c3c;
const double PI = acos(-1.0);

struct Point {
double x;
double y;
Point(double a = 0, double b = 0) : x(a), y(b) {}
friend bool operator < (const Point &l, const Point &r) {
return l.y < r.y || (l.y == r.y && l.x < r.x);
}
} p[MAXN], ch[MAXN * 2], tmp[MAXN];
// p, point   ch, convex hull

int f[MAXN][MAXN], c[MAXN][MAXN];		// rec and cast
int P;

double mult(const Point &a, const Point &b, const Point &o) {
return (a.x - o.x) * (b.y - o.y) >= (b.x - o.x) * (a.y - o.y);
}

double Graham(Point p[], int n, Point res[]) {
int top = 1;
sort(p, p + n);
if (n == 0) return 0;
res[0] = p[0];
if (n == 1) return 0;
res[1] = p[1];
if (n == 2) return dis(p[0], p[1]) * 2;
res[2] = p[2];
for (int i = 2; i < n; i++) {
while (top && (mult(p[i], res[top], res[top - 1])))
top--;
res[++top] = p[i];
}
int len = top;
res[++top] = p[n - 2];
for (int i = n - 3; i >= 0; i--) {
while (top != len && (mult(p[i], res[top], res[top - 1])))
top--;
res[++top] = p[i];
}
}

int calc(int i, int j) {
return (abs((int)ch[i].x + (int)ch[j].x) * abs((int)ch[i].y + (int)ch[j].y)) % P;
}

int dp(int l, int r) {
if (f[l][r]) return f[l][r];
if (r - l <= 2) return 0;
int ans = INF;
for (int i = l + 1; i < r; i++) {
ans = min(ans, dp(l, i) + dp(i, r) + c[l][i] + c[i][r]);
}
return f[l][r] = ans;
}

int main() {
int n;
while (~scanf("%d%d", &n, &P)) {
for (int i = 0; i < n; i++)
scanf("%lf%lf", &p[i].x, &p[i].y);
if (n <= 3) {
puts("0");
continue;
}
if (Graham(p, n, ch) < n)
puts("I can't cut.");
else {
memset(f, 0, sizeof(f));
for (int i = 0; i < n; i++)
for (int j = i + 2; j < n; j++)
c[i][j] = c[j][i] = calc(i, j);

printf("%d\n", dp(0, n - 1));
}
}
return 0;
}


• 本文已收录于以下专栏：

## zoj 3537 Cake 【凸包 + 区间dp】 【最优三角剖分】

Cake Time Limit: 1 Second      Memory Limit: 32768 KB You want to hold a party. Here's a pol...
• chenzhenyu123456
• 2015年12月23日 12:58
• 569

## zoj3537 Cake 【凸包+最优三角形划分】

• xinag578
• 2016年04月06日 11:29
• 564

## zoj3537cake【最优三角剖分+凸包】

• zhou_yujia
• 2016年02月19日 10:26
• 608

## zoj 3537 cake 切蛋糕 区间DP+凸包+递归 最优三角形剖分

• hcbbt
• 2013年12月05日 23:09
• 3083

## zoj 3537 cake （区间DP+凸包+ 最优三角形剖分模型）

• kalilili
• 2015年03月09日 20:54
• 624

## zoj 3537 Cake(区间DP+最优三角形剖分)待续

1、http://acm.zju.edu.cn/onlinejudge/showProblem.do?problemId=4472 2、题目大意; 给出一个多边形，将这个多边形用不相交的线段分割成一个...
• sdjzping
• 2014年02月10日 16:01
• 452

## Zoj 3537(区间DP)

//呜呜哈哈 我的区间DP的第一道题 // Graham这个凸包函数是直接复制网上的好像是求凸包的个数 //dp[i][j]表示从i到j这个多边形全部分成小三角形所花费的最小费用(不包括边i→j)...
• u011411189
• 2013年11月21日 13:57
• 625

## ZOJ3537 - Cake

Cake Time Limit: 1 Second      Memory Limit: 32768 KB You want to hold a party. Here's a polyg...
• mowayao
• 2014年04月20日 21:22
• 755

## Zoj 3537 Cake (DP最优三角形剖分)

• qq_21057881
• 2016年03月27日 20:58
• 176

## Zoj 3537 Cake (DP_最优三角形剖分)

• STILLxjy
• 2016年10月18日 15:24
• 167

举报原因： 您举报文章：zoj 3537 cake 切蛋糕 区间DP+凸包+递归 最优三角形剖分 色情 政治 抄袭 广告 招聘 骂人 其他 (最多只允许输入30个字)