头像

wt20

$\href{https://www.acwing.com/blog/content/14829/}{个人主页}$




离线:3小时前


最近来访(7080)
用户头像
Alinase
用户头像
Ethanyyc
用户头像
ㅤ煜
用户头像
月薪上W就开摆
用户头像
垫底抽風
用户头像
WA自由人
用户头像
LeRoi
用户头像
xinxixue6
用户头像
不再坠入星河
用户头像
Rickyの水果摊
用户头像
lscll
用户头像
東_7
用户头像
有机物
用户头像
天行_闻天
用户头像
Clock_Paradox
用户头像
zhyou
用户头像
种花家的兔兔
用户头像
zedchou
用户头像
PinkStrokes
用户头像
yxc的小迷弟

新鲜事 原文

wt20
18小时前
做线段树上树修改太难受,于是学了个树剖


活动打卡代码 AcWing 2568. 树链剖分

wt20
18小时前
#include <iostream>
#include <cstring>
#include <algorithm>

using namespace std;

typedef long long LL;
const int N = 1e5 + 10, M = N * 2;
int n, m;
int h[N], e[M], ne[M], idx, w[N];
int id[N], nw[N], cnt;
int dep[N], sz[N], top[N], fa[N], son[N];
struct Node{
    int l, r;
    LL add, sum;
}tr[N * 4];

void add(int a, int b)  {
    e[idx] = b, ne[idx] = h[a], h[a] = idx ++ ;
}

void pushup(int u) {
    tr[u].sum = tr[u << 1].sum + tr[u << 1 | 1].sum;
}

void pushdown(int u) {
    auto &root = tr[u], &left = tr[u << 1], &right = tr[u << 1 | 1];
    if (root.add) {
        left.add += root.add, right.add += root.add;
        left.sum += (left.r - left.l + 1) * root.add;
        right.sum += (right.r - right.l + 1) * root.add;
        root.add = 0;
    }
}

void build(int u, int l, int r) {
    tr[u] = {l, r, 0, nw[r]};
    if (l == r) return;
    int mid = l + r >> 1;
    build(u << 1, l, mid), build(u << 1 | 1, mid + 1, r);
    pushup(u);
}

void update(int u, int l, int r, int d) {
    if (tr[u].l >= l && tr[u].r <= r) {
        tr[u].add += d;
        tr[u].sum += 1ll * (tr[u].r - tr[u].l + 1) * d;
    }
    else {
        pushdown(u);
        int mid = tr[u].l + tr[u].r >> 1;
        if (l <= mid) update(u << 1, l, r, d);
        if (r > mid) update(u << 1 | 1, l, r, d);
        pushup(u);
    }
}

LL query(int u, int l, int r) {
    if (tr[u].l >= l && tr[u].r <= r) return tr[u].sum;
    else {
        pushdown(u);
        int mid = tr[u].l + tr[u].r >> 1;
        LL res = 0;
        if (l <= mid) res = query(u << 1, l, r);
        if (r > mid) res += query(u << 1 | 1, l, r);
        return res;
    }
}

void dfs1(int u, int father, int depth) {
    dep[u] = depth, fa[u] = father, sz[u] = 1;
    for (int i = h[u]; ~i; i = ne[i]) {
        int j = e[i];
        if (j == father) continue;
        dfs1(j, u, depth + 1);
        sz[u] += sz[j];
        if (sz[son[u]] < sz[j]) son[u] = j;
    }
}

void dfs2(int u, int t) {
    id[u] = ++ cnt, nw[cnt] = w[u], top[u] = t;
    if (!son[u]) return;
    dfs2(son[u], t);
    for (int i = h[u]; ~i; i = ne[i]) {
        int j = e[i];
        if (j == fa[u] || j == son[u]) continue;
        dfs2(j, j); 
    }
}

void update_path(int u, int v, int k) {
    while (top[u] != top[v]) {
        if (dep[top[u]] < dep[top[v]]) swap(u, v);
        update(1, id[top[u]], id[u], k);
        u = fa[top[u]];
    }
    if (dep[u] < dep[v]) swap(u, v);
    update(1, id[v], id[u], k);
}

void update_tree(int u, int k) {
    update(1, id[u], id[u] + sz[u] - 1, k);
}

LL query_path(int u, int v) {
    LL res = 0;
    while (top[u] != top[v]) {
        if (dep[top[u]] < dep[top[v]]) swap(u, v);
        res += query(1, id[top[u]], id[u]);
        u = fa[top[u]];
    }
    if (dep[u] < dep[v]) swap(u, v);
    res += query(1, id[v], id[u]);
    return res;
}

LL query_tree(int u) {
    return query(1, id[u], id[u] + sz[u] - 1);
}

int main() {
    cin >> n;
    for (int i = 1; i <= n; i ++ ) cin >> w[i];
    memset(h, -1, sizeof h);
    for (int i = 1; i < n; i ++ ) {
        int a, b;
        cin >> a >> b;
        add(a, b), add(b, a);
    }
    dfs1(1, -1, 1);
    dfs2(1, 1);
    build(1, 1, n);
    cin >> m;
    while (m -- ) {
        int t, u, v, k;
        cin >> t >> u;
        if (t == 1) {
            cin >> v >> k;
            update_path(u, v, k);
        }
        else if (t == 2) {
            cin >> k;
            update_tree(u, k);
        }
        else if (t == 3) {
            cin >> v;
            cout << query_path(u, v) << endl;
        }
        else cout << query_tree(u) << endl;
    }
    return 0;
}


新鲜事 原文

wt20
2天前
AcWing《蓝桥杯集训·每日一题》拼团优惠!https://www.acwing.com/activity/content/introduction/2869/group_buy/121897/ 找2人,速来!


新鲜事 原文

wt20
14天前
AcWing【集日历瓜分10000AC币活动】赠送4月日历!https://www.acwing.com/file_system/gui/taskbar/widgets/clock/calendar/send/49782/4/fe42a786/


新鲜事 原文

wt20
17天前
AcWing【集日历瓜分10000AC币活动】赠送10月日历!https://www.acwing.com/file_system/gui/taskbar/widgets/clock/calendar/send/49782/10/39c90bc1/ 速抢!!!!!


新鲜事 原文

wt20
18天前
还有人没有日历(除了10的,速速来找我,缺10的欢迎关注等待随机时间发10


活动打卡代码 AcWing 4729. 解密

wt20
18天前

除夕后补,今天要写第二分块(



新鲜事 原文

wt20
18天前
AcWing【集日历瓜分10000AC币活动】赠送10月日历!https://www.acwing.com/file_system/gui/taskbar/widgets/clock/calendar/send/49782/10/072e594b/


新鲜事 原文

wt20
18天前
除夕快乐!


新鲜事 原文

wt20
18天前
AcWing【集日历瓜分10000AC币活动】赠送10月日历!https://www.acwing.com/file_system/gui/taskbar/widgets/clock/calendar/send/49782/10/a48dc5c2/