|
| 1 | +#ifndef WEILYCODER_TARJAN_HPP |
| 2 | +#define WEILYCODER_TARJAN_HPP |
| 3 | + |
| 4 | +#include <cstddef> |
| 5 | +#include <stack> |
| 6 | +#include <utility> |
| 7 | +#include <vector> |
| 8 | + |
| 9 | +namespace weilycoder { |
| 10 | +template <typename ptr_t = size_t> struct TwoEdgeConnectedComponents { |
| 11 | + ptr_t dfs_time = 0, edge_time = 1; |
| 12 | + |
| 13 | + std::vector<bool> in_stack; |
| 14 | + std::stack<ptr_t> stk; |
| 15 | + std::vector<ptr_t> dfn, low; |
| 16 | + std::vector<std::vector<std::pair<ptr_t, ptr_t>>> graph; |
| 17 | + |
| 18 | + std::vector<std::vector<ptr_t>> eccs; |
| 19 | + |
| 20 | + TwoEdgeConnectedComponents(ptr_t n) |
| 21 | + : in_stack(n, false), dfn(n, 0), low(n, 0), graph(n) {} |
| 22 | + |
| 23 | + void add_edge(ptr_t u, ptr_t v) { |
| 24 | + graph[u].emplace_back(v, edge_time); |
| 25 | + graph[v].emplace_back(u, edge_time); |
| 26 | + edge_time++; |
| 27 | + } |
| 28 | + |
| 29 | + void tarjan(ptr_t u, ptr_t parent_edge) { |
| 30 | + dfn[u] = low[u] = ++dfs_time; |
| 31 | + stk.push(u); |
| 32 | + in_stack[u] = true; |
| 33 | + |
| 34 | + for (const auto &[v, edge_id] : graph[u]) { |
| 35 | + if (edge_id == parent_edge) |
| 36 | + continue; |
| 37 | + if (!dfn[v]) |
| 38 | + tarjan(v, edge_id), low[u] = std::min(low[u], low[v]); |
| 39 | + else if (in_stack[v]) |
| 40 | + low[u] = std::min(low[u], dfn[v]); |
| 41 | + } |
| 42 | + |
| 43 | + if (dfn[u] == low[u]) { |
| 44 | + eccs.emplace_back(); |
| 45 | + eccs.back().push_back(u); |
| 46 | + in_stack[u] = false; |
| 47 | + while (stk.top() != u) |
| 48 | + eccs.back().push_back(stk.top()), in_stack[stk.top()] = false, stk.pop(); |
| 49 | + stk.pop(); |
| 50 | + } |
| 51 | + } |
| 52 | + |
| 53 | + void solve() { |
| 54 | + for (size_t i = 0; i < graph.size(); ++i) |
| 55 | + if (!dfn[i]) |
| 56 | + tarjan(i, 0); |
| 57 | + } |
| 58 | +}; |
| 59 | +} // namespace weilycoder |
| 60 | + |
| 61 | +#endif |
0 commit comments