diff --git a/AcWing/2568/2568.cpp b/AcWing/2568/2568.cpp new file mode 100644 index 00000000..49ef499c --- /dev/null +++ b/AcWing/2568/2568.cpp @@ -0,0 +1,153 @@ +#include + +using namespace std; + +int n, m, a[100005], x, y, op, u, v, k; +vector g[100005]; +int cnt, dep[100005], fa[100005], siz[100005], son[100005], id[100005], top[100005], w[100005]; + +struct node { + int l, r; + long long s, d; + + node() + : l(0), r(0), s(0), d(0) {} + node(int _l, int _r) + : l(_l), r(_r), s(0), d(0) {} +} tr[100005 << 2]; + +void dfs1(int u, int father) { + dep[u] = dep[father] + 1; + fa[u] = father; + siz[u] = 1; + for (int v : g[u]) { + if (v == father) continue; + dfs1(v, u); + siz[u] += siz[v]; + if (siz[son[u]] < siz[v]) son[u] = v; + } +} + +void dfs2(int u, int t) { + id[u] = ++cnt; + w[cnt] = a[u]; + top[u] = t; + if (!son[u]) return; + dfs2(son[u], t); + for (int v : g[u]) { + if (v == fa[u]) continue; + if (v == son[u]) continue; + dfs2(v, v); + } +} + +void pushup(int u) { + tr[u].s = tr[u << 1].s + tr[u << 1 | 1].s; +} + +void pushdown(int u) { + tr[u << 1].s += (tr[u << 1].r - tr[u << 1].l + 1) * tr[u].d; + tr[u << 1].d += tr[u].d; + tr[u << 1 | 1].s += (tr[u << 1 | 1].r - tr[u << 1 | 1].l + 1) * tr[u].d; + tr[u << 1 | 1].d += tr[u].d; + tr[u].d = 0; +} + +void build(int u, int l, int r) { + tr[u] = node(l, r); + if (l == r) { + tr[u].s = w[r]; + return; + } + int mid = l + r >> 1; + build(u << 1, l, mid); + build(u << 1 | 1, mid + 1, r); + pushup(u); +} + +void modify(int u, int l, int r, int k) { + if (l <= tr[u].l && tr[u].r <= r) { + tr[u].s += (tr[u].r - tr[u].l + 1) * k; + tr[u].d += k; + return; + } + int mid = tr[u].l + tr[u].r >> 1; + pushdown(u); + if (l <= mid) modify(u << 1, l, r, k); + if (r > mid) modify(u << 1 | 1, l, r, k); + pushup(u); +} + +long long query(int u, int l, int r) { + if (l <= tr[u].l && tr[u].r <= r) { + return tr[u].s; + } + int mid = tr[u].l + tr[u].r >> 1; + long long sum = 0; + pushdown(u); + if (l <= mid) sum += query(u << 1, l, r); + if (r > mid) sum += query(u << 1 | 1, l, r); + return sum; +} + +void modify_path(int u, int v, int k) { + while (top[u] != top[v]) { + if (dep[top[u]] < dep[top[v]]) swap(u, v); + modify(1, id[top[u]], id[u], k); + u = fa[top[u]]; + } + if (dep[u] < dep[v]) swap(u, v); + modify(1, id[v], id[u], k); +} + +long long query_path(int u, int v) { + long long sum = 0; + while (top[u] != top[v]) { + if (dep[top[u]] < dep[top[v]]) swap(u, v); + sum += query(1, id[top[u]], id[u]); + u = fa[top[u]]; + } + if (dep[u] < dep[v]) swap(u, v); + sum += query(1, id[v], id[u]); + return sum; +} + +void modify_tree(int u, int k) { + modify(1, id[u], id[u] + siz[u] - 1, k); +} + +long long query_tree(int u) { + return query(1, id[u], id[u] + siz[u] - 1); +} + +int main() { + cin >> n; + for (int i = 1; i <= n; i++) { + cin >> a[i]; + } + for (int i = 1; i < n; i++) { + cin >> x >> y; + g[x].push_back(y); + g[y].push_back(x); + } + dfs1(1, 0); + dfs2(1, 1); + build(1, 1, n); + cin >> m; + for (int i = 1; i <= m; i++) { + cin >> op >> u; + if (op == 1) { + cin >> v >> k; + modify_path(u, v, k); + } else if (op == 2) { + cin >> k; + modify_tree(u, k); + } else if (op == 3) { + cin >> v; + cout << query_path(u, v) << endl; + } else { + cout << query_tree(u) << endl; + } + } + return 0; +}