#include #include #include #include #define MAX(a, b) ((a) > (b) ? (a) : (b)) #define MIN(a, b) ((a) < (b) ? (a) : (b)) #define IS_LEFT_CHILD(p, c) ((p) && (c) && ((p)->b_left == (c))) #define IS_RIGHT_CHILD(p, c) ((p) && (c) && ((p)->b_right == (c))) #define HAS_LEFT_CHILD(x) ((x) && ((x)->b_left)) #define HAS_RIGHT_CHILD(x) ((x) && ((x)->b_right)) #define HAS_NO_CHILDREN(x) ((x) && (!(x)->b_left) && (!(x)->b_right)) #define HAS_ONE_CHILD(x) ((HAS_LEFT_CHILD(x) && !HAS_RIGHT_CHILD(x)) || (!HAS_LEFT_CHILD(x) && HAS_RIGHT_CHILD(x))) #define HAS_TWO_CHILDREN(x) (HAS_LEFT_CHILD(x) && HAS_RIGHT_CHILD(x)) #define HEIGHT(x) ((x) ? (x)->b_height : 0) static inline void update_height(btree_node_t *x) { x->b_height = MAX(HEIGHT(x->b_left), HEIGHT((x->b_right))) + 1; } static inline int bf(btree_node_t *x) { int bf = 0; if (!x) { return bf; } if (x->b_right) { bf += x->b_right->b_height; } if (x->b_left) { bf -= x->b_left->b_height; } return bf; } static void rotate_left(btree_t *tree, btree_node_t *x) { btree_node_t *y = x->b_right; btree_node_t *p = x->b_parent; if (y->b_left) { y->b_left->b_parent = x; } x->b_right = y->b_left; if (!p) { tree->b_root = y; } else if (x == p->b_left) { p->b_left = y; } else { p->b_right = y; } x->b_parent = y; y->b_left = x; y->b_parent = p; p = x; while (p) { update_height(p); p = p->b_parent; } } static void rotate_right(btree_t *tree, btree_node_t *y) { btree_node_t *x = y->b_left; btree_node_t *p = y->b_parent; if (x->b_right) { x->b_right->b_parent = y; } y->b_left = x->b_right; if (!p) { tree->b_root = x; } else if (x == p->b_left) { p->b_left = x; } else { p->b_right = x; } y->b_parent = x; x->b_right = y; x->b_parent = p; p = y; while (p) { update_height(p); p = p->b_parent; } } static void rotate_double_left(btree_t *tree, btree_node_t *z) { btree_node_t *x = z->b_right; btree_node_t *y = x->b_left; rotate_right(tree, x); rotate_left(tree, z); update_height(z); update_height(x); update_height(y); btree_node_t *p = y->b_parent; while (p) { p->b_height--; p = p->b_parent; } } static void rotate_double_right(btree_t *tree, btree_node_t *z) { btree_node_t *x = z->b_left; btree_node_t *y = x->b_right; rotate_left(tree, x); rotate_right(tree, z); update_height(z); update_height(x); update_height(y); btree_node_t *p = y->b_parent; while (p) { p->b_height--; p = p->b_parent; } } static void fix_tree(btree_t *tree, btree_node_t *w) { btree_node_t *z = NULL, *y = NULL, *x = NULL; z = w; while (z) { if (bf(z) >= -1 && bf(z) <= 1) { goto next_ancestor; } if (IS_LEFT_CHILD(z, y)) { if (IS_LEFT_CHILD(y, x)) { rotate_right(tree, z); } else { rotate_double_right(tree, z); } } else { if (IS_LEFT_CHILD(y, x)) { rotate_double_left(tree, z); } else { rotate_left(tree, z); } } next_ancestor: x = y; y = z; z = z->b_parent; } } void btree_insert(btree_t *tree, btree_node_t *node) { if (!tree->b_root) { tree->b_root = node; node->b_parent = NULL; node->b_height = 1; return; } btree_node_t *cur = tree->b_root; while (1) { btree_node_t **nextp = NULL; if (node->b_key >= cur->b_key) { nextp = &cur->b_right; } else { nextp = &cur->b_left; } if (*nextp) { cur = *nextp; } else { node->b_parent = cur; *nextp = node; break; } } node->b_height = 0; cur = node; while (cur) { cur->b_height++; cur = cur->b_parent; } fix_tree(tree, node); } static btree_node_t *remove_node_with_no_children(btree_t *tree, btree_node_t *node) { printf("remove_node_with_no_children()\n"); btree_node_t *w = node->b_parent; btree_node_t *p = node->b_parent; node->b_parent = NULL; if (!p) { tree->b_root = NULL; } else if (IS_LEFT_CHILD(p, node)) { p->b_left = NULL; } else { p->b_right = NULL; } while (p) { update_height(p); p = p->b_parent; } return w; } static btree_node_t *replace_node_with_one_subtree(btree_t *tree, btree_node_t *node) { printf("remove_node_with_one_subtree()\n"); btree_node_t *p = node->b_parent; btree_node_t *z = NULL; if (HAS_LEFT_CHILD(node)) { z = node->b_left; } else { z = node->b_right; } btree_node_t *w = z; if (IS_LEFT_CHILD(p, node)) { p->b_left = z; } else if (IS_RIGHT_CHILD(p, node)) { p->b_right = z; } z->b_parent = p; node->b_parent = NULL; node->b_left = node->b_right = NULL; while (z) { update_height(z); z = z->b_parent; } return w; } static btree_node_t *replace_node_with_two_subtrees(btree_t *tree, btree_node_t *node) { printf("remove_node_with_two_subtrees()\n"); btree_node_t *cur = node->b_left; while (cur->b_right) { cur = cur->b_right; } btree_node_t *p = cur->b_parent; btree_node_t *k = p; if (IS_LEFT_CHILD(p, cur)) { p->b_left = NULL; } else { p->b_right = NULL; } p = node->b_parent; if (p) { if (IS_LEFT_CHILD(p, node)) { p->b_left = cur; } else { p->b_right = cur; } } cur->b_left = node->b_left; cur->b_right = node->b_right; cur->b_height = node->b_height; cur->b_parent = p; btree_node_t *w = cur; if (!p) { tree->b_root = cur; } if (tree->b_root->b_left) { tree->b_root->b_left->b_parent = tree->b_root; } if (tree->b_root->b_right) { tree->b_root->b_right->b_parent = tree->b_root; } while (k) { update_height(k); k = k->b_parent; } return w; } void btree_delete(btree_t *tree, btree_node_t *node) { btree_node_t *w = NULL; if (HAS_NO_CHILDREN(node)) { w = remove_node_with_no_children(tree, node); } else if (HAS_ONE_CHILD(node)) { w = replace_node_with_one_subtree(tree, node); } else if (HAS_TWO_CHILDREN(node)) { w = replace_node_with_two_subtrees(tree, node); } fix_tree(tree, w); }