hamiltonian.H

# ifndef HAMILTONIAN_H
# define HAMILTONIAN_H

# include <tpl_graph.H>

namespace Aleph 
{

template <class GT, 
          class SN = Dft_Show_Node<GT>,
          class SA = Dft_Show_Arc<GT> >
class Test_Hamiltonian_Sufficiency
{
  SN & sn;
  SA & sa;

  bool test_graph(GT & g)
  {
    I(not g.is_digraph());

    const size_t & n = g.get_num_nodes();

    for (Node_Iterator<GT, SN> i(g, sn); i.has_curr(); i.next())
      {
        typename GT::Node * src = i.get_curr();
        const size_t & nsrc = g.get_num_arcs(src);

        i.next();

        for (Node_Iterator<GT, SN> j(i); j.has_curr(); j.next())
          if (nsrc + g.get_num_arcs(j.get_curr()) < n)
            return false;
      }

    return true;
  }

  bool test_digraph(GT & g)
  {
    I(g.is_digraph());

    g.reset_counter_nodes();

      // una primera pasada sobre arcos para contar grados de entrada
    for (Arc_Iterator<GT, SA> it(g, sa); it.has_curr(); it.next())
      NODE_COUNTER(it.get_tgt_node())++;

    const size_t & n = g.get_num_nodes();

    for (Node_Iterator<GT, SN> i(g, sn); i.has_curr(); i.next())
      {
        typename GT::Node * src = i.get_curr();

        for (Node_Iterator<GT, SN> j(g, sn); j.has_curr(); j.next())
          {
            typename GT::Node * tgt = j.get_current();

            if (src == tgt)
              continue;

            if (g.get_num_arcs(src) + NODE_COUNTER(tgt) >= n)
              continue;

                // buscar si hay arco src-->tgt
            bool terminate = true;
            for (Node_Arc_Iterator<GT, SA> it(src, sa); it.has_curr(); 
                 it.next())
              if (src == it.get_tgt_node())
                {
                  terminate = false;
                  break;
                }
            
            if (terminate)
              return false;
          }
      } 

    return true;
  }

public:

  Test_Hamiltonian_Sufficiency(SN && __sn = SN(), SA && __sa = SA()) 
    : sn(__sn), sa(__sa)
  {
    // empty 
  }

      // retorna true si el grafo o digrafo satisface la prueba de
      // suficiencia para ser hamiltoniano.
  bool operator () (GT & g)
  {
    if (g.is_digraph())
      return test_digraph(g);
    else
      return test_graph(g);
  }
};



} // end namespace Aleph

# endif // HAMILTONIAN_H