Submission #839358

#TimeUsernameProblemLanguageResultExecution timeMemory
839358t6twotwoTwo Currencies (JOI23_currencies)C++17
100 / 100
2957 ms97844 KiB
#include <bits/stdc++.h>
using namespace std;
using ll = long long;
struct node {
    node *l, *r;
    ll s;
    node() : s(0) {
        l = r = nullptr;
    }
    node(ll v) : s(v) {
        l = r = nullptr;
    }
    node(node *u, node *v) : l(u), r(v) {
        s = l->s + r->s;
    }
};
node *bld(int l, int r) {
    if (l + 1 == r) {
        return new node();
    }
    int m = (l + r) / 2;
    return new node(bld(l, m), bld(m, r));
}
node *add(node *p, int l, int r, int i, int v) {
    if (l + 1 == r) {
        return new node(p->s + v);
    }
    int m = (l + r) / 2;
    if (i < m) {
        return new node(add(p->l, l, m, i, v), p->r);
    } else {
        return new node(p->l, add(p->r, m, r, i, v));
    }
}
ll qry(node *p, int l, int r, int L, int R) {
    if (R <= l || r <= L) {
        return 0;
    }
    if (L <= l && r <= R) {
        return p->s;
    }
    int m = (l + r) / 2;
    return qry(p->l, l, m, L, R) + qry(p->r, m, r, L, R);
}
int main() {
    ios::sync_with_stdio(false);
    cin.tie(nullptr);
    int N, M, Q;
    cin >> N >> M >> Q;
    vector<pair<int, int>> edges(N - 1);
    vector<vector<int>> adj(N);
    for (int i = 0; i < N - 1; i++) {
        int x, y;
        cin >> x >> y;
        x--, y--;
        edges[i] = {x, y};
        adj[x].push_back(i);
        adj[y].push_back(i);
    }
    vector<pair<int, int>> ch(M);
    for (int i = 0; i < M; i++) {
        int p, c;
        cin >> p >> c;
        ch[i] = {c, p - 1};
    }
    sort(ch.begin(), ch.end());
    vector<int> t(N - 1);
    {
        auto dfs = [&](auto dfs, int x, int p) -> void {
            for (int z : adj[x]) {
                auto [u, v] = edges[z];
                int y = x ^ u ^ v;
                if (y != p) {
                    t[z] = y;
                    dfs(dfs, y, x);
                }
            }
        };
        dfs(dfs, 0, -1);
    }
    for (int i = 0; i < N; i++) {
        for (int &j : adj[i]) {
            auto [x, y] = edges[j];
            j = i ^ x ^ y;
        }
    }
    vector<int> dep(N), par(N, -1), sz(N);
    auto init = [&](auto init, int x) -> void {
        if (par[x] != -1) {
            adj[x].erase(find(adj[x].begin(), adj[x].end(), par[x]));
        }
        sz[x] = 1;
        for (int &y : adj[x]) {
            par[y] = x;
            dep[y] = dep[x] + 1;
            init(init, y);
            sz[x] += sz[y];
            if (sz[y] > sz[adj[x][0]]) {
                swap(y, adj[x][0]);
            }
        }
    };
    init(init, 0);
    vector<int> top(N), in(N);
    int timer = 0;
    auto dfs = [&](auto dfs, int x) -> void {
        in[x] = timer++;
        for (int y : adj[x]) {
            top[y] = adj[x][0] == y ? top[x] : y;
            dfs(dfs, y);
        }
    };
    dfs(dfs, 0);
    vector<node*> q{bld(0, N)};
    for (int i = 0; i < M; i++) {
        auto [c, p] = ch[i];
        q.push_back(add(q.back(), 0, N, in[t[p]], c));
    }
    vector<int> S(Q), T(Q), X(Q);
    vector<ll> Y(Q);
    for (int i = 0; i < Q; i++) {
        cin >> S[i] >> T[i] >> X[i] >> Y[i];
        S[i]--, T[i]--;
    }
    vector<int> a(Q);
    for (int i = 0; i < Q; i++) {
        int lo = 0, hi = M;
        while (lo < hi) {
            int mi = (lo + hi + 1) / 2;
            int x = S[i], y = T[i]; ll s = 0;
            while (top[x] != top[y]) {
                if (dep[top[x]] < dep[top[y]]) {
                    swap(x, y);
                }
                s += qry(q[mi], 0, N, in[top[x]], in[x] + 1);
                x = par[top[x]];
            }
            if (dep[x] < dep[y]) {
                swap(x, y);
            }
            s += qry(q[mi], 0, N, in[y] + 1, in[x] + 1);
            if (s <= Y[i]) {
                lo = mi;
            } else {
                hi = mi - 1;
            }
        }
        a[i] = lo;
    }
    vector<ll> bit(N + 1);
    auto add = [&](int i, int v) {
        for (i++; i <= N; i += i & -i) {
            bit[i] += v;
        }
    };
    auto sum = [&](int i) {
        ll s = 0;
        for (; i; i -= i & -i) {
            s += bit[i];
        }
        return s;
    };
    auto get = [&](int x, int y) {
        ll s = 0;
        while (top[x] != top[y]) {
            if (dep[top[x]] < dep[top[y]]) {
                swap(x, y);
            }
            s += sum(in[x] + 1) - sum(in[top[x]]);
            x = par[top[x]];
        }
        if (dep[x] < dep[y]) {
            swap(x, y);
        }
        s += sum(in[x] + 1) - sum(in[y] + 1);
        return s;
    };
    vector<ll> ans(Q);
    vector<vector<int>> f(M + 1);
    for (int i = 0; i < Q; i++) {
        f[a[i]].push_back(i);
    }
    bit.assign(N + 1, 0);
    for (int i = M; i >= 0; i--) {
        if (i < M) {
            auto [c, p] = ch[i];
            add(in[t[p]], 1);
        }
        for (int j : f[i]) {
            ans[j] = max(-1LL, X[j] - get(S[j], T[j]));
        }
    }
    for (int i = 0; i < Q; i++) {
        cout << ans[i] << "\n";
    }
    return 6/22;
}
#Verdict Execution timeMemoryGrader output
Fetching results...
#Verdict Execution timeMemoryGrader output
Fetching results...
#Verdict Execution timeMemoryGrader output
Fetching results...
#Verdict Execution timeMemoryGrader output
Fetching results...