:heavy_check_mark: test/Tree/heavy_light_decomposition/yukicoder-399.cpp

Depends on

Code

// competitive-verifier: PROBLEM https://yukicoder.me/problems/no/399

#include "Tree/heavy_light_decomposition.hpp"
#include <atcoder/lazysegtree>
#include <bits/stdc++.h>
using namespace std;

struct S
{
    long long value;
    int size;
};
using F = long long;

S op(S a, S b) { return {a.value + b.value, a.size + b.size}; }
S e() { return {0, 0}; }
S mapping(F f, S x) { return {x.value + f * x.size, x.size}; }
F composition(F f, F g) { return f + g; }
F id() { return 0; }

int main()
{
    int N;
    cin >> N;

    HeavyLightDecomposition hld(N);

    for (int i = 0; i < N - 1; i++)
    {
        int u, v;
        cin >> u >> v;
        hld.add_edge(u - 1, v - 1);
    }

    atcoder::lazy_segtree<S, op, e, F, mapping, composition, id> seg(vector<S>(N, {1, 1}));

    int Q;
    cin >> Q;

    long long ans = 0;

    for (int i = 0; i < Q; i++)
    {
        int A, B;
        cin >> A >> B;
        A--, B--;

        auto ranges = hld.get_path_ranges(A, B);
        S sum = e();
        for (auto [l, r, _] : ranges)
        {
            sum = op(sum, seg.prod(l, r));
            seg.apply(l, r, 1);
        }

        ans += sum.value;
    }

    cout << ans << endl;

    return 0;
}
#line 1 "test/Tree/heavy_light_decomposition/yukicoder-399.cpp"
// competitive-verifier: PROBLEM https://yukicoder.me/problems/no/399

#line 1 "Tree/heavy_light_decomposition.hpp"
#include <bits/stdc++.h>
using namespace std;

class HeavyLightDecomposition
{
public:
    HeavyLightDecomposition(int n)
        : n(n), edges(n), parent(n, -1), depth(n, 0), heavy(n, -1), head(n), pos(n), sz(n), built(false)
    {
    }

    void add_edge(int u, int v)
    {
        assert(0 <= u && u < n);
        assert(0 <= v && v < n);
        edges[u].push_back(v);
        edges[v].push_back(u);
    }

    void build(int root = 0)
    {
        assert(0 <= root && root < n);
        if (built)
        {
            return;
        }

        parent[root] = -1;
        depth[root] = 0;
        dfs_size(root);

        int time = 0;
        decompose(root, root, time);
        built = true;
    }

    vector<tuple<int, int, bool>> get_path_ranges(int from, int to)
    {
        assert(0 <= from && from < n);
        assert(0 <= to && to < n);
        ensure_built();

        list<tuple<int, int, bool>> from_ranges, to_ranges;

        while (head[from] != head[to])
        {
            if (depth[head[to]] < depth[head[from]])
            {
                from_ranges.push_back({pos[head[from]], pos[from] + 1, true});
                from = parent[head[from]];
            }
            else
            {
                to_ranges.push_front({pos[head[to]], pos[to] + 1, false});
                to = parent[head[to]];
            }
        }

        if (depth[to] < depth[from])
        {
            to_ranges.push_front({pos[to], pos[from] + 1, true});
        }
        else
        {
            from_ranges.push_back({pos[from], pos[to] + 1, false});
        }

        vector<tuple<int, int, bool>> result;
        result.insert(result.end(), from_ranges.begin(), from_ranges.end());
        result.insert(result.end(), to_ranges.begin(), to_ranges.end());

        return result;
    }

    array<int, 2> get_subtree_range(int u)
    {
        assert(0 <= u && u < n);

        ensure_built();
        return {pos[u], pos[u] + sz[u]};
    }

    int lca(int u, int v)
    {
        assert(0 <= u && u < n);
        assert(0 <= v && v < n);

        ensure_built();
        while (head[u] != head[v])
        {
            if (depth[head[u]] > depth[head[v]])
            {
                swap(u, v);
            }
            v = parent[head[v]];
        }
        return depth[u] < depth[v] ? u : v;
    }

    int get_pos(int v)
    {
        assert(0 <= v && v < n);

        ensure_built();
        return pos[v];
    }

    int get_parent(int v)
    {
        assert(0 <= v && v < n);

        ensure_built();
        return parent[v];
    }

    int get_subtree_size(int u)
    {
        assert(0 <= u && u < n);

        ensure_built();
        return sz[u];
    }

private:
    int n;
    vector<vector<int>> edges;
    vector<int> parent, depth, heavy, head, pos, sz;
    bool built;

    void ensure_built()
    {
        if (!built)
        {
            build();
        }
    }

    int dfs_size(int v)
    {
        sz[v] = 1;
        int max_size = 0;
        for (int u : edges[v])
        {
            if (u == parent[v])
            {
                continue;
            }

            parent[u] = v;
            depth[u] = depth[v] + 1;
            int sub = dfs_size(u);
            sz[v] += sub;
            if (max_size < sub)
            {
                max_size = sub;
                heavy[v] = u;
            }
        }

        return sz[v];
    }

    void decompose(int v, int h, int &time)
    {
        head[v] = h;
        pos[v] = time++;
        if (heavy[v] != -1)
        {
            decompose(heavy[v], h, time);
        }

        for (int u : edges[v])
        {
            if (u == parent[v] || u == heavy[v])
            {
                continue;
            }

            decompose(u, u, time);
        }
    }
};
#line 4 "test/Tree/heavy_light_decomposition/yukicoder-399.cpp"
#include <atcoder/lazysegtree>
#line 6 "test/Tree/heavy_light_decomposition/yukicoder-399.cpp"
using namespace std;

struct S
{
    long long value;
    int size;
};
using F = long long;

S op(S a, S b) { return {a.value + b.value, a.size + b.size}; }
S e() { return {0, 0}; }
S mapping(F f, S x) { return {x.value + f * x.size, x.size}; }
F composition(F f, F g) { return f + g; }
F id() { return 0; }

int main()
{
    int N;
    cin >> N;

    HeavyLightDecomposition hld(N);

    for (int i = 0; i < N - 1; i++)
    {
        int u, v;
        cin >> u >> v;
        hld.add_edge(u - 1, v - 1);
    }

    atcoder::lazy_segtree<S, op, e, F, mapping, composition, id> seg(vector<S>(N, {1, 1}));

    int Q;
    cin >> Q;

    long long ans = 0;

    for (int i = 0; i < Q; i++)
    {
        int A, B;
        cin >> A >> B;
        A--, B--;

        auto ranges = hld.get_path_ranges(A, B);
        S sum = e();
        for (auto [l, r, _] : ranges)
        {
            sum = op(sum, seg.prod(l, r));
            seg.apply(l, r, 1);
        }

        ans += sum.value;
    }

    cout << ans << endl;

    return 0;
}

Test cases

Env Name Status Elapsed Memory
g++ 0_1.txt :heavy_check_mark: AC 5 ms 3 MB
g++ 0_2.txt :heavy_check_mark: AC 4 ms 3 MB
g++ 1_0.txt :heavy_check_mark: AC 4 ms 3 MB
g++ 1_1.txt :heavy_check_mark: AC 4 ms 3 MB
g++ 1_2.txt :heavy_check_mark: AC 6 ms 4 MB
g++ 1_3.txt :heavy_check_mark: AC 32 ms 5 MB
g++ 1_4_1.txt :heavy_check_mark: AC 428 ms 18 MB
g++ 1_4_2.txt :heavy_check_mark: AC 415 ms 18 MB
g++ 1_5_1.txt :heavy_check_mark: AC 431 ms 18 MB
g++ 1_5_2.txt :heavy_check_mark: AC 432 ms 18 MB
g++ 2_0_1.txt :heavy_check_mark: AC 8 ms 4 MB
g++ 2_1.txt :heavy_check_mark: AC 27 ms 5 MB
g++ 2_2_1.txt :heavy_check_mark: AC 334 ms 18 MB
g++ 2_3_1.txt :heavy_check_mark: AC 332 ms 18 MB
g++ 3_0_1.txt :heavy_check_mark: AC 124 ms 34 MB
g++ 3_1_1.txt :heavy_check_mark: AC 150 ms 34 MB
g++ 3_2_1.txt :heavy_check_mark: AC 216 ms 25 MB
g++ system_test1.txt :heavy_check_mark: AC 446 ms 18 MB
g++ system_test2.txt :heavy_check_mark: AC 423 ms 18 MB
Back to top page