bzoj1208【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
(abs(3-2) + abs(2-4)=3，最后一个领养者没有宠物可以领养)

HINT

Source

所以我们就可以用一颗平衡树维护没有主人的宠物或者没有宠物的主人，并且用一个标记记录现在平衡树维护的是主人还是宠物。

那么对于每次操作，如果平衡树为空或者该操作与标记相同，就把这个数插入到平衡树中，反之则删去前驱和后继中更接近它的一个数。

注意这道题输出结果要取模，为什么每次都看不见取模呢...=_=

#include<iostream> #include<cstdio> #include<cmath> #include<cstring> #include<algorithm> #include<cstdlib> #define F(i,j,n) for(int i=j;i<=n;i++) #define D(i,j,n) for(int i=j;i>=n;i--) #define LL long long #define pa pair<int,int> #define MAXN 80005 #define INF 1000000000000 #define mod 1000000 using namespace std; int rt=0,tot=0,l[MAXN],r[MAXN]; LL n,x,y,flag=0,sum=0,sz=0,v[MAXN],rnd[MAXN]; inline LL read() { LL x=0,f=1;char ch=getchar(); while (ch<'0'||ch>'9'){if (ch=='-') f=-1;ch=getchar();} while (ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();} return x*f; } inline void rturn(int &k){int tmp=l[k];l[k]=r[tmp];r[tmp]=k;k=tmp;} inline void lturn(int &k){int tmp=r[k];r[k]=l[tmp];l[tmp]=k;k=tmp;} inline void ins(int &k,LL x) { if (!k){k=++tot;l[k]=r[k]=0;v[k]=x;rnd[k]=rand();return;} if (x<v[k]){ins(l[k],x);if (rnd[l[k]]<rnd[k]) rturn(k);} else {ins(r[k],x);if (rnd[r[k]]<rnd[k]) lturn(k);} } inline void del(int &k,LL x) { if (v[k]==x) { if (!l[k]||!r[k]) k=l[k]+r[k]; else { if (rnd[l[k]]<rnd[r[k]]){rturn(k);del(k,x);} else {lturn(k);del(k,x);} } return; } else if (v[k]>x) del(l[k],x); else del(r[k],x); } inline LL pre(int k,LL x) { if (!k) return -INF; if (x==v[k]) return v[k]; if (x<v[k]) return pre(l[k],x); else { LL tmp=pre(r[k],x); return tmp==-INF?v[k]:tmp; } } inline LL suc(int k,LL x) { if (!k) return INF; if (x==v[k]) return v[k]; if (x>v[k]) return suc(r[k],x); else { LL tmp=suc(l[k],x); return tmp==INF?v[k]:tmp; } } inline void solve(LL x) { LL ans=pre(rt,x),tmp=suc(rt,x); if (abs(x-tmp)<abs(x-ans)) ans=tmp; sum=(sum+abs(x-ans))%mod; del(rt,ans); } int main() { n=read(); while(n--) { x=read();y=read(); if (!sz){ins(rt,y);flag=x;sz++;} else if (flag==x){ins(rt,y);sz++;} else {solve(y);sz--;} } printf("%lld\n",sum); }