CF1500E Subset Trick
XLVII.CF1500E Subset Trick
考虑对于每个集合大小 \(i\)，找到所有大小为 \(i\) 的集合中元素和最小的一个 \(l_i\) 与最大的一个 \(r_i\)。则，所有 \(x\in[l_i,r_i)\) 均不合法。
于是我们就要求 \(\Big|\bigcup\limits_{i=0}^{n}[l_i,r_i)\Big|\)。
直接求不容易；但是我们可以正难则反，考虑找出那些不被包含在其中的 \(x\)。虽然这一 \(x\) 的数量是无限的，但是发现若设所有数之和为 \(s\)，则所有 \(x\in[s,\infty)\) 显然均不被包含，可以不考虑。于是我们仅需考虑那些 \(\in[0,s)\) 且不被包含的数即可。
显然的是，\(l_i\) 必然是前 \(i\) 小的数之和，\(r_i\) 必然是前 \(i\) 大的数之和。
假如我们在平面直角坐标系中画出所有 \((i,l_i)\) 与 \((i,r_i)\) 的话，就会发现所有 \(l\) 连成的折线是凹的，所有 \(r\) 连成的这些是凸的，且这两条折线是中心对称的。
于是我们仅需考虑前一半中 \(x\) 的数量，然后乘二就得到了全体 \(x\) 的数量（当然，若 \(x\) 是奇数，还要额外考虑中间的那个位置，但是随便特判一下就出来了）。
现在考虑每个不被包含在其中的 \(x\) 的区间，发现其都是形如 \([r_{i-1},l_i)\) 的形式。
一个非常棒的观察是，因为 \(l\) 凹 \(r\) 凸且二者对称，所以在前一半中 \(r\) 的增量总是不小于 \(l\) 的增量，这也意味着上述区间若非空则所有这样的 \(i\) 是一段前缀。于是我们可以用什么数据结构维护一下 \(l\) 与 \(r\)，然后套个二分就能找到最大的这样的 \(i\)。
现在考虑用什么数据结构。平衡树显然是可行的，但是据说被卡常了。于是离散化后上线段树维护即可。
时间复杂度 \(O(n\log^2n)\)。那个二分不能放到线段树上实现，因为要维护的是两棵不同的线段树。
代码：
#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
int n,m,lim,sz;
ll a[200100],b[200100],sum;
vector<ll>v;
struct ST{
#define lson x<<1
#define rson x<<1|1
#define mid ((l+r)>>1)
struct SegTree{int sz;ll sum,tag;}seg[1600100];
void ADD(int x,ll y){seg[x].tag+=y,seg[x].sum+=seg[x].sz*y;}
void pushdown(int x){ADD(lson,seg[x].tag),ADD(rson,seg[x].tag),seg[x].tag=0;}
void pushup(int x){seg[x].sum=seg[lson].sum+seg[rson].sum,seg[x].sz=seg[lson].sz+seg[rson].sz;};
void turnon(int x,int l,int r,int P){if(l>P||r<P)return;if(l==r)seg[x].sum=seg[x].tag,seg[x].sz=1;else pushdown(x),turnon(lson,l,mid,P),turnon(rson,mid+1,r,P),pushup(x);}
void turnoff(int x,int l,int r,int P){if(l>P||r<P)return;if(l==r)seg[x].sum=0,seg[x].sz=0;else pushdown(x),turnoff(lson,l,mid,P),turnoff(rson,mid+1,r,P),pushup(x);}
void modify(int x,int l,int r,int L,int R,ll val){if(l>R||r<L)return;if(L<=l&&r<=R)ADD(x,val);else pushdown(x),modify(lson,l,mid,L,R,val),modify(rson,mid+1,r,L,R,val),pushup(x);}
ll kth(int x,int l,int r,int k){if(l==r)return seg[x].sum;pushdown(x);if(seg[lson].sz>=k)return kth(lson,l,mid,k);else return kth(rson,mid+1,r,k-seg[lson].sz);}
ll query(int x,int l,int r,int k){if(l==r)return seg[x].sum;pushdown(x);if(seg[lson].sz>=k)return query(lson,l,mid,k);return seg[lson].sum+query(rson,mid+1,r,k-seg[lson].sz);}
void iterate(int x,int l,int r){printf("%d[%d,%d]:%d,%d\n",x,l,r,seg[x].sum,seg[x].sz);if(l!=r)pushdown(x),iterate(lson,l,mid),iterate(rson,mid+1,r);}
}s1,s2;//s1 is minimum,s2 is maximum
void turnon(int P){s1.turnon(1,1,lim,P),s1.modify(1,1,lim,P,lim,v[P-1]);s2.turnon(1,1,lim,lim-P+1),s2.modify(1,1,lim,lim-P+1,lim,v[P-1]);sz++,sum+=v[P-1];}
void turnoff(int P){s1.turnoff(1,1,lim,P),s1.modify(1,1,lim,P,lim,-v[P-1]);s2.turnoff(1,1,lim,lim-P+1),s2.modify(1,1,lim,lim-P+1,lim,-v[P-1]);sz--,sum-=v[P-1];}
ll query(){
if(sz<=1)return 0;
int l=1,r=sz>>1;
while(l<r){
int md=(l+r+1)>>1;
if(s1.kth(1,1,lim,md)>=s2.kth(1,1,lim,md-1))l=md;
else r=md-1;
}
ll ret=s1.query(1,1,lim,l);if(l>1)ret-=s2.query(1,1,lim,l-1);
ret<<=1;
if(sz&1)ret+=max(0ll,s1.kth(1,1,lim,(sz+1)>>1)-s2.kth(1,1,lim,sz>>1));
// printf("(%d,%d)\n",sum,sz);
return sum-ret;
}
int main(){
scanf("%d%d",&n,&m);
for(int i=1;i<=n;i++)scanf("%lld",&a[i]),v.push_back(a[i]);
for(int i=1,x;i<=m;i++)scanf("%d%lld",&x,&b[i]),v.push_back(b[i]),b[i]=(x==1?b[i]:-b[i]);
sort(v.begin(),v.end()),v.resize(lim=unique(v.begin(),v.end())-v.begin());
for(int i=1;i<=n;i++)a[i]=lower_bound(v.begin(),v.end(),a[i])-v.begin()+1;
for(int i=1;i<=m;i++){
if(b[i]>0)b[i]=lower_bound(v.begin(),v.end(),b[i])-v.begin()+1;
else b[i]=-(lower_bound(v.begin(),v.end(),-b[i])-v.begin()+1);
}
for(int i=1;i<=n;i++)turnon(a[i]);
printf("%lld\n",query());
for(int i=1;i<=m;i++){
if(b[i]>0)turnon(b[i]);
else turnoff(-b[i]);
// s1.iterate(1,1,lim),puts("");
// s2.iterate(1,1,lim),puts("");
printf("%lld\n",query());
}
return 0;
}
