test/yosupo/tree/lca_link_cut.test.cpp

View this file on GitHub
Last update: 2025-10-29 12:51:11+07:00
Problem: https://judge.yosupo.jp/problem/lca

Depends on

Link Cut Tree (tree/link_cut.hpp)

Code

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

#include <bits/stdc++.h>
using namespace std;

#include "tree/link_cut.hpp"

struct node : lct_node<node> {
    void pull() {}
};

int main() {
    cin.tie(nullptr)->sync_with_stdio(false);
    int n, q;
    cin >> n >> q;
    LCT<node> lct(n);
    for (int i = 1; i < n; i++) {
        int p;
        cin >> p;
        lct.link(i, p);
    }
    lct.evert(0);
    while (q--) {
        int u, v;
        cin >> u >> v, cout << lct.lca(u, v) << "\n";
    }
}

#line 1 "test/yosupo/tree/lca_link_cut.test.cpp"
#define PROBLEM "https://judge.yosupo.jp/problem/lca"

#include <bits/stdc++.h>
using namespace std;

#line 2 "tree/link_cut.hpp"

template<class node> struct LCT {
    using ptr = node *;
    vector<node> nodes;

    LCT(int n = 0) : nodes(n) {}

    auto operator[](int i) { return &nodes[i]; }
    int id(ptr t) { return int(t - &nodes[0]); }

    void splay(int v) { splay(&nodes[v]); }
    void expose(int v) { expose(&nodes[v]); }
    void expose_path(int u, int v) { evert(u), expose(v); }
    void evert(int v) { evert(&nodes[v]); }
    void link(int v, int p) { link(&nodes[v], &nodes[p]); }
    void cut(int v, int p) { cut(&nodes[v], &nodes[p]); }
    int lca(int u, int v) {
        ptr l = lca(&nodes[u], &nodes[v]);
        return l ? id(l) : -1;
    }
    template<class... T> void set(int v, T &&...x) {
        expose(v), nodes[v].set(x...);
    }
    template<class... T> void add(int v, T &&...x) {
        expose(v), nodes[v].add(x...);
    }

    static void splay(ptr t) {
        for (t->_push(); state(t); rot(t)) {
            ptr p = t->p, g = p->p;
            if (g) g->_push();
            p->_push(), t->_push();
            if (state(p)) rot(state(t) == state(p) ? p : t);
        }
    }
    static ptr expose(ptr t) {
        ptr prv = 0;
        for (ptr cur = t; cur; cur = cur->p) {
            splay(cur);
            if (cur->r) cur->_add_light(cur->r);
            if (prv) cur->_sub_light(prv);
            attach(cur, 1, exchange(prv, cur));
        }
        splay(t);
        return prv;
    }
    static void evert(ptr t) { expose(t), t->_flip(); }
    static void link(ptr v, ptr p) {
        evert(v), expose(p);
        assert(!v->p && !p->r);
        attach(p, 1, v);
    }
    static void cut(ptr v, ptr p) {
        evert(p), expose(v);
        assert(!v->p && v->l == p);
        attach(v, 0, p->p = 0);
    }
    static ptr lca(ptr u, ptr v) {
        if (u == v) return u;
        expose(u);
        ptr l = expose(v);
        return u->p ? l : 0;
    }

    static ptr &ch(ptr t, bool d) { return d ? t->r : t->l; }
    static void attach(ptr p, int d, ptr c) {
        if (c) c->p = p;
        ch(p, d) = c, p->pull();
    }
    static int state(ptr t) {
        if (!t->p) return 0;
        return t == t->p->l ? -1 : t == t->p->r ? 1 : 0;
    }
    static void rot(ptr t) {
        ptr p = t->p, g = p->p;
        int d = state(t) == 1, dp = state(p);
        t->p = p->p;
        if (g) dp ? attach(g, dp == 1, t) : g->_change_light(p, t);
        attach(p, d, ch(t, !d));
        attach(t, !d, p);
    }
};

template<class node> struct lct_node {
    using ptr = node *;
    static ptr $(ptr t) {
        static node nil;
        return t ?: &nil;
    }
    ptr p = 0, l = 0, r = 0;
    bool rev = 0;

    node *as_derived() { return (node *)this; };
    void _push() {
        if (exchange(rev, 0)) $(l)->_flip(), $(r)->_flip();
        if constexpr (has_push<node>) as_derived()->push();
    }
    void _flip() {
        swap(l, r), rev ^= 1;
        if constexpr (has_flip<node>) as_derived()->flip();
    }

    void _add_light(ptr c) {
        if constexpr (has_add_light<node>) as_derived()->add_light(c);
    }
    void _sub_light(ptr c) {
        if constexpr (has_sub_light<node>) as_derived()->sub_light(c);
    }
    void _change_light(ptr prv, ptr cur) {
        if constexpr (has_change_light<node>)
            as_derived()->change_light(prv, cur);
    }

// rip compile time
#define CHECK(a) \
    template<class T, class = void> struct has_##a##_t : false_type {}; \
    template<class T> struct has_##a##_t<T, void_t<decltype(&T::a)>> : true_type {}; \
    template<class T> static constexpr bool has_##a = has_##a##_t<T>::value;
    CHECK(push) CHECK(flip)
    CHECK(add_light) CHECK(sub_light) CHECK(change_light)
#undef CHECK
};
#line 7 "test/yosupo/tree/lca_link_cut.test.cpp"

struct node : lct_node<node> {
    void pull() {}
};

int main() {
    cin.tie(nullptr)->sync_with_stdio(false);
    int n, q;
    cin >> n >> q;
    LCT<node> lct(n);
    for (int i = 1; i < n; i++) {
        int p;
        cin >> p;
        lct.link(i, p);
    }
    lct.evert(0);
    while (q--) {
        int u, v;
        cin >> u >> v, cout << lct.lca(u, v) << "\n";
    }
}

Test cases

Env	Name	Status	Elapsed	Memory
g++	almost_line_00	AC	488 ms	19 MB
g++	almost_line_01	AC	492 ms	19 MB
g++	binary_00	AC	748 ms	19 MB
g++	binary_01	AC	862 ms	19 MB
g++	binary_02	AC	816 ms	19 MB
g++	example_00	AC	5 ms	4 MB
g++	line_00	AC	375 ms	16 MB
g++	line_01	AC	394 ms	18 MB
g++	line_02	AC	272 ms	5 MB
g++	line_03	AC	59 ms	17 MB
g++	line_04	AC	104 ms	12 MB
g++	max_line_00	AC	485 ms	19 MB
g++	max_line_01	AC	474 ms	19 MB
g++	max_line_02	AC	477 ms	19 MB
g++	max_random_00	AC	568 ms	19 MB
g++	max_random_01	AC	548 ms	19 MB
g++	max_random_02	AC	632 ms	19 MB
g++	path_graph_root_centroid_00	AC	483 ms	19 MB
g++	path_graph_root_centroid_01	AC	501 ms	19 MB
g++	path_graph_root_centroid_02	AC	480 ms	19 MB
g++	random_00	AC	440 ms	16 MB
g++	random_01	AC	440 ms	18 MB
g++	random_02	AC	250 ms	5 MB
g++	random_03	AC	158 ms	17 MB
g++	random_04	AC	157 ms	12 MB
clang++	almost_line_00	AC	572 ms	19 MB
clang++	almost_line_01	AC	574 ms	19 MB
clang++	binary_00	AC	904 ms	19 MB
clang++	binary_01	AC	892 ms	19 MB
clang++	binary_02	AC	918 ms	19 MB
clang++	example_00	AC	5 ms	4 MB
clang++	line_00	AC	448 ms	16 MB
clang++	line_01	AC	457 ms	18 MB
clang++	line_02	AC	308 ms	5 MB
clang++	line_03	AC	63 ms	17 MB
clang++	line_04	AC	118 ms	12 MB
clang++	max_line_00	AC	558 ms	19 MB
clang++	max_line_01	AC	570 ms	19 MB
clang++	max_line_02	AC	594 ms	19 MB
clang++	max_random_00	AC	622 ms	19 MB
clang++	max_random_01	AC	631 ms	19 MB
clang++	max_random_02	AC	702 ms	19 MB
clang++	path_graph_root_centroid_00	AC	560 ms	19 MB
clang++	path_graph_root_centroid_01	AC	566 ms	19 MB
clang++	path_graph_root_centroid_02	AC	572 ms	19 MB
clang++	random_00	AC	513 ms	16 MB
clang++	random_01	AC	509 ms	18 MB
clang++	random_02	AC	280 ms	5 MB
clang++	random_03	AC	182 ms	17 MB
clang++	random_04	AC	182 ms	12 MB