doc/quicktour.dox
author deba
Sat, 17 Nov 2007 21:41:01 +0000
changeset 2515 caa640aa9a7e
parent 1640 9c7834ac5e64
child 2553 bfced05fa852
permissions -rw-r--r--
Changing max flow interface
alpar@2391
     1
/* -*- C++ -*-
alpar@2391
     2
 *
alpar@2391
     3
 * This file is a part of LEMON, a generic C++ optimization library
alpar@2391
     4
 *
alpar@2391
     5
 * Copyright (C) 2003-2007
alpar@2391
     6
 * Egervary Jeno Kombinatorikus Optimalizalasi Kutatocsoport
alpar@2391
     7
 * (Egervary Research Group on Combinatorial Optimization, EGRES).
alpar@2391
     8
 *
alpar@2391
     9
 * Permission to use, modify and distribute this software is granted
alpar@2391
    10
 * provided that this copyright notice appears in all copies. For
alpar@2391
    11
 * precise terms see the accompanying LICENSE file.
alpar@2391
    12
 *
alpar@2391
    13
 * This software is provided "AS IS" with no warranty of any kind,
alpar@2391
    14
 * express or implied, and with no claim as to its suitability for any
alpar@2391
    15
 * purpose.
alpar@2391
    16
 *
alpar@2391
    17
 */
alpar@2391
    18
athos@1169
    19
/**
athos@1169
    20
alpar@1170
    21
\page quicktour Quick Tour to LEMON
alpar@1170
    22
athos@1580
    23
Let us first answer the question <b>"What do I want to use LEMON for?"</b>. 
athos@1175
    24
LEMON is a C++ library, so you can use it if you want to write C++ 
athos@1175
    25
programs. What kind of tasks does the library LEMON help to solve? 
athos@1175
    26
It helps to write programs that solve optimization problems that arise
athos@1175
    27
frequently when <b>designing and testing certain networks</b>, for example
athos@1175
    28
in telecommunication, computer networks, and other areas that I cannot
athos@1175
    29
think of now. A very natural way of modelling these networks is by means
athos@1183
    30
of a <b> graph</b> (we will always mean a directed graph by that and say
athos@1183
    31
<b> undirected graph </b> otherwise). 
athos@1175
    32
So if you want to write a program that works with 
athos@1183
    33
graphs then you might find it useful to use our library LEMON. LEMON 
athos@1183
    34
defines various graph concepts depending on what you want to do with the 
athos@1183
    35
graph: a very good description can be found in the page
athos@1183
    36
about \ref graphs "graphs".
athos@1175
    37
athos@1514
    38
You will also want to assign data to the edges or nodes of the graph, for
athos@1514
    39
example a length or capacity function defined on the edges. You can do this in
athos@1514
    40
LEMON using so called \b 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 of any type. Read more about maps \ref maps-page "here".
athos@1175
    41
athos@1528
    42
In this quick tour we want to show you some facilities LEMON library can provide through examples (simple demo programs). The examples will only show part of the functionality, but links will always be given to reach complete details. 
athos@1528
    43
You will find links next to the code fragments that help to download full demo programs: save them on your computer and compile them according to the description in the page about \ref getstart "How to start using LEMON". 
athos@1528
    44
athos@1528
    45
Have fun!
athos@1175
    46
athos@1522
    47
<ul> <li> The first thing to discuss is the way one can create data structures
athos@1522
    48
like graphs and maps in a program using LEMON. 
athos@1522
    49
//There are more graph types
athos@1522
    50
//implemented in LEMON and you can implement your own graph type just as well:
athos@1522
    51
//read more about this in the already mentioned page on \ref graphs "graphs".
athos@1522
    52
athos@1522
    53
First we show how to add nodes and edges to a graph manually. We will also
athos@1522
    54
define a map on the edges of the graph. After this we show the way one can
athos@1522
    55
read a graph (and perhaps maps on it) from a stream (e.g. a file). Of course
athos@1522
    56
we also have routines that write a graph (and perhaps maps) to a stream
athos@1522
    57
(file): this will also be shown. LEMON supports the DIMACS file formats to
athos@1534
    58
read network optimization problems, but more importantly we also have our own
athos@1522
    59
file format that gives a more flexible way to store data related to network
athos@1522
    60
optimization.
athos@1522
    61
athos@1530
    62
<ol> <li>The following code shows how to build a graph from scratch
athos@1530
    63
and iterate on its nodes and edges.  This example also shows how to
athos@1530
    64
give a map on the edges of the graph.  The type Listgraph is one of
athos@1530
    65
the LEMON graph types: the typedefs in the beginning are for
athos@1530
    66
convenience and we will assume them later as well.
athos@1522
    67
alpar@1640
    68
\dontinclude hello_lemon.cc
alpar@1640
    69
\skip include
alpar@1640
    70
\until }
athos@1522
    71
athos@1530
    72
See the whole program in file \ref hello_lemon.cc in the \c demo subdir of
athos@1526
    73
LEMON package.
athos@1175
    74
athos@1526
    75
    If you want to read more on the LEMON graph structures and
athos@1526
    76
concepts, read the page about \ref graphs "graphs".
athos@1522
    77
athos@1530
    78
athos@1530
    79
<li>LEMON has an own file format for storing graphs, maps on edges/nodes and some other things. Instead of any explanation let us give a
athos@1530
    80
short example file in this format: read the detailed description of the LEMON
athos@1530
    81
graph file format and input-output routines here: \ref graph-io-page.
athos@1530
    82
athos@1530
    83
So here is a file describing a graph of 6 nodes (0 to 5), two nodemaps
athos@1530
    84
(called \c coordinates_x and \c coordinates_y), several edges, an edge map
athos@1530
    85
called \c capacity and two designated nodes (called \c source and \c target).
athos@1530
    86
athos@1541
    87
\verbatim
athos@1541
    88
@nodeset
athos@1541
    89
id      coordinates_x   coordinates_y
athos@1541
    90
5       796.398 208.035
athos@1541
    91
4       573.002 63.002
athos@1541
    92
3       568.549 401.748
athos@1541
    93
2       277.889 68.476
athos@1541
    94
1       288.248 397.327
athos@1541
    95
0       102.239 257.532
athos@1541
    96
@edgeset
athos@1541
    97
                id      capacity
athos@1541
    98
4       5       6       8
athos@1541
    99
3       5       5       8
athos@1541
   100
2       4       4       5
athos@1541
   101
1       4       3       8
athos@1541
   102
1       3       2       5
athos@1541
   103
0       2       1       10
athos@1541
   104
0       1       0       10
athos@1541
   105
#This is a comment here
athos@1541
   106
@nodes
athos@1541
   107
source 0
athos@1541
   108
target 5
athos@1541
   109
@edges 
athos@1541
   110
@attributes 
athos@1541
   111
author "Attila BERNATH"
athos@1541
   112
@end
athos@1541
   113
\endverbatim
athos@1530
   114
athos@1530
   115
Finally let us give a simple example that reads a graph from a file and writes
athos@1530
   116
it to the standard output.
athos@1530
   117
alpar@1640
   118
\dontinclude reader_writer_demo.cc
alpar@1640
   119
\skip include
alpar@1640
   120
\until return
alpar@1640
   121
\until }
athos@1530
   122
athos@1530
   123
See the whole program in file \ref reader_writer_demo.cc.
athos@1530
   124
athos@1526
   125
<li> The following code shows how to read a graph from a stream
athos@1526
   126
(e.g. a file) in the DIMACS file format (find the documentation of the
athos@1526
   127
DIMACS file formats on the web).
athos@1522
   128
athos@1181
   129
\code
athos@1181
   130
Graph g;
athos@1181
   131
std::ifstream f("graph.dim");
athos@1181
   132
readDimacs(f, g);
athos@1181
   133
\endcode
athos@1522
   134
athos@1526
   135
One can also store network (graph+capacity on the edges) instances and
athos@1526
   136
other things (minimum cost flow instances etc.) in DIMACS format and
athos@1534
   137
read these in LEMON: to see the details read the documentation of the
athos@1534
   138
\ref dimacs.h "Dimacs file format reader". 
athos@1522
   139
athos@1514
   140
</ol>
athos@1514
   141
<li> If you want to solve some transportation problems in a network then 
athos@1175
   142
you will want to find shortest paths between nodes of a graph. This is 
athos@1175
   143
usually solved using Dijkstra's algorithm. A utility
athos@1175
   144
that solves this is  the \ref lemon::Dijkstra "LEMON Dijkstra class".
athos@1522
   145
The following code is a simple program using the 
athos@1530
   146
\ref lemon::Dijkstra "LEMON Dijkstra class": it calculates the shortest path between node \c s and \c t in a graph \c g.
athos@1530
   147
We omit the part reading the graph  \c g and the length map \c len.
athos@1175
   148
athos@1528
   149
\dontinclude dijkstra_demo.cc
athos@1528
   150
\skip ListGraph
athos@1530
   151
\until Graph g
athos@1530
   152
...
athos@1530
   153
\skip Dijkstra algorithm
athos@1528
   154
\until std::cout << g.id(s)
athos@1175
   155
alpar@1287
   156
See the whole program in \ref dijkstra_demo.cc.
athos@1183
   157
athos@1530
   158
Some explanation: after instantiating a member of the Dijkstra class
athos@1530
   159
we run the Dijkstra algorithm from node \c s. After this we read some
athos@1530
   160
of the results.  You can do much more with the Dijkstra class, for
athos@1530
   161
example you can run it step by step and gain full control of the
athos@1530
   162
execution. For a detailed description, see the documentation of the
athos@1530
   163
\ref lemon::Dijkstra "LEMON Dijkstra class".
athos@1183
   164
athos@1183
   165
zsuzska@1578
   166
<li> If you want to design a network and want to minimize the total
zsuzska@1578
   167
length of wires then you might be looking for a <b>minimum spanning
zsuzska@1578
   168
tree</b> in an undirected graph. This can be found using the Kruskal
zsuzska@1584
   169
algorithm: the function \ref lemon::kruskal "LEMON Kruskal " does this
zsuzska@1584
   170
job for you.  
zsuzska@1584
   171
zsuzska@1584
   172
First make a graph \c g and a cost map \c
zsuzska@1584
   173
edge_cost_map, then make a bool edgemap \c tree_map or a vector \c
zsuzska@1584
   174
tree_edge_vec for the algorithm output. After calling the function it
zsuzska@1584
   175
gives back the weight of the minimum spanning tree and the \c tree_map or
zsuzska@1584
   176
the \c tree_edge_vec contains the edges of the tree.
zsuzska@1584
   177
zsuzska@1584
   178
If you want to store the edges in a bool edgemap, then use the
zsuzska@1584
   179
function as follows:
athos@1175
   180
zsuzska@1578
   181
\dontinclude kruskal_demo.cc
zsuzska@1584
   182
\skip Kruskal with boolmap; 
zsuzska@1584
   183
\until  std::endl
zsuzska@1578
   184
zsuzska@1584
   185
And if you rather use a vector instead of a bool map:
athos@1580
   186
zsuzska@1584
   187
\skip Kruskal with vector; 
zsuzska@1584
   188
\until std::endl
zsuzska@1578
   189
zsuzska@1578
   190
See the whole program in \ref kruskal_demo.cc.
zsuzska@1578
   191
zsuzska@1578
   192
athos@1511
   193
athos@1517
   194
<li>Many problems in network optimization can be formalized by means
athos@1517
   195
of a linear programming problem (LP problem, for short). In our
athos@1517
   196
library we decided not to write an LP solver, since such packages are
athos@1517
   197
available in the commercial world just as well as in the open source
athos@1517
   198
world, and it is also a difficult task to compete these. Instead we
athos@1517
   199
decided to develop an interface that makes it easier to use these
athos@1517
   200
solvers together with LEMON. The advantage of this approach is
athos@1517
   201
twofold. Firstly our C++ interface is more comfortable than the
athos@1517
   202
solvers' native interface. Secondly, changing the underlying solver in
athos@1517
   203
a certain software using LEMON's LP interface needs zero effort. So,
athos@1517
   204
for example, one may try his idea using a free solver, demonstrate its
athos@1517
   205
usability for a customer and if it works well, but the performance
athos@1517
   206
should be improved, then one may decide to purchase and use a better
athos@1517
   207
commercial solver.
athos@1517
   208
athos@1517
   209
So far we have an
athos@1526
   210
interface for the commercial LP solver software \b CPLEX (developed by ILOG)
athos@1514
   211
and for the open source solver \b GLPK (a shorthand for Gnu Linear Programming
athos@1517
   212
Toolkit).
athos@1514
   213
athos@1514
   214
We will show two examples, the first one shows how simple it is to formalize
athos@1514
   215
and solve an LP problem in LEMON, while the second one shows how LEMON
athos@1514
   216
facilitates solving network optimization problems using LP solvers.
athos@1514
   217
athos@1514
   218
<ol>
athos@1514
   219
<li>The following code shows how to solve an LP problem using the LEMON lp
athos@1517
   220
interface. The code together with the comments is self-explanatory.
athos@1511
   221
athos@1530
   222
\dontinclude lp_demo.cc
athos@1530
   223
\skip A default solver is taken
athos@1530
   224
\until End of LEMON style code
athos@1175
   225
athos@1514
   226
See the whole code in \ref lp_demo.cc.
athos@1514
   227
athos@1517
   228
<li>The second example shows how easy it is to formalize a max-flow
athos@1517
   229
problem as an LP problem using the LEMON LP interface: we are looking
athos@1517
   230
for a real valued function defined on the edges of the digraph
athos@1517
   231
satisfying the nonnegativity-, the capacity constraints and the
athos@1517
   232
flow-conservation constraints and giving the largest flow value
athos@1517
   233
between to designated nodes.
athos@1517
   234
athos@1517
   235
In the following code we suppose that we already have the graph \c g,
athos@1517
   236
the capacity map \c cap, the source node \c s and the target node \c t
athos@1517
   237
in the memory. We will also omit the typedefs.
athos@1517
   238
athos@1530
   239
\dontinclude lp_maxflow_demo.cc
athos@1530
   240
\skip Define a map on the edges for the variables of the LP problem
athos@1530
   241
\until lp.max();
athos@1530
   242
\skip Solve with the underlying solver
athos@1530
   243
\until lp.solve();
athos@1517
   244
athos@1517
   245
athos@1517
   246
The complete program can be found in file \ref lp_maxflow_demo.cc. After compiling run it in the form:
athos@1517
   247
athos@1528
   248
<tt>./lp_maxflow_demo < sample.lgf</tt>
athos@1517
   249
athos@1528
   250
where sample.lgf is a file in the lemon format containing a maxflow instance (designated "source", "target" nodes and "capacity" map on the edges).
athos@1517
   251
athos@1517
   252
athos@1514
   253
athos@1514
   254
</ol>
athos@1514
   255
</ul>
athos@1175
   256
athos@1175
   257
*/