Tomgao4116

导航

 

对顶堆维护区间中位数板子:每个操作复杂度: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;
	}
}
posted on 2025-10-08 16:05  Tomgao4116  阅读(2)  评论(0)    收藏  举报