[Educational Round 59][Codeforces 1107G. Vasya and Maximum Profit]

咸鱼了好久...出来冒个泡_(:з」∠)_

题目连接：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\)，再扫一遍就好了。

#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<<];
void up(int x,int mid)
{
    t[x].s=t[x*].s+t[x*+].s;
    t[x].w=max(t[x*].w,t[x*+].w);
    t[x].lw=max(t[x*].lw,t[x*].s+t[x*+].lw);
    t[x].rw=max(t[x*+].rw,t[x*+].s+t[x*].rw);
    t[x].w=max(t[x].w,t[x*].rw+t[x*+].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>>;
    Build(l,mid,x*);
    Build(mid+,r,x*+);
    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>>;
    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*);
    ask(ll,rr,mid+,r,x*+);
}
int main()
{
    scanf("%I64d%I64d",&n,&b);
    for(LL i=;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>=;i--)d[i]-=d[i-];
    Build(,n,);
    d[]=;
    l[]=,r[n]=n;
    for(LL i=;i<=n;i++)
      {
      LL _=i;
      while(_> && d[i]>=d[_-])
        _=l[_-];
      l[i]=_;
      }
    for(LL i=n-;i>=;i--)
      {
      LL _=i;
      while(_<n && d[i]>=d[_+])
        _=r[_+];
      r[i]=_;
      }
    for(LL i=;i<=n;i++)
      {
      S=;
      L=R=M=-(1e18);
      ask(l[i]-,r[i],,n,);
      ans=max(ans,M-d[i]*d[i]);
      }
    printf("%I64d\n",ans);
    return ;
}