連結成分の可換モノイド和を取得できる Union Find
(library/datastructure/union_find/union_find_component_sum.hpp)

• View this file on GitHub
• Last update: 2023-07-09 04:04:16+09:00
• Include: #include "library/datastructure/union_find/union_find_component_sum.hpp"

連結成分の可換モノイド和を取得できる Union Find

Depends on

• Union Find (library/datastructure/union_find/union_find.hpp)

Verified with

• test/src/datastructure/union_find/union_find_component_sum/abc183_f.test.cpp

Code

#ifndef SUISEN_UNION_FIND_COMPONENT_SUM
#define SUISEN_UNION_FIND_COMPONENT_SUM

#include "library/datastructure/union_find/union_find.hpp"

namespace suisen {
    template <typename T, void(*merge_data)(T&, T)>
    struct UnionFindComponentSum : UnionFind {
        UnionFindComponentSum() : UnionFindComponentSum(0) {}
        explicit UnionFindComponentSum(int n, const T &init_value = T{}) : UnionFindComponentSum(std::vector<T>(n, init_value)) {}
        explicit UnionFindComponentSum(const std::vector<T> &init_values) : UnionFind(init_values.size()), _sum(init_values) {}

        bool merge(int x, int y) {
            x = root(x), y = root(y);
            bool res = UnionFind::merge(x, y);
            if (res) {
                if (root(x) == y) std::swap(x, y);
                merge_data(_sum[x], std::move(_sum[y]));
            }
            return res;
        }
        const T& sum(int x) {
            return _sum[root(x)];
        }
    private:
        std::vector<T> _sum;
    };
} // namespace suisen


#endif // SUISEN_UNION_FIND_COMPONENT_SUM

#line 1 "library/datastructure/union_find/union_find_component_sum.hpp"



#line 1 "library/datastructure/union_find/union_find.hpp"



#include <algorithm>
#include <vector>

namespace suisen {
    struct UnionFind {
        UnionFind() = default;
        explicit UnionFind(int _n) : _n(_n), _dat(_n, -1) {}
        // Get the root of `x`. equivalent to `operator[](x)`
        int root(int x) {
            static std::vector<int> buf;
            while (_dat[x] >= 0) buf.push_back(x), x = _dat[x];
            while (buf.size()) _dat[buf.back()] = x, buf.pop_back();
            return x;
        }
        // Get the root of `x`. euivalent to `root(x)`
        int operator[](int x) {
            return root(x);
        }
        // Merge two vertices `x` and `y`.
        bool merge(int x, int y) {
            x = root(x), y = root(y);
            if (x == y) return false;
            if (_dat[x] > _dat[y]) std::swap(x, y);
            _dat[x] += _dat[y], _dat[y] = x;
            return true;
        }
        // Check if `x` and `y` belongs to the same connected component.
        bool same(int x, int y) {
            return root(x) == root(y);
        }
        // Get the size of connected componet to which `x` belongs.
        int size(int x) {
            return -_dat[root(x)];
        }
        // Get all of connected components.
        std::vector<std::vector<int>> groups() {
            std::vector<std::vector<int>> res(_n);
            for (int i = 0; i < _n; ++i) res[root(i)].push_back(i);
            res.erase(std::remove_if(res.begin(), res.end(), [](const auto& g) { return g.empty(); }), res.end());
            return res;
        }
    protected:
        int _n;
        std::vector<int> _dat;
    };
} // namespace suisen


#line 5 "library/datastructure/union_find/union_find_component_sum.hpp"

namespace suisen {
    template <typename T, void(*merge_data)(T&, T)>
    struct UnionFindComponentSum : UnionFind {
        UnionFindComponentSum() : UnionFindComponentSum(0) {}
        explicit UnionFindComponentSum(int n, const T &init_value = T{}) : UnionFindComponentSum(std::vector<T>(n, init_value)) {}
        explicit UnionFindComponentSum(const std::vector<T> &init_values) : UnionFind(init_values.size()), _sum(init_values) {}

        bool merge(int x, int y) {
            x = root(x), y = root(y);
            bool res = UnionFind::merge(x, y);
            if (res) {
                if (root(x) == y) std::swap(x, y);
                merge_data(_sum[x], std::move(_sum[y]));
            }
            return res;
        }
        const T& sum(int x) {
            return _sum[root(x)];
        }
    private:
        std::vector<T> _sum;
    };
} // namespace suisen

Back to top page