洛谷 P2286 [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以后的结果。

输入输出样例

输入样例#1:
5
0 2
0 4
1 3
1 2
1 5
输出样例#1:
3
注:abs(3-2) + abs(2-4)=3,
最后一个领养者没有宠物可以领养。

分析:
题目要求一个满足插入,删除,查询前继与后继。显然就是平衡树(Treap)。 注意那些没领着宠物的人不会走,直到有宠物可以领或者是直到输入完数据也没有宠物可以领,然后失望的回家。所以要开两棵treap啦,显然在过程上像fft那样打上(*T)就可以了呀。

代码:

#include <iostream>
#include <cstdio>
#include <cmath>
#include <ctime>
#include <cstdlib>

const int mod=1000000;
const int maxn=80007;
const int inf=2147483647;

using namespace std;

int ans;
int n,op,x,root1,root2,cnt,num;

struct node{
    int data,key;
    int l,r;
    int size;
}a[maxn],b[maxn];

void updata(int x,node *t)
{
    t[x].size=t[t[x].l].size+t[t[x].r].size+1;
}

void rttr(int &x,node *t)
{
    int y=t[x].l;
    t[x].l=t[y].r;
    t[y].r=x;
    updata(x,t); updata(y,t);
    x=y;
}

void rttl(int &x,node *t)
{
    int y=t[x].r;
    t[x].r=t[y].l;
    t[y].l=x;
    updata(x,t); updata(y,t);
    x=y;
}

void ins(int &x,int k,node *t)
{
    if (x==0)
    {
        x=++cnt;
        t[cnt].size=1;
        t[cnt].data=k;
        t[cnt].key=rand();
        return;
    }
    if (k<=t[x].data)
    {
        ins(t[x].l,k,t);
        if (t[t[x].l].key<t[x].key) rttr(x,t);
    }
    else
    {
        ins(t[x].r,k,t);
        if (t[t[x].r].key<t[x].key) rttl(x,t);
    }
    updata(x,t);
}

int get_pre(int x,int k,node *t)
{
    int ans=0;
    while (x)
    {
        if (k<t[x].data) x=t[x].l;
        else
        {
            ans=max(ans,t[x].data);
            x=t[x].r;
        }
    }
    return ans;
}

int get_next(int x,int k,node *t)
{
    int ans=inf;
    while (x)
    {
        if (k>t[x].data) x=t[x].r;
        else
        {
            ans=min(ans,t[x].data);
            x=t[x].l;
        }
    }
    return ans;
}

void del(int &x,int k,node *t)
{
    if (t[x].data==k)
    {
        if ((t[x].l) || (t[x].r))
        {
            if ((t[x].l) && ((!t[x].r) || (t[t[x].l].key<t[t[x].r].key)))
            {
                rttr(x,t);
                del(t[x].r,k,t);
            }
            else
            {
                rttl(x,t);
                del(t[x].l,k,t);
            }
            updata(x,t);
        }
        else x=0;
        return;
    }
    if (k<t[x].data) del(t[x].l,k,t);
    else del(t[x].r,k,t);
    updata(x,t);
    return;
}

int main()
{
    scanf("%d",&n);     
    for (int i=1;i<=n;i++)
    {       
        scanf("%d%d",&op,&x);
        if (op==0)
        {
            if (num>=0) ins(root1,x,a);
            else
            {
                int l=get_pre(root2,x,b),r=get_next(root2,x,b);
                if ((l==0) || (r-x<x-l))
                {
                    ans=(ans+r-x)%mod;
                    del(root2,r,b);
                }
                else 
                {
                    ans=(ans+x-l)%mod;
                    del(root2,l,b);
                }
            }
            num++;
        }
        else
        {
            if (num<=0) ins(root2,x,b);
            else
            {
                int l=get_pre(root1,x,a),r=get_next(root1,x,a);
                if ((l==0) || (r-x<x-l))
                {
                    ans=(ans+r-x)%mod;
                    del(root1,r,a);
                }
                else 
                {
                    ans=(ans+x-l)%mod;
                    del(root1,l,a);
                }
            }
            num--;
        }
    }
    printf("%d",ans);
}

猜你喜欢

转载自blog.csdn.net/liangzihao1/article/details/80412128
今日推荐