5.4 省选模拟赛 修改 线段树优化dp 线段树上二分

LINK：修改

题面就不放了 大致说一下做法。不愧是dls出的题 以前没见过这种类型的 不过还是自己dp的时候写丑了。

从这道题中得到一个结论 dp方程要写的优美一点 不过写的过丑 优化都优化不了。

容易想到 f[i][j]表示前i个数最大值为aj的最大收益。

那么有(j<=a_i,f[i][j]=f[i-1][k]-a_i+j+b_i.j>a_i,f[i][j]=f[i-1][j])

值得注意的是这个转移不完全 在第二个转移的式子中 决策不全面 强行利用f[i][j] 的单调性进行覆盖 那么第一个决策中k就可以变成j 且第二个转移也是正确的。

const ll MAXN=5010,N=40;
ll n,ans;
ll a[MAXN],b[MAXN];
ll f[MAXN][MAXN];
struct wy{int id,x,y;}t[MAXN];
inline int cmp(wy a,wy b){return a.x<b.x;}
signed main()
{
	//freopen("1.in","r",stdin);
	get(n);
	rep(1,n,i)get(a[i]),t[i].x=a[i],t[i].id=i;
	rep(1,n,i)get(b[i]),t[i].y=b[i];
	sort(t+1,t+1+n,cmp);
	rep(1,n,i)
	{
		rep(1,n,j)
		{
			f[i][j]=-INF;
			if(a[j]<=a[i])f[i][j]=max(f[i][j],f[i-1][j]-a[i]+a[j]+b[i]);
			else f[i][j]=max(f[i][j],f[i-1][j]);
		}
		rep(2,n,j)f[i][t[j].id]=max(f[i][t[j-1].id],f[i][t[j].id]);
	}
	rep(1,n,i)ans=max(ans,f[n][i]); 
	putl(ans);return 0;
}

值得注意的是 需要赋初值.

这个dp 可以发现 非常不可优化 所以我就弃疗了。

原因是 状态转移方程写丑了 少考虑了情况 导致f值不再单调 需要后来的 单调覆盖。

观察第一个式子 容易想到 如果出现覆盖的情况 那么说明大的j没有前面(b_i)的贡献才会这样。

考虑后面的式子 其实是可以加入这个决策转移的 (f[i][j]=max(f[i-1][j],f[i-1][i]+b_i))

这样f数组显然 单调 那么久不需要单调覆盖了。

具体观察两个dp式子:(j>a_i,f[i][j]=max(f[i-1][j],f[i-1][i]+b_i),j<=a_i,f[i][j]=f[i-1][k]-a_i+j+b_i)

第二个转移 显然是一个加上常数 加上一个端点值j的操作 第一个操作由于f值单调 所以这是一个区间赋值操作。

可以发现没有区间最值 所以线段树可以维护这个转移。至于区间赋值操作 显然可以在查到值之后进行线段树上二分寻找区间即可。

注意 每次修改后注意pushup某个区间的左端点的值 为线段树上二分做准备。

const int MAXN=500010,N=40;
int n,num,cnt;ll ww;
int a[MAXN],b[MAXN],c[MAXN];
struct wy
{
	int l,r;
	ll sum;
	ll tag,tag1,tag2;
}t[MAXN<<2];
inline void discrete()
{
	sort(c+1,c+1+n);
	rep(1,n,i)if(i==1||c[i]!=c[i-1])c[++cnt]=c[i];
	rep(1,n,i)a[i]=lower_bound(c+1,c+1+cnt,a[i])-c;
}
inline void build(int p,int l,int r)
{
	l(p)=l;r(p)=r;sum(p)=0;
	tag2(p)=-INF;
	if(l==r)return;
	int mid=(l+r)>>1;
	build(zz,l,mid);
	build(yy,mid+1,r);
}
inline void pushdown(int p)
{
	if(tag2(p)!=-INF)
	{
		sum(zz)=tag2(p);
		sum(yy)=tag2(p);
		tag2(zz)=tag2(p);
		tag2(yy)=tag2(p);
		tag1(zz)=tag1(yy)=0;
		tag(zz)=tag(yy)=0;
		tag2(p)=-INF;
		//return;
	}
	sum(zz)+=tag(p);
	sum(yy)+=tag(p);
	tag(zz)+=tag(p);
	tag(yy)+=tag(p);
	int mid=(l(p)+r(p))>>1;
	sum(zz)+=tag1(p)*c[l(p)];
	sum(yy)+=tag1(p)*c[mid+1];
	tag1(zz)+=tag1(p);
	tag1(yy)+=tag1(p);
	tag(p)=tag1(p)=0;
}
inline void change(int p,int l,int r,ll v)
{
	if(l<=l(p)&&r>=r(p))
	{
		sum(p)+=v;tag(p)+=v;
		sum(p)+=c[l(p)];
		++tag1(p);
		return;
	}
	pushdown(p);
	int mid=(l(p)+r(p))>>1;
	if(l<=mid)change(zz,l,r,v);
	if(r>mid)change(yy,l,r,v);
	sum(p)=sum(zz);
}
inline ll ask(int p,int x)
{
	if(l(p)==r(p))return sum(p);
	int mid=(l(p)+r(p))>>1;
	pushdown(p);
	if(x<=mid)return ask(zz,x);
	return ask(yy,x);
}
inline int ask1(int p,int l)
{
	if(l(p)>=l&&sum(p)>ww)return 0;
	if(l(p)==r(p))return l(p);
	int mid=(l(p)+r(p))>>1;
	pushdown(p);
	if(l>mid)return ask1(yy,l);
	int w=ask1(yy,l);
	if(!w)return ask1(zz,l);
	return w;
}
inline void modify(int p,int l,int r,ll v)
{
	if(l<=l(p)&&r>=r(p))
	{
		sum(p)=v;
		tag1(p)=tag(p)=0;
		tag2(p)=v;
		return;
	}
	pushdown(p);
	int mid=(l(p)+r(p))>>1;
	if(l<=mid)modify(zz,l,r,v);
	if(r>mid)modify(yy,l,r,v);
	sum(p)=sum(zz);
}
signed main()
{
	//freopen("1.in","r",stdin);
	get(n);
	rep(1,n,i)get(a[i]),c[i]=a[i];
	rep(1,n,i)get(b[i]);
	discrete();
	build(1,1,cnt);
	rep(1,n,i)
	{
		change(1,1,a[i],-c[a[i]]+b[i]);
		ww=ask(1,a[i]);
		int r=ask1(1,a[i]);
		//cout<<r<<endl;
		modify(1,a[i],r,ww);
	}
	ll ans=0;
	rep(1,cnt,i)ans=max(ans,ask(1,i));
	putl(ans);return 0;
}