正题

题目链接:https://uoj.ac/problem/749

---

题目大意

如果有序列$a$，你每次取走一个数字后然后这个序列最前面的数字会被别人取走，直到序列为空。此时$f(a)$表示你最大能取走的权值和。

给出一个长度为$n$的序列$a$，$q$次询问区间$[l,r]$，求$f(a_{l\sim r})$。

$1\le n,q\le 2\times 10^5,1\le a_i\le 10^9$

---

解题思路

考虑一下最优的取法，我们的限制其实相当于前$2i-1$个数之中不能取走超过$i$个数。

一个暴力的贪心想法是不停往前走，用一个堆维护现在取了的数，走到不是$2i-1$时我们考虑是否拿这个$a_i$替代堆里最小的数，走到$2i-1$时直接丢进堆里就行了。

或者反过来更简单，走到一个数就丢进堆里，遇到$2i-1$时直接取堆中最大的数就好了。

我们现在能处理左或右端点固定的所有答案了，我们考虑分治，问题在于怎么去合并两个区间的答案。

为了防止大量的分类讨论，对于长度为奇数的区间询问，我们可以考虑去掉区间的末尾，因为这个位置肯定会取到。

然后我们考虑怎么合并区间，显然是前面一些选了的数字变成不选，后面一些不选的数字变成选了，显然是前面拿小的，后面拿大的。

我们考虑在分治的时候对于区间$[L,mid,R]$，对于$[L,mid]$部分，我们用主席树维护选了的数字，对于$[mid+1,R]$部分，我们用主席树维护没选的数字，然后对于一个询问，二分一个$k$，然后左边的取出前$k$小，右边的取出前$k$大进行比较即可。

时间复杂度：$O(n{\log }^{2}n)$

---

code

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<queue>
#define ll long long
using namespace std;
const ll N=2e5+10,M=N<<5;
ll n,m,cnt,a[N],b[N],l[N],r[N],ans[N],s[N],rt[N];
priority_queue<ll> q;
struct SegTree{ll cnt,w[M],c[M],ls[M],rs[M];ll Change(ll x,ll L,ll R,ll pos,ll val){ll p=++cnt;w[p]=w[x]+val;c[p]=c[x]+val*b[pos];if(L==R)return p;ll mid=(L+R)>>1;if(pos<=mid)ls[p]=Change(ls[x],L,mid,pos,val),rs[p]=rs[x];else rs[p]=Change(rs[x],mid+1,R,pos,val),ls[p]=ls[x];return p;}ll Find(ll x,ll L,ll R,ll k){if(L==R)return L;ll mid=(L+R)>>1;if(w[ls[x]]>=k)return Find(ls[x],L,mid,k);return Find(rs[x],mid+1,R,k-w[ls[x]]);}ll Ask(ll x,ll L,ll R,ll k){if(!k)return 0;if(L==R)return b[L]*k;ll mid=(L+R)>>1;if(w[ls[x]]>=k)return Ask(ls[x],L,mid,k);return Ask(rs[x],mid+1,R,k-w[ls[x]])+c[ls[x]];}
}T;
bool cmp(ll x,ll y){return l[x]<l[y];}
void solve(ll L,ll R,vector<ll> &now){if(L==R)return;ll mid=(L+R)>>1;vector<ll> pl,pr,p;for(ll i=0;i<now.size();i++)if(l[now[i]]<=mid&&r[now[i]]>mid)p.push_back(now[i]);else if(l[now[i]]<=mid)pl.push_back(now[i]);else pr.push_back(now[i]);solve(L,mid,pl);solve(mid+1,R,pr);if(!p.size())return;for(ll g=0;g<2;g++){T.cnt=rt[mid]=s[mid]=0;for(ll i=mid+1;i<=R;i++){if(((i-mid)&1)==g)q.push(-a[i]),s[i]=s[i-1]+b[a[i]],rt[i]=rt[i-1];else{if(q.size()&&a[i]>-q.top()){s[i]=s[i-1]-b[-q.top()]+b[a[i]];rt[i]=T.Change(rt[i-1],1,cnt,-q.top(),1);q.pop();q.push(-a[i]);}else rt[i]=T.Change(rt[i-1],1,cnt,a[i],1),s[i]=s[i-1];}}while(!q.empty())q.pop();for(ll i=mid;i>=L;i--){q.push(a[i]);if(((mid-i)&1)==g){rt[i]=T.Change((i==mid)?0:rt[i+1],1,cnt,q.top(),1);s[i]=((i!=mid)?s[i+1]:0)+b[q.top()];q.pop();}else if(i!=mid)rt[i]=rt[i+1],s[i]=s[i+1];}while(!q.empty())q.pop();for(ll i=0;i<p.size();i++){if(((mid-l[p[i]])&1)!=g)continue;ll X=T.w[rt[l[p[i]]]],Y=T.w[rt[r[p[i]]]];ll _l=1,_r=min(X,Y);if(p[i]==2)i++,i--;while(_l<=_r){ll _mid=(_l+_r)>>1;if(T.Find(rt[l[p[i]]],1,cnt,_mid)<T.Find(rt[r[p[i]]],1,cnt,Y-_mid+1))_l=_mid+1;else _r=_mid-1;}ans[p[i]]+=s[l[p[i]]]-T.Ask(rt[l[p[i]]],1,cnt,_r);ans[p[i]]+=s[r[p[i]]]+T.Ask(rt[r[p[i]]],1,cnt,Y)-T.Ask(rt[r[p[i]]],1,cnt,Y-_r);}}return;
}
signed main()
{scanf("%lld%lld",&n,&m);for(ll i=1;i<=n;i++)scanf("%lld",&a[i]),b[i]=a[i];sort(b+1,b+1+n);cnt=unique(b+1,b+1+n)-b-1;for(ll i=1;i<=n;i++)a[i]=lower_bound(b+1,b+1+cnt,a[i])-b;vector<ll> v;for(ll i=1;i<=m;i++){scanf("%lld%lld",&l[i],&r[i]);v.push_back(i);if((r[i]-l[i]+1)&1)ans[i]+=b[a[r[i]]],r[i]--;}sort(v.begin(),v.end());solve(1,n,v);for(ll i=1;i<=m;i++)printf("%lld\n",ans[i]);return 0;
}