#include #include #define BLOCKSIZE 256 typedef int item_t; typedef int key_t; typedef struct uf_n_t { int rank; item_t *item; struct uf_n_t *up; } uf_node_t; typedef uf_node_t object_t; typedef struct tr_n_t { key_t key; struct tr_n_t *left; struct tr_n_t *right; int height; } tree_node_t; tree_node_t *currentblock = NULL; int size_left; tree_node_t *free_list = NULL; tree_node_t *get_node() { tree_node_t *tmp; if( free_list != NULL ) { tmp = free_list; free_list = free_list -> left; } else { if( currentblock == NULL || size_left == 0) { currentblock = (tree_node_t *) malloc( BLOCKSIZE * sizeof(tree_node_t) ); size_left = BLOCKSIZE; } tmp = currentblock++; size_left -= 1; } return( tmp ); } void return_node(tree_node_t *node) { node->left = free_list; free_list = node; } uf_node_t *get_uf_node() { return( (uf_node_t *) get_node() ); } void return_uf_node(uf_node_t *node) { return_node( (tree_node_t *) node ); } uf_node_t *uf_insert(item_t *new_item) { uf_node_t *new_node; new_node = get_uf_node(); new_node->item = new_item; new_node->rank = 0; new_node->up = NULL; return( new_node ); } int same_class( uf_node_t *node1, uf_node_t *node2 ) { uf_node_t *root1, *root2, *tmp; /* find both roots */ for( root1 = node1; root1->up != NULL; root1 = root1->up) ; /* follow path to root for node1 */ for( root2 = node2; root2->up != NULL; root2 = root2->up) ; /* follow path to root for node2 */ /* make both paths point directly to their respective roots */ tmp = node1->up; while( tmp != root1 && tmp != NULL ) { node1->up = root1; node1 = tmp; tmp = node1->up; } tmp = node2->up; while( tmp != root2 && tmp != NULL ) { node2->up = root2; node2 = tmp; tmp = node2->up; } /* return result */ return( root1 == root2 ); } void join( uf_node_t *node1, uf_node_t *node2 ) { uf_node_t *root1, *root2, *new_root, *tmp; /* find both roots */ for( root1 = node1; root1->up != NULL; root1 = root1->up) ; /* follow path to root for node1 */ for( root2 = node2; root2->up != NULL; root2 = root2->up) ; /* follow path to root for node2 */ /* perform union by rank */ if( root1->rank > root2->rank ) { new_root = root1; root2->up = new_root; } else if( root1->rank < root2->rank ) { new_root = root2; root1->up = new_root; } else /* same rank */ { new_root = root1; root2->up = new_root; new_root->rank += 1; } /* make both paths point directly to the new root */ tmp = node1->up; while( tmp != new_root && tmp != NULL ) { node1->up = new_root; node1 = tmp; tmp = node1->up; } tmp = node2->up; while( tmp != new_root && tmp != NULL ) { node2->up = new_root; node2 = tmp; tmp = node2->up; } } tree_node_t *create_tree(void) { tree_node_t *tmp_node; tmp_node = get_node(); tmp_node->left = NULL; return( tmp_node ); } void left_rotation(tree_node_t *n) { tree_node_t *tmp_node; key_t tmp_key; tmp_node = n->left; tmp_key = n->key; n->left = n->right; n->key = n->right->key; n->right = n->left->right; n->left->right = n->left->left; n->left->left = tmp_node; n->left->key = tmp_key; } void right_rotation(tree_node_t *n) { tree_node_t *tmp_node; key_t tmp_key; tmp_node = n->right; tmp_key = n->key; n->right = n->left; n->key = n->left->key; n->left = n->right->left; n->right->left = n->right->right; n->right->right = tmp_node; n->right->key = tmp_key; } object_t *find(tree_node_t *tree, key_t query_key) { tree_node_t *tmp_node; if( tree->left == NULL ) return(NULL); else { tmp_node = tree; while( tmp_node->right != NULL ) { if( query_key < tmp_node->key ) tmp_node = tmp_node->left; else tmp_node = tmp_node->right; } if( tmp_node->key == query_key ) return( (object_t *) tmp_node->left ); else return( NULL ); } } int insert(tree_node_t *tree, key_t new_key, object_t *new_object) { tree_node_t *tmp_node; int finished; if( tree->left == NULL ) { tree->left = (tree_node_t *) new_object; tree->key = new_key; tree->height = 0; tree->right = NULL; } else { tree_node_t * path_stack[100]; int path_st_p = 0; tmp_node = tree; while( tmp_node->right != NULL ) { path_stack[path_st_p++] = tmp_node; if( new_key < tmp_node->key ) tmp_node = tmp_node->left; else tmp_node = tmp_node->right; } /* found the candidate leaf. Test whether key distinct */ if( tmp_node->key == new_key ) return( -1 ); /* key is distinct, now perform the insert */ { tree_node_t *old_leaf, *new_leaf; old_leaf = get_node(); old_leaf->left = tmp_node->left; old_leaf->key = tmp_node->key; old_leaf->right = NULL; old_leaf->height = 0; new_leaf = get_node(); new_leaf->left = (tree_node_t *) new_object; new_leaf->key = new_key; new_leaf->right = NULL; new_leaf->height = 0; if( tmp_node->key < new_key ) { tmp_node->left = old_leaf; tmp_node->right = new_leaf; tmp_node->key = new_key; } else { tmp_node->left = new_leaf; tmp_node->right = old_leaf; } tmp_node->height = 1; } /* rebalance */ finished = 0; while( path_st_p > 0 && !finished ) { int tmp_height, old_height; tmp_node = path_stack[--path_st_p]; old_height= tmp_node->height; if( tmp_node->left->height - tmp_node->right->height == 2 ) { if( tmp_node->left->left->height - tmp_node->right->height == 1 ) { right_rotation( tmp_node ); tmp_node->right->height = tmp_node->right->left->height + 1; tmp_node->height = tmp_node->right->height + 1; } else { left_rotation( tmp_node->left ); right_rotation( tmp_node ); tmp_height = tmp_node->left->left->height; tmp_node->left->height = tmp_height + 1; tmp_node->right->height = tmp_height + 1; tmp_node->height = tmp_height + 2; } } else if ( tmp_node->left->height - tmp_node->right->height == -2 ) { if( tmp_node->right->right->height - tmp_node->left->height == 1 ) { left_rotation( tmp_node ); tmp_node->left->height = tmp_node->left->right->height + 1; tmp_node->height = tmp_node->left->height + 1; } else { right_rotation( tmp_node->right ); left_rotation( tmp_node ); tmp_height = tmp_node->right->right->height; tmp_node->left->height = tmp_height + 1; tmp_node->right->height = tmp_height + 1; tmp_node->height = tmp_height + 2; } } else /* update height even if there was no rotation */ { if( tmp_node->left->height > tmp_node->right->height ) tmp_node->height = tmp_node->left->height + 1; else tmp_node->height = tmp_node->right->height + 1; } if( tmp_node->height == old_height ) finished = 1; } } return( 0 ); } object_t *delete(tree_node_t *tree, key_t delete_key) { tree_node_t *tmp_node, *upper_node, *other_node; object_t *deleted_object; int finished; if( tree->left == NULL ) return( NULL ); else if( tree->right == NULL ) { if( tree->key == delete_key ) { deleted_object = (object_t *) tree->left; tree->left = NULL; return( deleted_object ); } else return( NULL ); } else { tree_node_t * path_stack[100]; int path_st_p = 0; tmp_node = tree; while( tmp_node->right != NULL ) { path_stack[path_st_p++] = tmp_node; upper_node = tmp_node; if( delete_key < tmp_node->key ) { tmp_node = upper_node->left; other_node = upper_node->right; } else { tmp_node = upper_node->right; other_node = upper_node->left; } } if( tmp_node->key != delete_key ) deleted_object = NULL; else { upper_node->key = other_node->key; upper_node->left = other_node->left; upper_node->right = other_node->right; upper_node->height = other_node->height; deleted_object = (object_t *) tmp_node->left; return_node( tmp_node ); return_node( other_node ); } /*start rebalance*/ finished = 0; path_st_p -= 1; while( path_st_p > 0 && !finished ) { int tmp_height, old_height; tmp_node = path_stack[--path_st_p]; old_height= tmp_node->height; if( tmp_node->left->height - tmp_node->right->height == 2 ) { if( tmp_node->left->left->height - tmp_node->right->height == 1 ) { right_rotation( tmp_node ); tmp_node->right->height = tmp_node->right->left->height + 1; tmp_node->height = tmp_node->right->height + 1; } else { left_rotation( tmp_node->left ); right_rotation( tmp_node ); tmp_height = tmp_node->left->left->height; tmp_node->left->height = tmp_height + 1; tmp_node->right->height = tmp_height + 1; tmp_node->height = tmp_height + 2; } } else if ( tmp_node->left->height - tmp_node->right->height == -2 ) { if( tmp_node->right->right->height - tmp_node->left->height == 1 ) { left_rotation( tmp_node ); tmp_node->left->height = tmp_node->left->right->height + 1; tmp_node->height = tmp_node->left->height + 1; } else { right_rotation( tmp_node->right ); left_rotation( tmp_node ); tmp_height = tmp_node->right->right->height; tmp_node->left->height = tmp_height + 1; tmp_node->right->height = tmp_height + 1; tmp_node->height = tmp_height + 2; } } else /* update height even if there was no rotation */ { if( tmp_node->left->height > tmp_node->right->height ) tmp_node->height = tmp_node->left->height + 1; else tmp_node->height = tmp_node->right->height + 1; } if( tmp_node->height == old_height ) finished = 1; } /*end rebalance*/ return( deleted_object ); } } void check_tree( tree_node_t *tr, int depth, int lower, int upper ) { if( tr->left == NULL ) { printf("Tree Empty\n"); return; } if( tr->key < lower || tr->key >= upper ) printf("Wrong Key Order \n"); if( tr->right == NULL ) { if( *( (int *) tr->left) == 10*tr->key + 2 ) printf("%d(%d) ", tr->key, depth ); else printf("Wrong Object \n"); } else { check_tree(tr->left, depth+1, lower, tr->key ); check_tree(tr->right, depth+1, tr->key, upper ); } } int main() { tree_node_t *searchtree; char nextop; searchtree = create_tree(); printf("Made Tree: Height-Balanced Tree\n"); while( (nextop = getchar())!= 'q' ) { if( nextop == 'i' ) { int *new_item; uf_node_t *new_uf_node; int success; new_item = (int *) malloc(sizeof(int)); scanf(" %d", new_item); new_uf_node = uf_insert(new_item); success = insert( searchtree, *new_item, new_uf_node ); if ( success == 0 ) printf(" inserted item %d successfully\n", *new_item); else printf(" insert failed, success = %d\n", success); } if( nextop == 's' ) { int item1, item2; uf_node_t *n1, *n2; scanf(" %d %d", &item1, &item2); n1 = find( searchtree, item1); n2 = find( searchtree, item2); if( same_class( n1, n2) ) printf(" items %d and %d are in the same class\n", item1, item2); else printf(" items %d and %d are in distinct classes\n", item1, item2); } if( nextop == 'j' ) { int item1, item2; uf_node_t *n1, *n2; scanf(" %d %d", &item1, &item2); n1 = find( searchtree, item1); n2 = find( searchtree, item2); join( n1, n2 ); if( same_class( n1, n2) ) printf(" joined classes of items %d and %d\n", item1, item2); else printf(" join of classes %d and %d failed\n", item1, item2); } if( nextop == '?' ) { printf(" Checking tree\n"); check_tree(searchtree,0,-1000,1000); printf("\n"); if( searchtree->left != NULL ) printf("key in root is %d, height of tree is %d\n", searchtree->key, searchtree->height ); printf(" Finished Checking tree\n"); } } return(0); }