lemon-1.2: lemon/maps.h@0821fea54d53 (annotated)

alpar@209	1	/* -- mode: C++; indent-tabs-mode: nil; --
alpar@25	2	*
alpar@209	3	* This file is a part of LEMON, a generic C++ optimization library.
alpar@25	4	*
alpar@39	5	* Copyright (C) 2003-2008
alpar@25	6	* Egervary Jeno Kombinatorikus Optimalizalasi Kutatocsoport
alpar@25	7	* (Egervary Research Group on Combinatorial Optimization, EGRES).
alpar@25	8	*
alpar@25	9	* Permission to use, modify and distribute this software is granted
alpar@25	10	* provided that this copyright notice appears in all copies. For
alpar@25	11	* precise terms see the accompanying LICENSE file.
alpar@25	12	*
alpar@25	13	* This software is provided "AS IS" with no warranty of any kind,
alpar@25	14	* express or implied, and with no claim as to its suitability for any
alpar@25	15	* purpose.
alpar@25	16	*
alpar@25	17	*/
alpar@25	18
alpar@25	19	#ifndef LEMON_MAPS_H
alpar@25	20	#define LEMON_MAPS_H
alpar@25	21
alpar@25	22	#include <iterator>
alpar@25	23	#include <functional>
alpar@25	24	#include <vector>
alpar@25	25
deba@220	26	#include <lemon/core.h>
alpar@25	27
alpar@25	28	///\file
alpar@25	29	///\ingroup maps
alpar@25	30	///\brief Miscellaneous property maps
kpeter@80	31
alpar@25	32	#include <map>
alpar@25	33
alpar@25	34	namespace lemon {
alpar@25	35
alpar@25	36	/// \addtogroup maps
alpar@25	37	/// @{
alpar@25	38
alpar@25	39	/// Base class of maps.
alpar@25	40
kpeter@80	41	/// Base class of maps. It provides the necessary type definitions
kpeter@80	42	/// required by the map %concepts.
kpeter@80	43	template<typename K, typename V>
alpar@25	44	class MapBase {
alpar@25	45	public:
kpeter@80	46	/// \biref The key type of the map.
alpar@25	47	typedef K Key;
kpeter@80	48	/// \brief The value type of the map.
kpeter@80	49	/// (The type of objects associated with the keys).
kpeter@80	50	typedef V Value;
alpar@25	51	};
alpar@25	52
kpeter@80	53
alpar@25	54	/// Null map. (a.k.a. DoNothingMap)
alpar@25	55
kpeter@29	56	/// This map can be used if you have to provide a map only for
kpeter@80	57	/// its type definitions, or if you have to provide a writable map,
kpeter@80	58	/// but data written to it is not required (i.e. it will be sent to
kpeter@29	59	/// <tt>/dev/null</tt>).
kpeter@80	60	/// It conforms the \ref concepts::ReadWriteMap "ReadWriteMap" concept.
kpeter@80	61	///
kpeter@80	62	/// \sa ConstMap
kpeter@80	63	template<typename K, typename V>
kpeter@80	64	class NullMap : public MapBase<K, V> {
alpar@25	65	public:
kpeter@80	66	typedef MapBase<K, V> Parent;
alpar@25	67	typedef typename Parent::Key Key;
alpar@25	68	typedef typename Parent::Value Value;
kpeter@80	69
alpar@25	70	/// Gives back a default constructed element.
kpeter@80	71	Value operator[](const Key&) const { return Value(); }
alpar@25	72	/// Absorbs the value.
kpeter@80	73	void set(const Key&, const Value&) {}
alpar@25	74	};
alpar@25	75
kpeter@80	76	/// Returns a \ref NullMap class
kpeter@29	77
kpeter@80	78	/// This function just returns a \ref NullMap class.
kpeter@80	79	/// \relates NullMap
kpeter@80	80	template <typename K, typename V>
alpar@25	81	NullMap<K, V> nullMap() {
alpar@25	82	return NullMap<K, V>();
alpar@25	83	}
alpar@25	84
alpar@25	85
alpar@25	86	/// Constant map.
alpar@25	87
kpeter@82	88	/// This \ref concepts::ReadMap "readable map" assigns a specified
kpeter@82	89	/// value to each key.
kpeter@80	90	///
kpeter@80	91	/// In other aspects it is equivalent to \ref NullMap.
kpeter@80	92	/// So it conforms the \ref concepts::ReadWriteMap "ReadWriteMap"
kpeter@80	93	/// concept, but it absorbs the data written to it.
kpeter@80	94	///
kpeter@80	95	/// The simplest way of using this map is through the constMap()
kpeter@80	96	/// function.
kpeter@80	97	///
kpeter@80	98	/// \sa NullMap
kpeter@80	99	/// \sa IdentityMap
kpeter@80	100	template<typename K, typename V>
kpeter@80	101	class ConstMap : public MapBase<K, V> {
alpar@25	102	private:
kpeter@80	103	V _value;
alpar@25	104	public:
kpeter@80	105	typedef MapBase<K, V> Parent;
alpar@25	106	typedef typename Parent::Key Key;
alpar@25	107	typedef typename Parent::Value Value;
alpar@25	108
alpar@25	109	/// Default constructor
alpar@25	110
kpeter@29	111	/// Default constructor.
kpeter@80	112	/// The value of the map will be default constructed.
alpar@25	113	ConstMap() {}
kpeter@80	114
kpeter@29	115	/// Constructor with specified initial value
alpar@25	116
kpeter@29	117	/// Constructor with specified initial value.
kpeter@123	118	/// \param v The initial value of the map.
kpeter@80	119	ConstMap(const Value &v) : _value(v) {}
alpar@25	120
kpeter@80	121	/// Gives back the specified value.
kpeter@80	122	Value operator[](const Key&) const { return _value; }
alpar@25	123
kpeter@80	124	/// Absorbs the value.
kpeter@80	125	void set(const Key&, const Value&) {}
kpeter@80	126
kpeter@80	127	/// Sets the value that is assigned to each key.
kpeter@80	128	void setAll(const Value &v) {
kpeter@80	129	_value = v;
kpeter@80	130	}
kpeter@80	131
kpeter@80	132	template<typename V1>
kpeter@80	133	ConstMap(const ConstMap<K, V1> &, const Value &v) : _value(v) {}
alpar@25	134	};
alpar@25	135
kpeter@80	136	/// Returns a \ref ConstMap class
alpar@25	137
kpeter@80	138	/// This function just returns a \ref ConstMap class.
kpeter@80	139	/// \relates ConstMap
kpeter@80	140	template<typename K, typename V>
alpar@25	141	inline ConstMap<K, V> constMap(const V &v) {
alpar@25	142	return ConstMap<K, V>(v);
alpar@25	143	}
alpar@25	144
kpeter@123	145	template<typename K, typename V>
kpeter@123	146	inline ConstMap<K, V> constMap() {
kpeter@123	147	return ConstMap<K, V>();
kpeter@123	148	}
kpeter@123	149
alpar@25	150
alpar@25	151	template<typename T, T v>
kpeter@80	152	struct Const {};
alpar@25	153
alpar@25	154	/// Constant map with inlined constant value.
alpar@25	155
kpeter@82	156	/// This \ref concepts::ReadMap "readable map" assigns a specified
kpeter@82	157	/// value to each key.
kpeter@80	158	///
kpeter@80	159	/// In other aspects it is equivalent to \ref NullMap.
kpeter@80	160	/// So it conforms the \ref concepts::ReadWriteMap "ReadWriteMap"
kpeter@80	161	/// concept, but it absorbs the data written to it.
kpeter@80	162	///
kpeter@80	163	/// The simplest way of using this map is through the constMap()
kpeter@80	164	/// function.
kpeter@80	165	///
kpeter@80	166	/// \sa NullMap
kpeter@80	167	/// \sa IdentityMap
alpar@25	168	template<typename K, typename V, V v>
alpar@25	169	class ConstMap<K, Const<V, v> > : public MapBase<K, V> {
alpar@25	170	public:
alpar@25	171	typedef MapBase<K, V> Parent;
alpar@25	172	typedef typename Parent::Key Key;
alpar@25	173	typedef typename Parent::Value Value;
alpar@25	174
kpeter@80	175	/// Constructor.
kpeter@80	176	ConstMap() {}
kpeter@80	177
kpeter@80	178	/// Gives back the specified value.
kpeter@80	179	Value operator[](const Key&) const { return v; }
kpeter@80	180
kpeter@80	181	/// Absorbs the value.
kpeter@80	182	void set(const Key&, const Value&) {}
alpar@25	183	};
alpar@25	184
kpeter@80	185	/// Returns a \ref ConstMap class with inlined constant value
alpar@25	186
kpeter@80	187	/// This function just returns a \ref ConstMap class with inlined
kpeter@80	188	/// constant value.
kpeter@80	189	/// \relates ConstMap
kpeter@80	190	template<typename K, typename V, V v>
alpar@25	191	inline ConstMap<K, Const<V, v> > constMap() {
alpar@25	192	return ConstMap<K, Const<V, v> >();
alpar@25	193	}
alpar@25	194
alpar@25	195
kpeter@82	196	/// Identity map.
kpeter@82	197
kpeter@82	198	/// This \ref concepts::ReadMap "read-only map" gives back the given
kpeter@82	199	/// key as value without any modification.
kpeter@80	200	///
kpeter@80	201	/// \sa ConstMap
kpeter@80	202	template <typename T>
kpeter@80	203	class IdentityMap : public MapBase<T, T> {
kpeter@80	204	public:
kpeter@80	205	typedef MapBase<T, T> Parent;
kpeter@80	206	typedef typename Parent::Key Key;
kpeter@80	207	typedef typename Parent::Value Value;
kpeter@80	208
kpeter@80	209	/// Gives back the given value without any modification.
kpeter@82	210	Value operator[](const Key &k) const {
kpeter@82	211	return k;
kpeter@80	212	}
kpeter@80	213	};
kpeter@80	214
kpeter@80	215	/// Returns an \ref IdentityMap class
kpeter@80	216
kpeter@80	217	/// This function just returns an \ref IdentityMap class.
kpeter@80	218	/// \relates IdentityMap
kpeter@80	219	template<typename T>
kpeter@80	220	inline IdentityMap<T> identityMap() {
kpeter@80	221	return IdentityMap<T>();
kpeter@80	222	}
kpeter@80	223
kpeter@80	224
kpeter@80	225	/// \brief Map for storing values for integer keys from the range
kpeter@80	226	/// <tt>[0..size-1]</tt>.
kpeter@80	227	///
kpeter@80	228	/// This map is essentially a wrapper for \c std::vector. It assigns
kpeter@80	229	/// values to integer keys from the range <tt>[0..size-1]</tt>.
kpeter@80	230	/// It can be used with some data structures, for example
kpeter@80	231	/// \ref UnionFind, \ref BinHeap, when the used items are small
kpeter@80	232	/// integers. This map conforms the \ref concepts::ReferenceMap
kpeter@80	233	/// "ReferenceMap" concept.
kpeter@80	234	///
kpeter@80	235	/// The simplest way of using this map is through the rangeMap()
kpeter@80	236	/// function.
kpeter@80	237	template <typename V>
kpeter@80	238	class RangeMap : public MapBase<int, V> {
kpeter@80	239	template <typename V1>
kpeter@80	240	friend class RangeMap;
kpeter@80	241	private:
kpeter@80	242
kpeter@80	243	typedef std::vector<V> Vector;
kpeter@80	244	Vector _vector;
kpeter@80	245
alpar@25	246	public:
alpar@25	247
kpeter@80	248	typedef MapBase<int, V> Parent;
kpeter@80	249	/// Key type
kpeter@45	250	typedef typename Parent::Key Key;
kpeter@80	251	/// Value type
kpeter@45	252	typedef typename Parent::Value Value;
kpeter@80	253	/// Reference type
kpeter@80	254	typedef typename Vector::reference Reference;
kpeter@80	255	/// Const reference type
kpeter@80	256	typedef typename Vector::const_reference ConstReference;
kpeter@80	257
kpeter@80	258	typedef True ReferenceMapTag;
kpeter@80	259
kpeter@80	260	public:
kpeter@80	261
kpeter@80	262	/// Constructor with specified default value.
kpeter@80	263	RangeMap(int size = 0, const Value &value = Value())
kpeter@80	264	: _vector(size, value) {}
kpeter@80	265
kpeter@80	266	/// Constructs the map from an appropriate \c std::vector.
kpeter@80	267	template <typename V1>
kpeter@80	268	RangeMap(const std::vector<V1>& vector)
kpeter@80	269	: _vector(vector.begin(), vector.end()) {}
kpeter@80	270
kpeter@80	271	/// Constructs the map from another \ref RangeMap.
kpeter@80	272	template <typename V1>
kpeter@80	273	RangeMap(const RangeMap<V1> &c)
kpeter@80	274	: _vector(c._vector.begin(), c._vector.end()) {}
kpeter@80	275
kpeter@80	276	/// Returns the size of the map.
kpeter@80	277	int size() {
kpeter@80	278	return _vector.size();
kpeter@80	279	}
kpeter@80	280
kpeter@80	281	/// Resizes the map.
kpeter@80	282
kpeter@80	283	/// Resizes the underlying \c std::vector container, so changes the
kpeter@80	284	/// keyset of the map.
kpeter@80	285	/// \param size The new size of the map. The new keyset will be the
kpeter@80	286	/// range <tt>[0..size-1]</tt>.
kpeter@80	287	/// \param value The default value to assign to the new keys.
kpeter@80	288	void resize(int size, const Value &value = Value()) {
kpeter@80	289	_vector.resize(size, value);
kpeter@80	290	}
kpeter@80	291
kpeter@80	292	private:
kpeter@80	293
kpeter@80	294	RangeMap& operator=(const RangeMap&);
kpeter@80	295
kpeter@80	296	public:
kpeter@80	297
kpeter@80	298	///\e
kpeter@80	299	Reference operator[](const Key &k) {
kpeter@80	300	return _vector[k];
kpeter@80	301	}
kpeter@80	302
kpeter@80	303	///\e
kpeter@80	304	ConstReference operator[](const Key &k) const {
kpeter@80	305	return _vector[k];
kpeter@80	306	}
kpeter@80	307
kpeter@80	308	///\e
kpeter@80	309	void set(const Key &k, const Value &v) {
kpeter@80	310	_vector[k] = v;
kpeter@80	311	}
kpeter@80	312	};
kpeter@80	313
kpeter@80	314	/// Returns a \ref RangeMap class
kpeter@80	315
kpeter@80	316	/// This function just returns a \ref RangeMap class.
kpeter@80	317	/// \relates RangeMap
kpeter@80	318	template<typename V>
kpeter@80	319	inline RangeMap<V> rangeMap(int size = 0, const V &value = V()) {
kpeter@80	320	return RangeMap<V>(size, value);
kpeter@80	321	}
kpeter@80	322
kpeter@80	323	/// \brief Returns a \ref RangeMap class created from an appropriate
kpeter@80	324	/// \c std::vector
kpeter@80	325
kpeter@80	326	/// This function just returns a \ref RangeMap class created from an
kpeter@80	327	/// appropriate \c std::vector.
kpeter@80	328	/// \relates RangeMap
kpeter@80	329	template<typename V>
kpeter@80	330	inline RangeMap<V> rangeMap(const std::vector<V> &vector) {
kpeter@80	331	return RangeMap<V>(vector);
kpeter@80	332	}
kpeter@80	333
kpeter@80	334
kpeter@80	335	/// Map type based on \c std::map
kpeter@80	336
kpeter@80	337	/// This map is essentially a wrapper for \c std::map with addition
kpeter@80	338	/// that you can specify a default value for the keys that are not
kpeter@80	339	/// stored actually. This value can be different from the default
kpeter@80	340	/// contructed value (i.e. \c %Value()).
kpeter@80	341	/// This type conforms the \ref concepts::ReferenceMap "ReferenceMap"
kpeter@80	342	/// concept.
kpeter@80	343	///
kpeter@80	344	/// This map is useful if a default value should be assigned to most of
kpeter@80	345	/// the keys and different values should be assigned only to a few
kpeter@80	346	/// keys (i.e. the map is "sparse").
kpeter@80	347	/// The name of this type also refers to this important usage.
kpeter@80	348	///
kpeter@80	349	/// Apart form that this map can be used in many other cases since it
kpeter@80	350	/// is based on \c std::map, which is a general associative container.
kpeter@80	351	/// However keep in mind that it is usually not as efficient as other
kpeter@80	352	/// maps.
kpeter@80	353	///
kpeter@80	354	/// The simplest way of using this map is through the sparseMap()
kpeter@80	355	/// function.
kpeter@80	356	template <typename K, typename V, typename Compare = std::less<K> >
kpeter@80	357	class SparseMap : public MapBase<K, V> {
kpeter@80	358	template <typename K1, typename V1, typename C1>
kpeter@80	359	friend class SparseMap;
kpeter@80	360	public:
kpeter@80	361
kpeter@80	362	typedef MapBase<K, V> Parent;
kpeter@80	363	/// Key type
kpeter@80	364	typedef typename Parent::Key Key;
kpeter@80	365	/// Value type
kpeter@80	366	typedef typename Parent::Value Value;
kpeter@80	367	/// Reference type
kpeter@80	368	typedef Value& Reference;
kpeter@80	369	/// Const reference type
kpeter@80	370	typedef const Value& ConstReference;
alpar@25	371
kpeter@45	372	typedef True ReferenceMapTag;
kpeter@45	373
alpar@25	374	private:
kpeter@80	375
kpeter@80	376	typedef std::map<K, V, Compare> Map;
kpeter@80	377	Map _map;
alpar@25	378	Value _value;
alpar@25	379
alpar@25	380	public:
alpar@25	381
kpeter@80	382	/// \brief Constructor with specified default value.
kpeter@80	383	SparseMap(const Value &value = Value()) : _value(value) {}
kpeter@80	384	/// \brief Constructs the map from an appropriate \c std::map, and
kpeter@47	385	/// explicitly specifies a default value.
kpeter@80	386	template <typename V1, typename Comp1>
kpeter@80	387	SparseMap(const std::map<Key, V1, Comp1> &map,
kpeter@80	388	const Value &value = Value())
alpar@25	389	: _map(map.begin(), map.end()), _value(value) {}
kpeter@80	390
kpeter@80	391	/// \brief Constructs the map from another \ref SparseMap.
kpeter@80	392	template<typename V1, typename Comp1>
kpeter@80	393	SparseMap(const SparseMap<Key, V1, Comp1> &c)
alpar@25	394	: _map(c._map.begin(), c._map.end()), _value(c._value) {}
alpar@25	395
alpar@25	396	private:
alpar@25	397
kpeter@80	398	SparseMap& operator=(const SparseMap&);
alpar@25	399
alpar@25	400	public:
alpar@25	401
alpar@25	402	///\e
alpar@25	403	Reference operator[](const Key &k) {
alpar@25	404	typename Map::iterator it = _map.lower_bound(k);
alpar@25	405	if (it != _map.end() && !_map.key_comp()(k, it->first))
alpar@209	406	return it->second;
alpar@25	407	else
alpar@209	408	return _map.insert(it, std::make_pair(k, _value))->second;
alpar@25	409	}
alpar@25	410
kpeter@80	411	///\e
alpar@25	412	ConstReference operator[](const Key &k) const {
alpar@25	413	typename Map::const_iterator it = _map.find(k);
alpar@25	414	if (it != _map.end())
alpar@209	415	return it->second;
alpar@25	416	else
alpar@209	417	return _value;
alpar@25	418	}
alpar@25	419
kpeter@80	420	///\e
kpeter@80	421	void set(const Key &k, const Value &v) {
alpar@25	422	typename Map::iterator it = _map.lower_bound(k);
alpar@25	423	if (it != _map.end() && !_map.key_comp()(k, it->first))
alpar@209	424	it->second = v;
alpar@25	425	else
alpar@209	426	_map.insert(it, std::make_pair(k, v));
alpar@25	427	}
alpar@25	428
kpeter@80	429	///\e
kpeter@80	430	void setAll(const Value &v) {
kpeter@80	431	_value = v;
alpar@25	432	_map.clear();
kpeter@80	433	}
kpeter@80	434	};
alpar@25	435
kpeter@80	436	/// Returns a \ref SparseMap class
kpeter@45	437
kpeter@80	438	/// This function just returns a \ref SparseMap class with specified
kpeter@80	439	/// default value.
kpeter@80	440	/// \relates SparseMap
kpeter@80	441	template<typename K, typename V, typename Compare>
kpeter@80	442	inline SparseMap<K, V, Compare> sparseMap(const V& value = V()) {
kpeter@80	443	return SparseMap<K, V, Compare>(value);
kpeter@54	444	}
kpeter@45	445
kpeter@80	446	template<typename K, typename V>
kpeter@80	447	inline SparseMap<K, V, std::less<K> > sparseMap(const V& value = V()) {
kpeter@80	448	return SparseMap<K, V, std::less<K> >(value);
kpeter@45	449	}
alpar@25	450
kpeter@80	451	/// \brief Returns a \ref SparseMap class created from an appropriate
kpeter@80	452	/// \c std::map
alpar@25	453
kpeter@80	454	/// This function just returns a \ref SparseMap class created from an
kpeter@80	455	/// appropriate \c std::map.
kpeter@80	456	/// \relates SparseMap
kpeter@80	457	template<typename K, typename V, typename Compare>
kpeter@80	458	inline SparseMap<K, V, Compare>
kpeter@80	459	sparseMap(const std::map<K, V, Compare> &map, const V& value = V())
kpeter@80	460	{
kpeter@80	461	return SparseMap<K, V, Compare>(map, value);
kpeter@45	462	}
alpar@25	463
alpar@25	464	/// @}
alpar@25	465
alpar@25	466	/// \addtogroup map_adaptors
alpar@25	467	/// @{
alpar@25	468
kpeter@80	469	/// Composition of two maps
kpeter@80	470
kpeter@82	471	/// This \ref concepts::ReadMap "read-only map" returns the
kpeter@80	472	/// composition of two given maps. That is to say, if \c m1 is of
kpeter@80	473	/// type \c M1 and \c m2 is of \c M2, then for
kpeter@80	474	/// \code
kpeter@80	475	/// ComposeMap<M1, M2> cm(m1,m2);
kpeter@80	476	/// \endcode
kpeter@80	477	/// <tt>cm[x]</tt> will be equal to <tt>m1[m2[x]]</tt>.
alpar@25	478	///
kpeter@80	479	/// The \c Key type of the map is inherited from \c M2 and the
kpeter@80	480	/// \c Value type is from \c M1.
kpeter@80	481	/// \c M2::Value must be convertible to \c M1::Key.
kpeter@80	482	///
kpeter@80	483	/// The simplest way of using this map is through the composeMap()
kpeter@80	484	/// function.
kpeter@80	485	///
kpeter@80	486	/// \sa CombineMap
kpeter@80	487	///
kpeter@80	488	/// \todo Check the requirements.
kpeter@80	489	template <typename M1, typename M2>
kpeter@80	490	class ComposeMap : public MapBase<typename M2::Key, typename M1::Value> {
kpeter@80	491	const M1 &_m1;
kpeter@80	492	const M2 &_m2;
alpar@25	493	public:
kpeter@80	494	typedef MapBase<typename M2::Key, typename M1::Value> Parent;
alpar@25	495	typedef typename Parent::Key Key;
alpar@25	496	typedef typename Parent::Value Value;
alpar@25	497
kpeter@80	498	/// Constructor
kpeter@80	499	ComposeMap(const M1 &m1, const M2 &m2) : _m1(m1), _m2(m2) {}
kpeter@80	500
alpar@25	501	/// \e
kpeter@80	502	typename MapTraits<M1>::ConstReturnValue
kpeter@80	503	operator[](const Key &k) const { return _m1[_m2[k]]; }
alpar@25	504	};
alpar@25	505
kpeter@80	506	/// Returns a \ref ComposeMap class
alpar@25	507
kpeter@80	508	/// This function just returns a \ref ComposeMap class.
kpeter@80	509	///
kpeter@80	510	/// If \c m1 and \c m2 are maps and the \c Value type of \c m2 is
kpeter@80	511	/// convertible to the \c Key of \c m1, then <tt>composeMap(m1,m2)[x]</tt>
kpeter@80	512	/// will be equal to <tt>m1[m2[x]]</tt>.
kpeter@80	513	///
kpeter@80	514	/// \relates ComposeMap
kpeter@80	515	template <typename M1, typename M2>
kpeter@80	516	inline ComposeMap<M1, M2> composeMap(const M1 &m1, const M2 &m2) {
kpeter@80	517	return ComposeMap<M1, M2>(m1, m2);
alpar@25	518	}
alpar@25	519
kpeter@80	520
kpeter@80	521	/// Combination of two maps using an STL (binary) functor.
kpeter@80	522
kpeter@82	523	/// This \ref concepts::ReadMap "read-only map" takes two maps and a
kpeter@80	524	/// binary functor and returns the combination of the two given maps
kpeter@80	525	/// using the functor.
kpeter@80	526	/// That is to say, if \c m1 is of type \c M1 and \c m2 is of \c M2
kpeter@80	527	/// and \c f is of \c F, then for
kpeter@80	528	/// \code
kpeter@80	529	/// CombineMap<M1,M2,F,V> cm(m1,m2,f);
kpeter@80	530	/// \endcode
kpeter@80	531	/// <tt>cm[x]</tt> will be equal to <tt>f(m1[x],m2[x])</tt>.
alpar@26	532	///
kpeter@80	533	/// The \c Key type of the map is inherited from \c M1 (\c M1::Key
kpeter@80	534	/// must be convertible to \c M2::Key) and the \c Value type is \c V.
kpeter@80	535	/// \c M2::Value and \c M1::Value must be convertible to the
kpeter@80	536	/// corresponding input parameter of \c F and the return type of \c F
kpeter@80	537	/// must be convertible to \c V.
kpeter@80	538	///
kpeter@80	539	/// The simplest way of using this map is through the combineMap()
kpeter@80	540	/// function.
kpeter@80	541	///
kpeter@80	542	/// \sa ComposeMap
kpeter@80	543	///
kpeter@80	544	/// \todo Check the requirements.
kpeter@80	545	template<typename M1, typename M2, typename F,
alpar@209	546	typename V = typename F::result_type>
kpeter@80	547	class CombineMap : public MapBase<typename M1::Key, V> {
kpeter@80	548	const M1 &_m1;
kpeter@80	549	const M2 &_m2;
kpeter@80	550	F _f;
alpar@25	551	public:
kpeter@80	552	typedef MapBase<typename M1::Key, V> Parent;
alpar@25	553	typedef typename Parent::Key Key;
alpar@25	554	typedef typename Parent::Value Value;
alpar@25	555
kpeter@80	556	/// Constructor
kpeter@80	557	CombineMap(const M1 &m1, const M2 &m2, const F &f = F())
kpeter@80	558	: _m1(m1), _m2(m2), _f(f) {}
kpeter@80	559	/// \e
kpeter@80	560	Value operator[](const Key &k) const { return _f(_m1[k],_m2[k]); }
kpeter@80	561	};
alpar@25	562
kpeter@80	563	/// Returns a \ref CombineMap class
alpar@25	564
kpeter@80	565	/// This function just returns a \ref CombineMap class.
kpeter@80	566	///
kpeter@80	567	/// For example, if \c m1 and \c m2 are both maps with \c double
kpeter@80	568	/// values, then
kpeter@80	569	/// \code
kpeter@80	570	/// combineMap(m1,m2,std::plus<double>())
kpeter@80	571	/// \endcode
kpeter@80	572	/// is equivalent to
kpeter@80	573	/// \code
kpeter@80	574	/// addMap(m1,m2)
kpeter@80	575	/// \endcode
kpeter@80	576	///
kpeter@80	577	/// This function is specialized for adaptable binary function
kpeter@80	578	/// classes and C++ functions.
kpeter@80	579	///
kpeter@80	580	/// \relates CombineMap
kpeter@80	581	template<typename M1, typename M2, typename F, typename V>
kpeter@80	582	inline CombineMap<M1, M2, F, V>
kpeter@80	583	combineMap(const M1 &m1, const M2 &m2, const F &f) {
kpeter@80	584	return CombineMap<M1, M2, F, V>(m1,m2,f);
alpar@25	585	}
alpar@25	586
kpeter@80	587	template<typename M1, typename M2, typename F>
kpeter@80	588	inline CombineMap<M1, M2, F, typename F::result_type>
kpeter@80	589	combineMap(const M1 &m1, const M2 &m2, const F &f) {
kpeter@80	590	return combineMap<M1, M2, F, typename F::result_type>(m1,m2,f);
kpeter@80	591	}
alpar@25	592
kpeter@80	593	template<typename M1, typename M2, typename K1, typename K2, typename V>
kpeter@80	594	inline CombineMap<M1, M2, V (*)(K1, K2), V>
kpeter@80	595	combineMap(const M1 &m1, const M2 &m2, V (*f)(K1, K2)) {
kpeter@80	596	return combineMap<M1, M2, V (*)(K1, K2), V>(m1,m2,f);
kpeter@80	597	}
kpeter@80	598
kpeter@80	599
kpeter@80	600	/// Converts an STL style (unary) functor to a map
kpeter@80	601
kpeter@82	602	/// This \ref concepts::ReadMap "read-only map" returns the value
kpeter@80	603	/// of a given functor. Actually, it just wraps the functor and
kpeter@80	604	/// provides the \c Key and \c Value typedefs.
alpar@26	605	///
kpeter@80	606	/// Template parameters \c K and \c V will become its \c Key and
kpeter@80	607	/// \c Value. In most cases they have to be given explicitly because
kpeter@80	608	/// a functor typically does not provide \c argument_type and
kpeter@80	609	/// \c result_type typedefs.
kpeter@80	610	/// Parameter \c F is the type of the used functor.
kpeter@29	611	///
kpeter@80	612	/// The simplest way of using this map is through the functorToMap()
kpeter@80	613	/// function.
kpeter@80	614	///
kpeter@80	615	/// \sa MapToFunctor
kpeter@80	616	template<typename F,
alpar@209	617	typename K = typename F::argument_type,
alpar@209	618	typename V = typename F::result_type>
kpeter@80	619	class FunctorToMap : public MapBase<K, V> {
kpeter@123	620	F _f;
kpeter@80	621	public:
kpeter@80	622	typedef MapBase<K, V> Parent;
kpeter@80	623	typedef typename Parent::Key Key;
kpeter@80	624	typedef typename Parent::Value Value;
alpar@25	625
kpeter@80	626	/// Constructor
kpeter@80	627	FunctorToMap(const F &f = F()) : _f(f) {}
kpeter@80	628	/// \e
kpeter@80	629	Value operator[](const Key &k) const { return _f(k); }
kpeter@80	630	};
kpeter@80	631
kpeter@80	632	/// Returns a \ref FunctorToMap class
kpeter@80	633
kpeter@80	634	/// This function just returns a \ref FunctorToMap class.
kpeter@80	635	///
kpeter@80	636	/// This function is specialized for adaptable binary function
kpeter@80	637	/// classes and C++ functions.
kpeter@80	638	///
kpeter@80	639	/// \relates FunctorToMap
kpeter@80	640	template<typename K, typename V, typename F>
kpeter@80	641	inline FunctorToMap<F, K, V> functorToMap(const F &f) {
kpeter@80	642	return FunctorToMap<F, K, V>(f);
kpeter@80	643	}
kpeter@80	644
kpeter@80	645	template <typename F>
kpeter@80	646	inline FunctorToMap<F, typename F::argument_type, typename F::result_type>
kpeter@80	647	functorToMap(const F &f)
kpeter@80	648	{
kpeter@80	649	return FunctorToMap<F, typename F::argument_type,
kpeter@80	650	typename F::result_type>(f);
kpeter@80	651	}
kpeter@80	652
kpeter@80	653	template <typename K, typename V>
kpeter@80	654	inline FunctorToMap<V ()(K), K, V> functorToMap(V (f)(K)) {
kpeter@80	655	return FunctorToMap<V (*)(K), K, V>(f);
kpeter@80	656	}
kpeter@80	657
kpeter@80	658
kpeter@80	659	/// Converts a map to an STL style (unary) functor
kpeter@80	660
kpeter@80	661	/// This class converts a map to an STL style (unary) functor.
kpeter@80	662	/// That is it provides an <tt>operator()</tt> to read its values.
kpeter@80	663	///
kpeter@80	664	/// For the sake of convenience it also works as a usual
kpeter@80	665	/// \ref concepts::ReadMap "readable map", i.e. <tt>operator[]</tt>
kpeter@80	666	/// and the \c Key and \c Value typedefs also exist.
kpeter@80	667	///
kpeter@80	668	/// The simplest way of using this map is through the mapToFunctor()
kpeter@80	669	/// function.
kpeter@80	670	///
kpeter@80	671	///\sa FunctorToMap
kpeter@80	672	template <typename M>
kpeter@80	673	class MapToFunctor : public MapBase<typename M::Key, typename M::Value> {
kpeter@80	674	const M &_m;
alpar@25	675	public:
alpar@25	676	typedef MapBase<typename M::Key, typename M::Value> Parent;
alpar@25	677	typedef typename Parent::Key Key;
alpar@25	678	typedef typename Parent::Value Value;
alpar@25	679
kpeter@80	680	typedef typename Parent::Key argument_type;
kpeter@80	681	typedef typename Parent::Value result_type;
kpeter@80	682
kpeter@80	683	/// Constructor
kpeter@80	684	MapToFunctor(const M &m) : _m(m) {}
kpeter@80	685	/// \e
kpeter@80	686	Value operator()(const Key &k) const { return _m[k]; }
kpeter@80	687	/// \e
kpeter@80	688	Value operator[](const Key &k) const { return _m[k]; }
alpar@25	689	};
kpeter@45	690
kpeter@80	691	/// Returns a \ref MapToFunctor class
kpeter@80	692
kpeter@80	693	/// This function just returns a \ref MapToFunctor class.
kpeter@80	694	/// \relates MapToFunctor
kpeter@45	695	template<typename M>
kpeter@80	696	inline MapToFunctor<M> mapToFunctor(const M &m) {
kpeter@80	697	return MapToFunctor<M>(m);
kpeter@45	698	}
alpar@25	699
alpar@25	700
kpeter@80	701	/// \brief Map adaptor to convert the \c Value type of a map to
kpeter@80	702	/// another type using the default conversion.
kpeter@80	703
kpeter@80	704	/// Map adaptor to convert the \c Value type of a \ref concepts::ReadMap
kpeter@80	705	/// "readable map" to another type using the default conversion.
kpeter@80	706	/// The \c Key type of it is inherited from \c M and the \c Value
kpeter@80	707	/// type is \c V.
kpeter@80	708	/// This type conforms the \ref concepts::ReadMap "ReadMap" concept.
alpar@26	709	///
kpeter@80	710	/// The simplest way of using this map is through the convertMap()
kpeter@80	711	/// function.
kpeter@80	712	template <typename M, typename V>
kpeter@80	713	class ConvertMap : public MapBase<typename M::Key, V> {
kpeter@80	714	const M &_m;
kpeter@80	715	public:
kpeter@80	716	typedef MapBase<typename M::Key, V> Parent;
kpeter@80	717	typedef typename Parent::Key Key;
kpeter@80	718	typedef typename Parent::Value Value;
kpeter@80	719
kpeter@80	720	/// Constructor
kpeter@80	721
kpeter@80	722	/// Constructor.
kpeter@80	723	/// \param m The underlying map.
kpeter@80	724	ConvertMap(const M &m) : _m(m) {}
kpeter@80	725
kpeter@80	726	/// \e
kpeter@80	727	Value operator[](const Key &k) const { return _m[k]; }
kpeter@80	728	};
kpeter@80	729
kpeter@80	730	/// Returns a \ref ConvertMap class
kpeter@80	731
kpeter@80	732	/// This function just returns a \ref ConvertMap class.
kpeter@80	733	/// \relates ConvertMap
kpeter@80	734	template<typename V, typename M>
kpeter@80	735	inline ConvertMap<M, V> convertMap(const M &map) {
kpeter@80	736	return ConvertMap<M, V>(map);
kpeter@80	737	}
kpeter@80	738
kpeter@80	739
kpeter@80	740	/// Applies all map setting operations to two maps
kpeter@80	741
kpeter@80	742	/// This map has two \ref concepts::WriteMap "writable map" parameters
kpeter@80	743	/// and each write request will be passed to both of them.
kpeter@80	744	/// If \c M1 is also \ref concepts::ReadMap "readable", then the read
kpeter@80	745	/// operations will return the corresponding values of \c M1.
kpeter@29	746	///
kpeter@80	747	/// The \c Key and \c Value types are inherited from \c M1.
kpeter@80	748	/// The \c Key and \c Value of \c M2 must be convertible from those
kpeter@80	749	/// of \c M1.
kpeter@80	750	///
kpeter@80	751	/// The simplest way of using this map is through the forkMap()
kpeter@80	752	/// function.
kpeter@80	753	template<typename M1, typename M2>
kpeter@80	754	class ForkMap : public MapBase<typename M1::Key, typename M1::Value> {
kpeter@80	755	M1 &_m1;
kpeter@80	756	M2 &_m2;
kpeter@80	757	public:
kpeter@80	758	typedef MapBase<typename M1::Key, typename M1::Value> Parent;
kpeter@80	759	typedef typename Parent::Key Key;
kpeter@80	760	typedef typename Parent::Value Value;
alpar@25	761
kpeter@80	762	/// Constructor
kpeter@80	763	ForkMap(M1 &m1, M2 &m2) : _m1(m1), _m2(m2) {}
kpeter@80	764	/// Returns the value associated with the given key in the first map.
kpeter@80	765	Value operator[](const Key &k) const { return _m1[k]; }
kpeter@80	766	/// Sets the value associated with the given key in both maps.
kpeter@80	767	void set(const Key &k, const Value &v) { _m1.set(k,v); _m2.set(k,v); }
kpeter@80	768	};
kpeter@80	769
kpeter@80	770	/// Returns a \ref ForkMap class
kpeter@80	771
kpeter@80	772	/// This function just returns a \ref ForkMap class.
kpeter@80	773	/// \relates ForkMap
kpeter@80	774	template <typename M1, typename M2>
kpeter@80	775	inline ForkMap<M1,M2> forkMap(M1 &m1, M2 &m2) {
kpeter@80	776	return ForkMap<M1,M2>(m1,m2);
kpeter@80	777	}
kpeter@80	778
kpeter@80	779
kpeter@80	780	/// Sum of two maps
kpeter@80	781
kpeter@82	782	/// This \ref concepts::ReadMap "read-only map" returns the sum
kpeter@80	783	/// of the values of the two given maps.
kpeter@80	784	/// Its \c Key and \c Value types are inherited from \c M1.
kpeter@80	785	/// The \c Key and \c Value of \c M2 must be convertible to those of
kpeter@80	786	/// \c M1.
kpeter@80	787	///
kpeter@80	788	/// If \c m1 is of type \c M1 and \c m2 is of \c M2, then for
kpeter@80	789	/// \code
kpeter@80	790	/// AddMap<M1,M2> am(m1,m2);
kpeter@80	791	/// \endcode
kpeter@80	792	/// <tt>am[x]</tt> will be equal to <tt>m1[x]+m2[x]</tt>.
kpeter@80	793	///
kpeter@80	794	/// The simplest way of using this map is through the addMap()
kpeter@80	795	/// function.
kpeter@80	796	///
kpeter@80	797	/// \sa SubMap, MulMap, DivMap
kpeter@80	798	/// \sa ShiftMap, ShiftWriteMap
kpeter@80	799	template<typename M1, typename M2>
alpar@25	800	class AddMap : public MapBase<typename M1::Key, typename M1::Value> {
kpeter@80	801	const M1 &_m1;
kpeter@80	802	const M2 &_m2;
alpar@25	803	public:
alpar@25	804	typedef MapBase<typename M1::Key, typename M1::Value> Parent;
alpar@25	805	typedef typename Parent::Key Key;
alpar@25	806	typedef typename Parent::Value Value;
alpar@25	807
kpeter@80	808	/// Constructor
kpeter@80	809	AddMap(const M1 &m1, const M2 &m2) : _m1(m1), _m2(m2) {}
kpeter@80	810	/// \e
kpeter@80	811	Value operator[](const Key &k) const { return _m1[k]+_m2[k]; }
alpar@25	812	};
alpar@25	813
kpeter@80	814	/// Returns an \ref AddMap class
kpeter@80	815
kpeter@80	816	/// This function just returns an \ref AddMap class.
alpar@25	817	///
kpeter@80	818	/// For example, if \c m1 and \c m2 are both maps with \c double
kpeter@80	819	/// values, then <tt>addMap(m1,m2)[x]</tt> will be equal to
kpeter@80	820	/// <tt>m1[x]+m2[x]</tt>.
kpeter@80	821	///
kpeter@80	822	/// \relates AddMap
kpeter@80	823	template<typename M1, typename M2>
kpeter@80	824	inline AddMap<M1, M2> addMap(const M1 &m1, const M2 &m2) {
alpar@25	825	return AddMap<M1, M2>(m1,m2);
alpar@25	826	}
alpar@25	827
alpar@25	828
kpeter@80	829	/// Difference of two maps
kpeter@80	830
kpeter@82	831	/// This \ref concepts::ReadMap "read-only map" returns the difference
kpeter@80	832	/// of the values of the two given maps.
kpeter@80	833	/// Its \c Key and \c Value types are inherited from \c M1.
kpeter@80	834	/// The \c Key and \c Value of \c M2 must be convertible to those of
kpeter@80	835	/// \c M1.
alpar@25	836	///
kpeter@80	837	/// If \c m1 is of type \c M1 and \c m2 is of \c M2, then for
kpeter@80	838	/// \code
kpeter@80	839	/// SubMap<M1,M2> sm(m1,m2);
kpeter@80	840	/// \endcode
kpeter@80	841	/// <tt>sm[x]</tt> will be equal to <tt>m1[x]-m2[x]</tt>.
kpeter@29	842	///
kpeter@80	843	/// The simplest way of using this map is through the subMap()
kpeter@80	844	/// function.
kpeter@80	845	///
kpeter@80	846	/// \sa AddMap, MulMap, DivMap
kpeter@80	847	template<typename M1, typename M2>
kpeter@80	848	class SubMap : public MapBase<typename M1::Key, typename M1::Value> {
kpeter@80	849	const M1 &_m1;
kpeter@80	850	const M2 &_m2;
kpeter@80	851	public:
kpeter@80	852	typedef MapBase<typename M1::Key, typename M1::Value> Parent;
kpeter@80	853	typedef typename Parent::Key Key;
kpeter@80	854	typedef typename Parent::Value Value;
kpeter@80	855
kpeter@80	856	/// Constructor
kpeter@80	857	SubMap(const M1 &m1, const M2 &m2) : _m1(m1), _m2(m2) {}
kpeter@80	858	/// \e
kpeter@80	859	Value operator[](const Key &k) const { return _m1[k]-_m2[k]; }
kpeter@80	860	};
kpeter@80	861
kpeter@80	862	/// Returns a \ref SubMap class
kpeter@80	863
kpeter@80	864	/// This function just returns a \ref SubMap class.
kpeter@80	865	///
kpeter@80	866	/// For example, if \c m1 and \c m2 are both maps with \c double
kpeter@80	867	/// values, then <tt>subMap(m1,m2)[x]</tt> will be equal to
kpeter@80	868	/// <tt>m1[x]-m2[x]</tt>.
kpeter@80	869	///
kpeter@80	870	/// \relates SubMap
kpeter@80	871	template<typename M1, typename M2>
kpeter@80	872	inline SubMap<M1, M2> subMap(const M1 &m1, const M2 &m2) {
kpeter@80	873	return SubMap<M1, M2>(m1,m2);
kpeter@80	874	}
kpeter@80	875
kpeter@80	876
kpeter@80	877	/// Product of two maps
kpeter@80	878
kpeter@82	879	/// This \ref concepts::ReadMap "read-only map" returns the product
kpeter@80	880	/// of the values of the two given maps.
kpeter@80	881	/// Its \c Key and \c Value types are inherited from \c M1.
kpeter@80	882	/// The \c Key and \c Value of \c M2 must be convertible to those of
kpeter@80	883	/// \c M1.
kpeter@80	884	///
kpeter@80	885	/// If \c m1 is of type \c M1 and \c m2 is of \c M2, then for
kpeter@80	886	/// \code
kpeter@80	887	/// MulMap<M1,M2> mm(m1,m2);
kpeter@80	888	/// \endcode
kpeter@80	889	/// <tt>mm[x]</tt> will be equal to <tt>m1[x]*m2[x]</tt>.
kpeter@80	890	///
kpeter@80	891	/// The simplest way of using this map is through the mulMap()
kpeter@80	892	/// function.
kpeter@80	893	///
kpeter@80	894	/// \sa AddMap, SubMap, DivMap
kpeter@80	895	/// \sa ScaleMap, ScaleWriteMap
kpeter@80	896	template<typename M1, typename M2>
kpeter@80	897	class MulMap : public MapBase<typename M1::Key, typename M1::Value> {
kpeter@80	898	const M1 &_m1;
kpeter@80	899	const M2 &_m2;
kpeter@80	900	public:
kpeter@80	901	typedef MapBase<typename M1::Key, typename M1::Value> Parent;
kpeter@80	902	typedef typename Parent::Key Key;
kpeter@80	903	typedef typename Parent::Value Value;
kpeter@80	904
kpeter@80	905	/// Constructor
kpeter@80	906	MulMap(const M1 &m1,const M2 &m2) : _m1(m1), _m2(m2) {}
kpeter@80	907	/// \e
kpeter@80	908	Value operator[](const Key &k) const { return _m1[k]*_m2[k]; }
kpeter@80	909	};
kpeter@80	910
kpeter@80	911	/// Returns a \ref MulMap class
kpeter@80	912
kpeter@80	913	/// This function just returns a \ref MulMap class.
kpeter@80	914	///
kpeter@80	915	/// For example, if \c m1 and \c m2 are both maps with \c double
kpeter@80	916	/// values, then <tt>mulMap(m1,m2)[x]</tt> will be equal to
kpeter@80	917	/// <tt>m1[x]*m2[x]</tt>.
kpeter@80	918	///
kpeter@80	919	/// \relates MulMap
kpeter@80	920	template<typename M1, typename M2>
kpeter@80	921	inline MulMap<M1, M2> mulMap(const M1 &m1,const M2 &m2) {
kpeter@80	922	return MulMap<M1, M2>(m1,m2);
kpeter@80	923	}
kpeter@80	924
kpeter@80	925
kpeter@80	926	/// Quotient of two maps
kpeter@80	927
kpeter@82	928	/// This \ref concepts::ReadMap "read-only map" returns the quotient
kpeter@80	929	/// of the values of the two given maps.
kpeter@80	930	/// Its \c Key and \c Value types are inherited from \c M1.
kpeter@80	931	/// The \c Key and \c Value of \c M2 must be convertible to those of
kpeter@80	932	/// \c M1.
kpeter@80	933	///
kpeter@80	934	/// If \c m1 is of type \c M1 and \c m2 is of \c M2, then for
kpeter@80	935	/// \code
kpeter@80	936	/// DivMap<M1,M2> dm(m1,m2);
kpeter@80	937	/// \endcode
kpeter@80	938	/// <tt>dm[x]</tt> will be equal to <tt>m1[x]/m2[x]</tt>.
kpeter@80	939	///
kpeter@80	940	/// The simplest way of using this map is through the divMap()
kpeter@80	941	/// function.
kpeter@80	942	///
kpeter@80	943	/// \sa AddMap, SubMap, MulMap
kpeter@80	944	template<typename M1, typename M2>
kpeter@80	945	class DivMap : public MapBase<typename M1::Key, typename M1::Value> {
kpeter@80	946	const M1 &_m1;
kpeter@80	947	const M2 &_m2;
kpeter@80	948	public:
kpeter@80	949	typedef MapBase<typename M1::Key, typename M1::Value> Parent;
kpeter@80	950	typedef typename Parent::Key Key;
kpeter@80	951	typedef typename Parent::Value Value;
kpeter@80	952
kpeter@80	953	/// Constructor
kpeter@80	954	DivMap(const M1 &m1,const M2 &m2) : _m1(m1), _m2(m2) {}
kpeter@80	955	/// \e
kpeter@80	956	Value operator[](const Key &k) const { return _m1[k]/_m2[k]; }
kpeter@80	957	};
kpeter@80	958
kpeter@80	959	/// Returns a \ref DivMap class
kpeter@80	960
kpeter@80	961	/// This function just returns a \ref DivMap class.
kpeter@80	962	///
kpeter@80	963	/// For example, if \c m1 and \c m2 are both maps with \c double
kpeter@80	964	/// values, then <tt>divMap(m1,m2)[x]</tt> will be equal to
kpeter@80	965	/// <tt>m1[x]/m2[x]</tt>.
kpeter@80	966	///
kpeter@80	967	/// \relates DivMap
kpeter@80	968	template<typename M1, typename M2>
kpeter@80	969	inline DivMap<M1, M2> divMap(const M1 &m1,const M2 &m2) {
kpeter@80	970	return DivMap<M1, M2>(m1,m2);
kpeter@80	971	}
kpeter@80	972
kpeter@80	973
kpeter@80	974	/// Shifts a map with a constant.
kpeter@80	975
kpeter@82	976	/// This \ref concepts::ReadMap "read-only map" returns the sum of
kpeter@80	977	/// the given map and a constant value (i.e. it shifts the map with
kpeter@80	978	/// the constant). Its \c Key and \c Value are inherited from \c M.
kpeter@80	979	///
kpeter@80	980	/// Actually,
kpeter@80	981	/// \code
kpeter@80	982	/// ShiftMap<M> sh(m,v);
kpeter@80	983	/// \endcode
kpeter@80	984	/// is equivalent to
kpeter@80	985	/// \code
kpeter@80	986	/// ConstMap<M::Key, M::Value> cm(v);
kpeter@80	987	/// AddMap<M, ConstMap<M::Key, M::Value> > sh(m,cm);
kpeter@80	988	/// \endcode
kpeter@80	989	///
kpeter@80	990	/// The simplest way of using this map is through the shiftMap()
kpeter@80	991	/// function.
kpeter@80	992	///
kpeter@80	993	/// \sa ShiftWriteMap
kpeter@80	994	template<typename M, typename C = typename M::Value>
alpar@25	995	class ShiftMap : public MapBase<typename M::Key, typename M::Value> {
kpeter@80	996	const M &_m;
kpeter@80	997	C _v;
alpar@25	998	public:
alpar@25	999	typedef MapBase<typename M::Key, typename M::Value> Parent;
alpar@25	1000	typedef typename Parent::Key Key;
alpar@25	1001	typedef typename Parent::Value Value;
alpar@25	1002
kpeter@80	1003	/// Constructor
alpar@25	1004
kpeter@80	1005	/// Constructor.
kpeter@80	1006	/// \param m The undelying map.
kpeter@80	1007	/// \param v The constant value.
kpeter@80	1008	ShiftMap(const M &m, const C &v) : _m(m), _v(v) {}
kpeter@80	1009	/// \e
kpeter@80	1010	Value operator[](const Key &k) const { return _m[k]+_v; }
alpar@25	1011	};
alpar@25	1012
kpeter@80	1013	/// Shifts a map with a constant (read-write version).
alpar@25	1014
kpeter@80	1015	/// This \ref concepts::ReadWriteMap "read-write map" returns the sum
kpeter@80	1016	/// of the given map and a constant value (i.e. it shifts the map with
kpeter@80	1017	/// the constant). Its \c Key and \c Value are inherited from \c M.
kpeter@80	1018	/// It makes also possible to write the map.
alpar@25	1019	///
kpeter@80	1020	/// The simplest way of using this map is through the shiftWriteMap()
kpeter@80	1021	/// function.
kpeter@80	1022	///
kpeter@80	1023	/// \sa ShiftMap
kpeter@80	1024	template<typename M, typename C = typename M::Value>
alpar@25	1025	class ShiftWriteMap : public MapBase<typename M::Key, typename M::Value> {
kpeter@80	1026	M &_m;
kpeter@80	1027	C _v;
alpar@25	1028	public:
alpar@25	1029	typedef MapBase<typename M::Key, typename M::Value> Parent;
alpar@25	1030	typedef typename Parent::Key Key;
alpar@25	1031	typedef typename Parent::Value Value;
alpar@25	1032
kpeter@80	1033	/// Constructor
alpar@25	1034
kpeter@80	1035	/// Constructor.
kpeter@80	1036	/// \param m The undelying map.
kpeter@80	1037	/// \param v The constant value.
kpeter@80	1038	ShiftWriteMap(M &m, const C &v) : _m(m), _v(v) {}
alpar@25	1039	/// \e
kpeter@80	1040	Value operator[](const Key &k) const { return _m[k]+_v; }
alpar@25	1041	/// \e
kpeter@80	1042	void set(const Key &k, const Value &v) { _m.set(k, v-_v); }
alpar@25	1043	};
alpar@25	1044
kpeter@80	1045	/// Returns a \ref ShiftMap class
kpeter@80	1046
kpeter@80	1047	/// This function just returns a \ref ShiftMap class.
kpeter@80	1048	///
kpeter@80	1049	/// For example, if \c m is a map with \c double values and \c v is
kpeter@80	1050	/// \c double, then <tt>shiftMap(m,v)[x]</tt> will be equal to
kpeter@80	1051	/// <tt>m[x]+v</tt>.
kpeter@80	1052	///
kpeter@80	1053	/// \relates ShiftMap
kpeter@80	1054	template<typename M, typename C>
kpeter@80	1055	inline ShiftMap<M, C> shiftMap(const M &m, const C &v) {
alpar@25	1056	return ShiftMap<M, C>(m,v);
alpar@25	1057	}
alpar@25	1058
kpeter@80	1059	/// Returns a \ref ShiftWriteMap class
kpeter@29	1060
kpeter@80	1061	/// This function just returns a \ref ShiftWriteMap class.
kpeter@80	1062	///
kpeter@80	1063	/// For example, if \c m is a map with \c double values and \c v is
kpeter@80	1064	/// \c double, then <tt>shiftWriteMap(m,v)[x]</tt> will be equal to
kpeter@80	1065	/// <tt>m[x]+v</tt>.
kpeter@80	1066	/// Moreover it makes also possible to write the map.
kpeter@80	1067	///
kpeter@80	1068	/// \relates ShiftWriteMap
kpeter@80	1069	template<typename M, typename C>
kpeter@80	1070	inline ShiftWriteMap<M, C> shiftWriteMap(M &m, const C &v) {
alpar@25	1071	return ShiftWriteMap<M, C>(m,v);
alpar@25	1072	}
alpar@25	1073
alpar@25	1074
kpeter@80	1075	/// Scales a map with a constant.
kpeter@80	1076
kpeter@82	1077	/// This \ref concepts::ReadMap "read-only map" returns the value of
kpeter@80	1078	/// the given map multiplied from the left side with a constant value.
kpeter@80	1079	/// Its \c Key and \c Value are inherited from \c M.
alpar@26	1080	///
kpeter@80	1081	/// Actually,
kpeter@80	1082	/// \code
kpeter@80	1083	/// ScaleMap<M> sc(m,v);
kpeter@80	1084	/// \endcode
kpeter@80	1085	/// is equivalent to
kpeter@80	1086	/// \code
kpeter@80	1087	/// ConstMap<M::Key, M::Value> cm(v);
kpeter@80	1088	/// MulMap<ConstMap<M::Key, M::Value>, M> sc(cm,m);
kpeter@80	1089	/// \endcode
alpar@25	1090	///
kpeter@80	1091	/// The simplest way of using this map is through the scaleMap()
kpeter@80	1092	/// function.
alpar@25	1093	///
kpeter@80	1094	/// \sa ScaleWriteMap
kpeter@80	1095	template<typename M, typename C = typename M::Value>
alpar@25	1096	class ScaleMap : public MapBase<typename M::Key, typename M::Value> {
kpeter@80	1097	const M &_m;
kpeter@80	1098	C _v;
alpar@25	1099	public:
alpar@25	1100	typedef MapBase<typename M::Key, typename M::Value> Parent;
alpar@25	1101	typedef typename Parent::Key Key;
alpar@25	1102	typedef typename Parent::Value Value;
alpar@25	1103
kpeter@80	1104	/// Constructor
alpar@25	1105
kpeter@80	1106	/// Constructor.
kpeter@80	1107	/// \param m The undelying map.
kpeter@80	1108	/// \param v The constant value.
kpeter@80	1109	ScaleMap(const M &m, const C &v) : _m(m), _v(v) {}
alpar@25	1110	/// \e
kpeter@80	1111	Value operator[](const Key &k) const { return _v*_m[k]; }
alpar@25	1112	};
alpar@25	1113
kpeter@80	1114	/// Scales a map with a constant (read-write version).
alpar@25	1115
kpeter@80	1116	/// This \ref concepts::ReadWriteMap "read-write map" returns the value of
kpeter@80	1117	/// the given map multiplied from the left side with a constant value.
kpeter@80	1118	/// Its \c Key and \c Value are inherited from \c M.
kpeter@80	1119	/// It can also be used as write map if the \c / operator is defined
kpeter@80	1120	/// between \c Value and \c C and the given multiplier is not zero.
kpeter@29	1121	///
kpeter@80	1122	/// The simplest way of using this map is through the scaleWriteMap()
kpeter@80	1123	/// function.
kpeter@80	1124	///
kpeter@80	1125	/// \sa ScaleMap
kpeter@80	1126	template<typename M, typename C = typename M::Value>
alpar@25	1127	class ScaleWriteMap : public MapBase<typename M::Key, typename M::Value> {
kpeter@80	1128	M &_m;
kpeter@80	1129	C _v;
alpar@25	1130	public:
alpar@25	1131	typedef MapBase<typename M::Key, typename M::Value> Parent;
alpar@25	1132	typedef typename Parent::Key Key;
alpar@25	1133	typedef typename Parent::Value Value;
alpar@25	1134
kpeter@80	1135	/// Constructor
alpar@25	1136
kpeter@80	1137	/// Constructor.
kpeter@80	1138	/// \param m The undelying map.
kpeter@80	1139	/// \param v The constant value.
kpeter@80	1140	ScaleWriteMap(M &m, const C &v) : _m(m), _v(v) {}
alpar@25	1141	/// \e
kpeter@80	1142	Value operator[](const Key &k) const { return _v*_m[k]; }
alpar@25	1143	/// \e
kpeter@80	1144	void set(const Key &k, const Value &v) { _m.set(k, v/_v); }
alpar@25	1145	};
alpar@25	1146
kpeter@80	1147	/// Returns a \ref ScaleMap class
kpeter@80	1148
kpeter@80	1149	/// This function just returns a \ref ScaleMap class.
kpeter@80	1150	///
kpeter@80	1151	/// For example, if \c m is a map with \c double values and \c v is
kpeter@80	1152	/// \c double, then <tt>scaleMap(m,v)[x]</tt> will be equal to
kpeter@80	1153	/// <tt>v*m[x]</tt>.
kpeter@80	1154	///
kpeter@80	1155	/// \relates ScaleMap
kpeter@80	1156	template<typename M, typename C>
kpeter@80	1157	inline ScaleMap<M, C> scaleMap(const M &m, const C &v) {
alpar@25	1158	return ScaleMap<M, C>(m,v);
alpar@25	1159	}
alpar@25	1160
kpeter@80	1161	/// Returns a \ref ScaleWriteMap class
kpeter@29	1162
kpeter@80	1163	/// This function just returns a \ref ScaleWriteMap class.
kpeter@80	1164	///
kpeter@80	1165	/// For example, if \c m is a map with \c double values and \c v is
kpeter@80	1166	/// \c double, then <tt>scaleWriteMap(m,v)[x]</tt> will be equal to
kpeter@80	1167	/// <tt>v*m[x]</tt>.
kpeter@80	1168	/// Moreover it makes also possible to write the map.
kpeter@80	1169	///
kpeter@80	1170	/// \relates ScaleWriteMap
kpeter@80	1171	template<typename M, typename C>
kpeter@80	1172	inline ScaleWriteMap<M, C> scaleWriteMap(M &m, const C &v) {
alpar@25	1173	return ScaleWriteMap<M, C>(m,v);
alpar@25	1174	}
alpar@25	1175
alpar@25	1176
kpeter@80	1177	/// Negative of a map
alpar@25	1178
kpeter@82	1179	/// This \ref concepts::ReadMap "read-only map" returns the negative
kpeter@80	1180	/// of the values of the given map (using the unary \c - operator).
kpeter@80	1181	/// Its \c Key and \c Value are inherited from \c M.
alpar@25	1182	///
kpeter@80	1183	/// If M::Value is \c int, \c double etc., then
kpeter@80	1184	/// \code
kpeter@80	1185	/// NegMap<M> neg(m);
kpeter@80	1186	/// \endcode
kpeter@80	1187	/// is equivalent to
kpeter@80	1188	/// \code
kpeter@80	1189	/// ScaleMap<M> neg(m,-1);
kpeter@80	1190	/// \endcode
kpeter@29	1191	///
kpeter@80	1192	/// The simplest way of using this map is through the negMap()
kpeter@80	1193	/// function.
kpeter@29	1194	///
kpeter@80	1195	/// \sa NegWriteMap
kpeter@80	1196	template<typename M>
alpar@25	1197	class NegMap : public MapBase<typename M::Key, typename M::Value> {
kpeter@80	1198	const M& _m;
alpar@25	1199	public:
alpar@25	1200	typedef MapBase<typename M::Key, typename M::Value> Parent;
alpar@25	1201	typedef typename Parent::Key Key;
alpar@25	1202	typedef typename Parent::Value Value;
alpar@25	1203
kpeter@80	1204	/// Constructor
kpeter@80	1205	NegMap(const M &m) : _m(m) {}
alpar@25	1206	/// \e
kpeter@80	1207	Value operator[](const Key &k) const { return -_m[k]; }
alpar@25	1208	};
alpar@25	1209
kpeter@80	1210	/// Negative of a map (read-write version)
kpeter@80	1211
kpeter@80	1212	/// This \ref concepts::ReadWriteMap "read-write map" returns the
kpeter@80	1213	/// negative of the values of the given map (using the unary \c -
kpeter@80	1214	/// operator).
kpeter@80	1215	/// Its \c Key and \c Value are inherited from \c M.
kpeter@80	1216	/// It makes also possible to write the map.
kpeter@80	1217	///
kpeter@80	1218	/// If M::Value is \c int, \c double etc., then
kpeter@80	1219	/// \code
kpeter@80	1220	/// NegWriteMap<M> neg(m);
kpeter@80	1221	/// \endcode
kpeter@80	1222	/// is equivalent to
kpeter@80	1223	/// \code
kpeter@80	1224	/// ScaleWriteMap<M> neg(m,-1);
kpeter@80	1225	/// \endcode
kpeter@80	1226	///
kpeter@80	1227	/// The simplest way of using this map is through the negWriteMap()
kpeter@80	1228	/// function.
kpeter@29	1229	///
kpeter@29	1230	/// \sa NegMap
kpeter@80	1231	template<typename M>
alpar@25	1232	class NegWriteMap : public MapBase<typename M::Key, typename M::Value> {
kpeter@80	1233	M &_m;
alpar@25	1234	public:
alpar@25	1235	typedef MapBase<typename M::Key, typename M::Value> Parent;
alpar@25	1236	typedef typename Parent::Key Key;
alpar@25	1237	typedef typename Parent::Value Value;
alpar@25	1238
kpeter@80	1239	/// Constructor
kpeter@80	1240	NegWriteMap(M &m) : _m(m) {}
alpar@25	1241	/// \e
kpeter@80	1242	Value operator[](const Key &k) const { return -_m[k]; }
alpar@25	1243	/// \e
kpeter@80	1244	void set(const Key &k, const Value &v) { _m.set(k, -v); }
alpar@25	1245	};
alpar@25	1246
kpeter@80	1247	/// Returns a \ref NegMap class
alpar@25	1248
kpeter@80	1249	/// This function just returns a \ref NegMap class.
kpeter@80	1250	///
kpeter@80	1251	/// For example, if \c m is a map with \c double values, then
kpeter@80	1252	/// <tt>negMap(m)[x]</tt> will be equal to <tt>-m[x]</tt>.
kpeter@80	1253	///
kpeter@80	1254	/// \relates NegMap
kpeter@80	1255	template <typename M>
alpar@25	1256	inline NegMap<M> negMap(const M &m) {
alpar@25	1257	return NegMap<M>(m);
alpar@25	1258	}
alpar@25	1259
kpeter@80	1260	/// Returns a \ref NegWriteMap class
kpeter@29	1261
kpeter@80	1262	/// This function just returns a \ref NegWriteMap class.
kpeter@80	1263	///
kpeter@80	1264	/// For example, if \c m is a map with \c double values, then
kpeter@80	1265	/// <tt>negWriteMap(m)[x]</tt> will be equal to <tt>-m[x]</tt>.
kpeter@80	1266	/// Moreover it makes also possible to write the map.
kpeter@80	1267	///
kpeter@80	1268	/// \relates NegWriteMap
kpeter@80	1269	template <typename M>
kpeter@80	1270	inline NegWriteMap<M> negWriteMap(M &m) {
alpar@25	1271	return NegWriteMap<M>(m);
alpar@25	1272	}
alpar@25	1273
alpar@25	1274
kpeter@80	1275	/// Absolute value of a map
kpeter@80	1276
kpeter@82	1277	/// This \ref concepts::ReadMap "read-only map" returns the absolute
kpeter@80	1278	/// value of the values of the given map.
kpeter@80	1279	/// Its \c Key and \c Value are inherited from \c M.
kpeter@80	1280	/// \c Value must be comparable to \c 0 and the unary \c -
kpeter@80	1281	/// operator must be defined for it, of course.
kpeter@80	1282	///
kpeter@80	1283	/// The simplest way of using this map is through the absMap()
kpeter@80	1284	/// function.
kpeter@80	1285	template<typename M>
alpar@25	1286	class AbsMap : public MapBase<typename M::Key, typename M::Value> {
kpeter@80	1287	const M &_m;
alpar@25	1288	public:
alpar@25	1289	typedef MapBase<typename M::Key, typename M::Value> Parent;
alpar@25	1290	typedef typename Parent::Key Key;
alpar@25	1291	typedef typename Parent::Value Value;
alpar@25	1292
kpeter@80	1293	/// Constructor
kpeter@80	1294	AbsMap(const M &m) : _m(m) {}
alpar@25	1295	/// \e
kpeter@80	1296	Value operator[](const Key &k) const {
kpeter@80	1297	Value tmp = _m[k];
alpar@25	1298	return tmp >= 0 ? tmp : -tmp;
alpar@25	1299	}
alpar@25	1300
alpar@25	1301	};
alpar@25	1302
kpeter@80	1303	/// Returns an \ref AbsMap class
kpeter@80	1304
kpeter@80	1305	/// This function just returns an \ref AbsMap class.
kpeter@80	1306	///
kpeter@80	1307	/// For example, if \c m is a map with \c double values, then
kpeter@80	1308	/// <tt>absMap(m)[x]</tt> will be equal to <tt>m[x]</tt> if
kpeter@80	1309	/// it is positive or zero and <tt>-m[x]</tt> if <tt>m[x]</tt> is
kpeter@80	1310	/// negative.
kpeter@80	1311	///
kpeter@80	1312	/// \relates AbsMap
kpeter@80	1313	template<typename M>
alpar@25	1314	inline AbsMap<M> absMap(const M &m) {
alpar@25	1315	return AbsMap<M>(m);
alpar@25	1316	}
alpar@25	1317
kpeter@82	1318	/// @}
alpar@209	1319
kpeter@82	1320	// Logical maps and map adaptors:
kpeter@82	1321
kpeter@82	1322	/// \addtogroup maps
kpeter@82	1323	/// @{
kpeter@82	1324
kpeter@82	1325	/// Constant \c true map.
kpeter@82	1326
kpeter@82	1327	/// This \ref concepts::ReadMap "read-only map" assigns \c true to
kpeter@82	1328	/// each key.
kpeter@82	1329	///
kpeter@82	1330	/// Note that
kpeter@82	1331	/// \code
kpeter@82	1332	/// TrueMap<K> tm;
kpeter@82	1333	/// \endcode
kpeter@82	1334	/// is equivalent to
kpeter@82	1335	/// \code
kpeter@82	1336	/// ConstMap<K,bool> tm(true);
kpeter@82	1337	/// \endcode
kpeter@82	1338	///
kpeter@82	1339	/// \sa FalseMap
kpeter@82	1340	/// \sa ConstMap
kpeter@82	1341	template <typename K>
kpeter@82	1342	class TrueMap : public MapBase<K, bool> {
kpeter@82	1343	public:
kpeter@82	1344	typedef MapBase<K, bool> Parent;
kpeter@82	1345	typedef typename Parent::Key Key;
kpeter@82	1346	typedef typename Parent::Value Value;
kpeter@82	1347
kpeter@82	1348	/// Gives back \c true.
kpeter@82	1349	Value operator[](const Key&) const { return true; }
kpeter@82	1350	};
kpeter@82	1351
kpeter@82	1352	/// Returns a \ref TrueMap class
kpeter@82	1353
kpeter@82	1354	/// This function just returns a \ref TrueMap class.
kpeter@82	1355	/// \relates TrueMap
kpeter@82	1356	template<typename K>
kpeter@82	1357	inline TrueMap<K> trueMap() {
kpeter@82	1358	return TrueMap<K>();
kpeter@82	1359	}
kpeter@82	1360
kpeter@82	1361
kpeter@82	1362	/// Constant \c false map.
kpeter@82	1363
kpeter@82	1364	/// This \ref concepts::ReadMap "read-only map" assigns \c false to
kpeter@82	1365	/// each key.
kpeter@82	1366	///
kpeter@82	1367	/// Note that
kpeter@82	1368	/// \code
kpeter@82	1369	/// FalseMap<K> fm;
kpeter@82	1370	/// \endcode
kpeter@82	1371	/// is equivalent to
kpeter@82	1372	/// \code
kpeter@82	1373	/// ConstMap<K,bool> fm(false);
kpeter@82	1374	/// \endcode
kpeter@82	1375	///
kpeter@82	1376	/// \sa TrueMap
kpeter@82	1377	/// \sa ConstMap
kpeter@82	1378	template <typename K>
kpeter@82	1379	class FalseMap : public MapBase<K, bool> {
kpeter@82	1380	public:
kpeter@82	1381	typedef MapBase<K, bool> Parent;
kpeter@82	1382	typedef typename Parent::Key Key;
kpeter@82	1383	typedef typename Parent::Value Value;
kpeter@82	1384
kpeter@82	1385	/// Gives back \c false.
kpeter@82	1386	Value operator[](const Key&) const { return false; }
kpeter@82	1387	};
kpeter@82	1388
kpeter@82	1389	/// Returns a \ref FalseMap class
kpeter@82	1390
kpeter@82	1391	/// This function just returns a \ref FalseMap class.
kpeter@82	1392	/// \relates FalseMap
kpeter@82	1393	template<typename K>
kpeter@82	1394	inline FalseMap<K> falseMap() {
kpeter@82	1395	return FalseMap<K>();
kpeter@82	1396	}
kpeter@82	1397
kpeter@82	1398	/// @}
kpeter@82	1399
kpeter@82	1400	/// \addtogroup map_adaptors
kpeter@82	1401	/// @{
kpeter@82	1402
kpeter@82	1403	/// Logical 'and' of two maps
kpeter@82	1404
kpeter@82	1405	/// This \ref concepts::ReadMap "read-only map" returns the logical
kpeter@82	1406	/// 'and' of the values of the two given maps.
kpeter@82	1407	/// Its \c Key type is inherited from \c M1 and its \c Value type is
kpeter@82	1408	/// \c bool. \c M2::Key must be convertible to \c M1::Key.
kpeter@82	1409	///
kpeter@82	1410	/// If \c m1 is of type \c M1 and \c m2 is of \c M2, then for
kpeter@82	1411	/// \code
kpeter@82	1412	/// AndMap<M1,M2> am(m1,m2);
kpeter@82	1413	/// \endcode
kpeter@82	1414	/// <tt>am[x]</tt> will be equal to <tt>m1[x]&&m2[x]</tt>.
kpeter@82	1415	///
kpeter@82	1416	/// The simplest way of using this map is through the andMap()
kpeter@82	1417	/// function.
kpeter@82	1418	///
kpeter@82	1419	/// \sa OrMap
kpeter@82	1420	/// \sa NotMap, NotWriteMap
kpeter@82	1421	template<typename M1, typename M2>
kpeter@82	1422	class AndMap : public MapBase<typename M1::Key, bool> {
kpeter@82	1423	const M1 &_m1;
kpeter@82	1424	const M2 &_m2;
kpeter@82	1425	public:
kpeter@82	1426	typedef MapBase<typename M1::Key, bool> Parent;
kpeter@82	1427	typedef typename Parent::Key Key;
kpeter@82	1428	typedef typename Parent::Value Value;
kpeter@82	1429
kpeter@82	1430	/// Constructor
kpeter@82	1431	AndMap(const M1 &m1, const M2 &m2) : _m1(m1), _m2(m2) {}
kpeter@82	1432	/// \e
kpeter@82	1433	Value operator[](const Key &k) const { return _m1[k]&&_m2[k]; }
kpeter@82	1434	};
kpeter@82	1435
kpeter@82	1436	/// Returns an \ref AndMap class
kpeter@82	1437
kpeter@82	1438	/// This function just returns an \ref AndMap class.
kpeter@82	1439	///
kpeter@82	1440	/// For example, if \c m1 and \c m2 are both maps with \c bool values,
kpeter@82	1441	/// then <tt>andMap(m1,m2)[x]</tt> will be equal to
kpeter@82	1442	/// <tt>m1[x]&&m2[x]</tt>.
kpeter@82	1443	///
kpeter@82	1444	/// \relates AndMap
kpeter@82	1445	template<typename M1, typename M2>
kpeter@82	1446	inline AndMap<M1, M2> andMap(const M1 &m1, const M2 &m2) {
kpeter@82	1447	return AndMap<M1, M2>(m1,m2);
kpeter@82	1448	}
kpeter@82	1449
kpeter@82	1450
kpeter@82	1451	/// Logical 'or' of two maps
kpeter@82	1452
kpeter@82	1453	/// This \ref concepts::ReadMap "read-only map" returns the logical
kpeter@82	1454	/// 'or' of the values of the two given maps.
kpeter@82	1455	/// Its \c Key type is inherited from \c M1 and its \c Value type is
kpeter@82	1456	/// \c bool. \c M2::Key must be convertible to \c M1::Key.
kpeter@82	1457	///
kpeter@82	1458	/// If \c m1 is of type \c M1 and \c m2 is of \c M2, then for
kpeter@82	1459	/// \code
kpeter@82	1460	/// OrMap<M1,M2> om(m1,m2);
kpeter@82	1461	/// \endcode
kpeter@82	1462	/// <tt>om[x]</tt> will be equal to <tt>m1[x]\|\|m2[x]</tt>.
kpeter@82	1463	///
kpeter@82	1464	/// The simplest way of using this map is through the orMap()
kpeter@82	1465	/// function.
kpeter@82	1466	///
kpeter@82	1467	/// \sa AndMap
kpeter@82	1468	/// \sa NotMap, NotWriteMap
kpeter@82	1469	template<typename M1, typename M2>
kpeter@82	1470	class OrMap : public MapBase<typename M1::Key, bool> {
kpeter@82	1471	const M1 &_m1;
kpeter@82	1472	const M2 &_m2;
kpeter@82	1473	public:
kpeter@82	1474	typedef MapBase<typename M1::Key, bool> Parent;
kpeter@82	1475	typedef typename Parent::Key Key;
kpeter@82	1476	typedef typename Parent::Value Value;
kpeter@82	1477
kpeter@82	1478	/// Constructor
kpeter@82	1479	OrMap(const M1 &m1, const M2 &m2) : _m1(m1), _m2(m2) {}
kpeter@82	1480	/// \e
kpeter@82	1481	Value operator[](const Key &k) const { return _m1[k]\|\|_m2[k]; }
kpeter@82	1482	};
kpeter@82	1483
kpeter@82	1484	/// Returns an \ref OrMap class
kpeter@82	1485
kpeter@82	1486	/// This function just returns an \ref OrMap class.
kpeter@82	1487	///
kpeter@82	1488	/// For example, if \c m1 and \c m2 are both maps with \c bool values,
kpeter@82	1489	/// then <tt>orMap(m1,m2)[x]</tt> will be equal to
kpeter@82	1490	/// <tt>m1[x]\|\|m2[x]</tt>.
kpeter@82	1491	///
kpeter@82	1492	/// \relates OrMap
kpeter@82	1493	template<typename M1, typename M2>
kpeter@82	1494	inline OrMap<M1, M2> orMap(const M1 &m1, const M2 &m2) {
kpeter@82	1495	return OrMap<M1, M2>(m1,m2);
kpeter@82	1496	}
kpeter@82	1497
alpar@25	1498
kpeter@80	1499	/// Logical 'not' of a map
kpeter@80	1500
kpeter@82	1501	/// This \ref concepts::ReadMap "read-only map" returns the logical
kpeter@80	1502	/// negation of the values of the given map.
kpeter@80	1503	/// Its \c Key is inherited from \c M and its \c Value is \c bool.
alpar@25	1504	///
kpeter@80	1505	/// The simplest way of using this map is through the notMap()
kpeter@80	1506	/// function.
alpar@25	1507	///
kpeter@80	1508	/// \sa NotWriteMap
kpeter@80	1509	template <typename M>
alpar@25	1510	class NotMap : public MapBase<typename M::Key, bool> {
kpeter@80	1511	const M &_m;
alpar@25	1512	public:
alpar@25	1513	typedef MapBase<typename M::Key, bool> Parent;
alpar@25	1514	typedef typename Parent::Key Key;
alpar@25	1515	typedef typename Parent::Value Value;
alpar@25	1516
alpar@25	1517	/// Constructor
kpeter@80	1518	NotMap(const M &m) : _m(m) {}
kpeter@80	1519	/// \e
kpeter@80	1520	Value operator[](const Key &k) const { return !_m[k]; }
alpar@25	1521	};
alpar@25	1522
kpeter@80	1523	/// Logical 'not' of a map (read-write version)
kpeter@80	1524
kpeter@80	1525	/// This \ref concepts::ReadWriteMap "read-write map" returns the
kpeter@80	1526	/// logical negation of the values of the given map.
kpeter@80	1527	/// Its \c Key is inherited from \c M and its \c Value is \c bool.
kpeter@80	1528	/// It makes also possible to write the map. When a value is set,
kpeter@80	1529	/// the opposite value is set to the original map.
kpeter@29	1530	///
kpeter@80	1531	/// The simplest way of using this map is through the notWriteMap()
kpeter@80	1532	/// function.
kpeter@80	1533	///
kpeter@80	1534	/// \sa NotMap
kpeter@80	1535	template <typename M>
alpar@25	1536	class NotWriteMap : public MapBase<typename M::Key, bool> {
kpeter@80	1537	M &_m;
alpar@25	1538	public:
alpar@25	1539	typedef MapBase<typename M::Key, bool> Parent;
alpar@25	1540	typedef typename Parent::Key Key;
alpar@25	1541	typedef typename Parent::Value Value;
alpar@25	1542
alpar@25	1543	/// Constructor
kpeter@80	1544	NotWriteMap(M &m) : _m(m) {}
kpeter@80	1545	/// \e
kpeter@80	1546	Value operator[](const Key &k) const { return !_m[k]; }
kpeter@80	1547	/// \e
kpeter@80	1548	void set(const Key &k, bool v) { _m.set(k, !v); }
alpar@25	1549	};
kpeter@80	1550
kpeter@80	1551	/// Returns a \ref NotMap class
kpeter@80	1552
kpeter@80	1553	/// This function just returns a \ref NotMap class.
kpeter@80	1554	///
kpeter@80	1555	/// For example, if \c m is a map with \c bool values, then
kpeter@80	1556	/// <tt>notMap(m)[x]</tt> will be equal to <tt>!m[x]</tt>.
kpeter@80	1557	///
kpeter@80	1558	/// \relates NotMap
kpeter@80	1559	template <typename M>
alpar@25	1560	inline NotMap<M> notMap(const M &m) {
alpar@25	1561	return NotMap<M>(m);
alpar@25	1562	}
kpeter@80	1563
kpeter@80	1564	/// Returns a \ref NotWriteMap class
kpeter@80	1565
kpeter@80	1566	/// This function just returns a \ref NotWriteMap class.
kpeter@80	1567	///
kpeter@80	1568	/// For example, if \c m is a map with \c bool values, then
kpeter@80	1569	/// <tt>notWriteMap(m)[x]</tt> will be equal to <tt>!m[x]</tt>.
kpeter@80	1570	/// Moreover it makes also possible to write the map.
kpeter@80	1571	///
kpeter@80	1572	/// \relates NotWriteMap
kpeter@80	1573	template <typename M>
kpeter@80	1574	inline NotWriteMap<M> notWriteMap(M &m) {
alpar@25	1575	return NotWriteMap<M>(m);
alpar@25	1576	}
alpar@25	1577
kpeter@82	1578
kpeter@82	1579	/// Combination of two maps using the \c == operator
kpeter@82	1580
kpeter@82	1581	/// This \ref concepts::ReadMap "read-only map" assigns \c true to
kpeter@82	1582	/// the keys for which the corresponding values of the two maps are
kpeter@82	1583	/// equal.
kpeter@82	1584	/// Its \c Key type is inherited from \c M1 and its \c Value type is
kpeter@82	1585	/// \c bool. \c M2::Key must be convertible to \c M1::Key.
kpeter@82	1586	///
kpeter@82	1587	/// If \c m1 is of type \c M1 and \c m2 is of \c M2, then for
kpeter@82	1588	/// \code
kpeter@82	1589	/// EqualMap<M1,M2> em(m1,m2);
kpeter@82	1590	/// \endcode
kpeter@82	1591	/// <tt>em[x]</tt> will be equal to <tt>m1[x]==m2[x]</tt>.
kpeter@82	1592	///
kpeter@82	1593	/// The simplest way of using this map is through the equalMap()
kpeter@82	1594	/// function.
kpeter@82	1595	///
kpeter@82	1596	/// \sa LessMap
kpeter@82	1597	template<typename M1, typename M2>
kpeter@82	1598	class EqualMap : public MapBase<typename M1::Key, bool> {
kpeter@82	1599	const M1 &_m1;
kpeter@82	1600	const M2 &_m2;
kpeter@82	1601	public:
kpeter@82	1602	typedef MapBase<typename M1::Key, bool> Parent;
kpeter@82	1603	typedef typename Parent::Key Key;
kpeter@82	1604	typedef typename Parent::Value Value;
kpeter@82	1605
kpeter@82	1606	/// Constructor
kpeter@82	1607	EqualMap(const M1 &m1, const M2 &m2) : _m1(m1), _m2(m2) {}
kpeter@82	1608	/// \e
kpeter@82	1609	Value operator[](const Key &k) const { return _m1[k]==_m2[k]; }
kpeter@82	1610	};
kpeter@82	1611
kpeter@82	1612	/// Returns an \ref EqualMap class
kpeter@82	1613
kpeter@82	1614	/// This function just returns an \ref EqualMap class.
kpeter@82	1615	///
kpeter@82	1616	/// For example, if \c m1 and \c m2 are maps with keys and values of
kpeter@82	1617	/// the same type, then <tt>equalMap(m1,m2)[x]</tt> will be equal to
kpeter@82	1618	/// <tt>m1[x]==m2[x]</tt>.
kpeter@82	1619	///
kpeter@82	1620	/// \relates EqualMap
kpeter@82	1621	template<typename M1, typename M2>
kpeter@82	1622	inline EqualMap<M1, M2> equalMap(const M1 &m1, const M2 &m2) {
kpeter@82	1623	return EqualMap<M1, M2>(m1,m2);
kpeter@82	1624	}
kpeter@82	1625
kpeter@82	1626
kpeter@82	1627	/// Combination of two maps using the \c < operator
kpeter@82	1628
kpeter@82	1629	/// This \ref concepts::ReadMap "read-only map" assigns \c true to
kpeter@82	1630	/// the keys for which the corresponding value of the first map is
kpeter@82	1631	/// less then the value of the second map.
kpeter@82	1632	/// Its \c Key type is inherited from \c M1 and its \c Value type is
kpeter@82	1633	/// \c bool. \c M2::Key must be convertible to \c M1::Key.
kpeter@82	1634	///
kpeter@82	1635	/// If \c m1 is of type \c M1 and \c m2 is of \c M2, then for
kpeter@82	1636	/// \code
kpeter@82	1637	/// LessMap<M1,M2> lm(m1,m2);
kpeter@82	1638	/// \endcode
kpeter@82	1639	/// <tt>lm[x]</tt> will be equal to <tt>m1[x]<m2[x]</tt>.
kpeter@82	1640	///
kpeter@82	1641	/// The simplest way of using this map is through the lessMap()
kpeter@82	1642	/// function.
kpeter@82	1643	///
kpeter@82	1644	/// \sa EqualMap
kpeter@82	1645	template<typename M1, typename M2>
kpeter@82	1646	class LessMap : public MapBase<typename M1::Key, bool> {
kpeter@82	1647	const M1 &_m1;
kpeter@82	1648	const M2 &_m2;
kpeter@82	1649	public:
kpeter@82	1650	typedef MapBase<typename M1::Key, bool> Parent;
kpeter@82	1651	typedef typename Parent::Key Key;
kpeter@82	1652	typedef typename Parent::Value Value;
kpeter@82	1653
kpeter@82	1654	/// Constructor
kpeter@82	1655	LessMap(const M1 &m1, const M2 &m2) : _m1(m1), _m2(m2) {}
kpeter@82	1656	/// \e
kpeter@82	1657	Value operator[](const Key &k) const { return _m1[k]<_m2[k]; }
kpeter@82	1658	};
kpeter@82	1659
kpeter@82	1660	/// Returns an \ref LessMap class
kpeter@82	1661
kpeter@82	1662	/// This function just returns an \ref LessMap class.
kpeter@82	1663	///
kpeter@82	1664	/// For example, if \c m1 and \c m2 are maps with keys and values of
kpeter@82	1665	/// the same type, then <tt>lessMap(m1,m2)[x]</tt> will be equal to
kpeter@82	1666	/// <tt>m1[x]<m2[x]</tt>.
kpeter@82	1667	///
kpeter@82	1668	/// \relates LessMap
kpeter@82	1669	template<typename M1, typename M2>
kpeter@82	1670	inline LessMap<M1, M2> lessMap(const M1 &m1, const M2 &m2) {
kpeter@82	1671	return LessMap<M1, M2>(m1,m2);
kpeter@82	1672	}
kpeter@82	1673
alpar@104	1674	namespace _maps_bits {
alpar@104	1675
alpar@104	1676	template <typename _Iterator, typename Enable = void>
alpar@104	1677	struct IteratorTraits {
alpar@104	1678	typedef typename std::iterator_traits<_Iterator>::value_type Value;
alpar@104	1679	};
alpar@104	1680
alpar@104	1681	template <typename _Iterator>
alpar@104	1682	struct IteratorTraits<_Iterator,
alpar@104	1683	typename exists<typename _Iterator::container_type>::type>
alpar@104	1684	{
alpar@104	1685	typedef typename _Iterator::container_type::value_type Value;
alpar@104	1686	};
alpar@104	1687
alpar@104	1688	}
alpar@104	1689
alpar@104	1690	/// \brief Writable bool map for logging each \c true assigned element
alpar@104	1691	///
kpeter@159	1692	/// A \ref concepts::WriteMap "writable" bool map for logging
alpar@104	1693	/// each \c true assigned element, i.e it copies subsequently each
alpar@104	1694	/// keys set to \c true to the given iterator.
kpeter@159	1695	/// The most important usage of it is storing certain nodes or arcs
kpeter@159	1696	/// that were marked \c true by an algorithm.
alpar@104	1697	///
kpeter@159	1698	/// There are several algorithms that provide solutions through bool
kpeter@159	1699	/// maps and most of them assign \c true at most once for each key.
kpeter@159	1700	/// In these cases it is a natural request to store each \c true
kpeter@159	1701	/// assigned elements (in order of the assignment), which can be
kpeter@167	1702	/// easily done with LoggerBoolMap.
kpeter@159	1703	///
kpeter@167	1704	/// The simplest way of using this map is through the loggerBoolMap()
kpeter@159	1705	/// function.
kpeter@159	1706	///
kpeter@159	1707	/// \tparam It The type of the iterator.
kpeter@159	1708	/// \tparam Ke The key type of the map. The default value set
kpeter@159	1709	/// according to the iterator type should work in most cases.
alpar@104	1710	///
alpar@104	1711	/// \note The container of the iterator must contain enough space
kpeter@159	1712	/// for the elements or the iterator should be an inserter iterator.
kpeter@159	1713	#ifdef DOXYGEN
kpeter@159	1714	template <typename It, typename Ke>
kpeter@159	1715	#else
alpar@104	1716	template <typename It,
alpar@209	1717	typename Ke=typename _maps_bits::IteratorTraits<It>::Value>
kpeter@159	1718	#endif
kpeter@167	1719	class LoggerBoolMap {
alpar@104	1720	public:
alpar@104	1721	typedef It Iterator;
alpar@104	1722
alpar@104	1723	typedef Ke Key;
alpar@104	1724	typedef bool Value;
alpar@104	1725
alpar@104	1726	/// Constructor
kpeter@167	1727	LoggerBoolMap(Iterator it)
alpar@104	1728	: _begin(it), _end(it) {}
alpar@104	1729
alpar@104	1730	/// Gives back the given iterator set for the first key
alpar@104	1731	Iterator begin() const {
alpar@104	1732	return _begin;
alpar@104	1733	}
alpar@104	1734
alpar@104	1735	/// Gives back the the 'after the last' iterator
alpar@104	1736	Iterator end() const {
alpar@104	1737	return _end;
alpar@104	1738	}
alpar@104	1739
alpar@104	1740	/// The set function of the map
kpeter@159	1741	void set(const Key& key, Value value) {
alpar@104	1742	if (value) {
alpar@209	1743	*_end++ = key;
alpar@104	1744	}
alpar@104	1745	}
alpar@104	1746
alpar@104	1747	private:
alpar@104	1748	Iterator _begin;
kpeter@159	1749	Iterator _end;
alpar@104	1750	};
alpar@209	1751
kpeter@167	1752	/// Returns a \ref LoggerBoolMap class
kpeter@159	1753
kpeter@167	1754	/// This function just returns a \ref LoggerBoolMap class.
kpeter@159	1755	///
kpeter@159	1756	/// The most important usage of it is storing certain nodes or arcs
kpeter@159	1757	/// that were marked \c true by an algorithm.
kpeter@159	1758	/// For example it makes easier to store the nodes in the processing
kpeter@159	1759	/// order of Dfs algorithm, as the following examples show.
kpeter@159	1760	/// \code
kpeter@159	1761	/// std::vector<Node> v;
kpeter@167	1762	/// dfs(g,s).processedMap(loggerBoolMap(std::back_inserter(v))).run();
kpeter@159	1763	/// \endcode
kpeter@159	1764	/// \code
kpeter@159	1765	/// std::vector<Node> v(countNodes(g));
kpeter@167	1766	/// dfs(g,s).processedMap(loggerBoolMap(v.begin())).run();
kpeter@159	1767	/// \endcode
kpeter@159	1768	///
kpeter@159	1769	/// \note The container of the iterator must contain enough space
kpeter@159	1770	/// for the elements or the iterator should be an inserter iterator.
kpeter@159	1771	///
kpeter@167	1772	/// \note LoggerBoolMap is just \ref concepts::WriteMap "writable", so
kpeter@159	1773	/// it cannot be used when a readable map is needed, for example as
kpeter@167	1774	/// \c ReachedMap for \ref Bfs, \ref Dfs and \ref Dijkstra algorithms.
kpeter@159	1775	///
kpeter@167	1776	/// \relates LoggerBoolMap
kpeter@159	1777	template<typename Iterator>
kpeter@167	1778	inline LoggerBoolMap<Iterator> loggerBoolMap(Iterator it) {
kpeter@167	1779	return LoggerBoolMap<Iterator>(it);
kpeter@159	1780	}
alpar@104	1781
deba@220	1782	/// Provides an immutable and unique id for each item in the graph.
deba@220	1783
deba@220	1784	/// The IdMap class provides a unique and immutable id for each item of the
deba@220	1785	/// same type (e.g. node) in the graph. This id is <ul><li>\b unique:
deba@220	1786	/// different items (nodes) get different ids <li>\b immutable: the id of an
deba@220	1787	/// item (node) does not change (even if you delete other nodes). </ul>
deba@220	1788	/// Through this map you get access (i.e. can read) the inner id values of
deba@220	1789	/// the items stored in the graph. This map can be inverted with its member
deba@220	1790	/// class \c InverseMap or with the \c operator() member.
deba@220	1791	///
deba@220	1792	template <typename _Graph, typename _Item>
deba@220	1793	class IdMap {
deba@220	1794	public:
deba@220	1795	typedef _Graph Graph;
deba@220	1796	typedef int Value;
deba@220	1797	typedef _Item Item;
deba@220	1798	typedef _Item Key;
deba@220	1799
deba@220	1800	/// \brief Constructor.
deba@220	1801	///
deba@220	1802	/// Constructor of the map.
deba@220	1803	explicit IdMap(const Graph& graph) : _graph(&graph) {}
deba@220	1804
deba@220	1805	/// \brief Gives back the \e id of the item.
deba@220	1806	///
deba@220	1807	/// Gives back the immutable and unique \e id of the item.
deba@220	1808	int operator[](const Item& item) const { return _graph->id(item);}
deba@220	1809
deba@220	1810	/// \brief Gives back the item by its id.
deba@220	1811	///
deba@220	1812	/// Gives back the item by its id.
deba@220	1813	Item operator()(int id) { return _graph->fromId(id, Item()); }
deba@220	1814
deba@220	1815	private:
deba@220	1816	const Graph* _graph;
deba@220	1817
deba@220	1818	public:
deba@220	1819
deba@220	1820	/// \brief The class represents the inverse of its owner (IdMap).
deba@220	1821	///
deba@220	1822	/// The class represents the inverse of its owner (IdMap).
deba@220	1823	/// \see inverse()
deba@220	1824	class InverseMap {
deba@220	1825	public:
deba@220	1826
deba@220	1827	/// \brief Constructor.
deba@220	1828	///
deba@220	1829	/// Constructor for creating an id-to-item map.
deba@220	1830	explicit InverseMap(const Graph& graph) : _graph(&graph) {}
deba@220	1831
deba@220	1832	/// \brief Constructor.
deba@220	1833	///
deba@220	1834	/// Constructor for creating an id-to-item map.
deba@220	1835	explicit InverseMap(const IdMap& map) : _graph(map._graph) {}
deba@220	1836
deba@220	1837	/// \brief Gives back the given item from its id.
deba@220	1838	///
deba@220	1839	/// Gives back the given item from its id.
deba@220	1840	///
deba@220	1841	Item operator[](int id) const { return _graph->fromId(id, Item());}
deba@220	1842
deba@220	1843	private:
deba@220	1844	const Graph* _graph;
deba@220	1845	};
deba@220	1846
deba@220	1847	/// \brief Gives back the inverse of the map.
deba@220	1848	///
deba@220	1849	/// Gives back the inverse of the IdMap.
deba@220	1850	InverseMap inverse() const { return InverseMap(*_graph);}
deba@220	1851
deba@220	1852	};
deba@220	1853
deba@220	1854
deba@220	1855	/// \brief General invertable graph-map type.
deba@220	1856
deba@220	1857	/// This type provides simple invertable graph-maps.
deba@220	1858	/// The InvertableMap wraps an arbitrary ReadWriteMap
deba@220	1859	/// and if a key is set to a new value then store it
deba@220	1860	/// in the inverse map.
deba@220	1861	///
deba@220	1862	/// The values of the map can be accessed
deba@220	1863	/// with stl compatible forward iterator.
deba@220	1864	///
deba@220	1865	/// \tparam _Graph The graph type.
deba@220	1866	/// \tparam _Item The item type of the graph.
deba@220	1867	/// \tparam _Value The value type of the map.
deba@220	1868	///
deba@220	1869	/// \see IterableValueMap
deba@220	1870	template <typename _Graph, typename _Item, typename _Value>
deba@220	1871	class InvertableMap
deba@220	1872	: protected ItemSetTraits<_Graph, _Item>::template Map<_Value>::Type {
deba@220	1873	private:
deba@220	1874
deba@220	1875	typedef typename ItemSetTraits<_Graph, _Item>::
deba@220	1876	template Map<_Value>::Type Map;
deba@220	1877	typedef _Graph Graph;
deba@220	1878
deba@220	1879	typedef std::map<_Value, _Item> Container;
deba@220	1880	Container _inv_map;
deba@220	1881
deba@220	1882	public:
deba@220	1883
deba@220	1884	/// The key type of InvertableMap (Node, Arc, Edge).
deba@220	1885	typedef typename Map::Key Key;
deba@220	1886	/// The value type of the InvertableMap.
deba@220	1887	typedef typename Map::Value Value;
deba@220	1888
deba@220	1889
deba@220	1890
deba@220	1891	/// \brief Constructor.
deba@220	1892	///
deba@220	1893	/// Construct a new InvertableMap for the graph.
deba@220	1894	///
deba@220	1895	explicit InvertableMap(const Graph& graph) : Map(graph) {}
deba@220	1896
deba@220	1897	/// \brief Forward iterator for values.
deba@220	1898	///
deba@220	1899	/// This iterator is an stl compatible forward
deba@220	1900	/// iterator on the values of the map. The values can
deba@220	1901	/// be accessed in the [beginValue, endValue) range.
deba@220	1902	///
deba@220	1903	class ValueIterator
deba@220	1904	: public std::iterator<std::forward_iterator_tag, Value> {
deba@220	1905	friend class InvertableMap;
deba@220	1906	private:
deba@220	1907	ValueIterator(typename Container::const_iterator _it)
deba@220	1908	: it(_it) {}
deba@220	1909	public:
deba@220	1910
deba@220	1911	ValueIterator() {}
deba@220	1912
deba@220	1913	ValueIterator& operator++() { ++it; return *this; }
deba@220	1914	ValueIterator operator++(int) {
deba@220	1915	ValueIterator tmp(*this);
deba@220	1916	operator++();
deba@220	1917	return tmp;
deba@220	1918	}
deba@220	1919
deba@220	1920	const Value& operator*() const { return it->first; }
deba@220	1921	const Value* operator->() const { return &(it->first); }
deba@220	1922
deba@220	1923	bool operator==(ValueIterator jt) const { return it == jt.it; }
deba@220	1924	bool operator!=(ValueIterator jt) const { return it != jt.it; }
deba@220	1925
deba@220	1926	private:
deba@220	1927	typename Container::const_iterator it;
deba@220	1928	};
deba@220	1929
deba@220	1930	/// \brief Returns an iterator to the first value.
deba@220	1931	///
deba@220	1932	/// Returns an stl compatible iterator to the
deba@220	1933	/// first value of the map. The values of the
deba@220	1934	/// map can be accessed in the [beginValue, endValue)
deba@220	1935	/// range.
deba@220	1936	ValueIterator beginValue() const {
deba@220	1937	return ValueIterator(_inv_map.begin());
deba@220	1938	}
deba@220	1939
deba@220	1940	/// \brief Returns an iterator after the last value.
deba@220	1941	///
deba@220	1942	/// Returns an stl compatible iterator after the
deba@220	1943	/// last value of the map. The values of the
deba@220	1944	/// map can be accessed in the [beginValue, endValue)
deba@220	1945	/// range.
deba@220	1946	ValueIterator endValue() const {
deba@220	1947	return ValueIterator(_inv_map.end());
deba@220	1948	}
deba@220	1949
deba@220	1950	/// \brief The setter function of the map.
deba@220	1951	///
deba@220	1952	/// Sets the mapped value.
deba@220	1953	void set(const Key& key, const Value& val) {
deba@220	1954	Value oldval = Map::operator[](key);
deba@220	1955	typename Container::iterator it = _inv_map.find(oldval);
deba@220	1956	if (it != _inv_map.end() && it->second == key) {
deba@220	1957	_inv_map.erase(it);
deba@220	1958	}
deba@220	1959	_inv_map.insert(make_pair(val, key));
deba@220	1960	Map::set(key, val);
deba@220	1961	}
deba@220	1962
deba@220	1963	/// \brief The getter function of the map.
deba@220	1964	///
deba@220	1965	/// It gives back the value associated with the key.
deba@220	1966	typename MapTraits<Map>::ConstReturnValue
deba@220	1967	operator[](const Key& key) const {
deba@220	1968	return Map::operator[](key);
deba@220	1969	}
deba@220	1970
deba@220	1971	/// \brief Gives back the item by its value.
deba@220	1972	///
deba@220	1973	/// Gives back the item by its value.
deba@220	1974	Key operator()(const Value& key) const {
deba@220	1975	typename Container::const_iterator it = _inv_map.find(key);
deba@220	1976	return it != _inv_map.end() ? it->second : INVALID;
deba@220	1977	}
deba@220	1978
deba@220	1979	protected:
deba@220	1980
deba@220	1981	/// \brief Erase the key from the map.
deba@220	1982	///
deba@220	1983	/// Erase the key to the map. It is called by the
deba@220	1984	/// \c AlterationNotifier.
deba@220	1985	virtual void erase(const Key& key) {
deba@220	1986	Value val = Map::operator[](key);
deba@220	1987	typename Container::iterator it = _inv_map.find(val);
deba@220	1988	if (it != _inv_map.end() && it->second == key) {
deba@220	1989	_inv_map.erase(it);
deba@220	1990	}
deba@220	1991	Map::erase(key);
deba@220	1992	}
deba@220	1993
deba@220	1994	/// \brief Erase more keys from the map.
deba@220	1995	///
deba@220	1996	/// Erase more keys from the map. It is called by the
deba@220	1997	/// \c AlterationNotifier.
deba@220	1998	virtual void erase(const std::vector<Key>& keys) {
deba@220	1999	for (int i = 0; i < int(keys.size()); ++i) {
deba@220	2000	Value val = Map::operator[](keys[i]);
deba@220	2001	typename Container::iterator it = _inv_map.find(val);
deba@220	2002	if (it != _inv_map.end() && it->second == keys[i]) {
deba@220	2003	_inv_map.erase(it);
deba@220	2004	}
deba@220	2005	}
deba@220	2006	Map::erase(keys);
deba@220	2007	}
deba@220	2008
deba@220	2009	/// \brief Clear the keys from the map and inverse map.
deba@220	2010	///
deba@220	2011	/// Clear the keys from the map and inverse map. It is called by the
deba@220	2012	/// \c AlterationNotifier.
deba@220	2013	virtual void clear() {
deba@220	2014	_inv_map.clear();
deba@220	2015	Map::clear();
deba@220	2016	}
deba@220	2017
deba@220	2018	public:
deba@220	2019
deba@220	2020	/// \brief The inverse map type.
deba@220	2021	///
deba@220	2022	/// The inverse of this map. The subscript operator of the map
deba@220	2023	/// gives back always the item what was last assigned to the value.
deba@220	2024	class InverseMap {
deba@220	2025	public:
deba@220	2026	/// \brief Constructor of the InverseMap.
deba@220	2027	///
deba@220	2028	/// Constructor of the InverseMap.
deba@220	2029	explicit InverseMap(const InvertableMap& inverted)
deba@220	2030	: _inverted(inverted) {}
deba@220	2031
deba@220	2032	/// The value type of the InverseMap.
deba@220	2033	typedef typename InvertableMap::Key Value;
deba@220	2034	/// The key type of the InverseMap.
deba@220	2035	typedef typename InvertableMap::Value Key;
deba@220	2036
deba@220	2037	/// \brief Subscript operator.
deba@220	2038	///
deba@220	2039	/// Subscript operator. It gives back always the item
deba@220	2040	/// what was last assigned to the value.
deba@220	2041	Value operator[](const Key& key) const {
deba@220	2042	return _inverted(key);
deba@220	2043	}
deba@220	2044
deba@220	2045	private:
deba@220	2046	const InvertableMap& _inverted;
deba@220	2047	};
deba@220	2048
deba@220	2049	/// \brief It gives back the just readable inverse map.
deba@220	2050	///
deba@220	2051	/// It gives back the just readable inverse map.
deba@220	2052	InverseMap inverse() const {
deba@220	2053	return InverseMap(*this);
deba@220	2054	}
deba@220	2055
deba@220	2056
deba@220	2057
deba@220	2058	};
deba@220	2059
deba@220	2060	/// \brief Provides a mutable, continuous and unique descriptor for each
deba@220	2061	/// item in the graph.
deba@220	2062	///
deba@220	2063	/// The DescriptorMap class provides a unique and continuous (but mutable)
deba@220	2064	/// descriptor (id) for each item of the same type (e.g. node) in the
deba@220	2065	/// graph. This id is <ul><li>\b unique: different items (nodes) get
deba@220	2066	/// different ids <li>\b continuous: the range of the ids is the set of
deba@220	2067	/// integers between 0 and \c n-1, where \c n is the number of the items of
deba@220	2068	/// this type (e.g. nodes) (so the id of a node can change if you delete an
deba@220	2069	/// other node, i.e. this id is mutable). </ul> This map can be inverted
deba@220	2070	/// with its member class \c InverseMap, or with the \c operator() member.
deba@220	2071	///
deba@220	2072	/// \tparam _Graph The graph class the \c DescriptorMap belongs to.
deba@220	2073	/// \tparam _Item The Item is the Key of the Map. It may be Node, Arc or
deba@220	2074	/// Edge.
deba@220	2075	template <typename _Graph, typename _Item>
deba@220	2076	class DescriptorMap
deba@220	2077	: protected ItemSetTraits<_Graph, _Item>::template Map<int>::Type {
deba@220	2078
deba@220	2079	typedef _Item Item;
deba@220	2080	typedef typename ItemSetTraits<_Graph, _Item>::template Map<int>::Type Map;
deba@220	2081
deba@220	2082	public:
deba@220	2083	/// The graph class of DescriptorMap.
deba@220	2084	typedef _Graph Graph;
deba@220	2085
deba@220	2086	/// The key type of DescriptorMap (Node, Arc, Edge).
deba@220	2087	typedef typename Map::Key Key;
deba@220	2088	/// The value type of DescriptorMap.
deba@220	2089	typedef typename Map::Value Value;
deba@220	2090
deba@220	2091	/// \brief Constructor.
deba@220	2092	///
deba@220	2093	/// Constructor for descriptor map.
deba@220	2094	explicit DescriptorMap(const Graph& _graph) : Map(_graph) {
deba@220	2095	Item it;
deba@220	2096	const typename Map::Notifier* nf = Map::notifier();
deba@220	2097	for (nf->first(it); it != INVALID; nf->next(it)) {
deba@220	2098	Map::set(it, _inv_map.size());
deba@220	2099	_inv_map.push_back(it);
deba@220	2100	}
deba@220	2101	}
deba@220	2102
deba@220	2103	protected:
deba@220	2104
deba@220	2105	/// \brief Add a new key to the map.
deba@220	2106	///
deba@220	2107	/// Add a new key to the map. It is called by the
deba@220	2108	/// \c AlterationNotifier.
deba@220	2109	virtual void add(const Item& item) {
deba@220	2110	Map::add(item);
deba@220	2111	Map::set(item, _inv_map.size());
deba@220	2112	_inv_map.push_back(item);
deba@220	2113	}
deba@220	2114
deba@220	2115	/// \brief Add more new keys to the map.
deba@220	2116	///
deba@220	2117	/// Add more new keys to the map. It is called by the
deba@220	2118	/// \c AlterationNotifier.
deba@220	2119	virtual void add(const std::vector<Item>& items) {
deba@220	2120	Map::add(items);
deba@220	2121	for (int i = 0; i < int(items.size()); ++i) {
deba@220	2122	Map::set(items[i], _inv_map.size());
deba@220	2123	_inv_map.push_back(items[i]);
deba@220	2124	}
deba@220	2125	}
deba@220	2126
deba@220	2127	/// \brief Erase the key from the map.
deba@220	2128	///
deba@220	2129	/// Erase the key from the map. It is called by the
deba@220	2130	/// \c AlterationNotifier.
deba@220	2131	virtual void erase(const Item& item) {
deba@220	2132	Map::set(_inv_map.back(), Map::operator[](item));
deba@220	2133	_inv_map[Map::operator[](item)] = _inv_map.back();
deba@220	2134	_inv_map.pop_back();
deba@220	2135	Map::erase(item);
deba@220	2136	}
deba@220	2137
deba@220	2138	/// \brief Erase more keys from the map.
deba@220	2139	///
deba@220	2140	/// Erase more keys from the map. It is called by the
deba@220	2141	/// \c AlterationNotifier.
deba@220	2142	virtual void erase(const std::vector<Item>& items) {
deba@220	2143	for (int i = 0; i < int(items.size()); ++i) {
deba@220	2144	Map::set(_inv_map.back(), Map::operator[](items[i]));
deba@220	2145	_inv_map[Map::operator[](items[i])] = _inv_map.back();
deba@220	2146	_inv_map.pop_back();
deba@220	2147	}
deba@220	2148	Map::erase(items);
deba@220	2149	}
deba@220	2150
deba@220	2151	/// \brief Build the unique map.
deba@220	2152	///
deba@220	2153	/// Build the unique map. It is called by the
deba@220	2154	/// \c AlterationNotifier.
deba@220	2155	virtual void build() {
deba@220	2156	Map::build();
deba@220	2157	Item it;
deba@220	2158	const typename Map::Notifier* nf = Map::notifier();
deba@220	2159	for (nf->first(it); it != INVALID; nf->next(it)) {
deba@220	2160	Map::set(it, _inv_map.size());
deba@220	2161	_inv_map.push_back(it);
deba@220	2162	}
deba@220	2163	}
deba@220	2164
deba@220	2165	/// \brief Clear the keys from the map.
deba@220	2166	///
deba@220	2167	/// Clear the keys from the map. It is called by the
deba@220	2168	/// \c AlterationNotifier.
deba@220	2169	virtual void clear() {
deba@220	2170	_inv_map.clear();
deba@220	2171	Map::clear();
deba@220	2172	}
deba@220	2173
deba@220	2174	public:
deba@220	2175
deba@220	2176	/// \brief Returns the maximal value plus one.
deba@220	2177	///
deba@220	2178	/// Returns the maximal value plus one in the map.
deba@220	2179	unsigned int size() const {
deba@220	2180	return _inv_map.size();
deba@220	2181	}
deba@220	2182
deba@220	2183	/// \brief Swaps the position of the two items in the map.
deba@220	2184	///
deba@220	2185	/// Swaps the position of the two items in the map.
deba@220	2186	void swap(const Item& p, const Item& q) {
deba@220	2187	int pi = Map::operator[](p);
deba@220	2188	int qi = Map::operator[](q);
deba@220	2189	Map::set(p, qi);
deba@220	2190	_inv_map[qi] = p;
deba@220	2191	Map::set(q, pi);
deba@220	2192	_inv_map[pi] = q;
deba@220	2193	}
deba@220	2194
deba@220	2195	/// \brief Gives back the \e descriptor of the item.
deba@220	2196	///
deba@220	2197	/// Gives back the mutable and unique \e descriptor of the map.
deba@220	2198	int operator[](const Item& item) const {
deba@220	2199	return Map::operator[](item);
deba@220	2200	}
deba@220	2201
deba@220	2202	/// \brief Gives back the item by its descriptor.
deba@220	2203	///
deba@220	2204	/// Gives back th item by its descriptor.
deba@220	2205	Item operator()(int id) const {
deba@220	2206	return _inv_map[id];
deba@220	2207	}
deba@220	2208
deba@220	2209	private:
deba@220	2210
deba@220	2211	typedef std::vector<Item> Container;
deba@220	2212	Container _inv_map;
deba@220	2213
deba@220	2214	public:
deba@220	2215	/// \brief The inverse map type of DescriptorMap.
deba@220	2216	///
deba@220	2217	/// The inverse map type of DescriptorMap.
deba@220	2218	class InverseMap {
deba@220	2219	public:
deba@220	2220	/// \brief Constructor of the InverseMap.
deba@220	2221	///
deba@220	2222	/// Constructor of the InverseMap.
deba@220	2223	explicit InverseMap(const DescriptorMap& inverted)
deba@220	2224	: _inverted(inverted) {}
deba@220	2225
deba@220	2226
deba@220	2227	/// The value type of the InverseMap.
deba@220	2228	typedef typename DescriptorMap::Key Value;
deba@220	2229	/// The key type of the InverseMap.
deba@220	2230	typedef typename DescriptorMap::Value Key;
deba@220	2231
deba@220	2232	/// \brief Subscript operator.
deba@220	2233	///
deba@220	2234	/// Subscript operator. It gives back the item
deba@220	2235	/// that the descriptor belongs to currently.
deba@220	2236	Value operator[](const Key& key) const {
deba@220	2237	return _inverted(key);
deba@220	2238	}
deba@220	2239
deba@220	2240	/// \brief Size of the map.
deba@220	2241	///
deba@220	2242	/// Returns the size of the map.
deba@220	2243	unsigned int size() const {
deba@220	2244	return _inverted.size();
deba@220	2245	}
deba@220	2246
deba@220	2247	private:
deba@220	2248	const DescriptorMap& _inverted;
deba@220	2249	};
deba@220	2250
deba@220	2251	/// \brief Gives back the inverse of the map.
deba@220	2252	///
deba@220	2253	/// Gives back the inverse of the map.
deba@220	2254	const InverseMap inverse() const {
deba@220	2255	return InverseMap(*this);
deba@220	2256	}
deba@220	2257	};
deba@220	2258
deba@220	2259	/// \brief Returns the source of the given arc.
deba@220	2260	///
deba@220	2261	/// The SourceMap gives back the source Node of the given arc.
deba@220	2262	/// \see TargetMap
deba@220	2263	template <typename Digraph>
deba@220	2264	class SourceMap {
deba@220	2265	public:
deba@220	2266
deba@220	2267	typedef typename Digraph::Node Value;
deba@220	2268	typedef typename Digraph::Arc Key;
deba@220	2269
deba@220	2270	/// \brief Constructor
deba@220	2271	///
deba@220	2272	/// Constructor
deba@220	2273	/// \param _digraph The digraph that the map belongs to.
deba@220	2274	explicit SourceMap(const Digraph& digraph) : _digraph(digraph) {}
deba@220	2275
deba@220	2276	/// \brief The subscript operator.
deba@220	2277	///
deba@220	2278	/// The subscript operator.
deba@220	2279	/// \param arc The arc
deba@220	2280	/// \return The source of the arc
deba@220	2281	Value operator[](const Key& arc) const {
deba@220	2282	return _digraph.source(arc);
deba@220	2283	}
deba@220	2284
deba@220	2285	private:
deba@220	2286	const Digraph& _digraph;
deba@220	2287	};
deba@220	2288
deba@220	2289	/// \brief Returns a \ref SourceMap class.
deba@220	2290	///
deba@220	2291	/// This function just returns an \ref SourceMap class.
deba@220	2292	/// \relates SourceMap
deba@220	2293	template <typename Digraph>
deba@220	2294	inline SourceMap<Digraph> sourceMap(const Digraph& digraph) {
deba@220	2295	return SourceMap<Digraph>(digraph);
deba@220	2296	}
deba@220	2297
deba@220	2298	/// \brief Returns the target of the given arc.
deba@220	2299	///
deba@220	2300	/// The TargetMap gives back the target Node of the given arc.
deba@220	2301	/// \see SourceMap
deba@220	2302	template <typename Digraph>
deba@220	2303	class TargetMap {
deba@220	2304	public:
deba@220	2305
deba@220	2306	typedef typename Digraph::Node Value;
deba@220	2307	typedef typename Digraph::Arc Key;
deba@220	2308
deba@220	2309	/// \brief Constructor
deba@220	2310	///
deba@220	2311	/// Constructor
deba@220	2312	/// \param _digraph The digraph that the map belongs to.
deba@220	2313	explicit TargetMap(const Digraph& digraph) : _digraph(digraph) {}
deba@220	2314
deba@220	2315	/// \brief The subscript operator.
deba@220	2316	///
deba@220	2317	/// The subscript operator.
deba@220	2318	/// \param e The arc
deba@220	2319	/// \return The target of the arc
deba@220	2320	Value operator[](const Key& e) const {
deba@220	2321	return _digraph.target(e);
deba@220	2322	}
deba@220	2323
deba@220	2324	private:
deba@220	2325	const Digraph& _digraph;
deba@220	2326	};
deba@220	2327
deba@220	2328	/// \brief Returns a \ref TargetMap class.
deba@220	2329	///
deba@220	2330	/// This function just returns a \ref TargetMap class.
deba@220	2331	/// \relates TargetMap
deba@220	2332	template <typename Digraph>
deba@220	2333	inline TargetMap<Digraph> targetMap(const Digraph& digraph) {
deba@220	2334	return TargetMap<Digraph>(digraph);
deba@220	2335	}
deba@220	2336
deba@220	2337	/// \brief Returns the "forward" directed arc view of an edge.
deba@220	2338	///
deba@220	2339	/// Returns the "forward" directed arc view of an edge.
deba@220	2340	/// \see BackwardMap
deba@220	2341	template <typename Graph>
deba@220	2342	class ForwardMap {
deba@220	2343	public:
deba@220	2344
deba@220	2345	typedef typename Graph::Arc Value;
deba@220	2346	typedef typename Graph::Edge Key;
deba@220	2347
deba@220	2348	/// \brief Constructor
deba@220	2349	///
deba@220	2350	/// Constructor
deba@220	2351	/// \param _graph The graph that the map belongs to.
deba@220	2352	explicit ForwardMap(const Graph& graph) : _graph(graph) {}
deba@220	2353
deba@220	2354	/// \brief The subscript operator.
deba@220	2355	///
deba@220	2356	/// The subscript operator.
deba@220	2357	/// \param key An edge
deba@220	2358	/// \return The "forward" directed arc view of edge
deba@220	2359	Value operator[](const Key& key) const {
deba@220	2360	return _graph.direct(key, true);
deba@220	2361	}
deba@220	2362
deba@220	2363	private:
deba@220	2364	const Graph& _graph;
deba@220	2365	};
deba@220	2366
deba@220	2367	/// \brief Returns a \ref ForwardMap class.
deba@220	2368	///
deba@220	2369	/// This function just returns an \ref ForwardMap class.
deba@220	2370	/// \relates ForwardMap
deba@220	2371	template <typename Graph>
deba@220	2372	inline ForwardMap<Graph> forwardMap(const Graph& graph) {
deba@220	2373	return ForwardMap<Graph>(graph);
deba@220	2374	}
deba@220	2375
deba@220	2376	/// \brief Returns the "backward" directed arc view of an edge.
deba@220	2377	///
deba@220	2378	/// Returns the "backward" directed arc view of an edge.
deba@220	2379	/// \see ForwardMap
deba@220	2380	template <typename Graph>
deba@220	2381	class BackwardMap {
deba@220	2382	public:
deba@220	2383
deba@220	2384	typedef typename Graph::Arc Value;
deba@220	2385	typedef typename Graph::Edge Key;
deba@220	2386
deba@220	2387	/// \brief Constructor
deba@220	2388	///
deba@220	2389	/// Constructor
deba@220	2390	/// \param _graph The graph that the map belongs to.
deba@220	2391	explicit BackwardMap(const Graph& graph) : _graph(graph) {}
deba@220	2392
deba@220	2393	/// \brief The subscript operator.
deba@220	2394	///
deba@220	2395	/// The subscript operator.
deba@220	2396	/// \param key An edge
deba@220	2397	/// \return The "backward" directed arc view of edge
deba@220	2398	Value operator[](const Key& key) const {
deba@220	2399	return _graph.direct(key, false);
deba@220	2400	}
deba@220	2401
deba@220	2402	private:
deba@220	2403	const Graph& _graph;
deba@220	2404	};
deba@220	2405
deba@220	2406	/// \brief Returns a \ref BackwardMap class
deba@220	2407
deba@220	2408	/// This function just returns a \ref BackwardMap class.
deba@220	2409	/// \relates BackwardMap
deba@220	2410	template <typename Graph>
deba@220	2411	inline BackwardMap<Graph> backwardMap(const Graph& graph) {
deba@220	2412	return BackwardMap<Graph>(graph);
deba@220	2413	}
deba@220	2414
deba@220	2415	/// \brief Potential difference map
deba@220	2416	///
deba@220	2417	/// If there is an potential map on the nodes then we
deba@220	2418	/// can get an arc map as we get the substraction of the
deba@220	2419	/// values of the target and source.
deba@220	2420	template <typename Digraph, typename NodeMap>
deba@220	2421	class PotentialDifferenceMap {
deba@220	2422	public:
deba@220	2423	typedef typename Digraph::Arc Key;
deba@220	2424	typedef typename NodeMap::Value Value;
deba@220	2425
deba@220	2426	/// \brief Constructor
deba@220	2427	///
deba@220	2428	/// Contructor of the map
deba@220	2429	explicit PotentialDifferenceMap(const Digraph& digraph,
deba@220	2430	const NodeMap& potential)
deba@220	2431	: _digraph(digraph), _potential(potential) {}
deba@220	2432
deba@220	2433	/// \brief Const subscription operator
deba@220	2434	///
deba@220	2435	/// Const subscription operator
deba@220	2436	Value operator[](const Key& arc) const {
deba@220	2437	return _potential[_digraph.target(arc)] -
deba@220	2438	_potential[_digraph.source(arc)];
deba@220	2439	}
deba@220	2440
deba@220	2441	private:
deba@220	2442	const Digraph& _digraph;
deba@220	2443	const NodeMap& _potential;
deba@220	2444	};
deba@220	2445
deba@220	2446	/// \brief Returns a PotentialDifferenceMap.
deba@220	2447	///
deba@220	2448	/// This function just returns a PotentialDifferenceMap.
deba@220	2449	/// \relates PotentialDifferenceMap
deba@220	2450	template <typename Digraph, typename NodeMap>
deba@220	2451	PotentialDifferenceMap<Digraph, NodeMap>
deba@220	2452	potentialDifferenceMap(const Digraph& digraph, const NodeMap& potential) {
deba@220	2453	return PotentialDifferenceMap<Digraph, NodeMap>(digraph, potential);
deba@220	2454	}
deba@220	2455
deba@220	2456	/// \brief Map of the node in-degrees.
deba@220	2457	///
deba@220	2458	/// This map returns the in-degree of a node. Once it is constructed,
deba@220	2459	/// the degrees are stored in a standard NodeMap, so each query is done
deba@220	2460	/// in constant time. On the other hand, the values are updated automatically
deba@220	2461	/// whenever the digraph changes.
deba@220	2462	///
deba@220	2463	/// \warning Besides addNode() and addArc(), a digraph structure may provide
deba@220	2464	/// alternative ways to modify the digraph. The correct behavior of InDegMap
deba@220	2465	/// is not guarantied if these additional features are used. For example
deba@220	2466	/// the functions \ref ListDigraph::changeSource() "changeSource()",
deba@220	2467	/// \ref ListDigraph::changeTarget() "changeTarget()" and
deba@220	2468	/// \ref ListDigraph::reverseArc() "reverseArc()"
deba@220	2469	/// of \ref ListDigraph will \e not update the degree values correctly.
deba@220	2470	///
deba@220	2471	/// \sa OutDegMap
deba@220	2472
deba@220	2473	template <typename _Digraph>
deba@220	2474	class InDegMap
deba@220	2475	: protected ItemSetTraits<_Digraph, typename _Digraph::Arc>
deba@220	2476	::ItemNotifier::ObserverBase {
deba@220	2477
deba@220	2478	public:
deba@220	2479
deba@220	2480	typedef _Digraph Digraph;
deba@220	2481	typedef int Value;
deba@220	2482	typedef typename Digraph::Node Key;
deba@220	2483
deba@220	2484	typedef typename ItemSetTraits<Digraph, typename Digraph::Arc>
deba@220	2485	::ItemNotifier::ObserverBase Parent;
deba@220	2486
deba@220	2487	private:
deba@220	2488
deba@220	2489	class AutoNodeMap
deba@220	2490	: public ItemSetTraits<Digraph, Key>::template Map<int>::Type {
deba@220	2491	public:
deba@220	2492
deba@220	2493	typedef typename ItemSetTraits<Digraph, Key>::
deba@220	2494	template Map<int>::Type Parent;
deba@220	2495
deba@220	2496	AutoNodeMap(const Digraph& digraph) : Parent(digraph, 0) {}
deba@220	2497
deba@220	2498	virtual void add(const Key& key) {
deba@220	2499	Parent::add(key);
deba@220	2500	Parent::set(key, 0);
deba@220	2501	}
deba@220	2502
deba@220	2503	virtual void add(const std::vector<Key>& keys) {
deba@220	2504	Parent::add(keys);
deba@220	2505	for (int i = 0; i < int(keys.size()); ++i) {
deba@220	2506	Parent::set(keys[i], 0);
deba@220	2507	}
deba@220	2508	}
deba@220	2509
deba@220	2510	virtual void build() {
deba@220	2511	Parent::build();
deba@220	2512	Key it;
deba@220	2513	typename Parent::Notifier* nf = Parent::notifier();
deba@220	2514	for (nf->first(it); it != INVALID; nf->next(it)) {
deba@220	2515	Parent::set(it, 0);
deba@220	2516	}
deba@220	2517	}
deba@220	2518	};
deba@220	2519
deba@220	2520	public:
deba@220	2521
deba@220	2522	/// \brief Constructor.
deba@220	2523	///
deba@220	2524	/// Constructor for creating in-degree map.
deba@220	2525	explicit InDegMap(const Digraph& digraph)
deba@220	2526	: _digraph(digraph), _deg(digraph) {
deba@220	2527	Parent::attach(_digraph.notifier(typename Digraph::Arc()));
deba@220	2528
deba@220	2529	for(typename Digraph::NodeIt it(_digraph); it != INVALID; ++it) {
deba@220	2530	_deg[it] = countInArcs(_digraph, it);
deba@220	2531	}
deba@220	2532	}
deba@220	2533
deba@220	2534	/// Gives back the in-degree of a Node.
deba@220	2535	int operator[](const Key& key) const {
deba@220	2536	return _deg[key];
deba@220	2537	}
deba@220	2538
deba@220	2539	protected:
deba@220	2540
deba@220	2541	typedef typename Digraph::Arc Arc;
deba@220	2542
deba@220	2543	virtual void add(const Arc& arc) {
deba@220	2544	++_deg[_digraph.target(arc)];
deba@220	2545	}
deba@220	2546
deba@220	2547	virtual void add(const std::vector<Arc>& arcs) {
deba@220	2548	for (int i = 0; i < int(arcs.size()); ++i) {
deba@220	2549	++_deg[_digraph.target(arcs[i])];
deba@220	2550	}
deba@220	2551	}
deba@220	2552
deba@220	2553	virtual void erase(const Arc& arc) {
deba@220	2554	--_deg[_digraph.target(arc)];
deba@220	2555	}
deba@220	2556
deba@220	2557	virtual void erase(const std::vector<Arc>& arcs) {
deba@220	2558	for (int i = 0; i < int(arcs.size()); ++i) {
deba@220	2559	--_deg[_digraph.target(arcs[i])];
deba@220	2560	}
deba@220	2561	}
deba@220	2562
deba@220	2563	virtual void build() {
deba@220	2564	for(typename Digraph::NodeIt it(_digraph); it != INVALID; ++it) {
deba@220	2565	_deg[it] = countInArcs(_digraph, it);
deba@220	2566	}
deba@220	2567	}
deba@220	2568
deba@220	2569	virtual void clear() {
deba@220	2570	for(typename Digraph::NodeIt it(_digraph); it != INVALID; ++it) {
deba@220	2571	_deg[it] = 0;
deba@220	2572	}
deba@220	2573	}
deba@220	2574	private:
deba@220	2575
deba@220	2576	const Digraph& _digraph;
deba@220	2577	AutoNodeMap _deg;
deba@220	2578	};
deba@220	2579
deba@220	2580	/// \brief Map of the node out-degrees.
deba@220	2581	///
deba@220	2582	/// This map returns the out-degree of a node. Once it is constructed,
deba@220	2583	/// the degrees are stored in a standard NodeMap, so each query is done
deba@220	2584	/// in constant time. On the other hand, the values are updated automatically
deba@220	2585	/// whenever the digraph changes.
deba@220	2586	///
deba@220	2587	/// \warning Besides addNode() and addArc(), a digraph structure may provide
deba@220	2588	/// alternative ways to modify the digraph. The correct behavior of OutDegMap
deba@220	2589	/// is not guarantied if these additional features are used. For example
deba@220	2590	/// the functions \ref ListDigraph::changeSource() "changeSource()",
deba@220	2591	/// \ref ListDigraph::changeTarget() "changeTarget()" and
deba@220	2592	/// \ref ListDigraph::reverseArc() "reverseArc()"
deba@220	2593	/// of \ref ListDigraph will \e not update the degree values correctly.
deba@220	2594	///
deba@220	2595	/// \sa InDegMap
deba@220	2596
deba@220	2597	template <typename _Digraph>
deba@220	2598	class OutDegMap
deba@220	2599	: protected ItemSetTraits<_Digraph, typename _Digraph::Arc>
deba@220	2600	::ItemNotifier::ObserverBase {
deba@220	2601
deba@220	2602	public:
deba@220	2603
deba@220	2604	typedef _Digraph Digraph;
deba@220	2605	typedef int Value;
deba@220	2606	typedef typename Digraph::Node Key;
deba@220	2607
deba@220	2608	typedef typename ItemSetTraits<Digraph, typename Digraph::Arc>
deba@220	2609	::ItemNotifier::ObserverBase Parent;
deba@220	2610
deba@220	2611	private:
deba@220	2612
deba@220	2613	class AutoNodeMap
deba@220	2614	: public ItemSetTraits<Digraph, Key>::template Map<int>::Type {
deba@220	2615	public:
deba@220	2616
deba@220	2617	typedef typename ItemSetTraits<Digraph, Key>::
deba@220	2618	template Map<int>::Type Parent;
deba@220	2619
deba@220	2620	AutoNodeMap(const Digraph& digraph) : Parent(digraph, 0) {}
deba@220	2621
deba@220	2622	virtual void add(const Key& key) {
deba@220	2623	Parent::add(key);
deba@220	2624	Parent::set(key, 0);
deba@220	2625	}
deba@220	2626	virtual void add(const std::vector<Key>& keys) {
deba@220	2627	Parent::add(keys);
deba@220	2628	for (int i = 0; i < int(keys.size()); ++i) {
deba@220	2629	Parent::set(keys[i], 0);
deba@220	2630	}
deba@220	2631	}
deba@220	2632	virtual void build() {
deba@220	2633	Parent::build();
deba@220	2634	Key it;
deba@220	2635	typename Parent::Notifier* nf = Parent::notifier();
deba@220	2636	for (nf->first(it); it != INVALID; nf->next(it)) {
deba@220	2637	Parent::set(it, 0);
deba@220	2638	}
deba@220	2639	}
deba@220	2640	};
deba@220	2641
deba@220	2642	public:
deba@220	2643
deba@220	2644	/// \brief Constructor.
deba@220	2645	///
deba@220	2646	/// Constructor for creating out-degree map.
deba@220	2647	explicit OutDegMap(const Digraph& digraph)
deba@220	2648	: _digraph(digraph), _deg(digraph) {
deba@220	2649	Parent::attach(_digraph.notifier(typename Digraph::Arc()));
deba@220	2650
deba@220	2651	for(typename Digraph::NodeIt it(_digraph); it != INVALID; ++it) {
deba@220	2652	_deg[it] = countOutArcs(_digraph, it);
deba@220	2653	}
deba@220	2654	}
deba@220	2655
deba@220	2656	/// Gives back the out-degree of a Node.
deba@220	2657	int operator[](const Key& key) const {
deba@220	2658	return _deg[key];
deba@220	2659	}
deba@220	2660
deba@220	2661	protected:
deba@220	2662
deba@220	2663	typedef typename Digraph::Arc Arc;
deba@220	2664
deba@220	2665	virtual void add(const Arc& arc) {
deba@220	2666	++_deg[_digraph.source(arc)];
deba@220	2667	}
deba@220	2668
deba@220	2669	virtual void add(const std::vector<Arc>& arcs) {
deba@220	2670	for (int i = 0; i < int(arcs.size()); ++i) {
deba@220	2671	++_deg[_digraph.source(arcs[i])];
deba@220	2672	}
deba@220	2673	}
deba@220	2674
deba@220	2675	virtual void erase(const Arc& arc) {
deba@220	2676	--_deg[_digraph.source(arc)];
deba@220	2677	}
deba@220	2678
deba@220	2679	virtual void erase(const std::vector<Arc>& arcs) {
deba@220	2680	for (int i = 0; i < int(arcs.size()); ++i) {
deba@220	2681	--_deg[_digraph.source(arcs[i])];
deba@220	2682	}
deba@220	2683	}
deba@220	2684
deba@220	2685	virtual void build() {
deba@220	2686	for(typename Digraph::NodeIt it(_digraph); it != INVALID; ++it) {
deba@220	2687	_deg[it] = countOutArcs(_digraph, it);
deba@220	2688	}
deba@220	2689	}
deba@220	2690
deba@220	2691	virtual void clear() {
deba@220	2692	for(typename Digraph::NodeIt it(_digraph); it != INVALID; ++it) {
deba@220	2693	_deg[it] = 0;
deba@220	2694	}
deba@220	2695	}
deba@220	2696	private:
deba@220	2697
deba@220	2698	const Digraph& _digraph;
deba@220	2699	AutoNodeMap _deg;
deba@220	2700	};
deba@220	2701
alpar@25	2702	/// @}
alpar@25	2703	}
alpar@25	2704
alpar@25	2705	#endif // LEMON_MAPS_H

author	Alpar Juttner <alpar@cs.elte.hu>
	Mon, 22 Sep 2008 08:12:44 +0100
changeset 275	0821fea54d53
parent 209	765619b7cbb2
child 280	e7f8647ce760
permissions	-rw-r--r--