Cake ZOJ - 3537

题目链接

本题也是区间dp三角剖分板子题，只不过加入了判断是否是凸包，计算顺序要用凸包顺序(凸包板)

#include<bits/stdc++.h>
using namespace std;
#define lowbit(x) ((x)&(-x))
typedef long long LL;
typedef pair<int,int> pii;

const double eps = 1e-;
const int INF = 0x3f3f3f3f;

int dp[][];
int f[][];

int dcmp(double x) {
    if(fabs(x) < eps) return ;
    return x < ?-:;
}

struct Point {
    double x, y;
    Point friend operator - (Point a, Point b) {
        return {a.x-b.x, a.y-b.y};
    }
} points[], s[];

double crossProduct(Point a, Point b) {
    return a.x*b.y - a.y*b.x;
}

double dis(Point a, Point b) {
    Point c = a-b;
    return sqrt(c.x*c.x+c.y*c.y);
}

bool isConvexHull(int n) {
    sort(points+, points++n, [&](Point& a, Point& b) {
        return a.y < b.y || (a.y == b.y && a.x < b.x);
    });
    sort(points+, points++n, [&](Point& a, Point& b) {
        double x = crossProduct(a-points[], b-points[]);
        if(dcmp(x) >  || (dcmp(x) ==  && dcmp(dis(a, points[])-dis(a, points[]))<)) return true;
        return false;
    });
    s[] = points[], s[] = points[];
    int t = ;
    for(int i = ; i <= n; ++i) {
        while(t >=  && dcmp(crossProduct(s[t] - s[t-], points[i] - s[t-])) <= ) t--;
        s[++t] = points[i];
    }
    return t == n;
}

int cost(Point a, Point b, int p) {
    return (abs(int(a.x+b.x)) * abs(int(a.y+b.y))) % p;
}

void run_case() {
    int n, p;
    while(cin >> n >> p) {
        for(int i = ; i <= n; ++i) cin >> points[i].x >> points[i].y;
        bool flag = isConvexHull(n);
        if(!flag) {
            cout << "I can't cut.\n";
        } else {
            for(int i = ; i <= n; i++){
                for(int j = i+; j <= n; j++){
                    f[i][j] = f[j][i] = cost(s[i],s[j], p);
                }
            }
            for(int i = n-; i >= ; --i) {
                for(int j = i+; j <= n; ++j) {
                    dp[i][j] = INF;
                    for(int k = i+; k < j; ++k)
                        dp[i][j] = min(dp[i][j], dp[i][k]+dp[k][j]+f[i][k]+f[k][j]);
                }
            }
            cout << dp[][n] << "\n";
        }
    }

}

int main() {
    ios::sync_with_stdio(false), cin.tie();
    cout.flags(ios::fixed);cout.precision();
    //int t; cin >> t;
    //while(t--)
    run_case();
    cout.flush();
    return ;
}