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

Depends on

Code

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

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

using S = int;
S op(S a, S b) { return a ^ b; }
S e() { return 0; }

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

    vector<int> C(N);
    for (int i = 0; i < N; i++)
    {
        cin >> C[i];
    }

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

    atcoder::segtree<S, op, e> seg(N);
    for (int i = 0; i < N; i++)
    {
        seg.set(hld.get_pos(i), C[i]);
    }

    for (int q = 0; q < Q; q++)
    {
        int T, x, y;
        cin >> T >> x >> y;
        x--;
        if (T == 1)
        {
            int xi = hld.get_pos(x);
            seg.set(xi, seg.get(xi) ^ y);
        }

        if (T == 2)
        {
            auto [l, r] = hld.get_subtree_range(x);
            cout << seg.prod(l, r) << endl;
        }
    }

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

#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-1641.cpp"
#include <atcoder/segtree>
#line 6 "test/Tree/heavy_light_decomposition/yukicoder-1641.cpp"
using namespace std;

using S = int;
S op(S a, S b) { return a ^ b; }
S e() { return 0; }

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

    vector<int> C(N);
    for (int i = 0; i < N; i++)
    {
        cin >> C[i];
    }

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

    atcoder::segtree<S, op, e> seg(N);
    for (int i = 0; i < N; i++)
    {
        seg.set(hld.get_pos(i), C[i]);
    }

    for (int q = 0; q < Q; q++)
    {
        int T, x, y;
        cin >> T >> x >> y;
        x--;
        if (T == 1)
        {
            int xi = hld.get_pos(x);
            seg.set(xi, seg.get(xi) ^ y);
        }

        if (T == 2)
        {
            auto [l, r] = hld.get_subtree_range(x);
            cout << seg.prod(l, r) << endl;
        }
    }

    return 0;
};

Test cases

Env Name Status Elapsed Memory
g++ 01_sample_01.txt :heavy_check_mark: AC 5 ms 3 MB
g++ 01_sample_02.txt :heavy_check_mark: AC 4 ms 3 MB
g++ 01_sample_03.txt :heavy_check_mark: AC 4 ms 3 MB
g++ 02_random_small_1.txt :heavy_check_mark: AC 5 ms 3 MB
g++ 02_random_small_10.txt :heavy_check_mark: AC 4 ms 3 MB
g++ 02_random_small_2.txt :heavy_check_mark: AC 5 ms 3 MB
g++ 02_random_small_3.txt :heavy_check_mark: AC 5 ms 3 MB
g++ 02_random_small_4.txt :heavy_check_mark: AC 5 ms 3 MB
g++ 02_random_small_5.txt :heavy_check_mark: AC 4 ms 3 MB
g++ 02_random_small_6.txt :heavy_check_mark: AC 5 ms 3 MB
g++ 02_random_small_7.txt :heavy_check_mark: AC 5 ms 3 MB
g++ 02_random_small_8.txt :heavy_check_mark: AC 5 ms 3 MB
g++ 02_random_small_9.txt :heavy_check_mark: AC 5 ms 3 MB
g++ 03_random_max_1.txt :heavy_check_mark: AC 289 ms 30 MB
g++ 03_random_max_2.txt :heavy_check_mark: AC 234 ms 30 MB
g++ 04_random_all_1.txt :heavy_check_mark: AC 8 ms 4 MB
g++ 04_random_all_2.txt :heavy_check_mark: AC 16 ms 4 MB
g++ 04_random_all_3.txt :heavy_check_mark: AC 18 ms 4 MB
g++ 04_random_all_4.txt :heavy_check_mark: AC 11 ms 4 MB
g++ 04_random_all_5.txt :heavy_check_mark: AC 11 ms 4 MB
g++ 05_killer_01.txt :heavy_check_mark: AC 166 ms 13 MB
Back to top page