nuclear@16: #include <stdio.h>
nuclear@16: #include <stdlib.h>
nuclear@16: #include <string.h>
nuclear@16: #include "inttypes.h"
nuclear@16: #include "rbtree.h"
nuclear@16: 
nuclear@16: #define INT2PTR(x)	((void*)(intptr_t)(x))
nuclear@16: #define PTR2INT(x)	((int)(intptr_t)(x))
nuclear@16: 
nuclear@16: struct rbtree {
nuclear@16: 	struct rbnode *root;
nuclear@16: 
nuclear@16: 	rb_alloc_func_t alloc;
nuclear@16: 	rb_free_func_t free;
nuclear@16: 
nuclear@16: 	rb_cmp_func_t cmp;
nuclear@16: 	rb_del_func_t del;
nuclear@16: 	void *del_cls;
nuclear@16: 
nuclear@16: 	struct rbnode *rstack, *iter;
nuclear@16: };
nuclear@16: 
nuclear@16: static int cmpaddr(const void *ap, const void *bp);
nuclear@16: static int cmpint(const void *ap, const void *bp);
nuclear@16: 
nuclear@16: static int count_nodes(struct rbnode *node);
nuclear@16: static void del_tree(struct rbnode *node, void (*delfunc)(struct rbnode*, void*), void *cls);
nuclear@16: static struct rbnode *insert(struct rbtree *rb, struct rbnode *tree, void *key, void *data);
nuclear@16: static struct rbnode *delete(struct rbtree *rb, struct rbnode *tree, void *key);
nuclear@16: /*static struct rbnode *find(struct rbtree *rb, struct rbnode *node, void *key);*/
nuclear@16: static void traverse(struct rbnode *node, void (*func)(struct rbnode*, void*), void *cls);
nuclear@16: 
nuclear@16: struct rbtree *rb_create(rb_cmp_func_t cmp_func)
nuclear@16: {
nuclear@16: 	struct rbtree *rb;
nuclear@16: 
nuclear@16: 	if(!(rb = malloc(sizeof *rb))) {
nuclear@16: 		return 0;
nuclear@16: 	}
nuclear@16: 	if(rb_init(rb, cmp_func) == -1) {
nuclear@16: 		free(rb);
nuclear@16: 		return 0;
nuclear@16: 	}
nuclear@16: 	return rb;
nuclear@16: }
nuclear@16: 
nuclear@16: void rb_free(struct rbtree *rb)
nuclear@16: {
nuclear@16: 	rb_destroy(rb);
nuclear@16: 	free(rb);
nuclear@16: }
nuclear@16: 
nuclear@16: 
nuclear@16: int rb_init(struct rbtree *rb, rb_cmp_func_t cmp_func)
nuclear@16: {
nuclear@16: 	memset(rb, 0, sizeof *rb);
nuclear@16: 
nuclear@16: 	if(cmp_func == RB_KEY_INT) {
nuclear@16: 		rb->cmp = cmpint;
nuclear@16: 	} else if(cmp_func == RB_KEY_STRING) {
nuclear@16: 		rb->cmp = (rb_cmp_func_t)strcmp;
nuclear@16: 	} else {
nuclear@16: 		rb->cmp = cmpaddr;
nuclear@16: 	}
nuclear@16: 
nuclear@16: 	rb->alloc = malloc;
nuclear@16: 	rb->free = free;
nuclear@16: 	return 0;
nuclear@16: }
nuclear@16: 
nuclear@16: void rb_destroy(struct rbtree *rb)
nuclear@16: {
nuclear@16: 	del_tree(rb->root, rb->del, rb->del_cls);
nuclear@16: }
nuclear@16: 
nuclear@16: void rb_set_allocator(struct rbtree *rb, rb_alloc_func_t alloc, rb_free_func_t free)
nuclear@16: {
nuclear@16: 	rb->alloc = alloc;
nuclear@16: 	rb->free = free;
nuclear@16: }
nuclear@16: 
nuclear@16: 
nuclear@16: void rb_set_compare_func(struct rbtree *rb, rb_cmp_func_t func)
nuclear@16: {
nuclear@16: 	rb->cmp = func;
nuclear@16: }
nuclear@16: 
nuclear@16: void rb_set_delete_func(struct rbtree *rb, rb_del_func_t func, void *cls)
nuclear@16: {
nuclear@16: 	rb->del = func;
nuclear@16: 	rb->del_cls = cls;
nuclear@16: }
nuclear@16: 
nuclear@16: 
nuclear@16: void rb_clear(struct rbtree *rb)
nuclear@16: {
nuclear@16: 	del_tree(rb->root, rb->del, rb->del_cls);
nuclear@16: 	rb->root = 0;
nuclear@16: }
nuclear@16: 
nuclear@16: int rb_copy(struct rbtree *dest, struct rbtree *src)
nuclear@16: {
nuclear@16: 	struct rbnode *node;
nuclear@16: 
nuclear@16: 	rb_clear(dest);
nuclear@16: 	rb_begin(src);
nuclear@16: 	while((node = rb_next(src))) {
nuclear@16: 		if(rb_insert(dest, node->key, node->data) == -1) {
nuclear@16: 			return -1;
nuclear@16: 		}
nuclear@16: 	}
nuclear@16: 	return 0;
nuclear@16: }
nuclear@16: 
nuclear@16: int rb_size(struct rbtree *rb)
nuclear@16: {
nuclear@16: 	return count_nodes(rb->root);
nuclear@16: }
nuclear@16: 
nuclear@16: int rb_insert(struct rbtree *rb, void *key, void *data)
nuclear@16: {
nuclear@16: 	rb->root = insert(rb, rb->root, key, data);
nuclear@16: 	rb->root->red = 0;
nuclear@16: 	return 0;
nuclear@16: }
nuclear@16: 
nuclear@16: int rb_inserti(struct rbtree *rb, int key, void *data)
nuclear@16: {
nuclear@16: 	rb->root = insert(rb, rb->root, INT2PTR(key), data);
nuclear@16: 	rb->root->red = 0;
nuclear@16: 	return 0;
nuclear@16: }
nuclear@16: 
nuclear@16: 
nuclear@16: int rb_delete(struct rbtree *rb, void *key)
nuclear@16: {
nuclear@16: 	rb->root = delete(rb, rb->root, key);
nuclear@16: 	rb->root->red = 0;
nuclear@16: 	return 0;
nuclear@16: }
nuclear@16: 
nuclear@16: int rb_deletei(struct rbtree *rb, int key)
nuclear@16: {
nuclear@16: 	rb->root = delete(rb, rb->root, INT2PTR(key));
nuclear@16: 	rb->root->red = 0;
nuclear@16: 	return 0;
nuclear@16: }
nuclear@16: 
nuclear@16: 
nuclear@16: struct rbnode *rb_find(struct rbtree *rb, void *key)
nuclear@16: {
nuclear@16: 	struct rbnode *node = rb->root;
nuclear@16: 
nuclear@16: 	while(node) {
nuclear@16: 		int cmp = rb->cmp(key, node->key);
nuclear@16: 		if(cmp == 0) {
nuclear@16: 			return node;
nuclear@16: 		}
nuclear@16: 		node = cmp < 0 ? node->left : node->right;
nuclear@16: 	}
nuclear@16: 	return 0;
nuclear@16: }
nuclear@16: 
nuclear@16: struct rbnode *rb_findi(struct rbtree *rb, int key)
nuclear@16: {
nuclear@16: 	return rb_find(rb, INT2PTR(key));
nuclear@16: }
nuclear@16: 
nuclear@16: 
nuclear@16: void rb_foreach(struct rbtree *rb, void (*func)(struct rbnode*, void*), void *cls)
nuclear@16: {
nuclear@16: 	traverse(rb->root, func, cls);
nuclear@16: }
nuclear@16: 
nuclear@16: 
nuclear@16: struct rbnode *rb_root(struct rbtree *rb)
nuclear@16: {
nuclear@16: 	return rb->root;
nuclear@16: }
nuclear@16: 
nuclear@16: void rb_begin(struct rbtree *rb)
nuclear@16: {
nuclear@16: 	rb->rstack = 0;
nuclear@16: 	rb->iter = rb->root;
nuclear@16: }
nuclear@16: 
nuclear@16: #define push(sp, x)		((x)->next = (sp), (sp) = (x))
nuclear@16: #define pop(sp)			((sp) = (sp)->next)
nuclear@16: #define top(sp)			(sp)
nuclear@16: 
nuclear@16: struct rbnode *rb_next(struct rbtree *rb)
nuclear@16: {
nuclear@16: 	struct rbnode *res = 0;
nuclear@16: 
nuclear@16: 	while(rb->rstack || rb->iter) {
nuclear@16: 		if(rb->iter) {
nuclear@16: 			push(rb->rstack, rb->iter);
nuclear@16: 			rb->iter = rb->iter->left;
nuclear@16: 		} else {
nuclear@16: 			rb->iter = top(rb->rstack);
nuclear@16: 			pop(rb->rstack);
nuclear@16: 			res = rb->iter;
nuclear@16: 			rb->iter = rb->iter->right;
nuclear@16: 			break;
nuclear@16: 		}
nuclear@16: 	}
nuclear@16: 	return res;
nuclear@16: }
nuclear@16: 
nuclear@16: void *rb_node_key(struct rbnode *node)
nuclear@16: {
nuclear@16: 	return node ? node->key : 0;
nuclear@16: }
nuclear@16: 
nuclear@16: int rb_node_keyi(struct rbnode *node)
nuclear@16: {
nuclear@16: 	return node ? PTR2INT(node->key) : 0;
nuclear@16: }
nuclear@16: 
nuclear@16: void *rb_node_data(struct rbnode *node)
nuclear@16: {
nuclear@16: 	return node ? node->data : 0;
nuclear@16: }
nuclear@16: 
nuclear@16: static int cmpaddr(const void *ap, const void *bp)
nuclear@16: {
nuclear@16: 	return ap < bp ? -1 : (ap > bp ? 1 : 0);
nuclear@16: }
nuclear@16: 
nuclear@16: static int cmpint(const void *ap, const void *bp)
nuclear@16: {
nuclear@16: 	return PTR2INT(ap) - PTR2INT(bp);
nuclear@16: }
nuclear@16: 
nuclear@16: 
nuclear@16: /* ---- left-leaning 2-3 red-black implementation ---- */
nuclear@16: 
nuclear@16: /* helper prototypes */
nuclear@16: static int is_red(struct rbnode *tree);
nuclear@16: static void color_flip(struct rbnode *tree);
nuclear@16: static struct rbnode *rot_left(struct rbnode *a);
nuclear@16: static struct rbnode *rot_right(struct rbnode *a);
nuclear@16: static struct rbnode *find_min(struct rbnode *tree);
nuclear@16: static struct rbnode *del_min(struct rbtree *rb, struct rbnode *tree);
nuclear@16: /*static struct rbnode *move_red_right(struct rbnode *tree);*/
nuclear@16: static struct rbnode *move_red_left(struct rbnode *tree);
nuclear@16: static struct rbnode *fix_up(struct rbnode *tree);
nuclear@16: 
nuclear@16: static int count_nodes(struct rbnode *node)
nuclear@16: {
nuclear@16: 	if(!node)
nuclear@16: 		return 0;
nuclear@16: 
nuclear@16: 	return 1 + count_nodes(node->left) + count_nodes(node->right);
nuclear@16: }
nuclear@16: 
nuclear@16: static void del_tree(struct rbnode *node, rb_del_func_t delfunc, void *cls)
nuclear@16: {
nuclear@16: 	if(!node)
nuclear@16: 		return;
nuclear@16: 
nuclear@16: 	del_tree(node->left, delfunc, cls);
nuclear@16: 	del_tree(node->right, delfunc, cls);
nuclear@16: 
nuclear@16: 	if(delfunc) {
nuclear@16: 		delfunc(node, cls);
nuclear@16: 	}
nuclear@16: 	free(node);
nuclear@16: }
nuclear@16: 
nuclear@16: static struct rbnode *insert(struct rbtree *rb, struct rbnode *tree, void *key, void *data)
nuclear@16: {
nuclear@16: 	int cmp;
nuclear@16: 
nuclear@16: 	if(!tree) {
nuclear@16: 		struct rbnode *node = rb->alloc(sizeof *node);
nuclear@16: 		node->red = 1;
nuclear@16: 		node->key = key;
nuclear@16: 		node->data = data;
nuclear@16: 		node->left = node->right = 0;
nuclear@16: 		return node;
nuclear@16: 	}
nuclear@16: 
nuclear@16: 	cmp = rb->cmp(key, tree->key);
nuclear@16: 
nuclear@16: 	if(cmp < 0) {
nuclear@16: 		tree->left = insert(rb, tree->left, key, data);
nuclear@16: 	} else if(cmp > 0) {
nuclear@16: 		tree->right = insert(rb, tree->right, key, data);
nuclear@16: 	} else {
nuclear@16: 		tree->data = data;
nuclear@16: 	}
nuclear@16: 
nuclear@16: 	/* fix right-leaning reds */
nuclear@16: 	if(is_red(tree->right)) {
nuclear@16: 		tree = rot_left(tree);
nuclear@16: 	}
nuclear@16: 	/* fix two reds in a row */
nuclear@16: 	if(is_red(tree->left) && is_red(tree->left->left)) {
nuclear@16: 		tree = rot_right(tree);
nuclear@16: 	}
nuclear@16: 
nuclear@16: 	/* if 4-node, split it by color inversion */
nuclear@16: 	if(is_red(tree->left) && is_red(tree->right)) {
nuclear@16: 		color_flip(tree);
nuclear@16: 	}
nuclear@16: 
nuclear@16: 	return tree;
nuclear@16: }
nuclear@16: 
nuclear@16: static struct rbnode *delete(struct rbtree *rb, struct rbnode *tree, void *key)
nuclear@16: {
nuclear@16: 	int cmp;
nuclear@16: 
nuclear@16: 	if(!tree) {
nuclear@16: 		return 0;
nuclear@16: 	}
nuclear@16: 
nuclear@16: 	cmp = rb->cmp(key, tree->key);
nuclear@16: 
nuclear@16: 	if(cmp < 0) {
nuclear@16: 		if(!is_red(tree->left) && !is_red(tree->left->left)) {
nuclear@16: 			tree = move_red_left(tree);
nuclear@16: 		}
nuclear@16: 		tree->left = delete(rb, tree->left, key);
nuclear@16: 	} else {
nuclear@16: 		/* need reds on the right */
nuclear@16: 		if(is_red(tree->left)) {
nuclear@16: 			tree = rot_right(tree);
nuclear@16: 		}
nuclear@16: 
nuclear@16: 		/* found it at the bottom (XXX what certifies left is null?) */
nuclear@16: 		if(cmp == 0 && !tree->right) {
nuclear@16: 			if(rb->del) {
nuclear@16: 				rb->del(tree, rb->del_cls);
nuclear@16: 			}
nuclear@16: 			rb->free(tree);
nuclear@16: 			return 0;
nuclear@16: 		}
nuclear@16: 
nuclear@16: 		if(!is_red(tree->right) && !is_red(tree->right->left)) {
nuclear@16: 			tree = move_red_left(tree);
nuclear@16: 		}
nuclear@16: 
nuclear@16: 		if(key == tree->key) {
nuclear@16: 			struct rbnode *rmin = find_min(tree->right);
nuclear@16: 			tree->key = rmin->key;
nuclear@16: 			tree->data = rmin->data;
nuclear@16: 			tree->right = del_min(rb, tree->right);
nuclear@16: 		} else {
nuclear@16: 			tree->right = delete(rb, tree->right, key);
nuclear@16: 		}
nuclear@16: 	}
nuclear@16: 
nuclear@16: 	return fix_up(tree);
nuclear@16: }
nuclear@16: 
nuclear@16: /*static struct rbnode *find(struct rbtree *rb, struct rbnode *node, void *key)
nuclear@16: {
nuclear@16: 	int cmp;
nuclear@16: 
nuclear@16: 	if(!node)
nuclear@16: 		return 0;
nuclear@16: 
nuclear@16: 	if((cmp = rb->cmp(key, node->key)) == 0) {
nuclear@16: 		return node;
nuclear@16: 	}
nuclear@16: 	return find(rb, cmp < 0 ? node->left : node->right, key);
nuclear@16: }*/
nuclear@16: 
nuclear@16: static void traverse(struct rbnode *node, void (*func)(struct rbnode*, void*), void *cls)
nuclear@16: {
nuclear@16: 	if(!node)
nuclear@16: 		return;
nuclear@16: 
nuclear@16: 	traverse(node->left, func, cls);
nuclear@16: 	func(node, cls);
nuclear@16: 	traverse(node->right, func, cls);
nuclear@16: }
nuclear@16: 
nuclear@16: /* helpers */
nuclear@16: 
nuclear@16: static int is_red(struct rbnode *tree)
nuclear@16: {
nuclear@16: 	return tree && tree->red;
nuclear@16: }
nuclear@16: 
nuclear@16: static void color_flip(struct rbnode *tree)
nuclear@16: {
nuclear@16: 	tree->red = !tree->red;
nuclear@16: 	tree->left->red = !tree->left->red;
nuclear@16: 	tree->right->red = !tree->right->red;
nuclear@16: }
nuclear@16: 
nuclear@16: static struct rbnode *rot_left(struct rbnode *a)
nuclear@16: {
nuclear@16: 	struct rbnode *b = a->right;
nuclear@16: 	a->right = b->left;
nuclear@16: 	b->left = a;
nuclear@16: 	b->red = a->red;
nuclear@16: 	a->red = 1;
nuclear@16: 	return b;
nuclear@16: }
nuclear@16: 
nuclear@16: static struct rbnode *rot_right(struct rbnode *a)
nuclear@16: {
nuclear@16: 	struct rbnode *b = a->left;
nuclear@16: 	a->left = b->right;
nuclear@16: 	b->right = a;
nuclear@16: 	b->red = a->red;
nuclear@16: 	a->red = 1;
nuclear@16: 	return b;
nuclear@16: }
nuclear@16: 
nuclear@16: static struct rbnode *find_min(struct rbnode *tree)
nuclear@16: {
nuclear@16: 	if(!tree)
nuclear@16: 		return 0;
nuclear@16: 
nuclear@16: 	while(tree->left) {
nuclear@16: 		tree = tree->left;
nuclear@16: 	}
nuclear@16: 	return tree;
nuclear@16: }
nuclear@16: 
nuclear@16: static struct rbnode *del_min(struct rbtree *rb, struct rbnode *tree)
nuclear@16: {
nuclear@16: 	if(!tree->left) {
nuclear@16: 		if(rb->del) {
nuclear@16: 			rb->del(tree->left, rb->del_cls);
nuclear@16: 		}
nuclear@16: 		rb->free(tree->left);
nuclear@16: 		return 0;
nuclear@16: 	}
nuclear@16: 
nuclear@16: 	/* make sure we've got red (3/4-nodes) at the left side so we can delete at the bottom */
nuclear@16: 	if(!is_red(tree->left) && !is_red(tree->left->left)) {
nuclear@16: 		tree = move_red_left(tree);
nuclear@16: 	}
nuclear@16: 	tree->left = del_min(rb, tree->left);
nuclear@16: 
nuclear@16: 	/* fix right-reds, red-reds, and split 4-nodes on the way up */
nuclear@16: 	return fix_up(tree);
nuclear@16: }
nuclear@16: 
nuclear@16: #if 0
nuclear@16: /* push a red link on this node to the right */
nuclear@16: static struct rbnode *move_red_right(struct rbnode *tree)
nuclear@16: {
nuclear@16: 	/* flipping it makes both children go red, so we have a red to the right */
nuclear@16: 	color_flip(tree);
nuclear@16: 
nuclear@16: 	/* if after the flip we've got a red-red situation to the left, fix it */
nuclear@16: 	if(is_red(tree->left->left)) {
nuclear@16: 		tree = rot_right(tree);
nuclear@16: 		color_flip(tree);
nuclear@16: 	}
nuclear@16: 	return tree;
nuclear@16: }
nuclear@16: #endif
nuclear@16: 
nuclear@16: /* push a red link on this node to the left */
nuclear@16: static struct rbnode *move_red_left(struct rbnode *tree)
nuclear@16: {
nuclear@16: 	/* flipping it makes both children go red, so we have a red to the left */
nuclear@16: 	color_flip(tree);
nuclear@16: 
nuclear@16: 	/* if after the flip we've got a red-red on the right-left, fix it */
nuclear@16: 	if(is_red(tree->right->left)) {
nuclear@16: 		tree->right = rot_right(tree->right);
nuclear@16: 		tree = rot_left(tree);
nuclear@16: 		color_flip(tree);
nuclear@16: 	}
nuclear@16: 	return tree;
nuclear@16: }
nuclear@16: 
nuclear@16: static struct rbnode *fix_up(struct rbnode *tree)
nuclear@16: {
nuclear@16: 	/* fix right-leaning */
nuclear@16: 	if(is_red(tree->right)) {
nuclear@16: 		tree = rot_left(tree);
nuclear@16: 	}
nuclear@16: 	/* change invalid red-red pairs into a proper 4-node */
nuclear@16: 	if(is_red(tree->left) && is_red(tree->left->left)) {
nuclear@16: 		tree = rot_right(tree);
nuclear@16: 	}
nuclear@16: 	/* split 4-nodes */
nuclear@16: 	if(is_red(tree->left) && is_red(tree->right)) {
nuclear@16: 		color_flip(tree);
nuclear@16: 	}
nuclear@16: 	return tree;
nuclear@16: }