开通博客赚积分
发布资源赚积分
分类

源码开发语言/平台
当前位置：首页 > 源码/资料 > 数值算法/人工智能 > 查看源码
planar_map.c
资源名称：leda.tar.gz　[点击查看]
上传用户：gzelex
上传日期：2007-01-07
资源大小：707k
文件大小：2k
源码类别：

数值算法/人工智能

开发平台：

MultiPlatform

planar_map.c：源码内容
  1. #include <LEDA/planar_map.h>
  2. #include <LEDA/graph_alg.h>
  4. void print_face(const PLANAR_MAP<int,int>& M, face f)
  5. { cout << string("F%d: ",M.inf(f));
  6.   list<node> L1 = M.adj_nodes(f);
  7.   node v;
  8.   forall(v,L1) M.print_node(v);
  9.   newline;
  10.  }
  13. main()
  15. {  GRAPH<int,int> G;
  16.    node v;
  17.    edge e;
  18.    face f;
  20.    test_graph(G);
  22.    edge_array<edge> rev(G);
  23.    Is_Bidirected(G,rev);
  25.    list<edge> L;
  26.    forall_edges(e,G) if (rev[e] == nil) L.append(e);
  28.    forall(e,L) G.new_edge(target(e),source(e));
  30.    if (!PLANAR(G,true))
  31.    { cerr << "G is not planar!n";
  32.      exit(1);
  33.     }
  35.    int count = 0;
  36.    forall_nodes(v,G) G[v] = count++;
  37.    count = 0;
  38.    forall_edges(e,G) G[e] = count++;
  41.    PLANAR_MAP<int,int>  M(G);
  43.    count = 0;
  44.    forall_faces(f,M) M[f] = count++;
  46.    newline;
  47.    forall_faces(f,M) print_face(M,f);
  49.    newline;
  50.    cout << "faces adjacent to nodes" << endl;
  52.    forall_nodes(v,M)
  53.    { M.print_node(v);
  54.      cout << " : ";
  55.      forall_adj_faces(f,v) cout << string("F%d ",M[f]);
  56.      newline;
  57.     }
  59. /*
  60.    M.triangulate();
  61.    forall_faces(f,M) print_face(M,f);
  62.    newline;
  63. */
  66.    newline;
  67.    while( Yes("split faces ? ") )
  68.    { list<face> FL= M.all_faces();
  69.       int i = 0;
  70.       forall(f,FL)  M.new_node(f,--i);
  71.    
  72.       forall_faces(f,M) print_face(M,f);
  73.       newline;
  74.    }
  75.    
  76.    
  77.    list<edge> E;
  78.    edge_array<bool> marked(M,false);
  80.    forall_edges(e,M)
  81.     if (!marked[e])
  82.     { E.append(e);
  83.       marked[M.reverse(e)] = true;
  84.      } 
  87.    forall(e,E)
  88.    { Yes("del_edge ");
  89.      M.del_edge(e);
  90.      forall_faces(f,M) print_face(M,f);
  91.      newline;
  93.   }
  95. }