线段树的一个好的实现，包括测度独立线段数的计算

#include  < stdio.h > 
#include  < stdlib.h > 
 
 const   int  maxn  =   5010 ;
 // 写一个线段树的过程 
 struct  Lines_tree
  ...{
    Lines_tree  *  lchild,  *  rchild;
     int  m;  // 测度 
      int  cnt;    // count 
      int  lines;  // 连续段数 
      int  lbd, rbd;  // 左右端点是否被覆盖  
      int  f, r;  // 左右端点 
 } ;
Lines_tree *  root;
 struct  rec  ...{ int  x, y, x1, y1;} r[maxn];
 struct  Line
  ...{
     int  x, y1, y2; int  sign;
    Line( int  a,  int  b,  int  c, int  d):x(a), y1(b), y2(c), sign(d)  ...{} 
     Line( void ):x( 0 ),y1( 0 ),y2( 0 ),sign( 0 )  ...{}  
} line[ 2 * maxn + 10 ];
 int  nr;
 int  ans;

 void  make_tree( int  a,  int  b, Lines_tree *  node)
  ...{
    node -> lines  =   0 ; node -> m  =   0 ; node -> cnt  =   0 ;
    node  ->  lbd  =   0 ; node  ->  rbd  =   0 ;
    node -> lchild  =  NULL; node -> rchild  =  NULL;
    node -> f  =  a; node -> r  =  b;
     if (b - a > 1 )
      ...{
        node -> lchild  =   new  Lines_tree;
        make_tree(a, (a + b) / 2 , node -> lchild);
        node -> rchild  =   new  Lines_tree;
        make_tree((a + b) / 2 , b, node -> rchild);
    } 
} 
 
 void  make( int  a,  int  b)
  ...{
     root  =   new  Lines_tree;
     make_tree(a, b, root);
} 
 
 void  update(Lines_tree  *  now)    // lbd, rbd, m的计算都在这个里面！ 
   ...{
     if (now -> cnt > 0 ) now -> m  =  now -> r - now -> f;
     else   if (now -> r == now -> f + 1 ) now -> m  =   0 ;
     else  now -> m  =  now -> lchild -> m  +  now -> rchild -> m;
} 
 
 void  update2(Lines_tree *  now)
  ...{
     if (now -> cnt > 0 )   ...{ now -> lbd  =   1 ; now -> rbd  =   1 ; now -> lines  =   1 ; } 
      else   if (now -> f + 1 == now -> r)   ...{now -> lbd  =   0 ; now -> rbd  =   0 ; now -> lines  =   0 ;} 
     else 
       ...{
        now -> lbd  =  now -> lchild -> lbd; now -> rbd  =  now -> rchild -> rbd;
        now -> lines  =  now -> lchild -> lines + now -> rchild -> lines  -  now -> lchild -> rbd * now -> rchild -> lbd;
    } 
} 
 
 void  insert( int  a,  int  b, Lines_tree  *  now)
  ...{
     if (a <= now -> f  &&  b >= now -> r)
        now -> cnt ++ ;
     if (now -> r - now -> f > 1 )
      ...{
         if (a < (now -> f + now -> r) / 2 )    insert(a, b, now -> lchild);
         if (b > (now -> f + now -> r) / 2 )    insert(a, b, now -> rchild);
    } 
    update(now);
    update2(now);
} 
 
 void  del( int  a,  int  b, Lines_tree  *  now)
  ...{
     if (a <= now -> f  &&  b >= now -> f)
      ...{
         if (a == now -> f) now -> lbd  =   0 ;
         if (b == now -> r) now -> rbd  =   0 ;
        now -> cnt -- ;
    } 
     if (now -> r - now -> f > 1 )
      ...{
         if (a < (now -> f + now -> r) / 2 )    del(a, b, now -> lchild);
         if (b > (now -> f + now -> r) / 2 )    del(a, b, now -> rchild);
    } 
    update(now);
    update2(now);
} 
 
 int  cmp( const   void   *  a,  const   void   *  b)
  ...{
     return  ( * (Line * )a).x  -  ( * (Line * )b).x;    // 这里不要写成-> 
 } 
 
 void  init()
  ...{
     // initiation
     // input 
      int  i;
    scanf( " %d " ,  & nr);
     for (i = 0 ; i < nr; i ++ )
      ...{
        scanf( " %d%d%d%d " ,  & r[i].x,  & r[i].y,  & r[i].x1,  & r[i].y1);
        line[ 2 * i]  =  Line(r[i].x, r[i].y, r[i].y1,  0 );
        line[ 2 * i + 1 ]  =  Line(r[i].x1, r[i].y, r[i].y1,  1 );
    }         
    qsort(line, nr * 2 ,  sizeof (line[ 0 ]), cmp);
     // pretreatment 
 } 
 
 void  work()
  ...{
     int  nowM  =   0 ;
     int  nowLine  =   0 ;
     int  lastM  =   0 ;
     int  lastLine  =   0 ;
     int  i;
     for (i = 0 ; i < nr * 2 ; i ++ )
      ...{
         if (line[i].sign == 0 )
            insert(line[i].y1, line[i].y2, root);
         else  del(line[i].y1, line[i].y2, root);
        nowM  =  root -> m;
        nowLine  =  root -> lines;
        ans  +=  lastLine  *   2   *  (line[i].x - line[i - 1 ].x);
        ans  +=  abs(nowM - lastM);
        lastM  =  nowM;
        lastLine  =  nowLine;
    } 
} 
 
 void  output()
  ...{
    printf( " %d  " , ans);
} 
 
 int  main()
  ...{
 //     freopen("t.in", "r", stdin); 
     make( - 10000 ,  10000 );
    init();
    work();
    output();
     return   0 ;
}