BZOJ1208[HNOI2004]宠物收养场——treap

凡凡开了一间宠物收养场。收养场提供两种服务:收养被主人遗弃的宠物和让新的主人领养这些宠物。

每个领养者都希望领养到自己满意的宠物,凡凡根据领养者的要求通过他自己发明的一个特殊的公式,得出该领养者希望领养的宠物的特点值a(a是一个正整数,a<2^31),而他也给每个处在收养场的宠物一个特点值。这样他就能够很方便的处理整个领养宠物的过程了,宠物收养场总是会有两种情况发生:被遗弃的宠物过多或者是想要收养宠物的人太多,而宠物太少。

被遗弃的宠物过多时,假若到来一个领养者,这个领养者希望领养的宠物的特点值为a,那么它将会领养一只目前未被领养的宠物中特点值最接近a的一只宠物。(任何两只宠物的特点值都不可能是相同的,任何两个领养者的希望领养宠物的特点值也不可能是一样的)如果有两只满足要求的宠物,即存在两只宠物他们的特点值分别为a-b和a+b,那么领养者将会领养特点值为a-b的那只宠物。

收养宠物的人过多,假若到来一只被收养的宠物,那么哪个领养者能够领养它呢?能够领养它的领养者,是那个希望被领养宠物的特点值最接近该宠物特点值的领养者,如果该宠物的特点值为a,存在两个领养者他们希望领养宠物的特点值分别为a-b和a+b,那么特点值为a-b的那个领养者将成功领养该宠物。

一个领养者领养了一个特点值为a的宠物,而它本身希望领养的宠物的特点值为b,那么这个领养者的不满意程度为abs(a-b)。

你得到了一年当中,领养者和被收养宠物到来收养所的情况,请你计算所有收养了宠物的领养者的不满意程度的总和。这一年初始时,收养所里面既没有宠物,也没有领养者。

输入格式: 

第一行为一个正整数n,n<=80000,表示一年当中来到收养场的宠物和领养者的总数。接下来的n行,按到来时间的先后顺序描述了一年当中来到收养场的宠物和领养者的情况。每行有两个正整数a, b,其中a=0表示宠物,a=1表示领养者,b表示宠物的特点值或是领养者希望领养宠物的特点值。(同一时间呆在收养所中的,要么全是宠物,要么全是领养者,这些宠物和领养者的个数不会超过10000个) 

输出格式:

仅有一个正整数,表示一年当中所有收养了宠物的领养者的不满意程度的总和mod 1000000以后的结果。

这是一道比较简单的平衡树练习题,因为在同一时间只有领养者或者宠物,所以可以只开一棵treap就可以(我为了求稳开了两棵)。

先讲一下两棵treap的做法。

一棵存领养者,另一棵存宠物。

首先,一共有两种操作,一种是来了领养者,另一种是来了宠物。

那么对于领养者有两种情况:

1.没有宠物;直接将领养者加到领养者的树上。

2.有宠物;在宠物的树上找到想要的宠物。

找宠物时要先找是否有符合的宠物,如果没有再找特点值最接近的(也就是找当前值的前驱和后继)。

对于宠物同样是这两种:

1.没有领养者;直接将宠物加到宠物的树上。

2.有领养者;在领养者的树上找到合适的领养者。

一棵树的做法和两棵树类似,只要每一次记录当前树上是宠物还是领养者就行了。

另外需要注意的是每次查找前驱和后继之后要把变量清零,还有要把根节点特点值设成最大值,因为一次查找时不一定前驱和后继都有。

#include<cstdio>
#include<algorithm>
#include<cmath>
#include<cstring>
#include<iostream>
#include<ctime>
using namespace std;
int n;
int opt;
int x;
int size[100010][2];
int rs[100010][2];
int ls[100010][2];
long long v[100010][2];
int w[100010][2];
int r[100010][2];
int tot[2];
int root[2];
long long ans;
int answer_back;
int answer_front;
int answer_front1;
int answer_back1;
void up(int x,int s)
{
    size[x][s]=size[rs[x][s]][s]+size[ls[x][s]][s]+w[x][s];
}
void rturn(int &x,int s)
{
    int t;
    t=ls[x][s];
    ls[x][s]=rs[t][s];
    rs[t][s]=x;
    size[t][s]=size[x][s];
    up(x,s);
    x=t;
}
void lturn(int &x,int s)
{
    int t;
    t=rs[x][s];
    rs[x][s]=ls[t][s];
    ls[t][s]=x;
    size[t][s]=size[x][s];
    up(x,s);
    x=t;
}
void insert_sum(int x,int &i,int s)
{
    if(!i)
    {
        i=++tot[s];
        w[i][s]=size[i][s]=1;
        v[i][s]=x;
        r[i][s]=rand();
        return ;
    }
    size[i][s]++;
    if(x==v[i][s])
    {
        w[i][s]++;
    }
    else if(x>v[i][s])
    {
        insert_sum(x,rs[i][s],s);
        if(r[rs[i][s]][s]<r[i][s])
        {
            lturn(i,s);
        }
    }
    else
    {
        insert_sum(x,ls[i][s],s);
        if(r[ls[i][s]][s]<r[i][s])    
        {
            rturn(i,s);
        }
    }
      
    return ;
}
void delete_sum(int x,int &i,int s)
{
    if(i==0)
    {
        return ;
    }
    if(v[i][s]==x)
    {
        if(w[i][s]>1)
        {
            w[i][s]--;
            size[i][s]--;
            return ;
        }
        if((ls[i][s]*rs[i][s])==0)
        {
            i=ls[i][s]+rs[i][s];
        }
        else if(r[ls[i][s]][s]<r[rs[i][s]][s])
        {
            rturn(i,s);
            delete_sum(x,i,s);
        }
        else
        {
            lturn(i,s);
            delete_sum(x,i,s);
        }
        return ;
    }
    size[i][s]--;
    if(v[i][s]<x)
    {
        delete_sum(x,rs[i][s],s);
    }
    else
    {
        delete_sum(x,ls[i][s],s);
    }
    return ;
}
int ask_num(int x,int i,int s)
{
    if(i==0)
    {
        return 0;
    }
    if(v[i][s]==x)
    {
        return size[ls[i][s]][s]+1;
    }
    if(v[i][s]<x)
    {
        return ask_num(x,rs[i][s],s)+size[ls[i][s]][s]+w[i][s];
    }
    return ask_num(x,ls[i][s],s);
}
void ask_front(int x,int i,int s)
{
    if(i==0)
    {
        return ;
    }
    if(v[i][s]<x)
    {
        answer_front1=i;
        ask_front(x,rs[i][s],s);
        return ;
    }
    else
    {
        ask_front(x,ls[i][s],s);
        return ;
    }
    return ;
}
void ask_back(int x,int i,int s)
{
    if(i==0)
    {
        return ;
    }
    if(v[i][s]>x)
    {
        answer_back1=i;
        ask_back(x,ls[i][s],s);
        return ;
    }
    else
    {
        ask_back(x,rs[i][s],s);
        return ;
    }
}
int main()
{
    srand(12378);
    scanf("%d",&n);
    v[0][0]=v[0][1]=1000000000;
    for(int i=1;i<=n;i++)
    {
        answer_front1=answer_back1=0;
        scanf("%d%d",&opt,&x);
        if(size[root[opt^1]][opt^1]==0)
        {
            insert_sum(x,root[opt],opt);
        }
        else
        {
            if(v[ask_num(x,root[opt^1],opt^1)][opt^1]==x)
            {
                delete_sum(x,root[opt^1],opt^1);
            }
            else
            {
                ask_front(x,root[opt^1],opt^1);
                ask_back(x,root[opt^1],opt^1);
                answer_front=v[answer_front1][opt^1];
                answer_back=v[answer_back1][opt^1];
                if(abs(answer_back-x)==abs(answer_front-x))
                {
                    ans+=abs(answer_front-x);
                    delete_sum(answer_front,root[opt^1],opt^1);
                }             
                else if(abs(answer_back-x)>abs(answer_front-x))
                {
                    ans+=abs(answer_front-x);
                    delete_sum(answer_front,root[opt^1],opt^1);
                }
                else
                {
                    ans+=abs(answer_back-x);
                    delete_sum(answer_back,root[opt^1],opt^1);
                } 
            }
        }
        ans%=1000000;
    }
    printf("%lld",ans%1000000);
    return 0;
}
原文地址:https://www.cnblogs.com/Khada-Jhin/p/8973804.html