bzoj 3594（树状数组优化dp）

传送门

题解：

设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;
}