树剖(重剖)RE求调(P3384)
  • 板块学术版
  • 楼主Eric_Shi
  • 当前回复11
  • 已保存回复11
  • 发布时间2023/8/3 10:42
  • 上次更新2023/11/3 06:12:45
查看原帖
树剖(重剖)RE求调(P3384)
813709
Eric_Shi楼主2023/8/3 10:42

测评记录

#include <iostream>
#include <algorithm>
using namespace std;

#define MOD p
#define MAX 200001
#define INF 0x7fffffff

int n, m, r, p, w[MAX], opt, x, y, z;
int fa[MAX], dep[MAX], size[MAX], son[MAX];
//节点的: 父亲 深度 子树的节点个数 重儿子
int top[MAX], dfn[MAX], rank[MAX], cnt = 0, wt[MAX];
//节点的: 所在重链顶部节点 Dfs序 Dfs序对应节点编号 cnt 点权

//链式前向星
struct Edge {
    int to, next;
} e[MAX];
int head[MAX], num = 0;
void add(int u, int v) {
    e[++num].to = v;
    e[num].next = head[u];
    head[u] = num;
}

//Dfs
void Dfs1(int u, int f, int deep) {
    fa[u] = f;
    dep[u] = deep;
    size[u] = 1;
    int maxson = -INF;
    for (int i = head[u]; i != 0; i = e[i].next) {
        int v = e[i].to;
        if (v == f)
            continue;
        Dfs1(v, u, deep + 1);
        size[u] += size[v];
        if (size[v] > maxson) {
            son[x] = v;
            maxson = size[v];
        }
    }
}
void Dfs2(int u, int topf) {
    dfn[u] = ++cnt;
    wt[cnt] = w[u];
    top[u] = topf;
    if (!son[u])
        return;
    Dfs2(son[u], topf);
    for (int i = head[u]; i != 0; i = e[i].next) {
        int v = e[i].to;
        if (v == fa[u] || v == son[u])
            continue;
        Dfs2(v, v);
    }
}

//线段树
int tree[MAX << 2], lazy[MAX << 2];
void Build(int rt, int l, int r) {
    if (l == r) {
        tree[rt] = wt[l];
        return;
    }
    int mid = (l + r) >> 1;
    Build(rt << 1, l, mid);
    Build(rt << 1 | 1, mid + 1, r);
    tree[rt] = tree[rt << 1] + tree[rt << 1 | 1];
    tree[rt] %= MOD;
}
void Update(int rt, int l, int r, int w) {
    tree[rt] += (r - l + 1) * w % MOD;
    lazy[rt] += w %= MOD;
} 
void Pushdown(int rt, int l, int r) {
    int mid = (l + r) >> 1;
    Update(rt << 1, l, mid, lazy[rt]);
    Update(rt << 1 | 1, mid + 1, r, lazy[rt]);
    lazy[rt] = 0;
}
void Modify(int rt, int l, int r, int s, int t, int w) {
    if (s <= l && r <= t) {
        Update(rt, l, r, w);
        return;
    }
    Pushdown(rt, l, r);
    int mid = (l + r) >> 1;
    if (s <= mid)
        Modify(rt << 1, l, mid, s, t, w);
    if (t > mid)
        Modify(rt << 1 | 1, mid + 1, r, s, t, w);
    tree[rt] = tree[rt << 1] + tree[rt << 1 | 1];
    tree[rt] %= MOD;
}
int Query(int rt, int l, int r, int s, int t) {
    if (s <= l && r <= t)
        return tree[rt] % MOD;
    Pushdown(rt, l, r);
    int mid = (l + r) >> 1;
    int lson = Query(rt << 1, l, mid, s, t) % MOD;
    int rson = Query(rt << 1 | 1, mid + 1, r, s, t) % MOD;
    if (t <= mid)
        return lson;
    else if (s > mid)
        return rson;
    else
        return lson + rson;
}

//Mode
void Mrange(int x, int y, int z) {
    z %= MOD;
    while (top[x] != top[y]) {
        if (dep[top[x]] < dep[top[y]])
            swap(x, y);
        Modify(1, 1, n, dfn[top[x]], dfn[x], z);
        x = fa[top[x]];
    }
    if (dep[x] > dep[y])
        swap(x, y);
    Modify (1, 1, n, dfn[x], dfn[y], z);
}
int Qrange(int x, int y) {
    int ans = 0;
    while (top[x] != top[y]) {
		if (dep[top[x]] < dep[top[y]])
            swap(x, y);
		ans += Query(1, 1, n, dfn[top[x]], dfn[x]);
		ans %= MOD;
		x = fa[top[x]];
	}
	if(dep[x] > dep[y])
        swap(x, y);
	ans += Query(1, 1, n, dfn[x], dfn[y]);
	ans %= MOD;
	return ans;
}
void Mtree(int x, int z) {
    Modify(1, 1, n, dfn[x], dfn[x]+size[x]-1, z);
}
int Qtree(int x) {
    return Query(1, 1, n, dfn[x], dfn[x] + size[x] - 1);
}

int main() {
    cin >> n >> m >> r >> p;
    for (int i = 1; i <= n; i++)
        cin >> w[i];
    for (int i = 1; i < n; i++) {
        cin >> x >> y;
        add(x, y);
        add(y, x);
    }
    Dfs1(r, 0, 1);
    Dfs2(r, r);
    Build(1, 1, n);
    while (m--) {
        cin >> opt;
        switch (opt) {
            case 1:
                cin >> x >> y >> z;
                Mrange(x, y, z);
                break;
            case 2:
                cin >> x >> y;
                cout << Qrange(x, y) << "\n";
                break;
            case 3:
                cin >> x >> z;
                Mtree(x, z);
                break;
            case 4:
                cin >> x;
                cout << Qtree(x) << "\n";
                break;
            default:
                break;
        }
    }

    return 0;
}

感谢大家!灌水区一个正经的都没有: (

2023/8/3 10:42
加载中...