bzoj3594&&方伯伯的玉米田

http://www.lydsy.com/JudgeOnline/problem.php?id=3594

考虑dp[i][j]表示以i结尾，拔高了j次的lis.

dp[i][j]=dp[k][l]+1枚举前一个最后以什么结尾，拔高了多少次。

这里需要满足a[i]+j-l>=a[k]

然后貌似就出来了：k<i,l<=j,a[k]+l<=a[i]+j. 活脱脱的一个三维偏序了。

但是这些值都不是定值（所以貌似CDQ不能做啊）

然后就上二维树状数组。

首先枚举了i保证了k<i，然后用树状数组维护(j,a[i]+j)这两维的lis.在前面找<=这个的，+1就是新的dp值。

然后树状数组下标不能为0，但是l和j可以为0，强制+1处理。

最后，貌似是由于要像背包那样只取一次计算dp值，所以内层循环要倒着枚举（不是很懂QAQ）

有点像dp[j][a[j]+j]=max(dp[j-xx][a[j]+j-xx]+1)        看成是带双条件的0,1背包，本来内两层都要倒着枚举，但是这里a[j]+j可以计算出来。

然后上代码：

#include<cstdio>
#include<iostream>
#include<cstring>
#include<cstdlib>
#include<algorithm>
#include<queue>
#include<cmath>
using namespace std;
const int maxn=10000+5;
int a[maxn];
int g[maxn+505][505];
int f[maxn+505][505];
int n,k,m;
int read()
{
int res=0;
int f=1;
char ch=getchar();
for(;ch>'9'||ch<'0';ch=getchar())
if(ch == '-') f = -1;
for(;ch>='0'&&ch<='9';ch=getchar())
res=res*10+ch-'0';
return res * f;
}
void updata(int x,int y,int val)
{
for(int i=x;i<=m;i+=i&-i)
{
for(int j=y;j<=k+1;j+=j&-j)
{
g[i][j]=max(g[i][j],val);
}
}
}
int query(int x,int y)
{
int ans=0;
for(int i=x;i;i-=i&-i)
{
for(int j=y;j;j-=j
4000
&-j)
{
ans=max(ans,g[i][j]);
}
}
return ans;
}
int main()
{
n=read(),k=read();
for(int i=1;i<=n;i++)
{
a[i]=read();
m=max(m,a[i]);
}
m+=k;
int ans=0;
for(int i=1;i<=n;i++)
{
for(int j=k;j>=0;j--)
{
f[i][j]=max(f[i][j],query(j+a[i],j+1)+1);
updata(j+a[i],j+1,f[i][j]);
ans=max(f[i][j],ans);
}
}
printf("%d\n",ans);
return 0;
}