0
3
0
... | ... |
@@ -184,97 +184,97 @@ |
184 | 184 |
typedef typename In::value_type::second_type Value; |
185 | 185 |
|
186 | 186 |
static Value kruskal(const Graph& graph, const In& in, Out& out) { |
187 | 187 |
return KruskalOutputSelector<Graph, In, Out>:: |
188 | 188 |
kruskal(graph, in, out); |
189 | 189 |
} |
190 | 190 |
|
191 | 191 |
}; |
192 | 192 |
|
193 | 193 |
template <typename Graph, typename In, typename Out> |
194 | 194 |
struct KruskalInputSelector<Graph, In, Out, |
195 | 195 |
typename enable_if<MapInputIndicator<In>, void>::type > |
196 | 196 |
{ |
197 | 197 |
typedef typename In::Value Value; |
198 | 198 |
static Value kruskal(const Graph& graph, const In& in, Out& out) { |
199 | 199 |
typedef typename In::Key MapArc; |
200 | 200 |
typedef typename In::Value Value; |
201 | 201 |
typedef typename ItemSetTraits<Graph, MapArc>::ItemIt MapArcIt; |
202 | 202 |
typedef std::vector<std::pair<MapArc, Value> > Sequence; |
203 | 203 |
Sequence seq; |
204 | 204 |
|
205 | 205 |
for (MapArcIt it(graph); it != INVALID; ++it) { |
206 | 206 |
seq.push_back(std::make_pair(it, in[it])); |
207 | 207 |
} |
208 | 208 |
|
209 | 209 |
std::sort(seq.begin(), seq.end(), PairComp<Sequence>()); |
210 | 210 |
return KruskalOutputSelector<Graph, Sequence, Out>:: |
211 | 211 |
kruskal(graph, seq, out); |
212 | 212 |
} |
213 | 213 |
}; |
214 | 214 |
|
215 | 215 |
template <typename T> |
216 | 216 |
struct RemoveConst { |
217 | 217 |
typedef T type; |
218 | 218 |
}; |
219 | 219 |
|
220 | 220 |
template <typename T> |
221 | 221 |
struct RemoveConst<const T> { |
222 | 222 |
typedef T type; |
223 | 223 |
}; |
224 | 224 |
|
225 | 225 |
template <typename Graph, typename In, typename Out> |
226 | 226 |
struct KruskalOutputSelector<Graph, In, Out, |
227 | 227 |
typename enable_if<SequenceOutputIndicator<Out>, void>::type > |
228 | 228 |
{ |
229 | 229 |
typedef typename In::value_type::second_type Value; |
230 | 230 |
|
231 | 231 |
static Value kruskal(const Graph& graph, const In& in, Out& out) { |
232 |
typedef |
|
232 |
typedef LoggerBoolMap<typename RemoveConst<Out>::type> Map; |
|
233 | 233 |
Map map(out); |
234 | 234 |
return _kruskal_bits::kruskal(graph, in, map); |
235 | 235 |
} |
236 | 236 |
|
237 | 237 |
}; |
238 | 238 |
|
239 | 239 |
template <typename Graph, typename In, typename Out> |
240 | 240 |
struct KruskalOutputSelector<Graph, In, Out, |
241 | 241 |
typename enable_if<MapOutputIndicator<Out>, void>::type > |
242 | 242 |
{ |
243 | 243 |
typedef typename In::value_type::second_type Value; |
244 | 244 |
|
245 | 245 |
static Value kruskal(const Graph& graph, const In& in, Out& out) { |
246 | 246 |
return _kruskal_bits::kruskal(graph, in, out); |
247 | 247 |
} |
248 | 248 |
}; |
249 | 249 |
|
250 | 250 |
} |
251 | 251 |
|
252 | 252 |
/// \ingroup spantree |
253 | 253 |
/// |
254 | 254 |
/// \brief Kruskal's algorithm to find a minimum cost tree of a graph. |
255 | 255 |
/// |
256 | 256 |
/// This function runs Kruskal's algorithm to find a minimum cost tree. |
257 | 257 |
/// Due to some C++ hacking, it accepts various input and output types. |
258 | 258 |
/// |
259 | 259 |
/// \param g The graph the algorithm runs on. |
260 | 260 |
/// It can be either \ref concepts::Digraph "directed" or |
261 | 261 |
/// \ref concepts::Graph "undirected". |
262 | 262 |
/// If the graph is directed, the algorithm consider it to be |
263 | 263 |
/// undirected by disregarding the direction of the arcs. |
264 | 264 |
/// |
265 | 265 |
/// \param in This object is used to describe the arc costs. It can be one |
266 | 266 |
/// of the following choices. |
267 | 267 |
/// - An STL compatible 'Forward Container' with |
268 | 268 |
/// <tt>std::pair<GR::Edge,X></tt> or |
269 | 269 |
/// <tt>std::pair<GR::Arc,X></tt> as its <tt>value_type</tt>, where |
270 | 270 |
/// \c X is the type of the costs. The pairs indicates the arcs |
271 | 271 |
/// along with the assigned cost. <em>They must be in a |
272 | 272 |
/// cost-ascending order.</em> |
273 | 273 |
/// - Any readable Arc map. The values of the map indicate the arc costs. |
274 | 274 |
/// |
275 | 275 |
/// \retval out Here we also have a choise. |
276 | 276 |
/// - It can be a writable \c bool arc map. After running the |
277 | 277 |
/// algorithm this will contain the found minimum cost spanning |
278 | 278 |
/// tree: the value of an arc will be set to \c true if it belongs |
279 | 279 |
/// to the tree, otherwise it will be set to \c false. The value of |
280 | 280 |
/// each arc will be set exactly once. |
... | ... |
@@ -1655,132 +1655,132 @@ |
1655 | 1655 |
/// Constructor |
1656 | 1656 |
LessMap(const M1 &m1, const M2 &m2) : _m1(m1), _m2(m2) {} |
1657 | 1657 |
/// \e |
1658 | 1658 |
Value operator[](const Key &k) const { return _m1[k]<_m2[k]; } |
1659 | 1659 |
}; |
1660 | 1660 |
|
1661 | 1661 |
/// Returns an \ref LessMap class |
1662 | 1662 |
|
1663 | 1663 |
/// This function just returns an \ref LessMap class. |
1664 | 1664 |
/// |
1665 | 1665 |
/// For example, if \c m1 and \c m2 are maps with keys and values of |
1666 | 1666 |
/// the same type, then <tt>lessMap(m1,m2)[x]</tt> will be equal to |
1667 | 1667 |
/// <tt>m1[x]<m2[x]</tt>. |
1668 | 1668 |
/// |
1669 | 1669 |
/// \relates LessMap |
1670 | 1670 |
template<typename M1, typename M2> |
1671 | 1671 |
inline LessMap<M1, M2> lessMap(const M1 &m1, const M2 &m2) { |
1672 | 1672 |
return LessMap<M1, M2>(m1,m2); |
1673 | 1673 |
} |
1674 | 1674 |
|
1675 | 1675 |
namespace _maps_bits { |
1676 | 1676 |
|
1677 | 1677 |
template <typename _Iterator, typename Enable = void> |
1678 | 1678 |
struct IteratorTraits { |
1679 | 1679 |
typedef typename std::iterator_traits<_Iterator>::value_type Value; |
1680 | 1680 |
}; |
1681 | 1681 |
|
1682 | 1682 |
template <typename _Iterator> |
1683 | 1683 |
struct IteratorTraits<_Iterator, |
1684 | 1684 |
typename exists<typename _Iterator::container_type>::type> |
1685 | 1685 |
{ |
1686 | 1686 |
typedef typename _Iterator::container_type::value_type Value; |
1687 | 1687 |
}; |
1688 | 1688 |
|
1689 | 1689 |
} |
1690 | 1690 |
|
1691 | 1691 |
/// \brief Writable bool map for logging each \c true assigned element |
1692 | 1692 |
/// |
1693 | 1693 |
/// A \ref concepts::WriteMap "writable" bool map for logging |
1694 | 1694 |
/// each \c true assigned element, i.e it copies subsequently each |
1695 | 1695 |
/// keys set to \c true to the given iterator. |
1696 | 1696 |
/// The most important usage of it is storing certain nodes or arcs |
1697 | 1697 |
/// that were marked \c true by an algorithm. |
1698 | 1698 |
/// |
1699 | 1699 |
/// There are several algorithms that provide solutions through bool |
1700 | 1700 |
/// maps and most of them assign \c true at most once for each key. |
1701 | 1701 |
/// In these cases it is a natural request to store each \c true |
1702 | 1702 |
/// assigned elements (in order of the assignment), which can be |
1703 |
/// easily done with |
|
1703 |
/// easily done with LoggerBoolMap. |
|
1704 | 1704 |
/// |
1705 |
/// The simplest way of using this map is through the |
|
1705 |
/// The simplest way of using this map is through the loggerBoolMap() |
|
1706 | 1706 |
/// function. |
1707 | 1707 |
/// |
1708 | 1708 |
/// \tparam It The type of the iterator. |
1709 | 1709 |
/// \tparam Ke The key type of the map. The default value set |
1710 | 1710 |
/// according to the iterator type should work in most cases. |
1711 | 1711 |
/// |
1712 | 1712 |
/// \note The container of the iterator must contain enough space |
1713 | 1713 |
/// for the elements or the iterator should be an inserter iterator. |
1714 | 1714 |
#ifdef DOXYGEN |
1715 | 1715 |
template <typename It, typename Ke> |
1716 | 1716 |
#else |
1717 | 1717 |
template <typename It, |
1718 | 1718 |
typename Ke=typename _maps_bits::IteratorTraits<It>::Value> |
1719 | 1719 |
#endif |
1720 |
class |
|
1720 |
class LoggerBoolMap { |
|
1721 | 1721 |
public: |
1722 | 1722 |
typedef It Iterator; |
1723 | 1723 |
|
1724 | 1724 |
typedef Ke Key; |
1725 | 1725 |
typedef bool Value; |
1726 | 1726 |
|
1727 | 1727 |
/// Constructor |
1728 |
|
|
1728 |
LoggerBoolMap(Iterator it) |
|
1729 | 1729 |
: _begin(it), _end(it) {} |
1730 | 1730 |
|
1731 | 1731 |
/// Gives back the given iterator set for the first key |
1732 | 1732 |
Iterator begin() const { |
1733 | 1733 |
return _begin; |
1734 | 1734 |
} |
1735 | 1735 |
|
1736 | 1736 |
/// Gives back the the 'after the last' iterator |
1737 | 1737 |
Iterator end() const { |
1738 | 1738 |
return _end; |
1739 | 1739 |
} |
1740 | 1740 |
|
1741 | 1741 |
/// The set function of the map |
1742 | 1742 |
void set(const Key& key, Value value) { |
1743 | 1743 |
if (value) { |
1744 | 1744 |
*_end++ = key; |
1745 | 1745 |
} |
1746 | 1746 |
} |
1747 | 1747 |
|
1748 | 1748 |
private: |
1749 | 1749 |
Iterator _begin; |
1750 | 1750 |
Iterator _end; |
1751 | 1751 |
}; |
1752 | 1752 |
|
1753 |
/// Returns a \ref |
|
1753 |
/// Returns a \ref LoggerBoolMap class |
|
1754 | 1754 |
|
1755 |
/// This function just returns a \ref |
|
1755 |
/// This function just returns a \ref LoggerBoolMap class. |
|
1756 | 1756 |
/// |
1757 | 1757 |
/// The most important usage of it is storing certain nodes or arcs |
1758 | 1758 |
/// that were marked \c true by an algorithm. |
1759 | 1759 |
/// For example it makes easier to store the nodes in the processing |
1760 | 1760 |
/// order of Dfs algorithm, as the following examples show. |
1761 | 1761 |
/// \code |
1762 | 1762 |
/// std::vector<Node> v; |
1763 |
/// dfs(g,s).processedMap( |
|
1763 |
/// dfs(g,s).processedMap(loggerBoolMap(std::back_inserter(v))).run(); |
|
1764 | 1764 |
/// \endcode |
1765 | 1765 |
/// \code |
1766 | 1766 |
/// std::vector<Node> v(countNodes(g)); |
1767 |
/// dfs(g,s).processedMap( |
|
1767 |
/// dfs(g,s).processedMap(loggerBoolMap(v.begin())).run(); |
|
1768 | 1768 |
/// \endcode |
1769 | 1769 |
/// |
1770 | 1770 |
/// \note The container of the iterator must contain enough space |
1771 | 1771 |
/// for the elements or the iterator should be an inserter iterator. |
1772 | 1772 |
/// |
1773 |
/// \note |
|
1773 |
/// \note LoggerBoolMap is just \ref concepts::WriteMap "writable", so |
|
1774 | 1774 |
/// it cannot be used when a readable map is needed, for example as |
1775 |
/// \c ReachedMap for Bfs, Dfs and Dijkstra algorithms. |
|
1775 |
/// \c ReachedMap for \ref Bfs, \ref Dfs and \ref Dijkstra algorithms. |
|
1776 | 1776 |
/// |
1777 |
/// \relates |
|
1777 |
/// \relates LoggerBoolMap |
|
1778 | 1778 |
template<typename Iterator> |
1779 |
inline StoreBoolMap<Iterator> storeBoolMap(Iterator it) { |
|
1780 |
return StoreBoolMap<Iterator>(it); |
|
1779 |
inline LoggerBoolMap<Iterator> loggerBoolMap(Iterator it) { |
|
1780 |
return LoggerBoolMap<Iterator>(it); |
|
1781 | 1781 |
} |
1782 | 1782 |
|
1783 | 1783 |
/// @} |
1784 | 1784 |
} |
1785 | 1785 |
|
1786 | 1786 |
#endif // LEMON_MAPS_H |
... | ... |
@@ -260,72 +260,72 @@ |
260 | 260 |
"Something is wrong with ShiftWriteMap"); |
261 | 261 |
check(scaleMap(id, 2.0)[1] == 2.0 && scaleMap(id, 2.0)[10] == 20.0, |
262 | 262 |
"Something is wrong with ScaleMap"); |
263 | 263 |
check(scaleWriteMap(id, 2.0)[1] == 2.0 && scaleWriteMap(id, 2.0)[10] == 20.0, |
264 | 264 |
"Something is wrong with ScaleWriteMap"); |
265 | 265 |
check(negMap(id)[1] == -1.0 && negMap(id)[-10] == 10.0, |
266 | 266 |
"Something is wrong with NegMap"); |
267 | 267 |
check(negWriteMap(id)[1] == -1.0 && negWriteMap(id)[-10] == 10.0, |
268 | 268 |
"Something is wrong with NegWriteMap"); |
269 | 269 |
check(absMap(id)[1] == 1.0 && absMap(id)[-10] == 10.0, |
270 | 270 |
"Something is wrong with AbsMap"); |
271 | 271 |
} |
272 | 272 |
|
273 | 273 |
// Logical maps: |
274 | 274 |
// - TrueMap, FalseMap |
275 | 275 |
// - AndMap, OrMap |
276 | 276 |
// - NotMap, NotWriteMap |
277 | 277 |
// - EqualMap, LessMap |
278 | 278 |
{ |
279 | 279 |
checkConcept<BoolMap, TrueMap<A> >(); |
280 | 280 |
checkConcept<BoolMap, FalseMap<A> >(); |
281 | 281 |
checkConcept<BoolMap, AndMap<BoolMap,BoolMap> >(); |
282 | 282 |
checkConcept<BoolMap, OrMap<BoolMap,BoolMap> >(); |
283 | 283 |
checkConcept<BoolMap, NotMap<BoolMap> >(); |
284 | 284 |
checkConcept<BoolWriteMap, NotWriteMap<BoolWriteMap> >(); |
285 | 285 |
checkConcept<BoolMap, EqualMap<DoubleMap,DoubleMap> >(); |
286 | 286 |
checkConcept<BoolMap, LessMap<DoubleMap,DoubleMap> >(); |
287 | 287 |
|
288 | 288 |
TrueMap<int> tm; |
289 | 289 |
FalseMap<int> fm; |
290 | 290 |
RangeMap<bool> rm(2); |
291 | 291 |
rm[0] = true; rm[1] = false; |
292 | 292 |
check(andMap(tm,rm)[0] && !andMap(tm,rm)[1] && !andMap(fm,rm)[0] && !andMap(fm,rm)[1], |
293 | 293 |
"Something is wrong with AndMap"); |
294 | 294 |
check(orMap(tm,rm)[0] && orMap(tm,rm)[1] && orMap(fm,rm)[0] && !orMap(fm,rm)[1], |
295 | 295 |
"Something is wrong with OrMap"); |
296 | 296 |
check(!notMap(rm)[0] && notMap(rm)[1], "Something is wrong with NotMap"); |
297 | 297 |
check(!notWriteMap(rm)[0] && notWriteMap(rm)[1], "Something is wrong with NotWriteMap"); |
298 | 298 |
|
299 | 299 |
ConstMap<int, double> cm(2.0); |
300 | 300 |
IdentityMap<int> im; |
301 | 301 |
ConvertMap<IdentityMap<int>, double> id(im); |
302 | 302 |
check(lessMap(id,cm)[1] && !lessMap(id,cm)[2] && !lessMap(id,cm)[3], |
303 | 303 |
"Something is wrong with LessMap"); |
304 | 304 |
check(!equalMap(id,cm)[1] && equalMap(id,cm)[2] && !equalMap(id,cm)[3], |
305 | 305 |
"Something is wrong with EqualMap"); |
306 | 306 |
} |
307 | 307 |
|
308 |
// |
|
308 |
// LoggerBoolMap |
|
309 | 309 |
{ |
310 | 310 |
typedef std::vector<int> vec; |
311 | 311 |
vec v1; |
312 | 312 |
vec v2(10); |
313 |
StoreBoolMap<std::back_insert_iterator<vec> > map1(std::back_inserter(v1)); |
|
314 |
StoreBoolMap<vec::iterator> map2(v2.begin()); |
|
313 |
LoggerBoolMap<std::back_insert_iterator<vec> > map1(std::back_inserter(v1)); |
|
314 |
LoggerBoolMap<vec::iterator> map2(v2.begin()); |
|
315 | 315 |
map1.set(10, false); |
316 | 316 |
map1.set(20, true); map2.set(20, true); |
317 | 317 |
map1.set(30, false); map2.set(40, false); |
318 | 318 |
map1.set(50, true); map2.set(50, true); |
319 | 319 |
map1.set(60, true); map2.set(60, true); |
320 | 320 |
check(v1.size() == 3 && v2.size() == 10 && |
321 | 321 |
v1[0]==20 && v1[1]==50 && v1[2]==60 && v2[0]==20 && v2[1]==50 && v2[2]==60, |
322 |
"Something is wrong with |
|
322 |
"Something is wrong with LoggerBoolMap"); |
|
323 | 323 |
|
324 | 324 |
int i = 0; |
325 |
for ( |
|
325 |
for ( LoggerBoolMap<vec::iterator>::Iterator it = map2.begin(); |
|
326 | 326 |
it != map2.end(); ++it ) |
327 |
check(v1[i++] == *it, "Something is wrong with |
|
327 |
check(v1[i++] == *it, "Something is wrong with LoggerBoolMap"); |
|
328 | 328 |
} |
329 | 329 |
|
330 | 330 |
return 0; |
331 | 331 |
} |
0 comments (0 inline)