Cake ZOJ - 3537 区间dp

最新推荐文章于 2021-09-11 20:52:38 发布

转载最新推荐文章于 2021-09-11 20:52:38 发布 · 248 阅读

区间dp 专栏收录该内容

14 篇文章

订阅专栏

本文介绍了一个关于凸包计算及区间动态规划的应用案例。通过分析一组点集，利用凸包算法筛选出边界上的点，并进一步运用区间DP解决点间的最小成本路径问题。代码中详细展示了如何实现这些算法步骤。

自己查错的时候的方向太逗比了。。

就是最后应该使用做完凸包后的res 点的坐标来弄。。

我却还是用原来的node

结果就会有点弄不到。。

我的debug应该要在关键位置才可以啊。

#include<iostream>
#include<string>
#include<cstring>
#include<algorithm>
#include<cmath>
#include<bits/stdc++.h>
using namespace std;
#define rep(x,n) for(ll (x)=0;(x)<(n);++(x))
#define rep1(x,n) for(ll (x)=1;(x)<=(n);++(x))
#define sf scanf
#define pf printf
#define mem(a,b) memset(a,b,sizeof(a));
const int maxn=305;
const int INF=1e9+7;
struct Point{
    double x,y;
}res[maxn],node[maxn];

double cross(Point a,Point b,Point c){
    return (a.x-c.x)*(b.y-c.y)-(a.y-c.y)*(b.x-c.x);
}
bool cmp(Point a,Point  b){
    if(a.x==b.x)return a.y<b.y;
    return a.x<b.x;
}
int convex(int n){
    int m=0;
    sort(node,node+n,cmp);
    for(int i=0;i<n;++i){
        while(m>1&&cross(res[m-1],node[i],res[m-2])<=0)
            m--;
        res[m++]=node[i];
    }
    int k=m;
    for(int i=n-2;i>=0;i--){
        while(m>k&&cross(res[m-1],node[i],res[m-2])<=0)
            m--;
        res[m++]=node[i];
    }
    if(n>1)m--;
    return m;
}

int n,p;
int dp[maxn][maxn];
int cost(int a,int  b){
    if(abs(a-b)<=1)return 0;
    return (abs((int)res[a].x+(int)res[b].x)*abs((int)res[a].y+(int)res[b].y))%p;
}
int main(){
    while(~sf("%d%d",&n,&p)){
        mem(dp,0);
        for(int i=0;i<n;++i){
            node[i].x=0;node[i].y=0;
            res[i].x=0;res[i].y=0;
        }
        int mini=0;
        double   miny=INF;
        for(int i=0;i<n;++i){
            sf("%lf%lf",&node[i].x,&node[i].y);
            if(node[i].y<miny){
                miny=node[i].y;
                mini=i;
            }
        }
        swap(node[0],node[mini]);
        sort(node,node+n,cmp);

        int m=convex(n);
        if(m<n)pf("I can't cut.\n");
        else{
            for(int step=3; step<n; step++) { // 区间dp，一般来说都可以用这种 增加步长 的方法
                for(int i=0; i+step<n; i++) {
                    int j = i+step;
                    dp[i][j]=INF;// 因为是求min，所以需要预设为 INF
                    for(int k=i+1; 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[0][n-1]);
            /*
            for(int s=2;s<n;++s){
                for(int i=0;i+s<n;++i){
                    int j=i+s;
                    dp[i][j]=INF;
                    for(int k=i+1;k<j;++k){
                        dp[i][j]=min(dp[i][j],dp[i][k]+dp[k][j]+cost(node[i],node[k])+cost(node[k],node[j]));
                    }
                }
            }*/
        }
    }
}