mango/sandbox/btree/btree.c

/*
  The Clear BSD License

  Copyright (c) 2023 Max Wash
  All rights reserved.

  Redistribution and use in source and binary forms, with or without
  modification, are permitted (subject to the limitations in the disclaimer
  below) provided that the following conditions are met:

  - Redistributions of source code must retain the above copyright notice,
    this list of conditions and the following disclaimer.

  - Redistributions in binary form must reproduce the above copyright
    notice, this list of conditions and the following disclaimer in the
    documentation and/or other materials provided with the distribution.
  
  - Neither the name of the copyright holder nor the names of its
    contributors may be used to endorse or promote products derived from this
    software without specific prior written permission.
 */

/* templated AVL binary tree implementation

   this file implements an extensible AVL binary tree data structure.

   the primary rule of an AVL binary tree is that for a given node N,
   the heights of N's left and right subtrees can differ by at most 1.

   the height of a subtree is the length of the longest path between
   the root of the subtree and a leaf node, including the root node itself.

   the height of a leaf node is 1.

   when a node is inserted into or deleted from the tree, this rule may
   be broken, in which the tree must be rotated to restore the balance.

   no more than one rotation is required for any insert operations,
   while multiple rotations may be required for a delete operation.

   there are four types of rotations that can be applied to a tree:
     - left rotation
     - right rotation
     - double left rotations
     - double right rotations

   by enforcing the balance rule, for a tree with n nodes, the worst-case
   performance for insert, delete, and search operations is guaranteed
   to be O(log n).

   this file intentionally excludes any kind of search function implementation.
   it is up to the programmer to implement their own tree node type
   using btree_node_t, and their own search function using btree_t.
   this allows the programmer to define their own node types with complex
   non-integer key types. btree.h contains a number of macros to help
   define these functions. the macros do all the work, you just have to
   provide a comparator function.
*/ 
   
#include <socks/btree.h>
#include <stddef.h>
#include <stdlib.h>
#include <stdio.h>
#include <assert.h>

#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;
}

/* perform a left rotation on a subtree

   if you have a tree like this:

        Z
       / \
      X   .
     / \
    .   Y
       / \
      .   .

   and you perform a left rotation on node X,
   you will get the following tree:

        Z
       / \
      Y   .
     / \
    X   .
   / \
  .   .
  
  note that this function does NOT update b_height for the rotated
  nodes. it is up to you to call update_height_to_root().
*/
static void rotate_left(btree_t *tree, btree_node_t *x)
{
	assert(x != NULL);

	btree_node_t *y = x->b_right;
	assert(y != NULL);

	assert(y == x->b_left || y == x->b_right);
	if (x->b_parent) {
		assert(x == x->b_parent->b_left || x == x->b_parent->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;
}

static void update_height_to_root(btree_node_t *x)
{
	while (x) {
		update_height(x);
		x = x->b_parent;
	}
}

/* perform a right rotation on a subtree

   if you have a tree like this:

        Z
       / \
      .   X
         / \
        Y   .
       / \
      .   .

   and you perform a right rotation on node X,
   you will get the following tree:

        Z
       / \
      .   Y
         / \
        .   X
           / \
	  .   .
	  
  note that this function does NOT update b_height for the rotated
  nodes. it is up to you to call update_height_to_root().
*/
static void rotate_right(btree_t *tree, btree_node_t *y)
{
	assert(y);

	btree_node_t *x = y->b_left;
	assert(x);

	assert(x == y->b_left || x == y->b_right);
	if (y->b_parent) {
		assert(y == y->b_parent->b_left || y == y->b_parent->b_right);
	}

	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 (y == p->b_left) {
		p->b_left = x;
	} else {
		p->b_right = x;
	}

	y->b_parent = x;
	x->b_right = y;
	x->b_parent = p;
}

/* for a given node Z, perform a right rotation on Z's right child,
   followed by a left rotation on Z itself.

   if you have a tree like this:

        Z
       / \
      .   X
         / \
        Y   .
       / \
      .   .

   and you perform a double-left rotation on node Z,
   you will get the following tree:

         Y
        / \ 
       /   \
      Z     X
     / \   / \
    .   . .   .
	  
  note that, unlike rotate_left and rotate_right, this function
  DOES update b_height for the rotated nodes (since it needs to be
  done in a certain order).
*/
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);

	while (y) {
		update_height(y);
		y = y->b_parent;
	}
}

/* for a given node Z, perform a left rotation on Z's left child,
   followed by a right rotation on Z itself.

   if you have a tree like this:

        Z
       / \
      X   .
     / \
    .   Y
       / \
      .   .

   and you perform a double-right rotation on node Z,
   you will get the following tree:

         Y
        / \ 
       /   \
      X     Z
     / \   / \
    .   . .   .
	  
  note that, unlike rotate_left and rotate_right, this function
  DOES update b_height for the rotated nodes (since it needs to be
  done in a certain order).
*/
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);

	while (y) {
		update_height(y);
		y = y->b_parent;
	}
}

/* run after an insert operation. checks that the balance factor
   of the local subtree is within the range -1 <= BF <= 1. if it
   is not, rotate the subtree to restore balance.

   note that at most one rotation should be required after a node
   is inserted into the tree.

   this function depends on all nodes in the tree having
   correct b_height values.

   @param w the node that was just inserted into the tree
*/
static void insert_fixup(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;
		}

		assert(x && y && z);
		assert(x == y->b_left || x == y->b_right);
		assert(y == z->b_left || y == z->b_right);

		if (IS_LEFT_CHILD(z, y)) {
			if (IS_LEFT_CHILD(y, x)) {
			    rotate_right(tree, z);
			    update_height_to_root(z);
			} else {
			    rotate_double_right(tree, z);
			}
		} else {
			if (IS_LEFT_CHILD(y, x)) {
			    rotate_double_left(tree, z);
			} else {
			    rotate_left(tree, z);
			    update_height_to_root(z);
			}
		}

next_ancestor:
		x = y;
		y = z;
		z = z->b_parent;
	}
}

/* run after a delete operation. checks that the balance factor
   of the local subtree is within the range -1 <= BF <= 1. if it
   is not, rotate the subtree to restore balance.

   note that, unlike insert_fixup, multiple rotations may be required
   to restore balance after a node is deleted.

   this function depends on all nodes in the tree having
   correct b_height values.

   @param w one of the following:
            - the parent of the node that was deleted if the node
	      had no children.
            - the parent of the node that replaced the deleted node
	      if the deleted node had two children.
	    - the node that replaced the node that was deleted, if
	      the node that was deleted had one child.
*/
static void delete_fixup(btree_t *tree, btree_node_t *w)
{
	btree_node_t *z = w;

	while (z) {
		if (bf(z) > 1) {
			if (bf(z->b_right) >= 0) {
				rotate_left(tree, z);
				update_height_to_root(z);
			} else {
				rotate_double_left(tree, z);
			}
		} else if (bf(z) < -1) {
			if (bf(z->b_left) <= 0) {
				rotate_right(tree, z);
				update_height_to_root(z);
			} else {
				rotate_double_right(tree, z);
			}
		}

		z = z->b_parent;
	}
}

/* updates b_height for all nodes between the inserted node and the root
   of the tree, and calls insert_fixup.

   @param node the node that was just inserted into the tree.
*/
void btree_insert_fixup(btree_t *tree, btree_node_t *node)
{
	node->b_height = 0;

	btree_node_t *cur = node;
	while (cur) {
		update_height(cur);
		cur = cur->b_parent;
	}

	insert_fixup(tree, node);
}

/* remove a node from a tree.

   this function assumes that `node` has no children, and therefore
   doesn't need to be replaced.

   updates b_height for all nodes between `node` and the tree root.

   @param node the node to delete.
*/
static btree_node_t *remove_node_with_no_children(btree_t *tree, btree_node_t *node)
{
	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;
}

/* remove a node from a tree.

   this function assumes that `node` has one child.
   the child of `node` is inherited by `node`'s parent, and `node` is removed.

   updates b_height for all nodes between the node that replaced
   `node` and the tree root.

   @param node the node to delete.
*/
static btree_node_t *replace_node_with_one_subtree(btree_t *tree, btree_node_t *node)
{
	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 (!p) {
		tree->b_root = z;
	} else 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;
}

/* remove a node from a tree.

   this function assumes that `node` has two children.
   find the in-order successor Y of `node` (the largest node in `node`'s left sub-tree),
   removes `node` from the tree and moves Y to where `node` used to be.

   if Y has a child (it will never have more than one), have Y's parent inherit
   Y's child.

   updates b_height for all nodes between the deepest node that was modified
   and the tree root.

   @param z the node to delete.
*/
static btree_node_t *replace_node_with_two_subtrees(btree_t *tree, btree_node_t *z)
{
	/* x will replace z */
	btree_node_t *x = z->b_left;

	while (x->b_right) {
		x = x->b_right;
	}

	/* y is the node that will replace x (if x has a left child) */
	btree_node_t *y = x->b_left;

	/* w is the starting point for the height update and fixup */
	btree_node_t *w = x;
	if (w->b_parent != z) {
		w = w->b_parent;
	}

	if (y) {
		w = y;
	}

	if (IS_LEFT_CHILD(x->b_parent, x)) {
		x->b_parent->b_left = y;
	} else if (IS_RIGHT_CHILD(x->b_parent, x)) {
		x->b_parent->b_right = y;
	}

	if (y) {
	    y->b_parent = x->b_parent;
	}

	if (IS_LEFT_CHILD(z->b_parent, z)) {
		z->b_parent->b_left = x;
	} else if (IS_RIGHT_CHILD(z->b_parent, z)) {
		z->b_parent->b_right = x;
	}

	x->b_parent = z->b_parent;
	x->b_left = z->b_left;
	x->b_right = z->b_right;

	if (x->b_left) {
		x->b_left->b_parent = x;
	}

	if (x->b_right) {
		x->b_right->b_parent = x;
	}

	if (!x->b_parent) {
		tree->b_root = x;
	}

	btree_node_t *cur = w;
	while (cur) {
		update_height(cur);
		cur = cur->b_parent;
	}

	return w;
}

/* delete a node from the tree and re-balance it afterwards */
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);
	}

	if (w) {
		delete_fixup(tree, w);
	}

	node->b_left = node->b_right = node->b_parent = NULL;
}

btree_node_t *btree_first(btree_t *tree)
{
	/* the first node in the tree is the node with the smallest key.
	   we keep moving left until we can't go any further */
	btree_node_t *cur = tree->b_root;
	if (!cur) {
		return NULL;
	}
	
	while (cur->b_left) {
		cur = cur->b_left;
	}

	return cur;
}

btree_node_t *btree_last(btree_t *tree)
{
	/* the first node in the tree is the node with the largest key.
	   we keep moving right until we can't go any further */
	btree_node_t *cur = tree->b_root;
	if (!cur) {
		return NULL;
	}
	
	while (cur->b_right) {
		cur = cur->b_right;
	}

	return cur;
}

btree_node_t *btree_next(btree_node_t *node)
{
	if (!node) {
		return NULL;
	}

	/* there are two possibilities for the next node:

	   1. if `node` has a right sub-tree, every node in this sub-tree is bigger
	      than node. the in-order successor of `node` is the smallest node in
	      this subtree.
	   2. if `node` has no right sub-tree, we've reached the largest node in
	      the sub-tree rooted at `node`. we need to go back to our parent
	      and continue the search elsewhere.
	*/
	if (node->b_right) {
		/* case 1: step into `node`'s right sub-tree and keep going
		   left to find the smallest node */
		btree_node_t *cur = node->b_right;
		while (cur->b_left) {
			cur = cur->b_left;
		}

		return cur;
	}

	/* case 2: keep stepping back up towards the root of the tree.
	   if we encounter a step where we are our parent's left child,
	   we've found a parent with a value larger than us. this parent
	   is the in-order successor of `node` */
	while (node->b_parent && node->b_parent->b_left != node) {
		node = node->b_parent;
	}
		
	return node->b_parent;
}	

btree_node_t *btree_prev(btree_node_t *node)
{
	if (!node) {
		return NULL;
	}

	/* there are two possibilities for the previous node:

	   1. if `node` has a left sub-tree, every node in this sub-tree is smaller
	      than `node`. the in-order predecessor of `node` is the largest node in
	      this subtree.
	   2. if `node` has no left sub-tree, we've reached the smallest node in
	      the sub-tree rooted at `node`. we need to go back to our parent
	      and continue the search elsewhere.
	*/
	if (node->b_left) {
		/* case 1: step into `node`'s left sub-tree and keep going
		   right to find the largest node */
		btree_node_t *cur = node->b_left;
		while (cur->b_right) {
			cur = cur->b_right;
		}

		return cur;
	}

	/* case 2: keep stepping back up towards the root of the tree.
	   if we encounter a step where we are our parent's right child,
	   we've found a parent with a value smaller than us. this parent
	   is the in-order predecessor of `node`. */
	while (node->b_parent && node->b_parent->b_right != node) {
		node = node->b_parent;
	}
		
	return node->b_parent;
}