【jzoj5347】【NOIP2017提高A组模拟9.5】【遥远的金字塔】【斜率优化动态规划】

description

这里写图片描述

solution

考虑动态规划,设f[i][j]表示前i个,用了j个矩阵的最大答案。 f[i][j]=maxk<=if[k][j1]+a[i](ik) ,其中a[i]表示第i段的长度。可以发现对于i的两个决策点 x,yx<y ,当 (f[y]f[x])/(yx)<a[i] 时x是比y更优的。这样就可以维护一个斜率单调递减的单调队列优化动态规划。

code

#include<cstdio>
#include<cmath>
#include<cstring>
#include<algorithm>
#define LF double
#define LL long long
#define ULL unsigned int
#define fo(i,j,k) for(int i=j;i<=k;i++)
#define fd(i,j,k) for(int i=j;i>=k;i--)
#define fr(i,j) for(int i=begin[j];i;i=next[i])
using namespace std;
int const mn=2*1e4+9;LL inf=1e18+7;
int t,n,K,q[mn];
LL f[mn],g[mn],a[mn];
int main(){
    freopen("pyramid.in","r",stdin);
    freopen("pyramid.out","w",stdout);
    scanf("%d%d",&n,&K);
    fo(i,1,n){
        int x,y;
        scanf("%d%d",&x,&y);
        a[i]=y-x+1;
    }
    fo(i,1,n)f[i]=-inf;
    fo(cas,1,K){
        int he=0,ti=0;
        q[0]=0;
        fo(i,1,n){
            while((he<ti)&&(f[q[he+1]]-f[q[he]]>a[i]*(q[he+1]-q[he])))he++;
            g[i]=f[q[he]]+a[i]*(i-q[he]);
            if(i!=n){
                while((he<ti)&&((f[q[ti]]-f[q[ti-1]])*(i-q[ti])<(f[i]-f[q[ti]])*(q[ti]-q[ti-1])))ti--;
                q[++ti]=i;
            }
        }
        fo(i,0,n)f[i]=g[i];
    }
    LL ans=0;
    fo(i,0,n)ans=max(ans,f[i]);
    printf("%lld",ans);
    return 0;
}
  • 0
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值