/*

   igraph library.

   Copyright (C) 2005-2020  The igraph development team

   This program is free software; you can redistribute it and/or modify

   it under the terms of the GNU General Public License as published by

   the Free Software Foundation; either version 2 of the License, or

   (at your option) any later version.

   This program is distributed in the hope that it will be useful,

   but WITHOUT ANY WARRANTY; without even the implied warranty of

   MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the

   GNU General Public License for more details.

   You should have received a copy of the GNU General Public License

   along with this program; if not, write to the Free Software

   Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA

   02110-1301 USA

*/

#include "igraph_eulerian.h"

#include "igraph_adjlist.h"

#include "igraph_bitset.h"

#include "igraph_interface.h"

#include "igraph_components.h"

#include "igraph_stack.h"

/**

 * \section about_eulerian

 *

 * <para>These functions calculate whether an Eulerian path or cycle exists

 * and if so, can find them.</para>

 */

/* solution adapted from https://www.geeksforgeeks.org/eulerian-path-and-circuit/

The function returns one of the following values

has_path is set to 1 if a path exists, 0 otherwise

has_cycle is set to 1 if a cycle exists, 0 otherwise

*/

static igraph_error_t igraph_i_is_eulerian_undirected(

        const igraph_t *graph, igraph_bool_t *has_path, igraph_bool_t *has_cycle, igraph_int_t *start_of_path) {

    igraph_int_t odd;

    igraph_vector_int_t degree;

    igraph_vector_int_t csize;

    /* boolean vector to mark singletons: */

    igraph_vector_int_t nonsingleton;

    igraph_int_t i, n, vsize;

    igraph_int_t cluster_count;

    /* number of self-looping singletons: */

    igraph_int_t es;

    /* will be set to 1 if there are non-isolated vertices, otherwise 0: */

    igraph_int_t ens;

    n = igraph_vcount(graph);

    if (igraph_ecount(graph) == 0 || n <= 1) {

        start_of_path = 0; /* in case the graph has one vertex with self-loops */

        *has_path = true;

        *has_cycle = true;

        return  IGRAPH_SUCCESS;

    }

    /* check for connectedness, but singletons are special since they affect

     * the Eulerian nature only if there is a self-loop AND another edge

     * somewhere else in the graph */

    IGRAPH_VECTOR_INT_INIT_FINALLY(&csize, 0);

    IGRAPH_CHECK(igraph_connected_components(graph, NULL, &csize, NULL, IGRAPH_WEAK));

    cluster_count = 0;

    vsize = igraph_vector_int_size(&csize);

    for (i = 0; i < vsize; i++) {

        if (VECTOR(csize)[i] > 1) {

            cluster_count++;

            if (cluster_count > 1) {

                /* disconnected edges, they'll never reach each other */

                *has_path = false;

                *has_cycle = false;

                igraph_vector_int_destroy(&csize);

                IGRAPH_FINALLY_CLEAN(1);

                return IGRAPH_SUCCESS;

            }

        }

    }

    igraph_vector_int_destroy(&csize);

    IGRAPH_FINALLY_CLEAN(1);

    /* the graph is connected except for singletons */

    /* find singletons (including those with self-loops) */

    IGRAPH_VECTOR_INT_INIT_FINALLY(&nonsingleton, 0);

    IGRAPH_CHECK(igraph_degree(graph, &nonsingleton, igraph_vss_all(), IGRAPH_ALL, IGRAPH_NO_LOOPS));

    /* check the degrees for odd/even:

     * - >= 2 odd means no cycle (1 odd is impossible)

     * - > 2 odd means no path

     * plus there are a few corner cases with singletons

     */

    IGRAPH_VECTOR_INT_INIT_FINALLY(&degree, 0);

    IGRAPH_CHECK(igraph_degree(graph, &degree, igraph_vss_all(), IGRAPH_ALL, IGRAPH_LOOPS));

    odd = 0;

    es = 0;

    ens = 0;

    for (i = 0; i < n; i++) {

        igraph_int_t deg = VECTOR(degree)[i];

        /* Eulerian is about edges, so skip free vertices */

        if (deg == 0) continue;

        if (!VECTOR(nonsingleton)[i]) {

            /* singleton with self loops */

            es++;

        } else {

            /* at least one non-singleton */

            ens = 1;

            /* note: self-loops count for two (in and out) */

            if (deg % 2) odd++;

        }

        if (es + ens > 1) {

            /* 2+ singletons with self loops or singleton with self-loops and

             * 1+ edges in the non-singleton part of the graph. */

            *has_path = false;

            *has_cycle = false;

            igraph_vector_int_destroy(&nonsingleton);

            igraph_vector_int_destroy(&degree);

            IGRAPH_FINALLY_CLEAN(2);

            return IGRAPH_SUCCESS;

        }

    }

    igraph_vector_int_destroy(&nonsingleton);

    IGRAPH_FINALLY_CLEAN(1);

    /* this is the usual algorithm on the connected part of the graph */

    if (odd > 2) {

        *has_path = false;

        *has_cycle = false;

    } else if (odd == 2) {

        *has_path = true;

        *has_cycle = false;

    } else {

        *has_path = true;

        *has_cycle = true;

    }

    /* set start of path if there is one but there is no cycle */

    /* note: we cannot do this in the previous loop because at that time we are

     * not sure yet if a path exists */

    for (i = 0; i < n; i++) {

        if ((*has_cycle && VECTOR(degree)[i] > 0) || (!*has_cycle && VECTOR(degree)[i] %2 == 1)) {

            *start_of_path = i;

            break;

        }

    }

    igraph_vector_int_destroy(&degree);

    IGRAPH_FINALLY_CLEAN(1);

    return IGRAPH_SUCCESS;

}

static igraph_error_t igraph_i_is_eulerian_directed(

        const igraph_t *graph, igraph_bool_t *has_path, igraph_bool_t *has_cycle, igraph_int_t *start_of_path) {

    igraph_int_t incoming_excess, outgoing_excess, n;

    igraph_int_t i, vsize;

    igraph_int_t cluster_count;

    igraph_vector_int_t out_degree, in_degree;

    igraph_vector_int_t csize;

    /* boolean vector to mark singletons: */

    igraph_vector_int_t nonsingleton;

    /* number of self-looping singletons: */

    igraph_int_t es;

    /* will be set to 1 if there are non-isolated vertices, otherwise 0: */

    igraph_int_t ens;

    n = igraph_vcount(graph);

    if (igraph_ecount(graph) == 0 || n <= 1) {

        start_of_path = 0; /* in case the graph has one vertex with self-loops */

        *has_path = true;

        *has_cycle = true;

        return IGRAPH_SUCCESS;

    }

    incoming_excess = 0;

    outgoing_excess = 0;

    /* check for weak connectedness, but singletons are special since they affect

     * the Eulerian nature only if there is a self-loop AND another edge

     * somewhere else in the graph */

    IGRAPH_VECTOR_INT_INIT_FINALLY(&csize, 0);

    IGRAPH_CHECK(igraph_connected_components(graph, NULL, &csize, NULL, IGRAPH_WEAK));

    cluster_count = 0;

    vsize = igraph_vector_int_size(&csize);

    for (i = 0; i < vsize; i++) {

        if (VECTOR(csize)[i] > 1) {

            cluster_count++;

            if (cluster_count > 1) {

                /* weakly disconnected edges, they'll never reach each other */

                *has_path = false;

                *has_cycle = false;

                igraph_vector_int_destroy(&csize);

                IGRAPH_FINALLY_CLEAN(1);

                return IGRAPH_SUCCESS;

            }

        }

    }

    igraph_vector_int_destroy(&csize);

    IGRAPH_FINALLY_CLEAN(1);

    /* the graph is weakly connected except for singletons */

    /* find the singletons (including those with self-loops) */

    IGRAPH_VECTOR_INT_INIT_FINALLY(&nonsingleton, 0);

    IGRAPH_CHECK(igraph_degree(graph, &nonsingleton, igraph_vss_all(), IGRAPH_ALL, IGRAPH_NO_LOOPS));

    /* checking if no. of incoming edges == outgoing edges

     * plus there are a few corner cases with singletons */

    IGRAPH_VECTOR_INT_INIT_FINALLY(&out_degree, 0);

    IGRAPH_CHECK(igraph_degree(graph, &out_degree, igraph_vss_all(), IGRAPH_OUT, IGRAPH_LOOPS));

    IGRAPH_VECTOR_INT_INIT_FINALLY(&in_degree, 0);

    IGRAPH_CHECK(igraph_degree(graph, &in_degree, igraph_vss_all(), IGRAPH_IN, IGRAPH_LOOPS));

    es = 0;

    ens = 0;

    *start_of_path = -1;

    for (i = 0; i < n; i++) {

        igraph_int_t degin = VECTOR(in_degree)[i];

        igraph_int_t degout = VECTOR(out_degree)[i];

        /* Eulerian is about edges, so skip free vertices */

        if (degin + degout == 0) continue;

        if (!VECTOR(nonsingleton)[i]) {

            /* singleton with self loops */

            es++;

            /* if we ever want a path, it has to be this self-loop */

            *start_of_path = i;

        } else {

            /* at least one non-singleton */

            ens = 1;

        }

        if (es + ens > 1) {

            /* 2+ singletons with self loops or singleton with self-loops and

             * 1+ edges in the non-singleton part of the graph. */

            *has_path = false;

            *has_cycle = false;

            igraph_vector_int_destroy(&nonsingleton);

            igraph_vector_int_destroy(&in_degree);

            igraph_vector_int_destroy(&out_degree);

            IGRAPH_FINALLY_CLEAN(3);

            return IGRAPH_SUCCESS;

        }

        /* as long as we have perfect balance, you can start

         * anywhere with an edge */

        if (*start_of_path == -1 && incoming_excess == 0 && outgoing_excess == 0) {

            *start_of_path = i;

        }

        /* same in and out (including self-loops, even in singletons) */

        if (degin == degout) {

            continue;

        }

        /* non-singleton, in != out */

        if (degin > degout) {

            incoming_excess += degin - degout;

        } else {

            outgoing_excess += degout - degin;

            if (outgoing_excess == 1) {

                *start_of_path = i;

            }

        }

        /* too much imbalance, either of the following:

         * 1. 1+ vertices have 2+ in/out

         * 2. 2+ nodes have 1+ in/out */

        if (incoming_excess > 1 || outgoing_excess > 1) {

            *has_path = false;

            *has_cycle = false;

            igraph_vector_int_destroy(&nonsingleton);

            igraph_vector_int_destroy(&in_degree);

            igraph_vector_int_destroy(&out_degree);

            IGRAPH_FINALLY_CLEAN(3);

            return IGRAPH_SUCCESS;

        }

    }

    *has_path = true;

    /* perfect edge balance -> strong connectivity */

    *has_cycle = (incoming_excess == 0) && (outgoing_excess == 0);

    /* either way, the start was set already */

    igraph_vector_int_destroy(&nonsingleton);

    igraph_vector_int_destroy(&in_degree);

    igraph_vector_int_destroy(&out_degree);

    IGRAPH_FINALLY_CLEAN(3);

    return IGRAPH_SUCCESS;

}

/**

 * \ingroup Eulerian

 * \function igraph_is_eulerian

 * \brief Checks whether an Eulerian path or cycle exists.

 *

 * An Eulerian path traverses each edge of the graph precisely once. A closed

 * Eulerian path is referred to as an Eulerian cycle.

 *

 * \param graph The graph object.

 * \param has_path Pointer to a Boolean, will be set to true if an Eulerian path exists.

 *    Must not be \c NULL.

 * \param has_cycle Pointer to a Boolean, will be set to true if an Eulerian cycle exists.

 *    Must not be \c NULL.

 * \return Error code:

 *         \c IGRAPH_ENOMEM, not enough memory for temporary data.

 *

 * Time complexity: O(|V|+|E|), the number of vertices plus the number of edges.

 *

 */

igraph_error_t igraph_is_eulerian(const igraph_t *graph, igraph_bool_t *has_path, igraph_bool_t *has_cycle) {

    igraph_int_t start_of_path = 0;

    if (igraph_is_directed(graph)) {

        IGRAPH_CHECK(igraph_i_is_eulerian_directed(graph, has_path, has_cycle, &start_of_path));

    } else {

        IGRAPH_CHECK(igraph_i_is_eulerian_undirected(graph, has_path, has_cycle, &start_of_path));

    }

    return IGRAPH_SUCCESS;

}

static igraph_error_t igraph_i_eulerian_path_undirected(

        const igraph_t *graph, igraph_vector_int_t *edge_res, igraph_vector_int_t *vertex_res,

        igraph_int_t start_of_path) {

    igraph_int_t curr;

    igraph_int_t n, m;

    igraph_inclist_t il;

    igraph_stack_int_t path, tracker, edge_tracker, edge_path;

    igraph_bitset_t visited_list;

    igraph_vector_int_t degree;

    n = igraph_vcount(graph);

    m = igraph_ecount(graph);

    if (edge_res) {

        igraph_vector_int_clear(edge_res);

    }

    if (vertex_res) {

        igraph_vector_int_clear(vertex_res);

    }

    if (m == 0 || n == 0) {

        return IGRAPH_SUCCESS;

    }

    IGRAPH_VECTOR_INT_INIT_FINALLY(&degree, 0);

    IGRAPH_CHECK(igraph_degree(graph, &degree, igraph_vss_all(), IGRAPH_ALL, IGRAPH_LOOPS));

    IGRAPH_STACK_INT_INIT_FINALLY(&path, n);

    IGRAPH_STACK_INT_INIT_FINALLY(&tracker, n);

    IGRAPH_STACK_INT_INIT_FINALLY(&edge_path, n);

    IGRAPH_STACK_INT_INIT_FINALLY(&edge_tracker, n);

    IGRAPH_BITSET_INIT_FINALLY(&visited_list, m);

    IGRAPH_CHECK(igraph_stack_int_push(&tracker, start_of_path));

    IGRAPH_CHECK(igraph_inclist_init(graph, &il, IGRAPH_OUT, IGRAPH_LOOPS_ONCE));

    IGRAPH_FINALLY(igraph_inclist_destroy, &il);

    curr = start_of_path;

    while (!igraph_stack_int_empty(&tracker)) {

        if (VECTOR(degree)[curr] != 0) {

            igraph_vector_int_t *incedges;

            igraph_int_t nc, edge = -1;

            igraph_int_t j, next;

            IGRAPH_CHECK(igraph_stack_int_push(&tracker, curr));

            incedges = igraph_inclist_get(&il, curr);

            nc = igraph_vector_int_size(incedges);

            IGRAPH_ASSERT(nc > 0);

            for (j = 0; j < nc; j++) {

                edge = VECTOR(*incedges)[j];

                if (!IGRAPH_BIT_TEST(visited_list, edge)) {

                    break;

                }

            }

            next = IGRAPH_OTHER(graph, edge, curr);

            IGRAPH_CHECK(igraph_stack_int_push(&edge_tracker, edge));

            /* remove edge here */

            VECTOR(degree)[curr]--;

            VECTOR(degree)[next]--;

            IGRAPH_BIT_SET(visited_list, edge);

            curr = next;

        } else { /* back track to find remaining circuit */

            igraph_int_t curr_e;

            IGRAPH_CHECK(igraph_stack_int_push(&path, curr));

            curr = igraph_stack_int_pop(&tracker);

            if (!igraph_stack_int_empty(&edge_tracker)) {

                curr_e = igraph_stack_int_pop(&edge_tracker);

                IGRAPH_CHECK(igraph_stack_int_push(&edge_path, curr_e));

            }

        }

    }

    if (edge_res) {

        IGRAPH_CHECK(igraph_vector_int_reserve(edge_res, m));

        while (!igraph_stack_int_empty(&edge_path)) {

            IGRAPH_CHECK(igraph_vector_int_push_back(edge_res, igraph_stack_int_pop(&edge_path)));

        }

    }

    if (vertex_res) {

        IGRAPH_CHECK(igraph_vector_int_reserve(vertex_res, m+1));

        while (!igraph_stack_int_empty(&path)) {

            IGRAPH_CHECK(igraph_vector_int_push_back(vertex_res, igraph_stack_int_pop(&path)));

        }

    }

    igraph_stack_int_destroy(&path);

    igraph_stack_int_destroy(&tracker);

    igraph_stack_int_destroy(&edge_path);

    igraph_stack_int_destroy(&edge_tracker);

    igraph_bitset_destroy(&visited_list);

    igraph_inclist_destroy(&il);

    igraph_vector_int_destroy(&degree);

    IGRAPH_FINALLY_CLEAN(7);

    return IGRAPH_SUCCESS;

}

/* solution adapted from https://www.geeksforgeeks.org/hierholzers-algorithm-directed-graph/ */

static igraph_error_t igraph_i_eulerian_path_directed(

        const igraph_t *graph, igraph_vector_int_t *edge_res, igraph_vector_int_t *vertex_res,

        igraph_int_t start_of_path) {

    igraph_int_t curr;

    igraph_int_t n, m;

    igraph_inclist_t il;

    igraph_stack_int_t path, tracker, edge_tracker, edge_path;

    igraph_bitset_t visited_list;

    igraph_vector_int_t remaining_out_edges;

    n = igraph_vcount(graph);

    m = igraph_ecount(graph);

    if (edge_res) {

        igraph_vector_int_clear(edge_res);

    }

    if (vertex_res) {

        igraph_vector_int_clear(vertex_res);

    }

    if (m == 0 || n == 0) {

        return IGRAPH_SUCCESS;

    }

    IGRAPH_STACK_INT_INIT_FINALLY(&path, n);

    IGRAPH_STACK_INT_INIT_FINALLY(&tracker, n);

    IGRAPH_STACK_INT_INIT_FINALLY(&edge_path, n);

    IGRAPH_STACK_INT_INIT_FINALLY(&edge_tracker, n);

    IGRAPH_BITSET_INIT_FINALLY(&visited_list, m);

    IGRAPH_CHECK(igraph_stack_int_push(&tracker, start_of_path));

    IGRAPH_CHECK(igraph_inclist_init(graph, &il, IGRAPH_OUT, IGRAPH_LOOPS_ONCE));

    IGRAPH_FINALLY(igraph_inclist_destroy, &il);

    IGRAPH_VECTOR_INT_INIT_FINALLY(&remaining_out_edges, 0);

    IGRAPH_CHECK(igraph_degree(graph, &remaining_out_edges, igraph_vss_all(), IGRAPH_OUT, IGRAPH_LOOPS));

    curr = start_of_path;

    while (!igraph_stack_int_empty(&tracker)) {

        if (VECTOR(remaining_out_edges)[curr] != 0) {

            igraph_vector_int_t *incedges;

            igraph_int_t nc, edge = -1;

            igraph_int_t j, next;

            IGRAPH_CHECK(igraph_stack_int_push(&tracker, curr));

            incedges = igraph_inclist_get(&il, curr);

            nc = igraph_vector_int_size(incedges);

            IGRAPH_ASSERT(nc > 0);

            for (j = 0; j < nc; j++) {

                edge = VECTOR(*incedges)[j];

                if (!IGRAPH_BIT_TEST(visited_list, edge)) {

                    break;

                }

            }

            next = IGRAPH_TO(graph, edge);

            IGRAPH_CHECK(igraph_stack_int_push(&edge_tracker, edge));

            /* remove edge here */

            VECTOR(remaining_out_edges)[curr]--;

            IGRAPH_BIT_SET(visited_list, edge);

            curr = next;

        } else { /* back track to find remaining circuit */

            igraph_int_t curr_e;

            IGRAPH_CHECK(igraph_stack_int_push(&path, curr));

            curr = igraph_stack_int_pop(&tracker);

            if (!igraph_stack_int_empty(&edge_tracker)) {

                curr_e = igraph_stack_int_pop(&edge_tracker);

                IGRAPH_CHECK(igraph_stack_int_push(&edge_path, curr_e));

            }

        }

    }

    if (edge_res) {

        IGRAPH_CHECK(igraph_vector_int_reserve(edge_res, m));

        while (!igraph_stack_int_empty(&edge_path)) {

            IGRAPH_CHECK(igraph_vector_int_push_back(edge_res, igraph_stack_int_pop(&edge_path)));

        }

    }

    if (vertex_res) {

        IGRAPH_CHECK(igraph_vector_int_reserve(vertex_res, m+1));

        while (!igraph_stack_int_empty(&path)) {

            IGRAPH_CHECK(igraph_vector_int_push_back(vertex_res, igraph_stack_int_pop(&path)));

        }

    }

    igraph_stack_int_destroy(&path);

    igraph_stack_int_destroy(&tracker);

    igraph_stack_int_destroy(&edge_path);

    igraph_stack_int_destroy(&edge_tracker);

    igraph_bitset_destroy(&visited_list);

    igraph_inclist_destroy(&il);

    igraph_vector_int_destroy(&remaining_out_edges);

    IGRAPH_FINALLY_CLEAN(7);

    return IGRAPH_SUCCESS;

}

/**

 * \ingroup Eulerian

 * \function igraph_eulerian_cycle

 * \brief Finds an Eulerian cycle.

 *

 * Finds an Eulerian cycle, if it exists. An Eulerian cycle is a closed path

 * that traverses each edge precisely once.

 *

 * </para><para>

 * If the graph has no edges, a zero-length cycle is returned.

 *

 * </para><para>

 * This function uses Hierholzer's algorithm.

 *

 * \param graph The graph object.

 * \param edge_res Pointer to an initialised vector. The indices of edges

 *                 belonging to the cycle will be stored here. May be \c NULL

 *                 if it is not needed by the caller.

 * \param vertex_res Pointer to an initialised vector. The indices of vertices

 *                   belonging to the cycle will be stored here. The first and

 *                   last vertex in the vector will be the same. May be \c NULL

 *                   if it is not needed by the caller.

 * \return Error code:

 *        \clist

 *        \cli IGRAPH_ENOMEM

 *           not enough memory for temporary data.

 *        \cli IGRAPH_ENOSOL

 *           graph does not have an Eulerian cycle.

 *        \endclist

 *

 * Time complexity: O(|V|+|E|), the number of vertices plus the number of edges.

 *

 */

igraph_error_t igraph_eulerian_cycle(

        const igraph_t *graph, igraph_vector_int_t *edge_res, igraph_vector_int_t *vertex_res) {

    igraph_bool_t has_cycle;

    igraph_bool_t has_path;

    igraph_int_t start_of_path = 0;

    if (igraph_is_directed(graph)) {

        IGRAPH_CHECK(igraph_i_is_eulerian_directed(graph, &has_path, &has_cycle, &start_of_path));

        if (!has_cycle) {

            IGRAPH_ERROR("The graph does not have an Eulerian cycle.", IGRAPH_ENOSOL);

        }

        IGRAPH_CHECK(igraph_i_eulerian_path_directed(graph, edge_res, vertex_res, start_of_path));

    } else {

        IGRAPH_CHECK(igraph_i_is_eulerian_undirected(graph, &has_path, &has_cycle, &start_of_path));

        if (!has_cycle) {

            IGRAPH_ERROR("The graph does not have an Eulerian cycle.", IGRAPH_ENOSOL);

        }

        IGRAPH_CHECK(igraph_i_eulerian_path_undirected(graph, edge_res, vertex_res, start_of_path));

    }

    return IGRAPH_SUCCESS;

}

/**

 * \ingroup Eulerian

 * \function igraph_eulerian_path

 * \brief Finds an Eulerian path.

 *

 * Finds an Eulerian path, if it exists. An Eulerian path traverses

 * each edge precisely once.

 *

 * </para><para>

 * If the graph has no edges, a zero-length path is returned.

 *

 * </para><para>

 * This function uses Hierholzer's algorithm.

 *

 * \param graph The graph object.

 * \param edge_res Pointer to an initialised vector. The indices of edges

 *                 belonging to the path will be stored here. May be \c NULL

 *                 if it is not needed by the caller.

 * \param vertex_res Pointer to an initialised vector. The indices of vertices

 *                   belonging to the path will be stored here. May be \c NULL

 *                   if it is not needed by the caller.

 * \return Error code:

 *        \clist

 *        \cli IGRAPH_ENOMEM

 *           not enough memory for temporary data.

 *        \cli IGRAPH_ENOSOL

 *           graph does not have an Eulerian path.

 *        \endclist

 *

 * Time complexity: O(|V|+|E|), the number of vertices plus the number of edges.

 *

 */

igraph_error_t igraph_eulerian_path(

        const igraph_t *graph, igraph_vector_int_t *edge_res, igraph_vector_int_t *vertex_res) {

    igraph_bool_t has_cycle;

    igraph_bool_t has_path;

    igraph_int_t start_of_path = 0;

    if (igraph_is_directed(graph)) {

        IGRAPH_CHECK(igraph_i_is_eulerian_directed(graph, &has_path, &has_cycle, &start_of_path));

        if (!has_path) {

            IGRAPH_ERROR("The graph does not have an Eulerian path.", IGRAPH_ENOSOL);

        }

        IGRAPH_CHECK(igraph_i_eulerian_path_directed(graph, edge_res, vertex_res, start_of_path));

    } else {

        IGRAPH_CHECK(igraph_i_is_eulerian_undirected(graph, &has_path, &has_cycle, &start_of_path));

        if (!has_path) {

            IGRAPH_ERROR("The graph does not have an Eulerian path.", IGRAPH_ENOSOL);

        }

        IGRAPH_CHECK(igraph_i_eulerian_path_undirected(graph, edge_res, vertex_res, start_of_path));

    }

    return IGRAPH_SUCCESS;

}