【HNOI2004】宠物收养所
P1218 - 【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,最后一个领养者没有宠物可以领养)
开两个平衡树,一个维护主人,一个维护宠物。然后有一个主人来了,若没有宠物,
把这个人丢到splay里,否则查询前驱和后继,找出最优的那个统计答案后删除。宠物来了同上。
1 #include<set> 2 #include<map> 3 #include<queue> 4 #include<stack> 5 #include<ctime> 6 #include<cmath> 7 #include<string> 8 #include<vector> 9 #include<cstdio> 10 #include<cstdlib> 11 #include<cstring> 12 #include<iostream> 13 #include<algorithm> 14 #define maxn 80010 15 #define mod 1000000 16 #define inf 1999999999 17 using namespace std; 18 int key[2][maxn],size[2][maxn],ch[2][maxn][2],pre[2][maxn],tot[2],root[2],Ans=0; 19 inline void Rotate(int id,int x,int kind){ 20 int y=pre[id][x]; 21 ch[id][y][!kind]=ch[id][x][kind]; 22 pre[id][ch[id][x][kind]]=y; 23 if(pre[y]) 24 ch[id][pre[id][y]][ch[id][pre[id][y]][1]==y]=x; 25 pre[id][x]=pre[id][y]; 26 ch[id][x][kind]=y; 27 pre[id][y]=x; 28 size[id][y]=size[id][ch[id][y][0]]+size[id][ch[id][y][1]]+1; 29 size[id][x]=size[id][ch[id][x][0]]+size[id][ch[id][x][1]]+1; 30 } 31 inline void Splay(int id,int r,int goal){ 32 while(pre[id][r]!=goal){ 33 if(pre[id][pre[id][r]]==goal) 34 Rotate(id,r,ch[id][pre[id][r]][0]==r); 35 else{ 36 int y=pre[id][r]; 37 int kind=ch[id][pre[id][y]][0]==y; 38 if(ch[id][y][kind]==r){ 39 Rotate(id,r,!kind); 40 Rotate(id,r,kind); 41 } 42 else{ 43 Rotate(id,y,kind); 44 Rotate(id,r,kind); 45 } 46 } 47 } 48 if(goal==0) root[id]=r; 49 } 50 inline void newnode(int id,int &x,int fa,int keyy){ 51 x=++tot[id]; 52 pre[id][x]=fa; 53 key[id][x]=keyy; 54 ch[id][x][1]=ch[id][x][0]=0; 55 size[id][x]=1; 56 } 57 inline void Insert(int id,int k){ 58 while(ch[id][root[id]][key[id][root[id]]<k]) 59 size[id][root[id]]++,root[id]=ch[id][root[id]][key[id][root[id]]<k]; 60 size[id][root[id]]++; 61 newnode(id,ch[id][root[id]][key[id][root[id]]<k],root[id],k); 62 Splay(id,ch[id][root[id]][key[id][root[id]]<k],0); 63 } 64 void get_pre(int id,int x,int k,int &ans){ 65 if(!x) return; 66 if(key[id][x]<=k) ans=key[id][x],get_pre(id,ch[id][x][1],k,ans); 67 else get_pre(id,ch[id][x][0],k,ans); 68 } 69 void get_nex(int id,int x,int k,int &ans){ 70 if(!x) return; 71 if(key[id][x]>=k) ans=key[id][x],get_nex(id,ch[id][x][0],k,ans); 72 else get_nex(id,ch[id][x][1],k,ans); 73 } 74 void find_pre(int id,int x,int k,int &ans){ 75 if(!x) return; 76 if(key[id][x]<k) ans=x,find_pre(id,ch[id][x][1],k,ans); 77 else find_pre(id,ch[id][x][0],k,ans); 78 } 79 void find_nex(int id,int x,int k,int &ans){ 80 if(!x) return; 81 if(key[id][x]>k) ans=x,find_nex(id,ch[id][x][0],k,ans); 82 else find_nex(id,ch[id][x][1],k,ans); 83 } 84 void Erase(int id,int k){ 85 int pre1=root[id]; 86 find_pre(id,root[id],k,pre1); 87 int nex1=root[id]; 88 find_nex(id,root[id],k,nex1); 89 Splay(id,pre1,0);Splay(id,nex1,pre1); 90 ch[id][nex1][0]=0; 91 size[id][nex1]--; 92 size[id][pre1]--; 93 } 94 inline void solve(int id,int k){ 95 int pre1=key[id][root[id]]; 96 get_pre(id,root[id],k,pre1); 97 int nex1=key[id][root[id]]; 98 get_nex(id,root[id],k,nex1); 99 if(abs(pre1-k)<=abs(nex1-k)) Ans+=abs(pre1-k),Erase(id,pre1); 100 else Ans+=abs(nex1-k),Erase(id,nex1); 101 Ans%=mod; 102 } 103 int main() 104 { 105 freopen("!.in","r",stdin); 106 freopen("!.out","w",stdout); 107 int n,type,x; 108 scanf("%d",&n); 109 newnode(0,root[0],0,-inf);Insert(0,inf); 110 newnode(1,root[1],0,-inf);Insert(1,inf); 111 for(int i=1;i<=n;i++){ 112 scanf("%d%d",&type,&x); 113 if(size[type^1][root[type^1]]>2) 114 solve(type^1,x); 115 else Insert(type,x); 116 } 117 printf("%d",Ans%mod); 118 }