P2120 [ZJOI2007]仓库建设(dp+斜率优化)

本文深入探讨了动态规划中的一种优化技巧——斜率优化。通过引入前缀和和斜率优化,将暴力DP的时间复杂度从O(n^2)降低到O(nlogn)。文章详细解析了如何通过化简不等式来确定决策的优劣,并提供了一段清晰的C++代码实现。

摘要生成于 C知道 ,由 DeepSeek-R1 满血版支持, 前往体验 >

思路

首先暴力DP显然,可以得20分

加上一个前缀和优化,可以得到40分

然后上斜率优化

\(sum_i\)\(\sum_{1}^iP_i\)\(sump_i\)\(\sum_{1}^{i}P_i\times Pos_i\)

则决策j优于决策i的条件可以表示为

\[ dp_j+C_t+pos_t \times(sum_t-sum_{j}) - (sump_t - sump_{j})<dp_i+C_t+pos_t \times(sum_t-sum_{i}) - (sump_t - sump_{i}) \]

化简即可

这里贴上我认为写的不错的一篇题解其实是懒得写式子了

dummy的题解

代码
#include <cstdio>
#include <cstring>
#include <algorithm>
#define int long long
using namespace std;
int dp[1000010],n,c[1000010],p[1000010],pos[1000010],sum[1000010],summ[1000010];
int queue[1000010],h,t;
int T(int x){
    return dp[x]+summ[x];
}
int S(int x){
    return sum[x];
}
signed main(){
    // freopen("6.in","r",stdin);
    memset(dp,0x3f,sizeof(dp));
    scanf("%lld",&n);
    for(int i=1;i<=n;i++){
        scanf("%lld %lld %lld",&pos[i],&p[i],&c[i]);
        sum[i]=sum[i-1]+p[i];
        summ[i]=summ[i-1]+p[i]*pos[i];
    }
    dp[0]=0;
    // dp[1]=c[1];
    for(int i=1;i<=n;i++){
        while(h<t&&1.0*(T(queue[h+1])-T(queue[h]))/((S(queue[h+1])-S(queue[h])))<pos[i])
            h++;
        dp[i]=dp[queue[h]]+pos[i]*(sum[i]-sum[queue[h]])-(summ[i]-summ[queue[h]])+c[i];
        while(h<t&&1.0*(T(i)-T(queue[t]))/((S(i)-S(queue[t])))<1.0*(T(queue[t-1])-T(queue[t]))/((S(queue[t-1])-S(queue[t]))))
            t--;
        queue[++t]=i;
    }
    printf("%lld\n",dp[n]);
    return 0;
}

转载于:https://www.cnblogs.com/dreagonm/p/10129500.html

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值