test/yosupo-judge/tree_diameter/main.test.cpp

• View this file on GitHub
• Last update: 2021-04-23 23:44:44+09:00
• Problem: https://judge.yosupo.jp/problem/tree_diameter
• Link: View error logs on GitHub Actions

Depends on

• Basic graph (Mylib/Graph/Template/graph.cpp)
• Tree diameter (Mylib/Graph/TreeUtils/tree_diameter.cpp)
• join function (Mylib/IO/join.cpp)

Code

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

#include <iostream>
#include "Mylib/Graph/Template/graph.cpp"
#include "Mylib/Graph/TreeUtils/tree_diameter.cpp"
#include "Mylib/IO/join.cpp"

namespace hl = haar_lib;

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

  int N;
  std::cin >> N;

  hl::tree<int64_t> tree(N);
  tree.read<0, false>(N - 1);

  auto [cost, path] = hl::tree_diameter(tree);

  std::cout
      << cost << " " << path.size() << "\n"
      << hl::join(path.begin(), path.end()) << "\n";

  return 0;
}

#line 1 "test/yosupo-judge/tree_diameter/main.test.cpp"
#define PROBLEM "https://judge.yosupo.jp/problem/tree_diameter"

#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/TreeUtils/tree_diameter.cpp"
#include <algorithm>
#include <stack>
#include <utility>
#line 7 "Mylib/Graph/TreeUtils/tree_diameter.cpp"

namespace haar_lib {
  template <typename T>
  std::pair<T, std::vector<int>> tree_diameter(const tree<T> &tr) {
    const int N = tr.size();

    std::vector<bool> check(N);
    std::vector<T> dp(N);
    std::stack<int> st;

    st.push(0);
    while (st.size()) {
      int i = st.top();
      st.pop();

      check[i] = true;

      for (const auto &e : tr[i]) {
        if (not check[e.to]) {
          dp[e.to] = dp[i] + e.cost;
          st.push(e.to);
        }
      }
    }

    const int u = std::max_element(dp.begin(), dp.end()) - dp.begin();

    dp.assign(N, 0);
    check.assign(N, false);

    std::vector<int> prev(N);

    st.push(u);
    while (st.size()) {
      int i = st.top();
      st.pop();

      check[i] = true;

      for (const auto &e : tr[i]) {
        if (not check[e.to]) {
          dp[e.to] = dp[i] + e.cost;
          st.push(e.to);
          prev[e.to] = i;
        }
      }
    }

    const int v = std::max_element(dp.begin(), dp.end()) - dp.begin();

    std::vector<int> ret;

    int cur = v;
    while (1) {
      ret.push_back(cur);
      if (cur == u) break;
      cur = prev[cur];
    }

    return std::make_pair(dp[v], ret);
  }
}  // namespace haar_lib
#line 3 "Mylib/IO/join.cpp"
#include <sstream>
#include <string>

namespace haar_lib {
  template <typename Iter>
  std::string join(Iter first, Iter last, std::string delim = " ") {
    std::stringstream s;

    for (auto it = first; it != last; ++it) {
      if (it != first) s << delim;
      s << *it;
    }

    return s.str();
  }
}  // namespace haar_lib
#line 7 "test/yosupo-judge/tree_diameter/main.test.cpp"

namespace hl = haar_lib;

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

  int N;
  std::cin >> N;

  hl::tree<int64_t> tree(N);
  tree.read<0, false>(N - 1);

  auto [cost, path] = hl::tree_diameter(tree);

  std::cout
      << cost << " " << path.size() << "\n"
      << hl::join(path.begin(), path.end()) << "\n";

  return 0;
}

Back to top page