/* -*- C++ -*-

 *

 * This file is a part of LEMON, a generic C++ optimization library

 *

 * Copyright (C) 2003-2006

 * Egervary Jeno Kombinatorikus Optimalizalasi Kutatocsoport

 * (Egervary Research Group on Combinatorial Optimization, EGRES).

 *

 * Permission to use, modify and distribute this software is granted

 * provided that this copyright notice appears in all copies. For

 * precise terms see the accompanying LICENSE file.

 *

 * This software is provided "AS IS" with no warranty of any kind,

 * express or implied, and with no claim as to its suitability for any

 * purpose.

 *

 */

#ifndef LEMON_UGRAPH_ADAPTOR_H

#define LEMON_UGRAPH_ADAPTOR_H

///\ingroup graph_adaptors

///\file

///\brief Several graph adaptors.

///

///This file contains several useful ugraph adaptor functions.

///

///\author Balazs Dezso

#include <lemon/invalid.h>

#include <lemon/maps.h>

#include <lemon/bits/graph_adaptor_extender.h>

#include <lemon/traits.h>

#include <iostream>

namespace lemon {

  /// \ingroup graph_adaptors

  ///

  /// \brief Base type for the Graph Adaptors

  ///

  /// \warning Graph adaptors are in even more experimental state than the 

  /// other parts of the lib. Use them at you own risk.

  ///

  /// This is the base type for most of LEMON graph adaptors. 

  /// This class implements a trivial graph adaptor i.e. it only wraps the 

  /// functions and types of the graph. The purpose of this class is to 

  /// make easier implementing graph adaptors. E.g. if an adaptor is 

  /// considered which differs from the wrapped graph only in some of its 

  /// functions or types, then it can be derived from GraphAdaptor, and only 

  /// the differences should be implemented.

  ///

  /// \author Balazs Dezso 

  template<typename _UGraph>

  class UGraphAdaptorBase {

  public:

    typedef _UGraph Graph;

    typedef Graph ParentGraph;

  protected:

    Graph* graph;

    UGraphAdaptorBase() : graph(0) {}

    void setGraph(Graph& _graph) { graph=&_graph; }

    Graph& getGraph() { return *graph; }

    const Graph& getGraph() const { return *graph; }

  public:

    UGraphAdaptorBase(Graph& _graph) : graph(&_graph) {}

    typedef typename Graph::Node Node;

    typedef typename Graph::Edge Edge;

    typedef typename Graph::UEdge UEdge;

    void first(Node& i) const { graph->first(i); }

    void first(Edge& i) const { graph->first(i); }

    void first(UEdge& i) const { graph->first(i); }

    void firstIn(Edge& i, const Node& n) const { graph->firstIn(i, n); }

    void firstOut(Edge& i, const Node& n ) const { graph->firstOut(i, n); }

    void firstInc(UEdge &i, bool &d, const Node &n) const {

      graph->firstInc(i, d, n);

    }

    void next(Node& i) const { graph->next(i); }

    void next(Edge& i) const { graph->next(i); }

    void next(UEdge& i) const { graph->next(i); }

    void nextIn(Edge& i) const { graph->nextIn(i); }

    void nextOut(Edge& i) const { graph->nextOut(i); }

    void nextInc(UEdge &i, bool &d) const { graph->nextInc(i, d); }

    Node source(const UEdge& e) const { return graph->source(e); }

    Node target(const UEdge& e) const { return graph->target(e); }

    Node source(const Edge& e) const { return graph->source(e); }

    Node target(const Edge& e) const { return graph->target(e); }

    typedef NodeNumTagIndicator<Graph> NodeNumTag;

    int nodeNum() const { return graph->nodeNum(); }

    typedef EdgeNumTagIndicator<Graph> EdgeNumTag;

    int edgeNum() const { return graph->edgeNum(); }

    typedef FindEdgeTagIndicator<Graph> FindEdgeTag;

    Edge findEdge(const Node& source, const Node& target, 

		  const Edge& prev = INVALID) {

      return graph->findEdge(source, target, prev);

    }

    UEdge findUEdge(const Node& source, const Node& target, 

                    const UEdge& prev = INVALID) {

      return graph->findUEdge(source, target, prev);

    }

    Node addNode() const { return graph->addNode(); }

    UEdge addEdge(const Node& source, const Node& target) const { 

      return graph->addEdge(source, target); 

    }

    void erase(const Node& i) const { graph->erase(i); }

    void erase(const Edge& i) const { graph->erase(i); }

    void clear() const { graph->clear(); }

    int id(const Node& v) const { return graph->id(v); }

    int id(const UEdge& e) const { return graph->id(e); }

    bool direction(const Edge& e) const { return graph->direction(e); }

    Edge direct(const UEdge& e, bool d) const { return graph->direct(e, d); }

    Edge direct(const UEdge& e, const Node& n) const { 

      return graph->direct(e, n); 

    }

    Node oppositeNode(const Node& n, const Edge& e) const { 

      return graph->oppositeNode(n, e); 

    }

    Edge oppositeEdge(const Edge& e) const { 

      return graph->oppositeEdge(e); 

    }

    template <typename _Value>

    class NodeMap : public Graph::template NodeMap<_Value> {

    public:

      typedef typename Graph::template NodeMap<_Value> Parent;

      explicit NodeMap(const UGraphAdaptorBase<Graph>& ga) 

	: Parent(*ga.graph) {}

      NodeMap(const UGraphAdaptorBase<Graph>& ga, const _Value& value)

	: Parent(*ga.graph, value) {}

      NodeMap& operator=(const NodeMap& cmap) {

	return operator=<NodeMap>(cmap);

      }

      template <typename CMap>

      NodeMap& operator=(const CMap& cmap) {

	checkConcept<concept::ReadMap<Node, _Value>, CMap>();

	const typename Parent::Graph* graph = Parent::getGraph();

	Node it;

	for (graph->first(it); it != INVALID; graph->next(it)) {

	  Parent::set(it, cmap[it]);

	}

	return *this;

      }

    };

    template <typename _Value>

    class EdgeMap : public Graph::template EdgeMap<_Value> {

    public:

      typedef typename Graph::template EdgeMap<_Value> Parent;

      explicit EdgeMap(const UGraphAdaptorBase<Graph>& ga) 

	: Parent(*ga.graph) {}

      EdgeMap(const UGraphAdaptorBase<Graph>& ga, const _Value& value)

	: Parent(*ga.graph, value) {}

      EdgeMap& operator=(const EdgeMap& cmap) {

	return operator=<EdgeMap>(cmap);

      }

      template <typename CMap>

      EdgeMap& operator=(const CMap& cmap) {

	checkConcept<concept::ReadMap<Edge, _Value>, CMap>();

	const typename Parent::Graph* graph = Parent::getGraph();

	Edge it;

	for (graph->first(it); it != INVALID; graph->next(it)) {

	  Parent::set(it, cmap[it]);

	}

	return *this;

      }

    };

    template <typename _Value>

    class UEdgeMap : public Graph::template UEdgeMap<_Value> {

    public:

      typedef typename Graph::template UEdgeMap<_Value> Parent;

      explicit UEdgeMap(const UGraphAdaptorBase<Graph>& ga) 

	: Parent(*ga.graph) {}

      UEdgeMap(const UGraphAdaptorBase<Graph>& ga, const _Value& value)

	: Parent(*ga.graph, value) {}

      UEdgeMap& operator=(const UEdgeMap& cmap) {

	return operator=<UEdgeMap>(cmap);

      }

      template <typename CMap>

      UEdgeMap& operator=(const CMap& cmap) {

	checkConcept<concept::ReadMap<UEdge, _Value>, CMap>();

	const typename Parent::Graph* graph = Parent::getGraph();

	UEdge it;

	for (graph->first(it); it != INVALID; graph->next(it)) {

	  Parent::set(it, cmap[it]);

	}

	return *this;

      }

    };

  };

  /// \ingroup graph_adaptors

  template <typename _UGraph>

  class UGraphAdaptor 

    : public UGraphAdaptorExtender< UGraphAdaptorBase<_UGraph> > { 

  public:

    typedef _UGraph Graph;

    typedef UGraphAdaptorExtender<UGraphAdaptorBase<_UGraph> > Parent;

  protected:

    UGraphAdaptor() : Parent() {}

  public:

    explicit UGraphAdaptor(Graph& _graph) { setGraph(_graph); }

  };

  template <typename _UGraph, typename NodeFilterMap, 

	    typename UEdgeFilterMap, bool checked = true>

  class SubUGraphAdaptorBase : public UGraphAdaptorBase<_UGraph> {

  public:

    typedef _UGraph Graph;

    typedef UGraphAdaptorBase<_UGraph> Parent;

  protected:

    NodeFilterMap* node_filter_map;

    UEdgeFilterMap* uedge_filter_map;

    SubUGraphAdaptorBase() 

      : Parent(), node_filter_map(0), uedge_filter_map(0) { }

    void setNodeFilterMap(NodeFilterMap& _node_filter_map) {

      node_filter_map=&_node_filter_map;

    }

    void setUEdgeFilterMap(UEdgeFilterMap& _uedge_filter_map) {

      uedge_filter_map=&_uedge_filter_map;

    }

  public:

    typedef typename Parent::Node Node;

    typedef typename Parent::Edge Edge;

    typedef typename Parent::UEdge UEdge;

    void first(Node& i) const { 

      Parent::first(i); 

      while (i!=INVALID && !(*node_filter_map)[i]) Parent::next(i); 

    }

    void first(Edge& i) const { 

      Parent::first(i); 

      while (i!=INVALID && (!(*uedge_filter_map)[i] 

	     || !(*node_filter_map)[Parent::source(i)]

	     || !(*node_filter_map)[Parent::target(i)])) Parent::next(i); 

    }

    void first(UEdge& i) const { 

      Parent::first(i); 

      while (i!=INVALID && (!(*uedge_filter_map)[i] 

	     || !(*node_filter_map)[Parent::source(i)]

	     || !(*node_filter_map)[Parent::target(i)])) Parent::next(i); 

    }

    void firstIn(Edge& i, const Node& n) const { 

      Parent::firstIn(i, n); 

      while (i!=INVALID && (!(*uedge_filter_map)[i] 

	     || !(*node_filter_map)[Parent::source(i)])) Parent::nextIn(i); 

    }

    void firstOut(Edge& i, const Node& n) const { 

      Parent::firstOut(i, n); 

      while (i!=INVALID && (!(*uedge_filter_map)[i] 

	     || !(*node_filter_map)[Parent::target(i)])) Parent::nextOut(i); 

    }

    void firstInc(UEdge& i, bool& d, const Node& n) const { 

      Parent::firstInc(i, d, n); 

      while (i!=INVALID && (!(*uedge_filter_map)[i] 

            || !(*node_filter_map)[Parent::target(i)])) Parent::nextInc(i, d); 

    }

    void next(Node& i) const { 

      Parent::next(i); 

      while (i!=INVALID && !(*node_filter_map)[i]) Parent::next(i); 

    }

    void next(Edge& i) const { 

      Parent::next(i); 

      while (i!=INVALID && (!(*uedge_filter_map)[i] 

	     || !(*node_filter_map)[Parent::source(i)]

	     || !(*node_filter_map)[Parent::target(i)])) Parent::next(i); 

    }

    void next(UEdge& i) const { 

      Parent::next(i); 

      while (i!=INVALID && (!(*uedge_filter_map)[i] 

	     || !(*node_filter_map)[Parent::source(i)]

	     || !(*node_filter_map)[Parent::target(i)])) Parent::next(i); 

    }

    void nextIn(Edge& i) const { 

      Parent::nextIn(i); 

      while (i!=INVALID && (!(*uedge_filter_map)[i] 

	     || !(*node_filter_map)[Parent::source(i)])) Parent::nextIn(i); 

    }

    void nextOut(Edge& i) const { 

      Parent::nextOut(i); 

      while (i!=INVALID && (!(*uedge_filter_map)[i] 

	     || !(*node_filter_map)[Parent::target(i)])) Parent::nextOut(i); 

    }

    void nextInc(UEdge& i, bool& d) const { 

      Parent::nextInc(i, d); 

      while (i!=INVALID && (!(*uedge_filter_map)[i] 

            || !(*node_filter_map)[Parent::source(i)])) Parent::nextInc(i, d); 

    }

    ///\e

    /// This function hides \c n in the graph, i.e. the iteration 

    /// jumps over it. This is done by simply setting the value of \c n  

    /// to be false in the corresponding node-map.

    void hide(const Node& n) const { node_filter_map->set(n, false); }

    ///\e

    /// This function hides \c e in the graph, i.e. the iteration 

    /// jumps over it. This is done by simply setting the value of \c e  

    /// to be false in the corresponding edge-map.

    void hide(const UEdge& e) const { uedge_filter_map->set(e, false); }

    ///\e

    /// The value of \c n is set to be true in the node-map which stores 

    /// hide information. If \c n was hidden previuosly, then it is shown 

    /// again

     void unHide(const Node& n) const { node_filter_map->set(n, true); }

    ///\e

    /// The value of \c e is set to be true in the edge-map which stores 

    /// hide information. If \c e was hidden previuosly, then it is shown 

    /// again

    void unHide(const UEdge& e) const { uedge_filter_map->set(e, true); }

    /// Returns true if \c n is hidden.

    ///\e

    ///

    bool hidden(const Node& n) const { return !(*node_filter_map)[n]; }

    /// Returns true if \c n is hidden.

    ///\e

    ///

    bool hidden(const UEdge& e) const { return !(*uedge_filter_map)[e]; }

    typedef False NodeNumTag;

    typedef False EdgeNumTag;

  };

  template <typename _UGraph, typename NodeFilterMap, typename UEdgeFilterMap>

  class SubUGraphAdaptorBase<_UGraph, NodeFilterMap, UEdgeFilterMap, false> 

    : public UGraphAdaptorBase<_UGraph> {

  public:

    typedef _UGraph Graph;

    typedef UGraphAdaptorBase<_UGraph> Parent;

  protected:

    NodeFilterMap* node_filter_map;

    UEdgeFilterMap* uedge_filter_map;

    SubUGraphAdaptorBase() : Parent(), 

			    node_filter_map(0), uedge_filter_map(0) { }

    void setNodeFilterMap(NodeFilterMap& _node_filter_map) {

      node_filter_map=&_node_filter_map;

    }

    void setUEdgeFilterMap(UEdgeFilterMap& _uedge_filter_map) {

      uedge_filter_map=&_uedge_filter_map;

    }

  public:

    typedef typename Parent::Node Node;

    typedef typename Parent::Edge Edge;

    typedef typename Parent::UEdge UEdge;

    void first(Node& i) const { 

      Parent::first(i); 

      while (i!=INVALID && !(*node_filter_map)[i]) Parent::next(i); 

    }

    void first(Edge& i) const { 

      Parent::first(i); 

      while (i!=INVALID && !(*uedge_filter_map)[i]) Parent::next(i); 

    }

    void first(UEdge& i) const { 

      Parent::first(i); 

      while (i!=INVALID && !(*uedge_filter_map)[i]) Parent::next(i); 

    }

    void firstIn(Edge& i, const Node& n) const { 

      Parent::firstIn(i, n); 

      while (i!=INVALID && !(*uedge_filter_map)[i]) Parent::nextIn(i); 

    }

    void firstOut(Edge& i, const Node& n) const { 

      Parent::firstOut(i, n); 

      while (i!=INVALID && !(*uedge_filter_map)[i]) Parent::nextOut(i); 

    }

    void firstInc(UEdge& i, bool& d, const Node& n) const { 

      Parent::firstInc(i, d, n); 

      while (i!=INVALID && !(*uedge_filter_map)[i]) Parent::nextInc(i, d); 

    }

    void next(Node& i) const { 

      Parent::next(i); 

      while (i!=INVALID && !(*node_filter_map)[i]) Parent::next(i); 

    }

    void next(Edge& i) const { 

      Parent::next(i); 

      while (i!=INVALID && !(*uedge_filter_map)[i]) Parent::next(i); 

    }

    void next(UEdge& i) const { 

      Parent::next(i); 

      while (i!=INVALID && !(*uedge_filter_map)[i]) Parent::next(i); 

    }

    void nextIn(Edge& i) const { 

      Parent::nextIn(i); 

      while (i!=INVALID && !(*uedge_filter_map)[i]) Parent::nextIn(i); 

    }

    void nextOut(Edge& i) const { 

      Parent::nextOut(i); 

      while (i!=INVALID && !(*uedge_filter_map)[i]) Parent::nextOut(i); 

    }

    void nextInc(UEdge& i, bool& d) const { 

      Parent::nextInc(i, d); 

      while (i!=INVALID && !(*uedge_filter_map)[i]) Parent::nextInc(i, d); 

    }

    ///\e

    /// This function hides \c n in the graph, i.e. the iteration 

    /// jumps over it. This is done by simply setting the value of \c n  

    /// to be false in the corresponding node-map.

    void hide(const Node& n) const { node_filter_map->set(n, false); }

    ///\e

    /// This function hides \c e in the graph, i.e. the iteration 

    /// jumps over it. This is done by simply setting the value of \c e  

    /// to be false in the corresponding edge-map.

    void hide(const UEdge& e) const { uedge_filter_map->set(e, false); }

    ///\e

    /// The value of \c n is set to be true in the node-map which stores 

    /// hide information. If \c n was hidden previuosly, then it is shown 

    /// again

     void unHide(const Node& n) const { node_filter_map->set(n, true); }

    ///\e

    /// The value of \c e is set to be true in the edge-map which stores 

    /// hide information. If \c e was hidden previuosly, then it is shown 

    /// again

    void unHide(const UEdge& e) const { uedge_filter_map->set(e, true); }

    /// Returns true if \c n is hidden.

    ///\e

    ///

    bool hidden(const Node& n) const { return !(*node_filter_map)[n]; }

    /// Returns true if \c n is hidden.

    ///\e

    ///

    bool hidden(const UEdge& e) const { return !(*uedge_filter_map)[e]; }

    typedef False NodeNumTag;

    typedef False EdgeNumTag;

  };

  /// \ingroup graph_adaptors

  ///

  /// \brief A graph adaptor for hiding nodes and edges from an undirected 

  /// graph.

  /// 

  /// SubUGraphAdaptor shows the undirected graph with filtered node-set and 

  /// edge-set. If the \c checked parameter is true then it filters the edgeset

  /// to do not get invalid edges without source or target.

  /// 

  /// If the \c checked template parameter is false then we have to note that 

  /// the node-iterator cares only the filter on the node-set, and the 

  /// edge-iterator cares only the filter on the edge-set.

  /// This way the edge-map

  /// should filter all edges which's source or target is filtered by the 

  /// node-filter.

  ///

  /// Note that \c n is of type \c SubGA::NodeIt, but it can be converted to

  /// \c Graph::Node that is why \c g.id(n) can be applied.

  /// 

  template<typename _UGraph, typename NodeFilterMap, 

	   typename UEdgeFilterMap, bool checked = true>

  class SubUGraphAdaptor : 

    public UGraphAdaptorExtender<

    SubUGraphAdaptorBase<_UGraph, NodeFilterMap, UEdgeFilterMap, checked> > {

  public:

    typedef _UGraph Graph;

    typedef UGraphAdaptorExtender<

      SubUGraphAdaptorBase<_UGraph, NodeFilterMap, UEdgeFilterMap> > Parent;

  protected:

    SubUGraphAdaptor() { }

  public:

    SubUGraphAdaptor(Graph& _graph, NodeFilterMap& _node_filter_map, 

		    UEdgeFilterMap& _uedge_filter_map) { 

      setGraph(_graph);

      setNodeFilterMap(_node_filter_map);

      setUEdgeFilterMap(_uedge_filter_map);

    }

  };

  template<typename UGraph, typename NodeFilterMap, typename EdgeFilterMap>

  SubUGraphAdaptor<const UGraph, NodeFilterMap, EdgeFilterMap>

  subUGraphAdaptor(const UGraph& graph, 

                   NodeFilterMap& nfm, EdgeFilterMap& efm) {

    return SubUGraphAdaptor<const UGraph, NodeFilterMap, EdgeFilterMap>

      (graph, nfm, efm);

  }

  template<typename UGraph, typename NodeFilterMap, typename EdgeFilterMap>

  SubUGraphAdaptor<const UGraph, const NodeFilterMap, EdgeFilterMap>

  subUGraphAdaptor(const UGraph& graph, 

                   NodeFilterMap& nfm, EdgeFilterMap& efm) {

    return SubUGraphAdaptor<const UGraph, const NodeFilterMap, EdgeFilterMap>

      (graph, nfm, efm);

  }

  template<typename UGraph, typename NodeFilterMap, typename EdgeFilterMap>

  SubUGraphAdaptor<const UGraph, NodeFilterMap, const EdgeFilterMap>

  subUGraphAdaptor(const UGraph& graph, 

                   NodeFilterMap& nfm, EdgeFilterMap& efm) {

    return SubUGraphAdaptor<const UGraph, NodeFilterMap, const EdgeFilterMap>

      (graph, nfm, efm);

  }

  template<typename UGraph, typename NodeFilterMap, typename EdgeFilterMap>

  SubUGraphAdaptor<const UGraph, const NodeFilterMap, const EdgeFilterMap>

  subUGraphAdaptor(const UGraph& graph, 

                   NodeFilterMap& nfm, EdgeFilterMap& efm) {

    return SubUGraphAdaptor<const UGraph, const NodeFilterMap, 

      const EdgeFilterMap>(graph, nfm, efm);

  }

  /// \ingroup graph_adaptors

  ///

  /// \brief An adaptor for hiding nodes from an undirected graph.

  ///

  ///

  /// An adaptor for hiding nodes from an undirected graph.

  /// This adaptor specializes SubUGraphAdaptor in the way that only

  /// the node-set 

  /// can be filtered. In usual case the checked parameter is true, we get the

  /// induced subgraph. But if the checked parameter is false then we can only

  /// filter only isolated nodes.

  template<typename _UGraph, typename NodeFilterMap, bool checked = true>

  class NodeSubUGraphAdaptor : 

    public SubUGraphAdaptor<_UGraph, NodeFilterMap, 

                            ConstMap<typename _UGraph::UEdge, bool>, checked> {

  public:

    typedef SubUGraphAdaptor<_UGraph, NodeFilterMap, 

                             ConstMap<typename _UGraph::UEdge, bool> > Parent;

    typedef _UGraph Graph;

  protected:

    ConstMap<typename _UGraph::UEdge, bool> const_true_map;

  public:

    NodeSubUGraphAdaptor(Graph& _graph, NodeFilterMap& _node_filter_map) : 

      Parent(), const_true_map(true) { 

      Parent::setGraph(_graph);

      Parent::setNodeFilterMap(_node_filter_map);

      Parent::setUEdgeFilterMap(const_true_map);

    }

  };

  template<typename UGraph, typename NodeFilterMap>

  NodeSubUGraphAdaptor<const UGraph, NodeFilterMap>

  nodeSubUGraphAdaptor(const UGraph& graph, NodeFilterMap& nfm) {

    return NodeSubUGraphAdaptor<const UGraph, NodeFilterMap>(graph, nfm);

  }

  template<typename UGraph, typename NodeFilterMap>

  NodeSubUGraphAdaptor<const UGraph, const NodeFilterMap>

  nodeSubUGraphAdaptor(const UGraph& graph, const NodeFilterMap& nfm) {

    return NodeSubUGraphAdaptor<const UGraph, const NodeFilterMap>(graph, nfm);

  }

  /// \brief An adaptor for hiding undirected edges from an undirected graph.

  ///

  /// \warning Graph adaptors are in even more experimental state

  /// than the other parts of the lib. Use them at you own risk.

  ///

  /// An adaptor for hiding undirected edges from an undirected graph.

  /// This adaptor specializes SubUGraphAdaptor in the way that

  /// only the edge-set 

  /// can be filtered.

  template<typename _UGraph, typename UEdgeFilterMap>

  class EdgeSubUGraphAdaptor : 

    public SubUGraphAdaptor<_UGraph, ConstMap<typename _UGraph::Node,bool>, 

                            UEdgeFilterMap, false> {

  public:

    typedef SubUGraphAdaptor<_UGraph, ConstMap<typename _UGraph::Node,bool>, 

                             UEdgeFilterMap, false> Parent;

    typedef _UGraph Graph;

  protected:

    ConstMap<typename Graph::Node, bool> const_true_map;

  public:

    EdgeSubUGraphAdaptor(Graph& _graph, UEdgeFilterMap& _uedge_filter_map) : 

      Parent(), const_true_map(true) { 

      Parent::setGraph(_graph);

      Parent::setNodeFilterMap(const_true_map);

      Parent::setUEdgeFilterMap(_uedge_filter_map);

    }

  };

  template<typename UGraph, typename EdgeFilterMap>

  EdgeSubUGraphAdaptor<const UGraph, EdgeFilterMap>

  edgeSubUGraphAdaptor(const UGraph& graph, EdgeFilterMap& efm) {

    return EdgeSubUGraphAdaptor<const UGraph, EdgeFilterMap>(graph, efm);

  }

  template<typename UGraph, typename EdgeFilterMap>

  EdgeSubUGraphAdaptor<const UGraph, const EdgeFilterMap>

  edgeSubUGraphAdaptor(const UGraph& graph, const EdgeFilterMap& efm) {

    return EdgeSubUGraphAdaptor<const UGraph, const EdgeFilterMap>(graph, efm);

  }

  template <typename _UGraph, typename _DirectionMap>

  class DirUGraphAdaptorBase {

  public:

    typedef _UGraph Graph;

    typedef _DirectionMap DirectionMap;

    typedef typename _UGraph::Node Node;

    typedef typename _UGraph::UEdge Edge;

    void first(Node& i) const { graph->first(i); }

    void first(Edge& i) const { graph->first(i); }

    void firstIn(Edge& i, const Node& n) const {

      bool d;

      graph->firstInc(i, d, n);

      while (i != INVALID && d == (*direction)[i]) graph->nextInc(i, d);

    }

    void firstOut(Edge& i, const Node& n ) const { 

      bool d;

      graph->firstInc(i, d, n);

      while (i != INVALID && d != (*direction)[i]) graph->nextInc(i, d);

    }

    void next(Node& i) const { graph->next(i); }

    void next(Edge& i) const { graph->next(i); }

    void nextIn(Edge& i) const {

      bool d = !(*direction)[i];

      graph->nextInc(i, d);

      while (i != INVALID && d == (*direction)[i]) graph->nextInc(i, d);

    }

    void nextOut(Edge& i) const {

      bool d = (*direction)[i];

      graph->nextInc(i, d);

      while (i != INVALID && d != (*direction)[i]) graph->nextInc(i, d);

    }

    Node source(const Edge& e) const { 

      return (*direction)[e] ? graph->source(e) : graph->target(e); 

    }

    Node target(const Edge& e) const { 

      return (*direction)[e] ? graph->target(e) : graph->source(e); 

    }

    typedef NodeNumTagIndicator<Graph> NodeNumTag;

    int nodeNum() const { return graph->nodeNum(); }

    typedef EdgeNumTagIndicator<Graph> EdgeNumTag;

    int edgeNum() const { return graph->uEdgeNum(); }

    typedef FindEdgeTagIndicator<Graph> FindEdgeTag;

    Edge findEdge(const Node& source, const Node& target, 

		  const Edge& prev = INVALID) {

      Edge edge = prev;

      bool d = edge == INVALID ? true : (*direction)[edge];

      if (d) {

        edge = graph->findUEdge(source, target, edge);

        while (edge != INVALID && !(*direction)[edge]) {

          graph->findUEdge(source, target, edge);

        }

        if (edge != INVALID) return edge;

      }

      graph->findUEdge(target, source, edge);

      while (edge != INVALID && (*direction)[edge]) {

        graph->findUEdge(source, target, edge);

      }

      return edge;

    }

    Node addNode() const { 

      return Node(graph->addNode()); 

    }

    Edge addEdge(const Node& source, const Node& target) const {

      Edge edge = graph->addEdge(source, target);

      (*direction)[edge] = graph->source(edge) == source;

      return edge; 

    }

    void erase(const Node& i) const { graph->erase(i); }

    void erase(const Edge& i) const { graph->erase(i); }

    void clear() const { graph->clear(); }

    int id(const Node& v) const { return graph->id(v); }

    int id(const Edge& e) const { return graph->id(e); }

    void reverseEdge(const Edge& edge) {

      direction->set(edge, !(*direction)[edge]);

    }

    template <typename _Value>

    class NodeMap : public _UGraph::template NodeMap<_Value> {

    public:

      typedef typename _UGraph::template NodeMap<_Value> Parent;

      explicit NodeMap(const DirUGraphAdaptorBase& ga) 

	: Parent(*ga.graph) { }

      NodeMap(const DirUGraphAdaptorBase& ga, const _Value& value)

	: Parent(*ga.graph, value) { }

    };

    template <typename _Value>

    class EdgeMap : public _UGraph::template UEdgeMap<_Value> {

    public:

      typedef typename _UGraph::template UEdgeMap<_Value> Parent;

      explicit EdgeMap(const DirUGraphAdaptorBase& ga) 

	: Parent(*ga.graph) { }

      EdgeMap(const DirUGraphAdaptorBase& ga, const _Value& value)

	: Parent(*ga.graph, value) { }

    };

  protected:

    Graph* graph;

    DirectionMap* direction;

    void setDirectionMap(DirectionMap& _direction) {

      direction = &_direction;

    }

    void setGraph(Graph& _graph) {

      graph = &_graph;

    }

  };

  /// \ingroup graph_adaptors

  /// \brief A directed graph is made from a undirected graph by an adaptor

  ///

  /// This adaptor gives a direction for each uedge in the undirected graph.

  /// The direction of the edges stored in the DirectionMap. This map is

  /// a bool map on the undirected edges. If the uedge is mapped to true

  /// then the direction of the directed edge will be the same as the

  /// default direction of the uedge. The edges can be easily reverted

  /// by the reverseEdge member in the adaptor.  

  template<typename _Graph, 

           typename DirectionMap = typename _Graph::template UEdgeMap<bool> > 

  class DirUGraphAdaptor : 

    public GraphAdaptorExtender<

    DirUGraphAdaptorBase<_Graph, DirectionMap> > {

  public:

    typedef _Graph Graph;

    typedef GraphAdaptorExtender<

      DirUGraphAdaptorBase<_Graph, DirectionMap> > Parent;

  protected:

    DirUGraphAdaptor() { }

  public:

    DirUGraphAdaptor(_Graph& _graph, DirectionMap& _direction_map) { 

      setGraph(_graph);

      setDirectionMap(_direction_map);

    }

  };

  template<typename UGraph, typename DirectionMap>

  DirUGraphAdaptor<const UGraph, DirectionMap>

  dirUGraphAdaptor(const UGraph& graph, DirectionMap& dm) {

    return DirUGraphAdaptor<const UGraph, DirectionMap>(graph, dm);

  }

  template<typename UGraph, typename DirectionMap>

  DirUGraphAdaptor<const UGraph, const DirectionMap>

  dirUGraphAdaptor(const UGraph& graph, const DirectionMap& dm) {

    return DirUGraphAdaptor<const UGraph, const DirectionMap>(graph, dm);

  }

}

#endif