HNOI2004 宠物收养所

探讨一种特殊算法,用于解决宠物收养所中宠物与领养者匹配问题,通过计算领养者的不满意程度,实现高效宠物分配。

Description

最近,阿Q开了一间宠物收养所。收养所提供两种服务:收养被主人遗弃的宠物和让新的主人领养这些宠物。 
  每个领养者都希望领养到自己满意的宠物,阿Q根据领养者的要求通过他自己发明的一个特殊的公式,得出该领养者希望领养的宠物的特点值a(a是一个正整数,a<2^31),而他也给每个处在收养所的宠物一个特点值。这样他就能够很方便的处理整个领养宠物的过程了,宠物收养所总是会有两种情况发生:被遗弃的宠物过多或者是想要收养宠物的人太多,而宠物太少。 
1、被遗弃的宠物过多时,假若到来一个领养者,这个领养者希望领养的宠物的特点值为a,那么它将会领养一只目前未被领养的宠物中特点值最接近a的一只宠物。(任何两只宠物的特点值都不可能是相同的,任何两个领养者的希望领养宠物的特点值也不可能是一样的)如果有两只满足要求的宠物,即存在两只宠物他们的特点值分别为a-b和a+b,那么领养者将会领养特点值为a-b的那只宠物。 
2、收养宠物的人过多,假若到来一只被收养的宠物,那么哪个领养者能够领养它呢?能够领养它的领养者,是那个希望被领养宠物的特点值最接近该宠物特点值的领养者,如果该宠物的特点值为a,存在两个领养者他们希望领养宠物的特点值分别为a-b和a+b,那么特点值为a-b的那个领养者将成功领养该宠物。 
  一个领养者领养了一个特点值为a的宠物,而它本身希望领养的宠物的特点值为b,那么这个领养者的不满意程度为abs(a-b)。 
  你的任务:你得到了一年当中,领养者和被收养宠物到来收养所的情况,希望你计算所有收养了宠物的领养者的不满意程度的总和。这一年初始时,收养所里面既没有宠物,也没有领养者。 

Input

  第一行为一个正整数n,n <= 80000,表示一年当中来到收养所的宠物和领养者的总数。 
  接下来的n行,按到来时间的先后顺序描述了一年当中来到收养所的宠物和领养者的情况。每行有两个正整数a, b,其中a=0表示宠物,a=1表示领养者,b表示宠物的特点值或是领养者希望领养宠物的特点值。(同一时间呆在收养所中的,要么全是宠物,要么全是领养者,这些宠物和领养者的个数不会超过10000个)

Output

  仅有一个正整数,表示一年当中所有收养了宠物的领养者的不满意程度的总和mod 1000000以后的结果。

Sample Input

5

0 2

0 4

1 3

1 2

1 5

Sample Output

3

Hint

样例说明:样例的计算方法:(abs(3-2) + abs(2-4)=3,最后一个领养者没有宠物可以领养)

没什么好讲的= =比普通平衡树还简单,不过这道题告诉了我们人是宠物,宠物是人o((⊙﹏⊙))o

#include<bits/stdc++.h>
using namespace std;
#define Inc(i,L,r) for(register int i=(L);i<=(r);++i)
#define Red(i,r,L) for(register int i=(r);i>=(L);--i)
const int N = 4e4+10;
#define Mod 1000000
long long sum;
struct Splay{
	int rt,cur,p[N],ch[N][2];
	int k[N],siz[N],cnt[N];
	#define Ls(v) ch[v][0]
	#define rs(v) ch[v][1]
	#define sum(v) siz[v]=siz[Ls(v)]+siz[rs(v)]+cnt[v]
	inline void rot(int x){
		int f=p[x],gf=p[f],type=ch[f][1]==x,son=ch[x][!type];
		ch[p[son]=f][type]=son,sum(f);
		ch[p[f]=x][!type]=f,sum(x);
		ch[p[x]=gf][ch[gf][1]==f]=x;
	}
	inline void splay(int x){
		while(p[x]){
			if(p[p[x]]&&((ch[p[p[x]]][1]==p[x])==(ch[p[x]][1]==x)))rot(p[x]);
			rot(x);
		}
		rt=x;
	}
	inline void insert(int &x,int fa,int v){
		if(k[x]==v)return ++siz[x],++cnt[x],splay(x),void();
		if(!x)return siz[x=++cur]=cnt[x]=1,p[x]=fa,k[x]=v,splay(x),void();
		insert(ch[x][k[x]<v],x,v);
	}
	inline int FindPre(int x){
		if(!Ls(x))return 0;
		x=Ls(x);
		while(rs(x))x=rs(x);
		return x;
	}
	inline int FindSub(int x){
		if(!rs(x))return 0;
		x=rs(x);
		while(Ls(x))x=Ls(x);
		return x;
	}
	inline void Delete(int x,int v){
		if(k[x]==v){
			if(cnt[x]>1)return --siz[x],--cnt[x],splay(x),void();
			splay(x);
			if(Ls(x)){
				int tmp=rs(x),New=FindPre(x);
				splay(New);
				rs(p[tmp]=New)=tmp;
			}else p[rs(x)]=0,rt=rs(x);
			return ;
		}
		Delete(ch[x][k[x]<v],v);
	}
	inline void solv(int ex){
		insert(rt,0,ex);
		int p1=FindPre(rt),p2=FindSub(rt),Ex1=0x3f3f3f3f,Ex2=0x3f3f3f3f;
		if(cnt[rt]>1){//本来就有期望值 
			Delete(rt,ex);
			Delete(rt,ex);
			return ;
		}
		if(p1)Ex1=k[p1];
		if(p2)Ex2=k[p2];
		if(abs(ex-Ex1)<=abs(ex-Ex2)){//第二种删法tle= =第一种跑得飞快= = 
			Delete(rt,ex);
			Delete(rt,Ex1);
//			Delete(rt,Ex1);
//			Delete(rt,ex);
			(sum+=abs(ex-Ex1))%=Mod;
		}else {
			Delete(rt,ex);
			Delete(rt,Ex2);
//			Delete(rt,Ex2);
//			Delete(rt,ex);
			(sum+=abs(ex-Ex2))%=Mod;
		}
	}
}sp[2];
int num[2];
int main(){
	int n;scanf("%d",&n);
	while(n--){
		int type,expect;scanf("%d%d",&type,&expect);
		if(num[!type])sp[!type].solv(expect),--num[!type];
		else sp[type].insert(sp[type].rt,0,expect),++num[type];
	}
	cout<<sum<<"\n";
	return 0;
}

 

评论
成就一亿技术人!
拼手气红包6.0元
还能输入1000个字符
 
红包 添加红包
表情包 插入表情
 条评论被折叠 查看
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值