ucx/avl.c

Tue, 19 May 2015 16:47:54 +0200

author
Mike Becker <universe@uap-core.de>
date
Tue, 19 May 2015 16:47:54 +0200
changeset 204
4477987d40cd
parent 203
3d999ea3f780
child 205
54a7ceb9151f
permissions
-rw-r--r--

better and better and better AVL API

universe@192 1 /*
universe@192 2 * DO NOT ALTER OR REMOVE COPYRIGHT NOTICES OR THIS HEADER.
universe@192 3 *
universe@192 4 * Copyright 2015 Olaf Wintermann. All rights reserved.
universe@192 5 *
universe@192 6 * Redistribution and use in source and binary forms, with or without
universe@192 7 * modification, are permitted provided that the following conditions are met:
universe@192 8 *
universe@192 9 * 1. Redistributions of source code must retain the above copyright
universe@192 10 * notice, this list of conditions and the following disclaimer.
universe@192 11 *
universe@192 12 * 2. Redistributions in binary form must reproduce the above copyright
universe@192 13 * notice, this list of conditions and the following disclaimer in the
universe@192 14 * documentation and/or other materials provided with the distribution.
universe@192 15 *
universe@192 16 * THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS IS"
universe@192 17 * AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
universe@192 18 * IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE
universe@192 19 * ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT HOLDER OR CONTRIBUTORS BE
universe@192 20 * LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR
universe@192 21 * CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF
universe@192 22 * SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS
universe@192 23 * INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN
universe@192 24 * CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE)
universe@192 25 * ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE
universe@192 26 * POSSIBILITY OF SUCH DAMAGE.
universe@192 27 */
universe@192 28
universe@192 29 #include "avl.h"
universe@192 30
universe@195 31 #define ptrcast(ptr) ((void*)(ptr))
universe@195 32
universe@195 33 static void ucx_avl_connect(UcxAVLTree *tree,
universe@195 34 UcxAVLNode *node, UcxAVLNode *child, intptr_t nullkey) {
universe@195 35 if (child) {
universe@195 36 child->parent = node;
universe@195 37 }
universe@195 38 // if child is NULL, nullkey decides if left or right pointer is cleared
universe@195 39 if (tree->cmpfunc(
universe@195 40 ptrcast(child ? child->key : nullkey),
universe@195 41 ptrcast(node->key), tree->userdata) > 0) {
universe@195 42 node->right = child;
universe@195 43 } else {
universe@195 44 node->left = child;
universe@195 45 }
universe@195 46 size_t lh = node->left ? node->left->height : 0;
universe@195 47 size_t rh = node->right ? node->right->height : 0;
universe@195 48 node->height = 1 + (lh > rh ? lh : rh);
universe@195 49 }
universe@195 50
universe@195 51 #define avlheight(node) ((node) ? (node)->height : 0)
universe@195 52
universe@195 53 static UcxAVLNode* avl_rotright(UcxAVLTree *tree, UcxAVLNode *l0) {
universe@195 54 UcxAVLNode *p = l0->parent;
universe@195 55 UcxAVLNode *l1 = l0->left;
universe@195 56 if (p) {
universe@195 57 ucx_avl_connect(tree, p, l1, 0);
universe@195 58 } else {
universe@195 59 l1->parent = NULL;
universe@195 60 }
universe@195 61 ucx_avl_connect(tree, l0, l1->right, l1->key);
universe@195 62 ucx_avl_connect(tree, l1, l0, 0);
universe@195 63 return l1;
universe@195 64 }
universe@195 65
universe@195 66 static UcxAVLNode* avl_rotleft(UcxAVLTree *tree, UcxAVLNode *l0) {
universe@195 67 UcxAVLNode *p = l0->parent;
universe@195 68 UcxAVLNode *l1 = l0->right;
universe@195 69 if (p) {
universe@195 70 ucx_avl_connect(tree, p, l1, 0);
universe@195 71 } else {
universe@195 72 l1->parent = NULL;
universe@195 73 }
universe@195 74 ucx_avl_connect(tree, l0, l1->left, l1->key);
universe@195 75 ucx_avl_connect(tree, l1, l0, 0);
universe@195 76 return l1;
universe@195 77 }
universe@195 78
universe@195 79 static void ucx_avl_balance(UcxAVLTree *tree, UcxAVLNode *n) {
universe@195 80 int lh = avlheight(n->left);
universe@195 81 int rh = avlheight(n->right);
universe@195 82 n->height = 1 + (lh > rh ? lh : rh);
universe@195 83
universe@195 84 if (lh - rh == 2) {
universe@195 85 UcxAVLNode *c = n->left;
universe@195 86 if (avlheight(c->right) - avlheight(c->left) == 1) {
universe@195 87 avl_rotleft(tree, c);
universe@195 88 }
universe@195 89 n = avl_rotright(tree, n);
universe@195 90 } else if (rh - lh == 2) {
universe@195 91 UcxAVLNode *c = n->right;
universe@195 92 if (avlheight(c->left) - avlheight(c->right) == 1) {
universe@195 93 avl_rotright(tree, c);
universe@195 94 }
universe@195 95 n = avl_rotleft(tree, n);
universe@195 96 }
universe@195 97
universe@195 98 if (n->parent) {
universe@195 99 ucx_avl_balance(tree, n->parent);
universe@195 100 } else {
universe@195 101 tree->root = n;
universe@195 102 }
universe@195 103 }
universe@195 104
universe@194 105 UcxAVLTree *ucx_avl_new(cmp_func cmpfunc) {
universe@194 106 return ucx_avl_new_a(cmpfunc, ucx_default_allocator());
universe@194 107 }
universe@194 108
universe@194 109 UcxAVLTree *ucx_avl_new_a(cmp_func cmpfunc, UcxAllocator *allocator) {
universe@196 110 UcxAVLTree *tree = almalloc(allocator, sizeof(UcxAVLTree));
universe@194 111 if (tree) {
universe@194 112 tree->allocator = allocator;
universe@194 113 tree->cmpfunc = cmpfunc;
universe@194 114 tree->root = NULL;
universe@194 115 tree->userdata = NULL;
universe@194 116 }
universe@194 117
universe@194 118 return tree;
universe@194 119 }
universe@194 120
universe@196 121 static void ucx_avl_free_node(UcxAllocator *al, UcxAVLNode *node) {
universe@196 122 if (node) {
universe@196 123 ucx_avl_free_node(al, node->left);
universe@196 124 ucx_avl_free_node(al, node->right);
universe@196 125 alfree(al, node);
universe@196 126 }
universe@196 127 }
universe@196 128
universe@196 129 void ucx_avl_free(UcxAVLTree *tree) {
universe@196 130 UcxAllocator *al = tree->allocator;
universe@196 131 ucx_avl_free_node(al, tree->root);
universe@196 132 alfree(al, tree);
universe@196 133 }
universe@196 134
universe@204 135 UcxAVLNode *ucx_avl_getn(UcxAVLTree *tree, intptr_t key) {
universe@195 136 UcxAVLNode *n = tree->root;
universe@195 137 int cmpresult;
universe@195 138 while (n && (cmpresult = tree->cmpfunc(
universe@195 139 ptrcast(key), ptrcast(n->key), tree->userdata))) {
universe@195 140 n = cmpresult > 0 ? n->right : n->left;
universe@195 141 }
universe@204 142 return n;
universe@204 143 }
universe@204 144
universe@204 145 void *ucx_avl_get(UcxAVLTree *tree, intptr_t key) {
universe@204 146 UcxAVLNode *n = ucx_avl_getn(tree, key);
universe@195 147 return n ? n->value : NULL;
universe@194 148 }
universe@194 149
universe@204 150 int ucx_avl_put(UcxAVLTree *tree, intptr_t key, void *value) {
universe@204 151 return ucx_avl_put_s(tree, key, value, NULL);
universe@204 152 }
universe@204 153
universe@204 154 int ucx_avl_put_s(UcxAVLTree *tree, intptr_t key, void *value,
universe@204 155 void **oldvalue) {
universe@195 156 if (tree->root) {
universe@195 157 UcxAVLNode *n = tree->root;
universe@195 158 int cmpresult;
universe@197 159 while ((cmpresult = tree->cmpfunc(
universe@197 160 ptrcast(key), ptrcast(n->key), tree->userdata))) {
universe@195 161 UcxAVLNode *m = cmpresult > 0 ? n->right : n->left;
universe@195 162 if (m) {
universe@195 163 n = m;
universe@195 164 } else {
universe@195 165 break;
universe@195 166 }
universe@195 167 }
universe@195 168
universe@195 169 if (cmpresult) {
universe@196 170 UcxAVLNode *e = almalloc(tree->allocator, sizeof(UcxAVLNode));
universe@204 171 if (e) {
universe@204 172 e->key = key; e->value = value; e->height = 1;
universe@204 173 e->parent = e->left = e->right = NULL;
universe@204 174 ucx_avl_connect(tree, n, e, 0);
universe@204 175 ucx_avl_balance(tree, n);
universe@204 176 return 0;
universe@204 177 } else {
universe@204 178 return 1;
universe@204 179 }
universe@195 180 } else {
universe@204 181 if (oldvalue) {
universe@204 182 *oldvalue = n->value;
universe@204 183 }
universe@195 184 n->value = value;
universe@204 185 return 0;
universe@195 186 }
universe@195 187 } else {
universe@196 188 tree->root = almalloc(tree->allocator, sizeof(UcxAVLNode));
universe@204 189 if (tree->root) {
universe@204 190 tree->root->key = key; tree->root->value = value;
universe@204 191 tree->root->height = 1;
universe@204 192 tree->root->parent = tree->root->left = tree->root->right = NULL;
universe@204 193
universe@204 194 if (oldvalue) {
universe@204 195 *oldvalue = NULL;
universe@204 196 }
universe@204 197
universe@204 198 return 0;
universe@204 199 } else {
universe@204 200 return 1;
universe@204 201 }
universe@195 202 }
universe@194 203 }
universe@194 204
universe@204 205 int ucx_avl_remove(UcxAVLTree *tree, intptr_t key) {
universe@204 206 return ucx_avl_remove_s(tree, key, NULL, NULL);
universe@204 207 }
universe@204 208
universe@204 209 int ucx_avl_remove_n(UcxAVLTree *tree, UcxAVLNode *node) {
universe@204 210 return ucx_avl_remove_s(tree, node->key, NULL, NULL);
universe@204 211 }
universe@204 212
universe@204 213 int ucx_avl_remove_s(UcxAVLTree *tree, intptr_t key,
universe@204 214 intptr_t *oldkey, void **oldvalue) {
universe@204 215
universe@195 216 UcxAVLNode *n = tree->root;
universe@195 217 int cmpresult;
universe@195 218 while (n && (cmpresult = tree->cmpfunc(
universe@195 219 ptrcast(key), ptrcast(n->key), tree->userdata))) {
universe@195 220 n = cmpresult > 0 ? n->right : n->left;
universe@195 221 }
universe@195 222 if (n) {
universe@204 223 if (oldkey) {
universe@204 224 *oldkey = n->key;
universe@204 225 }
universe@204 226 if (oldvalue) {
universe@204 227 *oldvalue = n->value;
universe@204 228 }
universe@204 229
universe@195 230 UcxAVLNode *p = n->parent;
universe@195 231 if (n->left && n->right) {
universe@195 232 UcxAVLNode *s = n->right;
universe@195 233 while (s->left) {
universe@195 234 s = s->left;
universe@195 235 }
universe@195 236 ucx_avl_connect(tree, s->parent, s->right, s->key);
universe@195 237 n->key = s->key; n->value = s->value;
universe@195 238 p = s->parent;
universe@196 239 alfree(tree->allocator, s);
universe@195 240 } else {
universe@195 241 if (p) {
universe@195 242 ucx_avl_connect(tree, p, n->right ? n->right:n->left, n->key);
universe@195 243 } else {
universe@195 244 tree->root = n->right ? n->right : n->left;
olaf@203 245 if (tree->root) {
olaf@202 246 tree->root->parent = NULL;
olaf@202 247 }
universe@195 248 }
universe@196 249 alfree(tree->allocator, n);
universe@195 250 }
universe@195 251
universe@195 252 if (p) {
universe@195 253 ucx_avl_balance(tree, p);
universe@195 254 }
universe@204 255
universe@204 256 return 0;
universe@195 257 } else {
universe@204 258 return 1;
universe@195 259 }
universe@194 260 }
universe@194 261
universe@199 262 static size_t ucx_avl_countn(UcxAVLNode *node) {
universe@199 263 if (node) {
universe@199 264 return 1 + ucx_avl_countn(node->left) + ucx_avl_countn(node->right);
universe@199 265 } else {
universe@199 266 return 0;
universe@199 267 }
universe@199 268 }
universe@199 269
universe@199 270 size_t ucx_avl_count(UcxAVLTree *tree) {
universe@199 271 return ucx_avl_countn(tree->root);
universe@199 272 }

mercurial