test/aoj/DPL_2_A/main.test.cpp

• View this file on GitHub
• Last update: 2021-04-23 23:44:44+09:00
• Problem: http://judge.u-aizu.ac.jp/onlinejudge/description.jsp?id=DPL_2_A
• Link: View error logs on GitHub Actions

Depends on

• Basic graph (Mylib/Graph/Template/graph.cpp)
• Travelling salesman problem (Mylib/Graph/travelling_salesman_problem.cpp)

Code

#define PROBLEM "http://judge.u-aizu.ac.jp/onlinejudge/description.jsp?id=DPL_2_A"

#include <iostream>
#include "Mylib/Graph/Template/graph.cpp"
#include "Mylib/Graph/travelling_salesman_problem.cpp"

namespace hl = haar_lib;

int main() {
  int V, E;
  std::cin >> V >> E;

  hl::graph<int> g(V);
  g.read<0>(E);

  std::cout << hl::travelling_salesman_problem(g, 0).value_or(-1) << std::endl;

  return 0;
}

#line 1 "test/aoj/DPL_2_A/main.test.cpp"
#define PROBLEM "http://judge.u-aizu.ac.jp/onlinejudge/description.jsp?id=DPL_2_A"

#include <iostream>
#line 3 "Mylib/Graph/Template/graph.cpp"
#include <vector>

namespace haar_lib {
  template <typename T>
  struct edge {
    int from, to;
    T cost;
    int index = -1;
    edge() {}
    edge(int from, int to, T cost) : from(from), to(to), cost(cost) {}
    edge(int from, int to, T cost, int index) : from(from), to(to), cost(cost), index(index) {}
  };

  template <typename T>
  struct graph {
    using weight_type = T;
    using edge_type   = edge<T>;

    std::vector<std::vector<edge<T>>> data;

    auto& operator[](size_t i) { return data[i]; }
    const auto& operator[](size_t i) const { return data[i]; }

    auto begin() const { return data.begin(); }
    auto end() const { return data.end(); }

    graph() {}
    graph(int N) : data(N) {}

    bool empty() const { return data.empty(); }
    int size() const { return data.size(); }

    void add_edge(int i, int j, T w, int index = -1) {
      data[i].emplace_back(i, j, w, index);
    }

    void add_undirected(int i, int j, T w, int index = -1) {
      add_edge(i, j, w, index);
      add_edge(j, i, w, index);
    }

    template <size_t I, bool DIRECTED = true, bool WEIGHTED = true>
    void read(int M) {
      for (int i = 0; i < M; ++i) {
        int u, v;
        std::cin >> u >> v;
        u -= I;
        v -= I;
        T w = 1;
        if (WEIGHTED) std::cin >> w;
        if (DIRECTED)
          add_edge(u, v, w, i);
        else
          add_undirected(u, v, w, i);
      }
    }
  };

  template <typename T>
  using tree = graph<T>;
}  // namespace haar_lib
#line 2 "Mylib/Graph/travelling_salesman_problem.cpp"
#include <algorithm>
#include <optional>
#line 6 "Mylib/Graph/travelling_salesman_problem.cpp"

namespace haar_lib {
  template <typename T>
  std::optional<T> travelling_salesman_problem(const graph<T> &g, int src) {
    const int n = g.size();

    std::vector<std::vector<std::optional<T>>> dp(n, std::vector<std::optional<T>>(1 << n));

    for (auto &e : g[src]) {
      if (not dp[e.to][1 << e.to]) {
        dp[e.to][1 << e.to] = e.cost;
      } else {
        dp[e.to][1 << e.to] = std::min(*dp[e.to][1 << e.to], e.cost);
      }
    }

    for (int s = 1; s < (1 << n); ++s) {
      for (int i = 0; i < n; ++i) {
        if (not(s & (1 << i))) continue;

        for (auto &e : g[i]) {
          if (s & (1 << e.to)) continue;

          if (dp[i][s]) {
            if (not dp[e.to][s | (1 << e.to)]) {
              dp[e.to][s | (1 << e.to)] = *dp[i][s] + e.cost;
            } else {
              dp[e.to][s | (1 << e.to)] = std::min(*dp[e.to][s | (1 << e.to)], *dp[i][s] + e.cost);
            }
          }
        }
      }
    }

    return dp[src][(1 << n) - 1];
  }
}  // namespace haar_lib
#line 6 "test/aoj/DPL_2_A/main.test.cpp"

namespace hl = haar_lib;

int main() {
  int V, E;
  std::cin >> V >> E;

  hl::graph<int> g(V);
  g.read<0>(E);

  std::cout << hl::travelling_salesman_problem(g, 0).value_or(-1) << std::endl;

  return 0;
}

Back to top page