test/src/graph/directed_eulerian_graph/eulerian_trail_directed.test.cpp

• View this file on GitHub
• Last update: 2024-01-30 19:31:48+09:00
• Problem: https://judge.yosupo.jp/problem/eulerian_trail_directed

Depends on

• Directed Eulerian Graph (library/graph/directed_eulerian_graph.hpp)

Code

#define PROBLEM "https://judge.yosupo.jp/problem/eulerian_trail_directed"

#include <map>
#include <iostream>
#include "library/graph/directed_eulerian_graph.hpp"

int main() {
    std::ios::sync_with_stdio(false);
    std::cin.tie(nullptr);

    int t;
    std::cin >> t;
    for (int tc = 0; tc < t; ++tc) {
        int n, m;
        std::cin >> n >> m;
        std::vector<std::map<int, std::vector<int>>> mps(n);
        suisen::DirectedEulerianGraph g(n);
        for (int e = 0; e < m; ++e) {
            int u, v;
            std::cin >> u >> v;
            g.add_edge(u, v);
            mps[u][v].push_back(e);
        }
        auto p = g.eulerian_trail();
        if (p) {
            std::cout << "Yes\n";
            for (int i = 0; i < m; ++i) {
                std::cout << (*p)[i] << ' ';
            }
            std::cout << p->back() << '\n';
            for (int i = 0; i < m; ++i) {
                int u = (*p)[i], v = (*p)[i + 1];
                auto &es = mps[u][v];
                if (i) std::cout << ' ';
                std::cout << es.back();
                es.pop_back();
            }
            std::cout << '\n';
        } else {
            std::cout << "No\n";
        }
    }
}

#line 1 "test/src/graph/directed_eulerian_graph/eulerian_trail_directed.test.cpp"
#define PROBLEM "https://judge.yosupo.jp/problem/eulerian_trail_directed"

#include <map>
#include <iostream>
#line 1 "library/graph/directed_eulerian_graph.hpp"



#include <algorithm>
#include <cassert>
#include <optional>
#include <vector>

namespace suisen {
    struct DirectedEulerianGraph {
        DirectedEulerianGraph() = default;
        DirectedEulerianGraph(int n) : n(n), g(n), in_deg(n, 0) {}

        void add_edge(int u, int v) {
            g[u].push_back(v);
            ++in_deg[v];
        }

        std::optional<std::vector<int>> eulerian_circuit(int start = -1) {
            std::size_t edge_num = 0;
            std::vector<std::vector<bool>> used(n);
            for (int i = 0; i < n; ++i) {
                const int sz = g[i].size();
                edge_num += sz;
                used[i].resize(sz, false);
                if (in_deg[i] != sz) return std::nullopt;
            }
            if (start < 0) {
                start = 0;
                for (int i = 0; i < n; ++i) if (in_deg[i]) start = i;
            }
            std::vector<int> res;
            std::vector<std::size_t> iter(n);
            auto dfs = [&](auto dfs, int u) -> void {
                for (std::size_t& i = iter[u]; i < g[u].size(); ++i) {
                    if (used[u][i]) continue;
                    const int v = g[u][i];
                    used[u][i] = true;
                    dfs(dfs, v);
                }
                res.push_back(u);
            };
            dfs(dfs, start);
            std::reverse(res.begin(), res.end());
            if (res.size() != edge_num + 1) return std::nullopt;
            return res;
        }
        std::optional<std::vector<int>> eulerian_trail(bool different_endpoints = false) {
            int s = -1, t = -1, invalid = -1;
            for (int i = 0; i < n; ++i) {
                int out_deg = g[i].size();
                if (out_deg == in_deg[i] + 1) {
                    (s < 0 ? s : invalid) = i;
                } else if (out_deg == in_deg[i] - 1) {
                    (t < 0 ? t : invalid) = i;
                } else if (out_deg != in_deg[i]) {
                    invalid = i;
                }
            }
            if (not different_endpoints and s < 0 and t < 0 and invalid < 0) {
                return eulerian_circuit();
            }
            if (s < 0 or t < 0 or invalid >= 0) return std::nullopt;
            add_edge(t, s);
            auto opt_res = eulerian_circuit(s);
            if (not opt_res) return std::nullopt;
            auto &res = *opt_res;
            res.pop_back();
            // remove edge (t, s)
            g[t].pop_back();
            --in_deg[s];
            if (res.back() == t) return res;
            auto is_ts_edge = [&](int u, int v) {
                return u == t and v == s;
            };
            std::rotate(res.begin(), std::adjacent_find(res.begin(), res.end(), is_ts_edge) + 1, res.end());
            return std::move(res);
        }

        const std::vector<int>& operator[](int u) const {
            return g[u];
        }
    private:
        int n;
        std::vector<std::vector<int>> g;
        std::vector<int> in_deg;
    };
}


#line 6 "test/src/graph/directed_eulerian_graph/eulerian_trail_directed.test.cpp"

int main() {
    std::ios::sync_with_stdio(false);
    std::cin.tie(nullptr);

    int t;
    std::cin >> t;
    for (int tc = 0; tc < t; ++tc) {
        int n, m;
        std::cin >> n >> m;
        std::vector<std::map<int, std::vector<int>>> mps(n);
        suisen::DirectedEulerianGraph g(n);
        for (int e = 0; e < m; ++e) {
            int u, v;
            std::cin >> u >> v;
            g.add_edge(u, v);
            mps[u][v].push_back(e);
        }
        auto p = g.eulerian_trail();
        if (p) {
            std::cout << "Yes\n";
            for (int i = 0; i < m; ++i) {
                std::cout << (*p)[i] << ' ';
            }
            std::cout << p->back() << '\n';
            for (int i = 0; i < m; ++i) {
                int u = (*p)[i], v = (*p)[i + 1];
                auto &es = mps[u][v];
                if (i) std::cout << ' ';
                std::cout << es.back();
                es.pop_back();
            }
            std::cout << '\n';
        } else {
            std::cout << "No\n";
        }
    }
}

Back to top page