BZOJ 1208: [HNOI2004]宠物收养所

BZOJ 1208: [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

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

---

题解

平衡树

因为，同一时间只有宠物或人，所以只要写一个平衡树就可以了

---

代码

#include<cstdio>
#define tt 1000000
using namespace std;

int n,sz,rt,ans,a,b,ret,con,t1,t2;
int son[80005][2],f[80005],key[80005];

int readln()
{
    int x=0;
    char ch=getchar();
    while (ch<'0'||ch>'9') ch=getchar();
    while ('0'<=ch&&ch<='9') x=x*10+ch-48,ch=getchar();
    return x;
}

bool get(int x) {return son[f[x]][1]==x;}

void rotate(int x)
{
    int fa=f[x],gfa=f[fa],p=get(x);
    son[fa][p]=son[x][p^1];f[son[fa][p]]=fa;
    son[x][p^1]=fa;f[fa]=x;
    f[x]=gfa;
    if (gfa) son[gfa][son[gfa][1]==fa]=x;
}

void splay(int x)
{
    for (int fa;fa=f[x];rotate(x)) if (f[fa]) rotate((get(x)==get(fa)?fa:x));
    rt=x;
}

void insert(int x)
{
    if (rt==0) {
        sz++;
        son[sz][0]=son[sz][1]=f[sz]=0;
        rt=sz;
        key[sz]=x;
        return;
    }
    int now=rt,fa=0;
    while (1)
    {
        if (x==key[now]){
            splay(now);break;
        }
        fa=now;
        now=son[now][key[now]<x];
        if (now==0)
        {
            sz++;
            son[sz][0]=son[sz][1]=0;
            f[sz]=fa;
            son[fa][key[fa]<x]=sz;
            key[sz]=x;
            splay(sz);
            break;
        }
    }
}

void pre(int now,int x)
{
    if (now==0) return;
    if (key[now]<=x) {t1=now;pre(son[now][1],x);}
    else pre(son[now][0],x);
}

void nxt(int now,int x)
{
    if (now==0) return;
    if (key[now]>=x) {t2=now;nxt(son[now][0],x);}
    else nxt(son[now][1],x);
}

void del(int x)
{
    splay(x);
    if (son[x][0]*son[x][1]==0) rt=son[x][0]+son[x][1];
    else {
        int now=son[x][1];
        while (son[now][0]) now=son[now][0];
        son[now][0]=son[x][0];f[son[x][0]]=now;
        rt=son[x][1];
    }
    f[rt]=0;
}

int main()
{
    ans=0;
    n=readln();
    while (n--)
    {
        a=readln();b=readln();
        if (!rt) con=a;
        if (con==a) insert(b);
        else {
            t1=t2=-1;
            pre(rt,b);nxt(rt,b);
            if (t1==-1) {ans+=(key[t2]-b)%tt;ans%=tt;del(t2);}
            else if (t2==-1) {ans+=(b-key[t1])%tt;ans%=tt;del(t1);}
            else {
                if (b-key[t1]>key[t2]-b) {ans+=(key[t2]-b)%tt;ans%=tt;del(t2);}
                else {ans+=(b-key[t1])%tt;ans%=tt;del(t1);}
            }
        }
    }
    printf("%d",ans);
    return 0;
}