bzoj 3594 [Scoi2014]方伯伯的玉米田

http://www.elijahqi.win/archives/3683
Description
方伯伯在自己的农田边散步，他突然发现田里的一排玉米非常的不美。
这排玉米一共有N株，它们的高度参差不齐。
方伯伯认为单调不下降序列很美，所以他决定先把一些玉米拔高，再把破坏美感的玉米拔除掉，使得剩下的玉米的高度构成一个单调不下降序列。
方伯伯可以选择一个区间，把这个区间的玉米全部拔高1单位高度，他可以进行最多K次这样的操作。拔玉米则可以随意选择一个集合的玉米拔掉。
问能最多剩多少株玉米，来构成一排美丽的玉米。

Input
第1行包含2个整数n，K，分别表示这排玉米的数目以及最多可进行多少次操作。
第2行包含n个整数，第i个数表示这排玉米，从左到右第i株玉米的高度ai。

Output
输出1个整数，最多剩下的玉米数。

Sample Input
3 1
2 1 3
Sample Output
3
HINT
1 < N < 10000，1 < K ≤ 500，1 ≤ ai ≤5000

Source
By 佚名提供

考虑一个i性质 就是每次拔高的区域一定是从区间右端点开始

因为选择对一段区域拔高之后 左区域+这段的最长不降可能增加

而这段+右区域 不增且可能减少

那么设dp[i][j]表示考虑完前i个现在在做第i个 拔高j次之后能得到的最长不降子序列

那么用二维树状数组维护 高度<=h[i]+j且拔高次数<=j 的最大值即可 然后每次二维区域选最大更新

#include<cstdio>
#include<cctype>
#include<algorithm>
using namespace std;
inline char gc(){
    static char now[1<<16],*S,*T;
    if(T==S){T=(S=now)+fread(now,1,1<<16,stdin);if (T==S) return EOF;}
    return *S++;
}
inline int read(){
    int x=0,f=1;char ch=gc();
    while(!isdigit(ch)) {if (ch=='-') f=-1;ch=gc();}
    while(isdigit(ch)) x=x*10+ch-'0',ch=gc();
    return x*f;
}
const int N=11000;
int mx,n,k,h[N],s[N][510];
inline void add(int x,int y,int v){
    for (int i=x;i<=mx+k;i+=i&-i)
        for (int j=y;j<=k+1;j+=j&-j) s[i][j]=max(s[i][j],v);
}
inline int query(int x,int y){
    static int tmp;tmp=0;
    for (int i=x;i;i-=i&-i)
        for (int j=y;j;j-=j&-j) tmp=max(tmp,s[i][j]);return tmp;
}
int main(){
    freopen("bzoj3594.in","r",stdin);
    n=read();k=read();int ans=0;
    for (int i=1;i<=n;++i) h[i]=read(),mx=max(h[i],mx);
    for (int i=1;i<=n;++i){
        for (int j=k;~j;--j){
            int p=query(h[i]+j,j+1);
            ans=max(ans,p+1);
            add(h[i]+j,j+1,p+1);
        }
    }printf("%d\n",ans);
    return 0;
}