Cross the Wall(斜率 dp,模板)

Cross the Wall

#include <cstring>
#include <iostream>
#include <cstdio>
#include <algorithm>

using namespace std;
#define ll long long
#define maxn 50010
ll dp[maxn][110], q[maxn];
struct node {
    ll w, h;
} peo[maxn];

bool cmp(node a, node b) {
    if (a.h == b.h) {
        return a.w > b.w;
    } else {
        return a.h > b.h;
    }
}

ll gety(int i, int j, int k) {
    return dp[j][k] - dp[i][k];
}

ll getx(int x1, int x2) {
    return peo[x1 + 1].h - peo[x2 + 1].h;
}

int main() {
    int n, k;
    while (scanf("%d%d", &n, &k) != EOF) {
        for (int i = 1; i <= n; i++) {
            scanf("%lld%lld", &peo[i].w, &peo[i].h);
        }
        sort(peo + 1, peo + n + 1, cmp);
        int num = 1;
        for (int i = 1; i <= n; i++) {
            if (peo[i].w > peo[num].w) {
                peo[++num] = peo[i];
            }
        }
        n = num;
        k = min(k, n);
        for (int i = 1; i <= n; i++) {
            dp[i][1] = peo[1].h * peo[i].w;
        }
        for (int j = 2; j <= k; j++) {
            int head = 0, tail = 0;
            q[tail++] = 0;
            for (int i = 1; i <= n; i++) {
                while (head + 1 < tail && gety(q[head], q[head + 1], j - 1) <= peo[i].w * getx(q[head], q[head + 1])) {
                    head++;
                }
                dp[i][j] = dp[q[head]][j - 1] + peo[q[head] + 1].h * peo[i].w;
                while (head + 1 < tail && gety(q[tail - 1], i, j - 1) * getx(q[tail - 2], q[tail - 1]) <=
                                          gety(q[tail - 2], q[tail - 1], j - 1) * getx(q[tail - 1], i)) {
                    tail--;
                }
                q[tail++] = i;
            }
        }
        ll ans = dp[n][1];
        for (int i = 2; i <= k; i++) {
            ans = min(ans, dp[n][i]);
        }
        printf("%lld\n", ans);
    }
    return 0;
}
  • 0
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值