BZOJ 1208: [HNOI2004]宠物收养所
1208: [HNOI2004]宠物收养所
Time Limit: 10 Sec Memory Limit: 162 MBSubmit: 11960 Solved: 4874
[Submit][Status][Discuss]
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
0 2
0 4
1 3
1 2
1 5
Sample Output
(abs(3-2) + abs(2-4)=3,最后一个领养者没有宠物可以领养)
HINT
Source
题解:
用平衡树看似需要两颗,其实只需要一颗并打上标记是人还是宠物即可(因为题目上说同一时刻任何宠物不能同时存在);
对于每次操作,只要找到前驱元素和后继元素,然后取差值小的即可,别忘了删掉使用的元素;
参考代码:
1 /************************************************************** 2 Problem: 1208 3 User: SongHL 4 Language: C++ 5 Result: Accepted 6 Time:104 ms 7 Memory:2544 kb 8 ****************************************************************/ 9 10 #include<bits/stdc++.h> 11 using namespace std; 12 #define mod 1000000 13 const int INF=0x3f3f3f3f; 14 inline int read() 15 { 16 int x=0,f=1;char ch=getchar(); 17 while(ch<'0'||ch>'9'){if(ch=='-') f=-1;ch=getchar();} 18 while(ch>='0'&&ch<='9'){x=(x<<3)+(x<<1)+ch-'0';ch=getchar();} 19 return x*f; 20 } 21 int tr[80001][2],fa[80001],num[80001]; 22 int n,t1,t2,rt,ans,size,kind,val,flag; 23 void rotate(int x,int &k) 24 { 25 int y=fa[x],z=fa[y],l,r; 26 if(tr[y][0]==x)l=0;else l=1;r=l^1; 27 if(y==k)k=x; 28 else{if(tr[z][0]==y)tr[z][0]=x;else tr[z][1]=x;} 29 fa[x]=z;fa[y]=x;fa[tr[x][r]]=y; 30 tr[y][l]=tr[x][r];tr[x][r]=y; 31 } 32 void splay(int x,int &k) 33 { 34 int y,z; 35 while(x!=k) 36 { 37 y=fa[x],z=fa[y]; 38 if(y!=k) 39 { 40 if((tr[y][0]==x)^(tr[z][0]==y))rotate(x,k); 41 else rotate(y,k); 42 } 43 rotate(x,k); 44 } 45 } 46 void Insert(int &k,int x,int last) 47 { 48 if(k==0){size++;k=size;num[k]=x;fa[k]=last;splay(k,rt);return;} 49 if(x<num[k]) Insert(tr[k][0],x,k); 50 else Insert(tr[k][1],x,k); 51 } 52 inline void del(int x) 53 { 54 splay(x,rt); 55 if(tr[x][0]*tr[x][1]==0){rt=tr[x][0]+tr[x][1];} 56 else 57 { 58 int k=tr[x][1]; 59 while(tr[k][0]) k=tr[k][0]; 60 tr[k][0]=tr[x][0];fa[tr[x][0]]=k; 61 rt=tr[x][1]; 62 } 63 fa[rt]=0; 64 } 65 void ask_before(int k,int x) 66 { 67 if(k==0)return; 68 if(num[k]<=x){t1=k;ask_before(tr[k][1],x);} 69 else ask_before(tr[k][0],x); 70 } 71 void ask_after(int k,int x) 72 { 73 if(k==0)return; 74 if(num[k]>=x){t2=k;ask_after(tr[k][0],x);} 75 else ask_after(tr[k][1],x); 76 } 77 int main() 78 { 79 n=read();rt=size=ans=0; 80 for(int i=1;i<=n;++i) 81 { 82 kind=read();val=read(); 83 if(!rt){flag=kind;Insert(rt,val,0);} 84 else if(kind==flag) Insert(rt,val,0); 85 else 86 { 87 t1=t2=-1; 88 ask_before(rt,val);ask_after(rt,val); 89 if(t1==-1) {ans=(ans+num[t2]-val)%mod;del(t2);} 90 else if(t2==-1) {ans=(ans+val-num[t1])%mod;del(t1);} 91 else 92 { 93 if(val-num[t1]>num[t2]-val) {ans=(ans+num[t2]-val)%mod;del(t2);} 94 else {ans=(ans+val-num[t1])%mod;del(t1);} 95 } 96 } 97 } 98 printf("%d\n",ans); 99 return 0; 100 }