BZOJ 1208: [HNOI2004]宠物收养所

1208: [HNOI2004]宠物收养所

Time Limit: 10 Sec Memory Limit: 162 MB
Submit: 7514 Solved: 2982
[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

5
0 2
0 4
1 3
1 2
1 5

Sample Output

3
(abs(3-2) + abs(2-4)=3，最后一个领养者没有宠物可以领养)

分析：

Splay模板题...其实可以用set水过去...然而第一次写Splay就拿来练手了...(QAQ为什么我的Splay这么长...)

代码：

  1 #include<algorithm>
  2 #include<iostream>
  3 #include<cstring>
  4 #include<cstdio>
  5 #include<cmath>
  6 //by NeighThorn
  7 #define inf 0x3f3f3f3f
  8 using namespace std;
  9 
 10 const int maxn=80000+5,MOD=1000000;
 11 
 12 int n,a,b,ans,cnt;
 13 int root,tot,w[maxn],fa[maxn],ls[maxn],rs[maxn];
 14 
 15 inline void zig(int x){
 16     int y=fa[x];
 17     if(rs[x])
 18         fa[rs[x]]=y,ls[y]=rs[x];
 19     else
 20         ls[y]=0;
 21     fa[x]=fa[y];
 22     if(fa[x]){
 23         if(ls[fa[x]]==y)
 24             ls[fa[x]]=x;
 25         else
 26             rs[fa[x]]=x;
 27     }
 28     fa[y]=x;rs[x]=y;
 29 }
 30 
 31 inline void zag(int x){
 32     int y=fa[x];
 33     if(ls[x])
 34         fa[ls[x]]=y,rs[y]=ls[x];
 35     else
 36         rs[y]=0;
 37     fa[x]=fa[y];
 38     if(fa[x]){
 39         if(ls[fa[x]]==y)
 40             ls[fa[x]]=x;
 41         else
 42             rs[fa[x]]=x;
 43     }
 44     fa[y]=x,ls[x]=y;
 45 }
 46 
 47 inline void splay(int x,int z){
 48     while(fa[x]!=z){
 49         int y=fa[x];
 50         if(fa[y]==z){
 51             if(ls[y]==x)
 52                 zig(x);
 53             else
 54                 zag(x);
 55         }
 56         else{
 57             if(ls[fa[y]]==y){
 58                 if(ls[y]==x)
 59                     zig(y),zig(x);
 60                 else
 61                     zag(x),zig(x);
 62             }
 63             else{
 64                 if(rs[y]==x)
 65                     zag(y),zag(x);
 66                 else
 67                     zig(x),zag(x);
 68             }
 69         }
 70     }
 71     if(!z)
 72         root=x;
 73 }
 74 
 75 inline void ins(int rt,int x){
 76     if(rt==0){
 77         w[++tot]=x;root=tot;
 78         return;
 79     }
 80     if(w[rt]==x)
 81         return;
 82     else if(x<w[rt]){
 83         if(!ls[rt])
 84             w[++tot]=x,ls[rt]=tot,fa[tot]=rt,splay(tot,0);
 85         else
 86             ins(ls[rt],x);
 87     }
 88     else{
 89         if(!rs[rt])
 90             w[++tot]=x,rs[rt]=tot,fa[tot]=rt,splay(tot,0);
 91         else
 92             ins(rs[rt],x);
 93     }
 94 }
 95 
 96 inline void del(int rt,int x){
 97     if(w[rt]==x){
 98         splay(rt,0);
 99         int l=ls[rt],r=rs[rt];
100         if(!l)
101             root=r,fa[r]=0;
102         else{
103             while(rs[l])
104                 l=rs[l];
105             splay(l,rt);
106             root=l,fa[l]=0,rs[l]=rs[rt],fa[rs[rt]]=l;
107         }
108     }
109     else if(x<w[rt])
110         del(ls[rt],x);
111     else
112         del(rs[rt],x);
113 }
114 
115 inline int pre(int rt,int x){
116     int res=-inf;
117     while(rt){
118         if(x==w[rt])
119             return x;
120         else{
121             if(w[rt]<x)
122                 res=max(res,w[rt]),rt=rs[rt];
123             else
124                 rt=ls[rt];
125         }
126     }
127     return res;
128 }
129 
130 inline int suf(int rt,int x){
131     int res=inf;
132     while(rt){
133         if(x==w[rt])
134             return x;
135         else{
136             if(w[rt]>x)
137                 res=min(res,w[rt]),rt=ls[rt];
138             else
139                 rt=rs[rt];
140         }
141     }
142     return res;
143 }
144 
145 signed main(void){
146     memset(ls,0,sizeof(ls));
147     memset(rs,0,sizeof(rs));
148     memset(fa,0,sizeof(fa));
149     scanf("%d",&n);tot=root=ans=cnt=0;
150     for(int i=1;i<=n;i++){
151         scanf("%d%d",&a,&b);
152         if(a==0){
153             if(cnt>=0)
154                 ins(root,b);
155             else{
156                 int lala1=pre(root,b),lala2=suf(root,b);
157                 if(abs(lala1-b)<=abs(lala2-b))
158                     (ans+=abs(lala1-b))%=MOD,del(root,lala1);
159                 else
160                     (ans+=abs(lala2-b))%=MOD,del(root,lala2);
161             }
162             cnt++;
163         }
164         else{
165             if(cnt<=0)
166                 ins(root,b);
167             else{
168                 int lala1=pre(root,b),lala2=suf(root,b);
169                 if(abs(lala1-b)<=abs(lala2-b))
170                     (ans+=abs(lala1-b))%=MOD,del(root,lala1);
171                 else
172                     (ans+=abs(lala2-b))%=MOD,del(root,lala2);
173             }
174             cnt--;
175         }
176     }
177     printf("%d
",ans);
178     return 0;
179 }

View Code

by NeighThorn