[BZOJ 1208][HNOI2004] 宠物收养所【Treap】
Problem:
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,最后一个领养者没有宠物可以领养)
Solution:
这是一道练习平衡树的好题,直接套用非旋 Treap 模板即可。
我们只需要记录当前 Treap 的类型,即存储的是宠物还是领养者。
如果新来的类型和 Treap 的类型相同,则直接插入到 Treap 中。否则在 Treap 中利用前驱/后继即可找出与当前特点值最近的两个值。另外,当 Treap 为空(即根节点为 0)时需要特判。
Code: O(nlogn) [2864K, 196MS]
#include<cstdio> #include<cstdlib> #include<cstring> #include<cmath> #include<iostream> #include<algorithm> #define MAXN 80005 #define INFINT 0x3f3f3f3f #define MOD 1000000 using namespace std; typedef pair<int, int> PII; inline void getint(int &num){ char ch; while(!isdigit(ch = getchar())); num = ch - '0'; while(isdigit(ch = getchar())) num = num * 10 + ch - '0'; } int n, tp, v, treap_tp, unpl = 0; struct Treap{ int root, val[MAXN], wt[MAXN], size[MAXN], lc[MAXN], rc[MAXN], topp; inline void update(int u) {size[u] = size[lc[u]] + size[rc[u]] + 1;} inline PII split(int u, int pos){ if(!pos) return make_pair(0, u); int l = lc[u], r = rc[u]; if(pos == size[l]) return lc[u] = 0, update(u), make_pair(l, u); else if(pos == size[l] + 1) return rc[u] = 0, update(u), make_pair(u, r); else if(pos < size[l]){ PII roots = split(l, pos); return lc[u] = roots.second, update(u), make_pair(roots.first, u); } else{ PII roots = split(r, pos - size[l] - 1); return rc[u] = roots.first, update(u), make_pair(u, roots.second); } } inline int merge(int u, int v){ if(!u || !v) return u | v; if(wt[u] < wt[v]) return rc[u] = merge(rc[u], v), update(u), u; else return lc[v] = merge(u, lc[v]), update(v), v; } inline int rnk(int x){ int u = root, delta = 0, res = INFINT; while(u){ if(x == val[u]) res = min(res, delta + size[lc[u]] + 1); if(x > val[u]) delta += size[lc[u]] + 1, u = rc[u]; else u = lc[u]; } return res < INFINT ? res : delta; } inline int kth(int k){ int u = root; while(u){ if(k == size[lc[u]] + 1) return val[u]; if(k < size[lc[u]] + 1) u = lc[u]; else k -= size[lc[u]] + 1, u = rc[u]; } } inline int pred(int x){ int u = root, res = -INFINT; while(u){ if(val[u] < x) res = max(res, val[u]), u = rc[u]; else u = lc[u]; } return res; } inline int succ(int x){ int u = root, res = INFINT; while(u){ if(val[u] > x) res = min(res, val[u]), u = lc[u]; else u = rc[u]; } return res; } inline void insert(int x){ int k = rnk(x); PII roots = split(root, k); val[++topp] = x, wt[topp] = rand(), size[topp] = 1; root = merge(merge(roots.first, topp), roots.second); } inline void erase(int x){ int k = rnk(x); PII rootr = split(root, k), rootl = split(rootr.first, k - 1); root = merge(rootl.first, rootr.second); } } treap; int main(){ getint(n); for(register int i = 1; i <= n; i++){ getint(tp), getint(v); if(!treap.root) treap_tp = tp, treap.insert(v); else if(tp == treap_tp) treap.insert(v); else{ int low = treap.pred(v + 1), up = treap.succ(v - 1); if(v - low <= up - v) unpl = (unpl + v - low) % MOD, treap.erase(low); else unpl = (unpl + up - v) % MOD, treap.erase(up); } } printf("%d\n", unpl); return 0; }
发表评论