bzoj1208[HDOI2004] 宠物收养所

Bzoj1208 [HNOI2004]宠物收养所 解题报告

Time Limit: 10 Sec   Memory Limit: 162 MB
Submit: 6443   Solved: 2510
[ Submit][ Status][ Discuss]

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

Analysis

      算法一：Splay或Treap

      建两棵平衡树，一棵是狗树，一棵是人树。当来了一只新的狗时，看看人树是不是空的，如果是空的，那就将狗插入狗树；如果不是空的，那就查询狗特点值的前驱和后继，然后判断那个应该领取它，删掉那个节点。人树同理。

      算法二：范浩强Treap

      当来了一只狗时，将人树根据狗的特点值split成两半，在找到左树中最大的点，右树中最小的点，比较该删哪个，如果该删的点在左边，就merge(split_l(左树, 特点值-1),右树)，右边对称。如果人树是空的，就将狗insert到狗树里去，人树对称。

      感觉fhq Treap和其他平衡树完全不在一个频道。。。

My Code

//bzoj1208[HDOI2004]宠物收养所 fhq Treap
#include <cstdio>
#include <algorithm>
#include <cstdlib>
#define inf 0x7fffffff
using namespace std;
struct NODE
{
	int k, w;NODE *l, *r;
	NODE(int kk, int ww, NODE* ll, NODE* rr){k=kk,w=ww;l=ll,r=rr;}
}*ROOT[2];
int N;
long long ans;
NODE* merge(NODE *a, NODE *b)
{
	return (!a||!b) ? (a?a:b) :
		a->w < b->w ? new NODE(a->k,a->w,a->l,merge(a->r,b))
					: new NODE(b->k,b->w,merge(a,b->l),b->r);
}
NODE* split_l(NODE *x, int k)
{
	return !x ? 0 :
		x->k <= k ? new NODE(x->k,x->w,x->l,split_l(x->r,k))
			  : split_l(x->l,k);
}
NODE* split_r(NODE *x, int k)
{
	return !x ? 0 :
		x->k > k ? new NODE(x->k,x->w,split_r(x->l,k),x->r)
			  : split_r(x->r,k);
}
NODE* insert(NODE *x, int k, int w)
{
	return (!x||x->w>w) ? new NODE(k,w,split_l(x,k),split_r(x,k)) :
		x->k < k ? new NODE(x->k,x->w,x->l,insert(x->r,k,w))
				 : new NODE(x->k,x->w,insert(x->l,k,w),x->r);
}
int find(NODE *x, bool f)
{
	return !x ? -1 :
		f ? (x->r?find(x->r,f):x->k)
		  : (x->l?find(x->l,f):x->k);
}
int main()
{
	int i, a, b, d1, d2, x, y;
	NODE *p1, *p2, *p;
	while(scanf("%d",&N)!=EOF)
	{
		ans=0;
		ROOT[0]=new NODE(0,0,0,0);
		ROOT[1]=new NODE(0,0,0,0);
		for(i=1;i<=N;i++)
		{
			scanf("%d%d",&a,&b);
			if((p=ROOT[!a])->l)
			{
				p1=split_l(p->l,b), p2=split_r(p->l,b);
				x=find(p1,1), y=find(p2,0);
				d1=x==-1?inf:abs(x-b);
				d2=y==-1?inf:abs(y-b);
				if(d2<d1)ans=(ans+d2)%1000000, p->l=merge(p1,split_r(p2,y));
				else ans=(ans+d1)%1000000, p->l=merge(split_l(p1,x-1),p2);
			}
			else ROOT[a]->l=insert(ROOT[a]->l,b,rand());
		}
		printf("%d\n",ans%1000000);
	}
	return 0;
}