Location: LEMON/LEMON-official/lemon/concepts/heap.h - annotation
Load file history
Support negative costs and bounds in NetworkSimplex (#270)
* The interface is reworked to support negative costs and bounds.
- ProblemType and problemType() are renamed to
SupplyType and supplyType(), see also #234.
- ProblemType type is introduced similarly to the LP interface.
- 'bool run()' is replaced by 'ProblemType run()' to handle
unbounded problem instances, as well.
- Add INF public member constant similarly to the LP interface.
* Remove capacityMap() and boundMaps(), see also #266.
* Update the problem definition in the MCF module.
* Remove the usage of Circulation (and adaptors) for checking feasibility.
Check feasibility by examining the artifical arcs instead (after solving
the problem).
* Additional check for unbounded negative cycles found during the
algorithm (it is possible now, since negative costs are allowed).
* Fix in the constructor (the value types needn't be integer any more),
see also #254.
* Improve and extend the doc.
* Rework the test file and add test cases for negative costs and bounds.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136 137 138 139 140 141 142 143 144 145 146 147 148 149 150 151 152 153 154 155 156 157 158 159 160 161 162 163 164 165 166 167 168 169 170 171 172 173 174 175 176 177 178 179 180 181 182 183 184 185 186 187 188 189 190 191 192 193 194 195 196 197 198 199 200 201 202 203 204 205 206 207 208 209 210 211 212 213 214 215 216 217 218 219 220 221 222 223 224 225 226 227 228 229 230 231 232 233 234 235 236 237 238 239 240 241 242 243 244 245 246 247 248 249 250 251 252 253 254 255 256 257 258 259 260 261 | r209:765619b7cbb2 r100:4f754b4cf82b r209:765619b7cbb2 r100:4f754b4cf82b r463:88ed40ad0d4f r100:4f754b4cf82b r100:4f754b4cf82b r100:4f754b4cf82b r100:4f754b4cf82b r100:4f754b4cf82b r100:4f754b4cf82b r100:4f754b4cf82b r100:4f754b4cf82b r100:4f754b4cf82b r100:4f754b4cf82b r100:4f754b4cf82b r100:4f754b4cf82b r100:4f754b4cf82b r100:4f754b4cf82b r100:4f754b4cf82b r113:18a7ee8fa56e r100:4f754b4cf82b r576:f5bc148f7e1f r576:f5bc148f7e1f r100:4f754b4cf82b r220:a5d8c039f218 r566:c786cd201266 r100:4f754b4cf82b r100:4f754b4cf82b r113:18a7ee8fa56e r100:4f754b4cf82b r113:18a7ee8fa56e r100:4f754b4cf82b r100:4f754b4cf82b r100:4f754b4cf82b r113:18a7ee8fa56e r100:4f754b4cf82b r606:c5fd2d996909 r606:c5fd2d996909 r606:c5fd2d996909 r606:c5fd2d996909 r606:c5fd2d996909 r606:c5fd2d996909 r606:c5fd2d996909 r606:c5fd2d996909 r606:c5fd2d996909 r606:c5fd2d996909 r606:c5fd2d996909 r606:c5fd2d996909 r606:c5fd2d996909 r606:c5fd2d996909 r606:c5fd2d996909 r606:c5fd2d996909 r100:4f754b4cf82b r100:4f754b4cf82b r100:4f754b4cf82b r606:c5fd2d996909 r606:c5fd2d996909 r606:c5fd2d996909 r606:c5fd2d996909 r113:18a7ee8fa56e r113:18a7ee8fa56e r100:4f754b4cf82b r113:18a7ee8fa56e r100:4f754b4cf82b r113:18a7ee8fa56e r113:18a7ee8fa56e r113:18a7ee8fa56e r113:18a7ee8fa56e r100:4f754b4cf82b r606:c5fd2d996909 r606:c5fd2d996909 r100:4f754b4cf82b r631:33c6b6e755cd r631:33c6b6e755cd r631:33c6b6e755cd r100:4f754b4cf82b r209:765619b7cbb2 r100:4f754b4cf82b r100:4f754b4cf82b r100:4f754b4cf82b r113:18a7ee8fa56e r113:18a7ee8fa56e r113:18a7ee8fa56e r113:18a7ee8fa56e r113:18a7ee8fa56e r100:4f754b4cf82b r100:4f754b4cf82b r100:4f754b4cf82b r100:4f754b4cf82b r100:4f754b4cf82b r100:4f754b4cf82b r113:18a7ee8fa56e r100:4f754b4cf82b r113:18a7ee8fa56e r100:4f754b4cf82b r100:4f754b4cf82b r113:18a7ee8fa56e r100:4f754b4cf82b r113:18a7ee8fa56e r100:4f754b4cf82b r100:4f754b4cf82b r113:18a7ee8fa56e r209:765619b7cbb2 r209:765619b7cbb2 r100:4f754b4cf82b r100:4f754b4cf82b r100:4f754b4cf82b r100:4f754b4cf82b r113:18a7ee8fa56e r100:4f754b4cf82b r113:18a7ee8fa56e r113:18a7ee8fa56e r100:4f754b4cf82b r100:4f754b4cf82b r113:18a7ee8fa56e r100:4f754b4cf82b r113:18a7ee8fa56e r113:18a7ee8fa56e r100:4f754b4cf82b r100:4f754b4cf82b r113:18a7ee8fa56e r100:4f754b4cf82b r113:18a7ee8fa56e r113:18a7ee8fa56e r100:4f754b4cf82b r100:4f754b4cf82b r113:18a7ee8fa56e r100:4f754b4cf82b r113:18a7ee8fa56e r209:765619b7cbb2 r100:4f754b4cf82b r100:4f754b4cf82b r113:18a7ee8fa56e r100:4f754b4cf82b r209:765619b7cbb2 r606:c5fd2d996909 r100:4f754b4cf82b r100:4f754b4cf82b r100:4f754b4cf82b r113:18a7ee8fa56e r113:18a7ee8fa56e r100:4f754b4cf82b r113:18a7ee8fa56e r113:18a7ee8fa56e r113:18a7ee8fa56e r113:18a7ee8fa56e r100:4f754b4cf82b r100:4f754b4cf82b r100:4f754b4cf82b r209:765619b7cbb2 r113:18a7ee8fa56e r100:4f754b4cf82b r113:18a7ee8fa56e r100:4f754b4cf82b r100:4f754b4cf82b r606:c5fd2d996909 r100:4f754b4cf82b r100:4f754b4cf82b r113:18a7ee8fa56e r100:4f754b4cf82b r113:18a7ee8fa56e r100:4f754b4cf82b r100:4f754b4cf82b r606:c5fd2d996909 r100:4f754b4cf82b r100:4f754b4cf82b r113:18a7ee8fa56e r100:4f754b4cf82b r100:4f754b4cf82b r113:18a7ee8fa56e r113:18a7ee8fa56e r113:18a7ee8fa56e r113:18a7ee8fa56e r113:18a7ee8fa56e r100:4f754b4cf82b r100:4f754b4cf82b r100:4f754b4cf82b r113:18a7ee8fa56e r100:4f754b4cf82b r113:18a7ee8fa56e r113:18a7ee8fa56e r100:4f754b4cf82b r100:4f754b4cf82b r113:18a7ee8fa56e r100:4f754b4cf82b r100:4f754b4cf82b r100:4f754b4cf82b r100:4f754b4cf82b r100:4f754b4cf82b r100:4f754b4cf82b r209:765619b7cbb2 r209:765619b7cbb2 r209:765619b7cbb2 r209:765619b7cbb2 r113:18a7ee8fa56e r209:765619b7cbb2 r209:765619b7cbb2 r209:765619b7cbb2 r209:765619b7cbb2 r209:765619b7cbb2 r209:765619b7cbb2 r100:4f754b4cf82b r209:765619b7cbb2 r209:765619b7cbb2 r209:765619b7cbb2 r209:765619b7cbb2 r209:765619b7cbb2 r209:765619b7cbb2 r209:765619b7cbb2 r209:765619b7cbb2 r100:4f754b4cf82b r209:765619b7cbb2 r209:765619b7cbb2 r209:765619b7cbb2 r209:765619b7cbb2 r100:4f754b4cf82b r209:765619b7cbb2 r209:765619b7cbb2 r209:765619b7cbb2 r209:765619b7cbb2 r100:4f754b4cf82b r209:765619b7cbb2 r209:765619b7cbb2 r209:765619b7cbb2 r209:765619b7cbb2 r209:765619b7cbb2 r209:765619b7cbb2 r100:4f754b4cf82b r209:765619b7cbb2 r209:765619b7cbb2 r209:765619b7cbb2 r100:4f754b4cf82b r209:765619b7cbb2 r209:765619b7cbb2 r209:765619b7cbb2 r209:765619b7cbb2 r209:765619b7cbb2 r209:765619b7cbb2 r100:4f754b4cf82b r209:765619b7cbb2 r209:765619b7cbb2 r209:765619b7cbb2 r100:4f754b4cf82b r209:765619b7cbb2 r209:765619b7cbb2 r100:4f754b4cf82b r209:765619b7cbb2 r209:765619b7cbb2 r209:765619b7cbb2 r209:765619b7cbb2 r209:765619b7cbb2 r209:765619b7cbb2 r209:765619b7cbb2 r100:4f754b4cf82b r100:4f754b4cf82b r100:4f754b4cf82b r100:4f754b4cf82b r100:4f754b4cf82b r100:4f754b4cf82b r576:f5bc148f7e1f | /* -*- mode: C++; indent-tabs-mode: nil; -*-
*
* This file is a part of LEMON, a generic C++ optimization library.
*
* Copyright (C) 2003-2009
* Egervary Jeno Kombinatorikus Optimalizalasi Kutatocsoport
* (Egervary Research Group on Combinatorial Optimization, EGRES).
*
* Permission to use, modify and distribute this software is granted
* provided that this copyright notice appears in all copies. For
* precise terms see the accompanying LICENSE file.
*
* This software is provided "AS IS" with no warranty of any kind,
* express or implied, and with no claim as to its suitability for any
* purpose.
*
*/
///\ingroup concept
///\file
///\brief The concept of heaps.
#ifndef LEMON_CONCEPTS_HEAP_H
#define LEMON_CONCEPTS_HEAP_H
#include <lemon/core.h>
#include <lemon/concept_check.h>
namespace lemon {
namespace concepts {
/// \addtogroup concept
/// @{
/// \brief The heap concept.
///
/// Concept class describing the main interface of heaps. A \e heap
/// is a data structure for storing items with specified values called
/// \e priorities in such a way that finding the item with minimum
/// priority is efficient. In a heap one can change the priority of an
/// item, add or erase an item, etc.
///
/// \tparam PR Type of the priority of the items.
/// \tparam IM A read and writable item map with int values, used
/// internally to handle the cross references.
/// \tparam Comp A functor class for the ordering of the priorities.
/// The default is \c std::less<PR>.
#ifdef DOXYGEN
template <typename PR, typename IM, typename Comp = std::less<PR> >
#else
template <typename PR, typename IM>
#endif
class Heap {
public:
/// Type of the item-int map.
typedef IM ItemIntMap;
/// Type of the priorities.
typedef PR Prio;
/// Type of the items stored in the heap.
typedef typename ItemIntMap::Key Item;
/// \brief Type to represent the states of the items.
///
/// Each item has a state associated to it. It can be "in heap",
/// "pre heap" or "post heap". The later two are indifferent
/// from the point of view of the heap, but may be useful for
/// the user.
///
/// The item-int map must be initialized in such way that it assigns
/// \c PRE_HEAP (<tt>-1</tt>) to any element to be put in the heap.
enum State {
IN_HEAP = 0, ///< = 0. The "in heap" state constant.
PRE_HEAP = -1, ///< = -1. The "pre heap" state constant.
POST_HEAP = -2 ///< = -2. The "post heap" state constant.
};
/// \brief The constructor.
///
/// The constructor.
/// \param map A map that assigns \c int values to keys of type
/// \c Item. It is used internally by the heap implementations to
/// handle the cross references. The assigned value must be
/// \c PRE_HEAP (<tt>-1</tt>) for every item.
explicit Heap(ItemIntMap &map) {}
/// \brief The number of items stored in the heap.
///
/// Returns the number of items stored in the heap.
int size() const { return 0; }
/// \brief Checks if the heap is empty.
///
/// Returns \c true if the heap is empty.
bool empty() const { return false; }
/// \brief Makes the heap empty.
///
/// Makes the heap empty.
void clear();
/// \brief Inserts an item into the heap with the given priority.
///
/// Inserts the given item into the heap with the given priority.
/// \param i The item to insert.
/// \param p The priority of the item.
void push(const Item &i, const Prio &p) {}
/// \brief Returns the item having minimum priority.
///
/// Returns the item having minimum priority.
/// \pre The heap must be non-empty.
Item top() const {}
/// \brief The minimum priority.
///
/// Returns the minimum priority.
/// \pre The heap must be non-empty.
Prio prio() const {}
/// \brief Removes the item having minimum priority.
///
/// Removes the item having minimum priority.
/// \pre The heap must be non-empty.
void pop() {}
/// \brief Removes an item from the heap.
///
/// Removes the given item from the heap if it is already stored.
/// \param i The item to delete.
void erase(const Item &i) {}
/// \brief The priority of an item.
///
/// Returns the priority of the given item.
/// \param i The item.
/// \pre \c i must be in the heap.
Prio operator[](const Item &i) const {}
/// \brief Sets the priority of an item or inserts it, if it is
/// not stored in the heap.
///
/// This method sets the priority of the given item if it is
/// already stored in the heap.
/// Otherwise it inserts the given item with the given priority.
///
/// \param i The item.
/// \param p The priority.
void set(const Item &i, const Prio &p) {}
/// \brief Decreases the priority of an item to the given value.
///
/// Decreases the priority of an item to the given value.
/// \param i The item.
/// \param p The priority.
/// \pre \c i must be stored in the heap with priority at least \c p.
void decrease(const Item &i, const Prio &p) {}
/// \brief Increases the priority of an item to the given value.
///
/// Increases the priority of an item to the given value.
/// \param i The item.
/// \param p The priority.
/// \pre \c i must be stored in the heap with priority at most \c p.
void increase(const Item &i, const Prio &p) {}
/// \brief Returns if an item is in, has already been in, or has
/// never been in the heap.
///
/// This method returns \c PRE_HEAP if the given item has never
/// been in the heap, \c IN_HEAP if it is in the heap at the moment,
/// and \c POST_HEAP otherwise.
/// In the latter case it is possible that the item will get back
/// to the heap again.
/// \param i The item.
State state(const Item &i) const {}
/// \brief Sets the state of an item in the heap.
///
/// Sets the state of the given item in the heap. It can be used
/// to manually clear the heap when it is important to achive the
/// better time complexity.
/// \param i The item.
/// \param st The state. It should not be \c IN_HEAP.
void state(const Item& i, State st) {}
template <typename _Heap>
struct Constraints {
public:
void constraints() {
typedef typename _Heap::Item OwnItem;
typedef typename _Heap::Prio OwnPrio;
typedef typename _Heap::State OwnState;
Item item;
Prio prio;
item=Item();
prio=Prio();
ignore_unused_variable_warning(item);
ignore_unused_variable_warning(prio);
OwnItem own_item;
OwnPrio own_prio;
OwnState own_state;
own_item=Item();
own_prio=Prio();
ignore_unused_variable_warning(own_item);
ignore_unused_variable_warning(own_prio);
ignore_unused_variable_warning(own_state);
_Heap heap1(map);
_Heap heap2 = heap1;
ignore_unused_variable_warning(heap1);
ignore_unused_variable_warning(heap2);
int s = heap.size();
ignore_unused_variable_warning(s);
bool e = heap.empty();
ignore_unused_variable_warning(e);
prio = heap.prio();
item = heap.top();
prio = heap[item];
own_prio = heap.prio();
own_item = heap.top();
own_prio = heap[own_item];
heap.push(item, prio);
heap.push(own_item, own_prio);
heap.pop();
heap.set(item, prio);
heap.decrease(item, prio);
heap.increase(item, prio);
heap.set(own_item, own_prio);
heap.decrease(own_item, own_prio);
heap.increase(own_item, own_prio);
heap.erase(item);
heap.erase(own_item);
heap.clear();
own_state = heap.state(own_item);
heap.state(own_item, own_state);
own_state = _Heap::PRE_HEAP;
own_state = _Heap::IN_HEAP;
own_state = _Heap::POST_HEAP;
}
_Heap& heap;
ItemIntMap& map;
};
};
/// @}
} // namespace lemon
}
#endif
|