[BZOJ 1208]宠物收养所

解析HNOI2004宠物收养所问题,介绍如何使用Splay树处理宠物与领养者匹配,计算领养过程中的不满意程度总和。

62. [HNOI2004] 宠物收养所

★★★   输入文件:pet.in   输出文件:pet.out   简单对比
时间限制:1 s   内存限制:128 MB

最近,阿Q开了一间宠物收养所。收养所提供两种服务:收养被主人遗弃的宠物和让新的主人领养这些宠物。每个领养者都希望领养到自己满意的宠物,阿Q根据领养者的要求通过他自己发明的一个特殊的公式,得出该领养者希望领养的宠物的特点值aa是一个正整数,a<2^31),而他也给每个处在收养所的宠物一个特点值。这样他就能够很方便的处理整个领养宠物的过程了,宠物收养所总是会有两种情况发生:被遗弃的宠物过多或者是想要收养宠物的人太多,而宠物太少。

1.被遗弃的宠物过多时,假若到来一个领养者,这个领养者希望领养的宠物的特点值为a,那么它将会领养一只目前未被领养的宠物中特点值最接近a的一只宠物。(任何两只宠物的特点值都不可能是相同的,任何两个领养者的希望领养宠物的特点值也不可能是一样的)如果有两只满足要求的宠物,即存在两只宠物他们的特点值分别为a-ba+b,那么领养者将会领养特点值为a-b的那只宠物。
2.收养宠物的人过多,假若到来一只被收养的宠物,那么哪个领养者能够领养它呢?能够领养它的领养者,是那个希望被领养宠物的特点值最接近该宠物特点值的领养者,如果该宠物的特点值为a,存在两个领养者他们希望领养宠物的特点值分别为a-ba+b,那么特点值为a-b的那个领养者将成功领养该宠物。
一个领养者领养了一个特点值为a的宠物,而它本身希望领养的宠物的特点值为b,那么这个领养者的不满意程度为abs(a-b)。
【任务描述】
你得到了一年当中,领养者和被收养宠物到来收养所的情况,希望你计算所有收养了宠物的领养者的不满意程度的总和。这一年初始时,收养所里面既没有宠物,也没有领养者。
【输入格式】
你将从文件当中读入数据。文件的第一行为一个正整数nn<=80000,表示一年当中来到收养所的宠物和领养者的总数。接下来的n行,按到来时间的先后顺序描述了一年当中来到收养所的宠物和领养者的情况。每行有两个正整数a, b,其中a=0表示宠物,a=1表示领养者,b表示宠物的特点值或是领养者希望领养宠物的特点值。(同一时间呆在收养所中的,要么全是宠物,要么全是领养者,这些宠物和领养者的个数不会超过10000个)
【输入样例】

5
0 2
0 4
1 3
1 2
1 5

【输出格式】
输出文件中仅有一个正整数,表示一年当中所有收养了宠物的领养者的不满意程度的总和mod 1000000以后的结果。
【输出样例】

3

(abs(3-2) + abs(2-4)=3,最后一个领养者没有宠物可以领养)
【运行限制】
运行时限:1秒钟
【评分方法】
本题目一共有十个测试点,每个测试点的分数为总分数的10%。对于每个测试点来说,如果你给出的答案正确,那么你将得到该测试点全部的分数,否则得0分。

 分裂时要注意左子树为空的时候
还有要把父边,子边都断掉再Splay 
嗯,断边很重要,mod很重要。



#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#define maxn 100010
#define mod 1000000
using namespace std;

int n,type,val,modo,now_num,ans,root,tot;
struct Tree{
	int ch[2],val,fa;
	#define val(x) tree[x].val
	#define fa(x) tree[x].fa
	#define c(x,y) tree[x].ch[y]
}tree[maxn];
void Rotate(int p,int x){
	int mark= p==c(x,1),y=c(p,mark^1),z=fa(x);
	if(x==c(z,0))c(z,0)=p;
	if(x==c(z,1))c(z,1)=p;
	if(y)fa(y)=x;
	fa(p)=z;c(p,mark^1)=x;fa(x)=p;c(x,mark)=y;
}
void Splay(int p){
	while(fa(p)){
		int x=fa(p),y=fa(x);
		if(!y)Rotate(p,x);
		else if(x==c(y,0)^p==c(x,0))Rotate(p,x),Rotate(p,y);
		else Rotate(x,y),Rotate(p,x);
	}root=p;
}
void Insert(int t){
	int now=root,fa=root;
	while(now){
		fa=now;
		if(t<=val(now))now=c(now,0);
		else now=c(now,1);
	}
	++tot;val(tot)=t;fa(tot)=fa;
	if(fa)c(fa,t>val(fa))=tot;
	Splay(tot);
}
int Pre(int t){
	int now=root,ans=-0x7ffffff,k=root;
	while(now){
		if(t>=val(now)){
			if(ans<=val(now)){
				ans=val(now);
				k=now;
			}now=c(now,1);
		}else now=c(now,0);
	}return k;
}
int Suc(int t){
	int now=root,ans=0x7ffffff,k=root;
	while(now){
		if(t<=val(now)){
			if(ans>=val(now)){
				ans=val(now);
				k=now;
			}now=c(now,0);
		}else now=c(now,1);
	}return k;
}
void Del(int t){
	int p=Pre(t),s=Suc(t);
	int kp=abs(val(p)-t),ks=abs(val(s)-t);
	if(kp<=ks)Splay(p),ans+=kp;
	else Splay(s),ans+=ks;
	ans%=mod;
	int r=c(root,1),l=c(root,0);
	fa(l)=fa(r)=0;c(root,0)=c(root,1)=0;
	if(!l){root=r;return;}
	while(c(l,1))l=c(l,1);
	Splay(l);c(root,1)=r;fa(r)=root;
}

int main(){
	scanf("%d",&n);
	for(int i=1;i<=n;i++){
		scanf("%d%d",&type,&val);
		if(!type)type=-1;
		if(!now_num)Insert(val);
		else if((now_num>0)^(type>0))Del(val);
		else Insert(val);
		now_num+=type;
		
	}printf("%d",ans%mod);
	return 0;	
}


评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值