对顶堆维护区间中位数板子:每个操作复杂度:O(logn)
点击查看代码
//加入a[p]的值:sum1,sum2为统计>mid与<=mid数之和分别是多少,l,r为取的数的区间范围
void add(ll p){
if(a[p]>mid){
sum2+=a[p];
q2.insert({a[p],p});
}else{
sum1+=a[p];
q1.insert({a[p],p});
}
if((r-l+1)%2){
if(q1.size()<=q2.size()){
pll v=*q2.begin();
sum2-=v.first;
sum1+=v.first;
q2.erase(q2.begin());
q1.insert(v);
}
if(q1.size()>q2.size()+1){
pll v=*q1.begin();
sum1-=v.first;
sum2+=v.first;
q1.erase(q1.begin());
q2.insert(v);
}
}else{
if(q1.size()<q2.size()){
pll v=*q2.begin();
sum2-=v.first;
sum1+=v.first;
q2.erase(q2.begin());
q1.insert(v);
}
if(q1.size()>=q2.size()+1){
pll v=*q1.begin();
sum1-=v.first;
sum2+=v.first;
q1.erase(q1.begin());
q2.insert(v);
}
}
pll x=*q1.begin();
mid=x.first;
val=(r-l+1)/2;
if((r-l+1)%2)val++;
cnt=val*mid-sum1+sum2-(r-l+1-val)*mid;
}
//删除a[p]的值:
void del(ll p){
if(q1.count({a[p],p})>0){
sum1-=a[p];
q1.erase({a[p],p});
}else{
sum2-=a[p];
q2.erase({a[p],p});
}
if((r-l)%2){
if(q1.size()<=q2.size()){
pll v=*q2.begin();
sum2-=v.first;
sum1+=v.first;
q2.erase(q2.begin());
q1.insert(v);
}
if(q1.size()>q2.size()+1){
pll v=*q1.begin();
sum1-=v.first;
sum2+=v.first;
q1.erase(q1.begin());
q2.insert(v);
}
}else{
if(q1.size()<q2.size()){
pll v=*q2.begin();
sum2-=v.first;
sum1+=v.first;
q2.erase(q2.begin());
q1.insert(v);
}
if(q1.size()>=q2.size()+1){
pll v=*q1.begin();
sum1-=v.first;
sum2+=v.first;
q1.erase(q1.begin());
q2.insert(v);
}
}
if(q1.size()){
pll x=*q1.begin();
mid=x.first,p=x.second;
val=(r-l)/2;
if((r-l)%2)val++;
cnt=val*mid-sum1+sum2-(r-l-val)*mid;
}else{
mid=0,cnt=0;
}
}