Program to implement a weighted sum of intervals tree.


#include <stdio.h>
#include <stdlib.h>

#define NEGINFTY -1000
#define BLOCKSIZE 256

typedef int object_t;
typedef int number_t;
typedef int key_t;
typedef struct tr_n_t { key_t        key;
                    struct tr_n_t  *left;
                    struct tr_n_t *right;
                    number_t      summand;
                    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;
}


tree_node_t *create_tree(void)
{  tree_node_t *tree;
   tree = get_node();
   tree->left = NULL;
   tree->summand = 0;
   /* need key -infty, use tree as non-NULL object ptr */
   insert( tree, NEGINFTY, (object_t *) tree );
   return( tree );
}

void left_rotation(tree_node_t *n)
{  tree_node_t *tmp_node;
   key_t        tmp_key; /* push down summand from n */
   n->left->summand += n->summand; 
   n->right->summand += n->summand;
   n->summand = 0;
   tmp_node = n->right; /* push down summand from n->right */
   tmp_node->left->summand += tmp_node->summand;
   tmp_node->right->summand += tmp_node->summand;
   tmp_node->summand = 0;
   tmp_node = n->left;  /* perform normal left rotation */
   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; /* push down summand from n */
   n->left->summand += n->summand;
   n->right->summand += n->summand;
   n->summand = 0;
   tmp_node = n->left; /* push down summand from n->left */
   tmp_node->left->summand += tmp_node->summand;
   tmp_node->right->summand += tmp_node->summand;
   tmp_node->summand = 0;
   tmp_node = n->right;/* perform normal right rotation */ 
   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 );
   }
}


number_t evaluate_sum(tree_node_t *tree, key_t query_key)
{  tree_node_t *current_tree_node;
   number_t sum;
   if( tree->left == NULL )
      return(0); 
   else
   {  current_tree_node = tree;
      sum = tree->summand;
      while( current_tree_node->right != NULL )
      {   if( query_key < current_tree_node->key )
              current_tree_node = current_tree_node->left;
          else
              current_tree_node = current_tree_node->right;
          sum += current_tree_node->summand;
      }
      return( sum );
   }
}


void insert_interval(tree_node_t *tree, 
                     key_t a, key_t b, number_t w)
{  tree_node_t *tmp_node;
   if ( find(tree, a) == NULL )
   {  insert( tree, a, (object_t *) tree );
   }  /* used treenode itself as non-NULL object pointer*/
   if ( find(tree, b) == NULL )
   {  insert( tree, b, (object_t *) tree );
   }
   tmp_node = tree; /* follow search path for a,*/
   while( tmp_node->right != NULL )
   {  /* add w to everything right of path */  
      if( a < tmp_node->key )
      {  tmp_node->right->summand += w;
         tmp_node = tmp_node->left;
      }
      else
         tmp_node = tmp_node->right;
   }
   tmp_node->summand += w; /* leaf with key a */
   tmp_node = tree; /* follow search path for b, */
   while( tmp_node->right != NULL )
   {  /* subtract w from everything right of path */  
      if( b < tmp_node->key )
      {  tmp_node->right->summand -= w;
         tmp_node = tmp_node->left;
      }
      else
         tmp_node = tmp_node->right;
   }
   tmp_node->summand -= w; /* leaf with key b */
}


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->summand = 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;
         old_leaf->summand = 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; 
         new_leaf->summand = 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 )
   {  printf("%d(%d)  ", tr->key, depth );
   }
   else
   {  check_tree(tr->left, depth+1, lower, tr->key ); 
      check_tree(tr->right, depth+1, tr->key, upper ); 
   }
}

int main()
{  tree_node_t *ws_tree;
   char nextop;
   ws_tree = create_tree();
   printf("Made Tree: weighted sum of intervals tree\n");
   while( (nextop = getchar())!= 'q' )
   { if( nextop == 'i' )
     { int ins_a, ins_b, ins_w, success;
       scanf(" [%d,%d[ %d", &ins_a, &ins_b, &ins_w);
       insert_interval( ws_tree, ins_a, ins_b, ins_w );
       printf("  inserted [%d,%d[with weight %d. Height now %d\n",
              ins_a, ins_b, ins_w, ws_tree->height );
     }  
     if( nextop == 'e' )
     { int eval_key, value;
       scanf(" %d", &eval_key);
       value = evaluate_sum( ws_tree, eval_key);
       printf(" the sum at key %d is %d\n", eval_key, value);
     }
     /*     if( nextop == 'd' )
     { int delkey, *delobj;
       scanf(" %d", &delkey);
       delobj = delete( searchtree, delkey);
       if( delobj == NULL )
         printf("  delete failed for key %d\n", delkey);
       else
         printf("  delete successful, deleted object %d, height is now %d\n", 
             *delobj, searchtree->height);
        }*/
     if( nextop == '?' )
     {  printf("  Checking tree\n"); 
        check_tree(ws_tree,0,-10000,10000);
        printf("\n");
        if( ws_tree->left != NULL )
      printf("key in root is %d, height of tree is %d\n", 
         ws_tree->key, ws_tree->height );
        printf("  Finished Checking tree\n"); 
     }
   }
   return(0);
}
naveed08st
 

Click Here to Leave a Comment Below 0 comments

Leave a Reply: