【BZOJ3594】方伯伯的玉米田(SCOI2014)-DP+二维树状数组

测试地址:方伯伯的玉米田
做法:本题需要用到DP+二维树状数组。
首先,我们发现每次拔高的区间都是一个后缀。这个自己画一画就大概能证出来了。
那么我们就有了一个状态定义:令 f(i,j) f ( i , j ) 为前 i i 个后缀被拔高了j次能获得的最长不下降子序列长度,有状态转移方程:
f(i,j)=1+max{f(p,q)|p<i,qj,apaijq} f ( i , j ) = 1 + max { f ( p , q ) | p < i , q ≤ j , a p − a i ≤ j − q }
我们发现 apaijq a p − a i ≤ j − q 可以写成 ap+qai+j a p + q ≤ a i + j ,那么一个状态 f(i,j) f ( i , j ) 就有了一个二维坐标 (j,ai+j) ( j , a i + j ) ,而我们要做的就是从小到大枚举 i i ,每次状态更新时只要找到一个点左下角所有状态的最大值即可,这个显然可以用二维树状数组维护,于是我们就完成了这一题,时间复杂度为O(nklognlogk)
以下是本人代码:

#include <bits/stdc++.h>
using namespace std;
int n,k,mxx,x[10010],mx[510][6010]={0};

int lowbit(int x)
{
    return x&(-x);
}

void modify(int x,int y,int c)
{
    for(int i=x;i<=k+1;i+=lowbit(i))
        for(int j=y;j<=mxx+k;j+=lowbit(j))
            mx[i][j]=max(mx[i][j],c);
}

int calc_max(int x,int y)
{
    int ans=0;
    for(int i=x;i;i-=lowbit(i))
        for(int j=y;j;j-=lowbit(j))
            ans=max(ans,mx[i][j]);
    return ans;
}

int main()
{
    scanf("%d%d",&n,&k);
    mxx=0;
    for(int i=1;i<=n;i++)
    {
        scanf("%d",&x[i]);
        mxx=max(mxx,x[i]);
    }

    for(int i=1;i<=n;i++)
        for(int j=k;j>=0;j--)
        {
            int now=calc_max(j+1,x[i]+j)+1;
            modify(j+1,x[i]+j,now);
        }
    printf("%d",calc_max(k+1,mxx+k));

    return 0;
}
  • 0
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

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

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

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值