8 programs. What kind of tasks does the library LEMON help to solve? |
8 programs. What kind of tasks does the library LEMON help to solve? |
9 It helps to write programs that solve optimization problems that arise |
9 It helps to write programs that solve optimization problems that arise |
10 frequently when <b>designing and testing certain networks</b>, for example |
10 frequently when <b>designing and testing certain networks</b>, for example |
11 in telecommunication, computer networks, and other areas that I cannot |
11 in telecommunication, computer networks, and other areas that I cannot |
12 think of now. A very natural way of modelling these networks is by means |
12 think of now. A very natural way of modelling these networks is by means |
13 of a <b> graph</b> (we will always mean a directed graph by that). |
13 of a <b> graph</b> (we will always mean a directed graph by that and say |
|
14 <b> undirected graph </b> otherwise). |
14 So if you want to write a program that works with |
15 So if you want to write a program that works with |
15 graphs then you might find it useful to use our library LEMON. |
16 graphs then you might find it useful to use our library LEMON. LEMON |
|
17 defines various graph concepts depending on what you want to do with the |
|
18 graph: a very good description can be found in the page |
|
19 about \ref graphs "graphs". |
16 |
20 |
17 |
21 You will also want to assign data to the edges or nodes of the graph, for example a length or capacity function defined on the edges. You can do this in LEMON using so called \ref maps "maps". You can define a map on the nodes or on the edges of the graph and the value of the map (the range of the function) can be practically almost any type. Read more about maps \ref maps-page "here". |
18 |
22 |
19 Some examples are the following (you will find links next to the code fragments that help to download full demo programs): |
23 Some examples are the following (you will find links next to the code fragments that help to download full demo programs): |
20 |
24 |
21 - First we give two examples that show how to instantiate a graph. The |
25 - First we give two examples that show how to instantiate a graph. The |
22 first one shows the methods that add nodes and edges, but one will |
26 first one shows the methods that add nodes and edges, but one will |
42 \endcode |
46 \endcode |
43 |
47 |
44 If you want to read more on the LEMON graph structures and concepts, read the page about \ref graphs "graphs". |
48 If you want to read more on the LEMON graph structures and concepts, read the page about \ref graphs "graphs". |
45 |
49 |
46 -# The following code shows how to read a graph from a stream (e.g. a file). LEMON supports the DIMACS file format: it can read a graph instance from a file |
50 -# The following code shows how to read a graph from a stream (e.g. a file). LEMON supports the DIMACS file format: it can read a graph instance from a file |
47 in that format (find the documentation of the format on the web). |
51 in that format (find the documentation of the DIMECS file format on the web). |
48 \code |
52 \code |
49 Graph g; |
53 Graph g; |
50 std::ifstream f("graph.dim"); |
54 std::ifstream f("graph.dim"); |
51 readDimacs(f, g); |
55 readDimacs(f, g); |
52 \endcode |
56 \endcode |
53 One can also store network (graph+capacity on the edges) instances and other things in DIMACS format: to see the details read the documentation of the \ref dimacs.h "Dimacs file format reader". |
57 One can also store network (graph+capacity on the edges) instances and other things in DIMACS format and use these in LEMON: to see the details read the documentation of the \ref dimacs.h "Dimacs file format reader". |
54 |
58 |
55 |
59 |
56 - If you want to solve some transportation problems in a network then |
60 - If you want to solve some transportation problems in a network then |
57 you will want to find shortest paths between nodes of a graph. This is |
61 you will want to find shortest paths between nodes of a graph. This is |
58 usually solved using Dijkstra's algorithm. A utility |
62 usually solved using Dijkstra's algorithm. A utility |
59 that solves this is the \ref lemon::Dijkstra "LEMON Dijkstra class". |
63 that solves this is the \ref lemon::Dijkstra "LEMON Dijkstra class". |
60 A simple program using the \ref lemon::Dijkstra "LEMON Dijkstra class" is |
64 The following code is a simple program using the \ref lemon::Dijkstra "LEMON |
61 as follows (we do not include the part that instantiates the graph and the length function): |
65 Dijkstra class" and it also shows how to define a map on the edges (the length |
|
66 function): |
62 |
67 |
63 \code |
68 \code |
64 typedef Graph::EdgeMap<int> LengthMap; |
69 |
65 Graph G; |
70 typedef ListGraph Graph; |
66 Node s, t; |
71 typedef Graph::Node Node; |
67 LengthMap cap(G); |
72 typedef Graph::Edge Edge; |
68 ... |
73 typedef Graph::EdgeMap<int> LengthMap; |
69 Dijkstra<Graph, LengthMap> |
74 |
70 dijkstra_test(G, cap); |
75 Graph g; |
71 dijkstra_test.run(s); |
76 |
|
77 //An example from Ahuja's book |
|
78 |
|
79 Node s=g.addNode(); |
|
80 Node v2=g.addNode(); |
|
81 Node v3=g.addNode(); |
|
82 Node v4=g.addNode(); |
|
83 Node v5=g.addNode(); |
|
84 Node t=g.addNode(); |
|
85 |
|
86 Edge s_v2=g.addEdge(s, v2); |
|
87 Edge s_v3=g.addEdge(s, v3); |
|
88 Edge v2_v4=g.addEdge(v2, v4); |
|
89 Edge v2_v5=g.addEdge(v2, v5); |
|
90 Edge v3_v5=g.addEdge(v3, v5); |
|
91 Edge v4_t=g.addEdge(v4, t); |
|
92 Edge v5_t=g.addEdge(v5, t); |
|
93 |
|
94 LengthMap len(g); |
|
95 |
|
96 len.set(s_v2, 10); |
|
97 len.set(s_v3, 10); |
|
98 len.set(v2_v4, 5); |
|
99 len.set(v2_v5, 8); |
|
100 len.set(v3_v5, 5); |
|
101 len.set(v4_t, 8); |
|
102 len.set(v5_t, 8); |
|
103 |
|
104 std::cout << "The id of s is " << g.id(s)<< ", the id of t is " << g.id(t)<<"."<<std::endl; |
|
105 |
|
106 std::cout << "Dijkstra algorithm test..." << std::endl; |
|
107 |
|
108 Dijkstra<Graph, LengthMap> dijkstra_test(g,len); |
|
109 |
|
110 dijkstra_test.run(s); |
|
111 |
|
112 |
|
113 std::cout << "The distance of node t from node s: " << dijkstra_test.dist(t)<<std::endl; |
|
114 |
|
115 std::cout << "The shortest path from s to t goes through the following nodes (the first one is t, the last one is s): "<<std::endl; |
|
116 |
|
117 for (Node v=t;v != s; v=dijkstra_test.predNode(v)){ |
|
118 std::cout << g.id(v) << "<-"; |
|
119 } |
|
120 std::cout << g.id(s) << std::endl; |
72 \endcode |
121 \endcode |
|
122 |
|
123 See the whole program in \file dijkstra_demo.cc. |
|
124 |
|
125 The first part of the code is self-explanatory: we build the graph and set the |
|
126 length values of the edges. Then we instantiate a member of the Dijkstra class |
|
127 and run the Dijkstra algorithm from node \c s. After this we read some of the |
|
128 results. |
|
129 You can do much more with the Dijkstra class, for example you can run it step |
|
130 by step and gain full control of the execution. For a detailed description, see the documentation of the \ref lemon::Dijkstra "LEMON Dijkstra class". |
|
131 |
73 |
132 |
74 - If you want to design a network and want to minimize the total length |
133 - If you want to design a network and want to minimize the total length |
75 of wires then you might be looking for a <b>minimum spanning tree</b> in |
134 of wires then you might be looking for a <b>minimum spanning tree</b> in |
76 an undirected graph. This can be found using the Kruskal algorithm: the |
135 an undirected graph. This can be found using the Kruskal algorithm: the |
77 class \ref lemon::Kruskal "LEMON Kruskal class" does this job for you. |
136 class \ref lemon::Kruskal "LEMON Kruskal class" does this job for you. |