博客
关于我
强烈建议你试试无所不能的chatGPT,快点击我
BZOJ 1208 [HNOI 2004] 宠物收养所 Splay
阅读量:4506 次
发布时间:2019-06-08

本文共 5665 字,大约阅读时间需要 18 分钟。

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,最后一个领养者没有宠物可以领养)

如标题所说是一颗普通的平衡树,只需插入一个数,查询前驱、后继,然后删除。

#include 
#include
#include
using namespace std;typedef long long LL;const int SZ = 1000010;const int INF = 1000000010;const int MOD = 1000000;struct node{ node *ch[2], *f; int sz, cnt, v; void maintain() { sz = cnt + ch[0] -> sz + ch[1] -> sz; } int cmp(int x) { if(x == v) return -1; return x < v ? 0 : 1; } int dir() { return f -> ch[1] == this; } void setc(node *x, int d) { (ch[d] = x) -> f= this; }}T[SZ], *root, *null;int Tcnt;node *newnode(int x, node *f){ node *k = T + (++ Tcnt); k -> v = x; k -> ch[0] = k -> ch[1] = null; k -> sz = k -> cnt = 1; k -> f = f; return k;}void rotate(node *p){ node *fa = p -> f; int d = p -> dir(); fa -> f -> setc(p, fa -> dir()); fa -> setc(p -> ch[d ^ 1], d); fa -> maintain(); p -> setc(fa, d ^ 1); p -> maintain(); if(fa == root) root = p;}void splay(node *p, node *rt = null){ while(p -> f != rt) { if(p ->f -> f == rt) rotate(p); else { if(p -> dir() == p -> f -> dir()) rotate(p -> f), rotate(p); else rotate(p), rotate(p); } } p -> maintain();}void insert(node *p, int x){ if(root == null) { root = newnode(x, null); return ; } while(p != null) { p -> sz ++; int d = p -> cmp(x); if(d == -1) { p -> cnt ++; return ; } if(p -> ch[d] == null) { p -> ch[d] = newnode(x, p); p = p -> ch[d]; break; } p = p -> ch[d]; } splay(p);}void erase(node *p, int x){ while(p != null) { p -> sz --; int d = p -> cmp(x); if(d == -1) { p -> cnt --; break; } p = p -> ch[d]; } if(p -> cnt) return ; splay(p); if(p -> ch[0] == null) { root = p -> ch[1]; root -> f = null; return ; } if(p -> ch[1] == null) { root = p -> ch[0]; root -> f = null; return ; } p = p -> ch[0]; while(p -> ch[1] != null) p = p -> ch[1]; splay(p, root); p -> ch[1] = root -> ch[1]; p -> ch[1] -> f = p; p -> f = null; p -> maintain(); root = p;}void init(){ null = newnode(-INF, null); null -> sz = null -> cnt = 0; root = null;}int ask_pre(node *p, int x){ int ans = INF; while(p != null) { if(p -> v <= x) ans = p -> v, p = p -> ch[1]; else p = p -> ch[0]; } return ans;}int ask_suf(node *p, int x){ int ans = INF; while(p != null) { if(p -> v >= x) ans = p -> v, p = p -> ch[0]; else p = p -> ch[1]; } return ans;}int main(){ init(); int n; LL ans = 0; bool haha = 0; int dq = 1; scanf("%d", &n); while(n --) { int kind, num; scanf("%d%d", &kind, &num); if(kind) { if(!root -> sz || haha) { haha = 1; insert(root, num); } else if(root -> sz && !haha) { int pre = ask_pre(root, num); int suf = ask_suf(root, num); if(abs(num - pre) <= abs(suf - num)) { ans = ((ans % MOD) + (abs(num - pre)) % MOD) % MOD; erase(root, pre); } else if(abs(num - pre) > abs(suf - num)) { ans = ((ans % MOD) + (abs(suf - num)) % MOD) % MOD; erase(root, suf); } } } else { if(!root -> sz || !haha) { insert(root, num); haha = 0; } else if(root -> sz && haha) { int pre = ask_pre(root, num); int suf = ask_suf(root, num); if(abs(num - pre) <= abs(suf - num)) { ans = ((ans % MOD) + (abs(num - pre) % MOD)) % MOD; erase(root, pre); } else if(abs(num - pre) > abs(num - suf)) { ans = ((ans % MOD) + (abs(num - suf) % MOD)) % MOD; erase(root, suf); } } } } printf("%lld", ans % MOD); return 0;}

转载于:https://www.cnblogs.com/Loi-Vampire/p/6017058.html

你可能感兴趣的文章
我答"编程为什么不用中文?": 中文API的意义和探索
查看>>
Swoft 2.0.5 更新,新增高效秒级定时任务、异常管理组件
查看>>
版本对比 两个版本号 一样长度
查看>>
计算百分比
查看>>
js 乘除算法 浮点 精度解决办法
查看>>
sql server 2012 链接服务器不能链接sql server 2000的解决方案 ,
查看>>
sqlserver2005版本的mdf文件,还没有log文件,
查看>>
错误“该伙伴事务管理器已经禁止了它对远程/网络事务的支持”解决方案
查看>>
System x 服务器制作ServerGuide U盘安装Windows Server 2008 操作系统 --不格式化盘
查看>>
java面试
查看>>
前端常见跨域解决方案(全)
查看>>
单点登录(Single Sign On)解决方案
查看>>
umi---className设置多个样式
查看>>
网页包抓取工具Fiddler工具简单设置
查看>>
周总结报告
查看>>
Selecting Courses POJ - 2239(我是沙雕吧 按时间点建边 || 匹配水题)
查看>>
Win+R指令(2)
查看>>
codeforces 578c - weekness and poorness - 三分
查看>>
数值微分方程
查看>>
动态规划--电路布线(circuit layout)
查看>>