BZOJ1208 HNOI2004 宠物收养所 【非旋转Treap】

114 篇文章 0 订阅
33 篇文章 0 订阅

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

HINT

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


脑洞清奇的非旋Treap做法,看到题目中描述了两个时间,送宠物和领宠物,我们可以很容易的发现我们需要维护一个可以支持插入(送宠物),删除(领宠物),查询前驱+后继(匹配领哪一只)的数据结构,果断敲非旋Treap,思路是维护两棵Treap,一棵记录当前的宠物有哪些,另一棵记录当前没有满足的领宠物需求有哪些,容易发现,两棵Treap在同一时刻只有一个不为空,所以我们只需要在询问的时候判断一下需要访问的那一棵树是不是空,如果为空就在相应的树中插入,否则在另一棵树中查询答案后删除

写的时候因为把节点标号误认为是权值GG掉了。。。

大爱非旋Treap


#include<bits/stdc++.h>
using namespace std;
#define LL long long
LL read(){
    LL ans=0,w=1;char c=getchar();
    while(!isdigit(c)&&c!='-')c=getchar();
    if(c=='-')w=-1,c=getchar();
    while(isdigit(c))ans=ans*10+c-'0',c=getchar();
    return ans*w;
}
#define N 80010
typedef pair<LL,LL> pi;
const LL Mod=1e6;
struct Treap{
    LL siz[N],val[N],key[N],son[N][2],tot,root;
    void clear(){tot=root=0;}
    LL new_treap_point(LL vl){
        LL newt=++tot;
        siz[newt]=1;
        son[newt][0]=son[newt][1]=0;
        val[newt]=vl;
        key[newt]=rand();
        return newt;
    }
    void update(LL o){
        siz[o]=siz[son[o][0]]+siz[son[o][1]]+1;
    }
    LL merge(LL a,LL b){
        if(!a)return b;
        if(!b)return a;
        if(key[a]<key[b]){
            son[a][1]=merge(son[a][1],b);
            update(a);return a;
        }else{
            son[b][0]=merge(a,son[b][0]);
            update(b);return b;
        }
    }
    pi split(LL o,LL k){
        if(!o)return pi(0,0);
        pi ans;
        if(siz[son[o][0]]>=k){
            pi tmp=split(son[o][0],k);
            son[o][0]=tmp.second;
            update(o);
            ans.first=tmp.first;
            ans.second=o;
        }else{
            pi tmp=split(son[o][1],k-siz[son[o][0]]-1);
            son[o][1]=tmp.first;
            update(o);
            ans.first=o;
            ans.second=tmp.second;
        }
        return ans;
    }
    LL getrank(LL o,LL vl){
        if(!o)return 0;
        if(val[o]==vl)return siz[son[o][0]]+1;
        else if(val[o]>vl)return getrank(son[o][0],vl);
        else return getrank(son[o][1],vl)+siz[son[o][0]]+1;
    }
    void insert(LL vl){
        LL k=getrank(root,vl);
        pi x=split(root,k);
        LL o=new_treap_point(vl);
        root=merge(merge(x.first,o),x.second);
    }
    void erase(LL vl){
        LL k=getrank(root,vl);
        pi x=split(root,k);
        root=merge(split(x.first,k-1).first,x.second);
    }
    LL get_pre(int vl){
        int k=getrank(root,vl);
        pi x=split(root,k);
        pi y=split(x.first,k-1);
        LL ans=val[y.second];
        root=merge(merge(y.first,y.second),x.second);
        return ans;
    }
    LL get_next(int vl){
        int k=getrank(root,vl);
        pi x=split(root,k);
        pi y=split(x.second,1);
        LL ans=val[y.first];
        root=merge(x.first,merge(y.first,y.second));
        return ans;
    }
    LL query(LL vl){
        LL ls=get_pre(vl),rs=get_next(vl);
        if(!ls)return rs;
        if(!rs)return ls;
        if(rs-vl<vl-ls)return rs;
        else return ls;
    }
}A,B;
LL n,ans=0;
int main(){
    A.clear();B.clear();
    n=read();
    for(LL i=1;i<=n;i++){
        LL op=read(),x=read();
        switch(op){
            case 0:{
                if(!B.siz[B.root])A.insert(x);
                else{
                    LL t=B.query(x);B.erase(t);
                    ans=(ans+abs(t-x))%Mod;
                }
                break;
            }
            case 1:{
                if(!A.siz[A.root])B.insert(x);
                else{
                    LL t=A.query(x);A.erase(t);
                    ans=(ans+abs(t-x))%Mod;
                }
                break;
            }
        }
    }
    printf("%lld",ans);
    return 0;
}
  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值