咸鱼了好久...出来冒个泡_(:з」∠)_
题目连接:1107G - Vasya and Maximum Profit
题目大意:给出\(n,a\)以及长度为\(n\)的数组\(c_i\)和长度为\(n\)的严格单调上升数组\(d_i\),求\(\max\limits_{1 \le l \le r \le n} (a\cdot(r-l+1)-\sum_{i=l}^{r}c_i-gap(l,r))\),其中\(gap(l, r) = \max\limits_{l \le i < r} (d_{i + 1} - d_i)^2\)
题解:首先将所有的\(c_i\)转换为\(a-c_i\),这样就变成了求\(\max\limits_{1 \le l \le r \le n} (\sum_{i=l}^{r}c_i-gap(l,r))\)。如果\(l,r\)确定的话,我们就能通过求前缀和以及区间内最大值来算出该区间对应的答案,但我们还需要进一步的优化。
考虑每一个\(d_{i + 1} - d_i\)能成为\(gap(l,r)\)的范围,即在区间\([L,R]\)中,\(\forall L \le l \le r \le R,gap(l,r)\le d_{i + 1} - d_i\)。这样我们只需要用线段树查询区间\([L,R]\)的最大子段和就能求出当\(gap(l,r) \le d_{i + 1} - d_i\)时的答案。先预处理所有的\(L,R\),再扫一遍就好了。
![](https://i-blog.csdnimg.cn/blog_migrate/8f900a89c6347c561fdf2122f13be562.gif)
![](https://i-blog.csdnimg.cn/blog_migrate/961ddebeb323a10fe0623af514929fc1.gif)
#include<bits/stdc++.h> using namespace std; #define N 300001 #define LL long long LL n,b,l[N],r[N],a[N],d[N],L,R,M,S,ans; struct rua{LL l,r,w,s,lw,rw;}t[N<<2]; void up(int x,int mid) { t[x].s=t[x*2].s+t[x*2+1].s; t[x].w=max(t[x*2].w,t[x*2+1].w); t[x].lw=max(t[x*2].lw,t[x*2].s+t[x*2+1].lw); t[x].rw=max(t[x*2+1].rw,t[x*2+1].s+t[x*2].rw); t[x].w=max(t[x].w,t[x*2].rw+t[x*2+1].lw); } void Build(int l,int r,int x) { t[x].l=l,t[x].r=r; if(l==r){t[x].w=t[x].lw=t[x].rw=t[x].s=a[l];return;} int mid=l+r>>1; Build(l,mid,x*2); Build(mid+1,r,x*2+1); up(x,mid); } void ask(int ll,int rr,int l,int r,int x) { if(ll>r || l>rr)return; int mid=l+r>>1; if(ll<=l && r<=rr) { M=max(M,max(t[x].w,R+t[x].lw)); L=max(L,S+t[x].lw); R=max(R+t[x].s,t[x].rw); M=max(M,max(L,R)); S+=t[x].s; return; } ask(ll,rr,l,mid,x*2); ask(ll,rr,mid+1,r,x*2+1); } int main() { scanf("%I64d%I64d",&n,&b); for(LL i=1;i<=n;i++) { scanf("%I64d%I64d",&d[i],&a[i]); a[i]=b-a[i],ans=max(ans,a[i]); } for(LL i=n;i>=1;i--)d[i]-=d[i-1]; Build(1,n,1); d[1]=0; l[2]=2,r[n]=n; for(LL i=3;i<=n;i++) { LL _=i; while(_>2 && d[i]>=d[_-1]) _=l[_-1]; l[i]=_; } for(LL i=n-1;i>=2;i--) { LL _=i; while(_<n && d[i]>=d[_+1]) _=r[_+1]; r[i]=_; } for(LL i=2;i<=n;i++) { S=0; L=R=M=-(1e18); ask(l[i]-1,r[i],1,n,1); ans=max(ans,M-d[i]*d[i]); } printf("%I64d\n",ans); return 0; }