题解:
设dp[i][j]表示前i个玉米用j次操作的LIS长度。
dp[i][j]=max{dp[x][y]}+1(x<i,y≤j,a[x]+y≤a[i]+j)
维护一个关于(j,max{a[i]}+j)的树状数组,每次倒序枚举j防止dp[i]内部自己更新自己,如果正序枚举可能发现这种情况:dp[i][k]更新到后来dp[i][j](k<j)。
由dp转移式知dp[i][...]只能被dp[比i小的数][...]更新到。
#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
using namespace std;
int n,m,mx=0,ans=0;
int dp[10004][504],c[5504][504],a[10004];
inline int read() {
int x=0;char c=getchar();
while (c<'0'||c>'9') c=getchar();
while (c>='0'&&c<='9') x=x*10+c-'0',c=getchar();
return x;
}
inline void modify(int x,int y,int v) {
for (int i=x;i<=mx+m;i+=i&-i)
for (int j=y;j<=m+1;j+=j&-j)
c[i][j]=max(c[i][j],v);
}
inline int query(int x,int y) {
int ret=0;
for (int i=x;i;i-=i&-i)
for (int j=y;j;j-=j&-j)
ret=max(ret,c[i][j]);
return ret;
}
int main() {
n=read(),m=read();
for (int i=1;i<=n;++i) a[i]=read(),mx=max(mx,a[i]);
for (int i=1;i<=n;++i) {
// for (int j=0;j<=m;++j) {
for (int j=m;~j;--j) {//因为j可能为0,所以BIT定义域强制加1
dp[i][j]=query(a[i]+j,j+1)+1;
ans=max(ans,dp[i][j]);
modify(a[i]+j,j+1,dp[i][j]);
}
}
printf("%d\n",ans);
return 0;
}