結果

問題 No.2296 Union Path Query (Hard)
コンテスト
ユーザー t98slider
提出日時 2023-03-20 20:30:57
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 596 ms / 7,000 ms
コード長 6,857 bytes
コンパイル時間 2,246 ms
コンパイル使用メモリ 202,036 KB
実行使用メモリ 42,384 KB
最終ジャッジ日時 2024-11-23 04:50:18
合計ジャッジ時間 22,052 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge2
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 4
other AC * 45
権限があれば一括ダウンロードができます

ソースコード

diff # 

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

class LinkCutTree {
    struct Node {
        Node *l, *r, *p;
        int key;
        long long sum;
        bool rev;

        explicit Node() : key(0), sum(0), rev(false),
                          l(nullptr), r(nullptr), p(nullptr) {}

        bool is_root() const { return not p or (p->l != this and p->r != this); }
    };

    void update(Node* t) {
        t->sum = t->key;
        if(t->l) t->sum += t->l->sum;
        if(t->r) t->sum += t->r->sum;
    }

    void rotate_right(Node* t) {
        Node *x = t->p, *y = x->p;
        if((x->l = t->r)) t->r->p = x;
        t->r = x, x->p = t;
        update(x), update(t);
        if((t->p = y)) {
            if(y->l == x) y->l = t;
            if(y->r == x) y->r = t;
            update(y);
        }
    }

    void rotate_left(Node* t) {
        Node *x = t->p, *y = x->p;
        if((x->r = t->l)) t->l->p = x;
        t->l = x, x->p = t;
        update(x), update(t);
        if((t->p = y)) {
            if(y->l == x) y->l = t;
            if(y->r == x) y->r = t;
            update(y);
        }
    }

    void toggle(Node* t) {
        std::swap(t->l, t->r);
        t->rev ^= true;
    }

    void push(Node* t) {
        if(t->rev) {
            if(t->l) toggle(t->l);
            if(t->r) toggle(t->r);
            t->rev = false;
        }
    }

    void splay(Node *t) {
        push(t);
        while(not t->is_root()) {
            Node* q = t->p;
            if(q->is_root()) {
                push(q), push(t);
                q->l == t ? rotate_right(t) : rotate_left(t);
            } else {
                Node* r = q->p;
                push(r), push(q), push(t);
                if(r->l == q) {
                if(q->l == t) rotate_right(q), rotate_right(t);
                else rotate_left(t), rotate_right(t);
                } else {
                if(q->r == t) rotate_left(q), rotate_left(t);
                else rotate_right(t), rotate_left(t);
                }
            }
        }
    }

    std::vector<Node*> Node_info;

    public:
    int N;
    LinkCutTree(int n) : N(n), Node_info(n) {
        for(int i = 0; i < N; i++){
            Node_info[i] = new Node();
        }
    }

    //vと根との間を Heavy-edge で繋げてvを木の根にする
    Node* expose(int v) {
        Node* t = Node_info[v];
        Node* rp = nullptr;
        for(Node* cur = t; cur; cur = cur->p) {
            splay(cur);
            cur->r = rp;
            update(cur);
            rp = cur;
        }
        splay(t);
        return rp;
    }

    //vをもともの木の根にする
    void evert(int v) {
        expose(v);
        toggle(Node_info[v]);
        push(Node_info[v]);
    }

    void link(int _child, int _parent) {
        Node *child = Node_info[_child], *parent = Node_info[_parent];
        evert(_child), expose(_parent);
        child->p = parent;
        parent->r = child;
        update(parent);
    }

    void cut(int _child) {
        Node *child = Node_info[_child];
        expose(_child);
        auto *parent = child->l;
        child->l = nullptr;
        parent->p = nullptr;
    }

    bool is_connected(int u, int v) {
        expose(u), expose(v);
        return Node_info[u] == Node_info[v] or Node_info[u]->p;
    }

    const long long &query(int u, int v) {
        evert(u);
        expose(v);
        return Node_info[v]->sum;
    }

    void set(int p, int v) {
        expose(p);
        Node_info[p]->key = v;
        update(Node_info[p]);
    }
};

struct dsu {
    public:
    int csz;
    dsu() : _n(0) {}
    dsu(int n) : _n(n), csz(n), parent_or_size(n, -1) {}

    int merge(int a, int b) {
        assert(0 <= a && a < _n);
        assert(0 <= b && b < _n);
        int x = leader(a), y = leader(b);
        if (x == y) return x;
        if (-parent_or_size[x] < -parent_or_size[y]) std::swap(x, y);
        csz--;
        parent_or_size[x] += parent_or_size[y];
        parent_or_size[y] = x;
        return x;
    }

    bool same(int a, int b) {
        assert(0 <= a && a < _n);
        assert(0 <= b && b < _n);
        return leader(a) == leader(b);
    }

    int leader(int a) {
        assert(0 <= a && a < _n);
        if (parent_or_size[a] < 0) return a;
        return parent_or_size[a] = leader(parent_or_size[a]);
    }

    int size(int a) {
        assert(0 <= a && a < _n);
        return -parent_or_size[leader(a)];
    }

    std::vector<std::vector<int>> groups() {
        std::vector<int> leader_buf(_n), group_size(_n);
        for (int i = 0; i < _n; i++) {
            leader_buf[i] = leader(i);
            group_size[leader_buf[i]]++;
        }
        std::vector<std::vector<int>> result(_n);
        for (int i = 0; i < _n; i++) {
            result[i].reserve(group_size[i]);
        }
        for (int i = 0; i < _n; i++) {
            result[leader_buf[i]].push_back(i);
        }
        result.erase(
            std::remove_if(result.begin(), result.end(),
                           [&](const std::vector<int>& v) { return v.empty(); }),
            result.end());
        return result;
    }

    private:
    int _n;
    // root node: -1 * component size
    // otherwise: parent
    std::vector<int> parent_or_size;
};

int main() {
    ios::sync_with_stdio(false);
    cin.tie(0);
    int N, Q;
    ll X;
    cin >> N >> X >> Q;
    LinkCutTree LCT(N + Q);

    dsu uf(N);
    vector<ll> dia(N);
    vector<array<int,2>> sv(N);
    for(int i = 0; i < N; i++) sv[i] = {i, i};

    auto update_diameter = [&](int x, int y){
        ll mx;
        int u, v;
        if(dia[x] >= dia[y]){
            mx = dia[x], u = sv[x][0], v = sv[x][1];
        }else{
            mx = dia[y], u = sv[y][0], v = sv[y][1];
        }
        for(int i = 0; i < 2; i++){
            for(int j = 0; j < 2; j++){
                ll d = LCT.query(sv[x][i], sv[y][j]);
                if(d > mx){
                    mx = d, u = sv[x][i], v = sv[y][j];
                }
            }
        }
        int r = uf.merge(x, y);
        dia[r] = mx;
        sv[r] = {u, v};
    };

    int type, v, u, w;
    for(int i = 0; i < Q; i++){
        cin >> type;
        if(type == 1){
            cin >> v >> w;
            LCT.link(v, i + N);
            LCT.link(X, i + N);
            LCT.set(i + N, w);
            update_diameter(uf.leader(v), uf.leader(X));
        }else if(type == 2){
            cin >> v >> u;
            if(!uf.same(v, u)){
                cout << -1 << '\n';
                continue;
            }
            ll d = LCT.query(v, u);
            (X += d) %= N;
            cout << d << '\n';
        }else if(type == 3){
            cin >> v;
            cout << dia[uf.leader(v)] << '\n';
        }else{
            cin >> v;
            X += v;
            if(X >= N) X -= N;
        }
    }
}
0