JZOJ 5451 Genocide

Cenocide

Description

给出一个长度为$n$的序列$a$，$a_i$表示选择第$i$个数的代价，设$c_i$表示是否选择第$i$个数，（$0$为不选，$1$为选择），对于某种选择方案，定义总收益为

接下来有$m$个道具，第$i$个道具能够将$a_{P_i}$修改成$X_i$。
对于每个道具，询问仅使用这个道具的情况下，最大的总收益为多少。

Data Constraint

Solution

考虑$m$=$1$时怎么做。
预处理出$a$的前缀和序列$q$。
设$f_i$表示$1$~$i$的数已处理完的最大收益，考虑使用斜率优化。
若由$j$转移到$i$比由$k$转移到$i$更优（$j$ < $k$），设仅跟$j$相关的项的和为$F_j^{'}$($F_j^{'}$=$f_j$+${j*(j-1)}\over 2$+$q_j$)，则可以得出下面的式子:${{F_j^{'}-F_k^{'}}\over{j-k}}<i$
用一个单调栈维护斜率单调递减即可。
由于第$i$个数可以不选择，因而$f_i$要与$f_{i-1}$取$min$。

现在考虑$m$<=$3*10^5$时怎么做。
同样的，我们设$g_j$表示$j$~$n$的数已处理完的最大收益，用上述同样的方法可以求得。
对于每一个$P_i$，若不选择这一个数，那么答案显然为$f_{P_i-1}$+$g_{P_i+1}$。
设$h_i$表示必须选择$i$时整个序列的最大收益，那么选择$P_i$这个数的收益即为$h_{P_i}$+$a_{P_i}$-$X_i$。
现在考虑如何求出$h$。
对于一个$h_x$，$h_x$=$max$($f_i$+$g_j$-$q_{j-1}$+$q_i$+${(j-i)*(j-i-1)}\over 2$)（$i$ < $x$ < $j$）

令F_i=$f_i$+${i*(i+1)}\over 2$+$q_i$,$G_i$=$g_i$+${i*(i-1)}\over 2$-$q_{i-1}$
则$h_x$=$max$($F_i$+$G_j$-$i*j$)（$i$ < $x$ < $j$）

考虑分治，对于分治区间（$l$ ,$r$ ），考虑$j$在$mid$右边，$i$在$mid$左边的情况来更新$x$在($l$,$mid$)内的$h$。
设$W_i$=$F_i$+$max$（$G_j$-$i$*$j$）,那么$h_x$每一次都可以由$W_i$更新到（$l$<$i$<$x$），做个前缀最大即可。

那现在考虑如何求出$W_i$，可以发现随着$i$的递增，决策点$j$越来越前，考虑对右区间维护一个斜率单调递减的单调栈（这里就不推式子来证明了，自行手推吧），随着$i$的增加，不断将栈顶元素退掉（当栈顶的斜率<$i$时便可退掉栈顶，这个自己推一下便可以发现了）,如此这般便能求出对于每个$i$的最优决策点了。
仔细观察便能发现这样做并不能让所有的情况更新到$h_x$，将$f$和$g$对调后再翻转，$a$也翻转，之后再做一遍分治所有的情况便都能更新到每一个$h_x$了。

Code

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>

#define fo(i,j,l) for(int i=j;i<=l;i++)
#define fd(i,j,l) for(int i=j;i>=l;i--)

using namespace std;
typedef long long ll;
typedef double db;
const ll N=66e4,maxn=1e12;
ll a[N],b[N],f[N],g[N],h[N],t[N],dd[N];
int top,n,m,j,k,l,i,c; 
ll F[N],G[N],d[N],bb;

double xl(ll a,ll b)
{return (F[a]-F[b])/(a-b*1.0);}

double xx(ll a,ll b)
{return (G[a]-G[b])/(a-b*1.0);}

ll max(ll a,ll b)
{return a>b?a:b;}

void ef(int l,int r,int f)
{
    if(l==r)return;
    int k=0,mid=(l+r)/2;
    fo(i,mid+1,r){
        while(k-1&&xx(d[k],i)>=xx(d[k-1],d[k]))k--;
        d[++k]=i;
    }
    t[l-1]=-maxn;
    fo(i,l,mid){
        while(k-1&&xx(d[k-1],d[k])<=i)k--;
        t[i]=max(t[i-1],F[i]+G[d[k]]-i*d[k]);
    }
    if(f==0)fo(i,l,mid)h[i]=max(h[i],t[i-1]);
    else fo(i,l,mid)h[n-i+1]=max(h[n-i+1],t[i-1]);
    ef(l,mid,f); ef(mid+1,r,f);
}

int main()
{
    cin>>n;
    fo(i,1,n)scanf("%lld",&a[i]),b[i]=b[i-1]+a[i];
    int t;
    d[t=1]=0;
    fo(i,1,n){
        while(t-1&&xl(d[t-1],d[t])<i)t--;
        f[i]=max(f[i-1],f[d[t]]+(ll)(i-d[t])*(i-d[t]+1)/2-b[i]+b[d[t]]);
        F[i]=f[i]+(ll)(i-1)*(ll)i/2+b[i];
        while(t-1&&xl(d[t],i)>=xl(d[t-1],d[t]))t--;
        d[++t]=i;
    }
    d[t=1]=0;
    fo(i,1,n)b[i]=b[i-1]+a[n-i+1];
    fo(i,1,n){
        while(t-1&&xl(d[t-1],d[t])<=i)t--;
        g[i]=max(g[i-1],g[d[t]]+(ll)(i-d[t])*(i-d[t]+1)/2-b[i]+b[d[t]]);
        F[i]=g[i]+(ll)(i-1)*(ll)i/2+b[i];
        while(t-1&&xl(d[t],i)>=xl(d[t-1],d[t]))t--;
        d[++t]=i;
    }
    fo(i,1,n/2)swap(g[i],g[n-i+1]);
    fo(i,1,n+1)dd[i]=f[i-1]+g[i+1],b[i]=b[i-1]+a[i];
    fo(i,0,n+1)F[i]=f[i]+(ll)i*(i+1)/2+b[i],G[i]=g[i]+(ll)i*(i-1)/2-b[i-1];
    fo(i,1,n)h[i]=-maxn;
    ef(0,n+1,0);
    fo(i,1,n+1)swap(f[i],g[i]),b[i]=b[i-1]+a[n-i+1];
    fo(i,1,n/2)swap(f[i],f[n-i+1]),swap(g[i],g[n-i+1]);
    fo(i,0,n+1)F[i]=f[i]+(ll)i*(i+1)/2+b[i],G[i]=g[i]+(ll)i*(i-1)/2-b[i-1];
    ef(0,n+1,1); 
    cin>>m;
    fo(i,1,m){
        scanf("%d%lld",&c,&bb);
        ll ans=max(dd[c],h[c]-(bb-a[c]));
        printf("%lld\n",ans);
    }
}