如果没有年龄区间

画图发现ans=\(\sum_{i=1}^n dis_i + n * dis_u - 2 * \sum_{i=1}^{n} dis_{lca(i, u)}\)

对\(\sum_{i=1}^{n} dis_{lca(i, u)}\)用树链剖分，对于每个点，都向上走到根，记录每条路经过次数，询问时从u向上跳，每次加边权*覆盖次数即可

现在有年龄区间，加个前缀和，计算点权在[L,R]内的点到点u的距离，用主席树

然而主席树怎么pushdown，加上标记永久化

# include <bits/stdc++.h>

# define IL inline

# define RG register

# define Fill(a, b) memset(a, b, sizeof(a))

using namespace std;

typedef long long ll;

const int _(3e6 + 10), __(1e7 + 10);

IL ll Read(){

    RG char c = getchar(); RG ll x = 0, z = 1;

    for(; c < '0' || c > '9'; c = getchar()) z = c == '-' ? -1 : 1;

    for(; c >= '0' && c <= '9'; c = getchar()) x = (x << 1) + (x << 3) + (c ^ 48);

    return x * z;

}

int n, m, Q, A, id[_], o[_];

int fst[_], nxt[_], to[_], w[_], size[_], fa[_], dfn[_], top[_], cnt, son[_], pos[_];

int rt[__], ls[__], rs[__], cov[__];

ll dis[_], sdis[_], sum[__];

IL bool Cmp(RG int x, RG int y){  return o[x] < o[y];  }

IL void Add(RG int u, RG int v, RG int f){  w[cnt] = f; to[cnt] = v; nxt[cnt] = fst[u]; fst[u] = cnt++;  }

IL void Dfs1(RG int u){

	size[u] = 1;

	for(RG int e = fst[u]; e != -1; e = nxt[e]){

		if(size[to[e]]) continue;

		dis[to[e]] = dis[u] + w[e]; fa[to[e]] = u;

		Dfs1(to[e]);

		size[u] += size[to[e]];

		if(size[to[e]] > size[son[u]]) son[u] = to[e];

	}

}

IL void Dfs2(RG int u, RG int Top){

	dfn[u] = ++cnt; top[u] = Top; pos[cnt] = u;

	if(son[u]) Dfs2(son[u], Top);

	for(RG int e = fst[u]; e != -1; e = nxt[e])

		if(!dfn[to[e]]) Dfs2(to[e], to[e]);

}

IL void Build(RG int &x, RG int l, RG int r){

	x = ++cnt;

	if(l == r) return;

	RG int mid = (l + r) >> 1;

	Build(ls[x], l, mid); Build(rs[x], mid + 1, r);

}

IL void Modify(RG int &x, RG int l, RG int r, RG int L, RG int R){

	sum[++cnt] = sum[x]; ls[cnt] = ls[x]; rs[cnt] = rs[x]; cov[cnt] = cov[x]; x = cnt;

	sum[x] += dis[pos[R]] - dis[fa[pos[L]]];

	if(L <= l && R >= r){  cov[x]++; return;  }

	RG int mid = (l + r) >> 1;

	if(R <= mid) Modify(ls[x], l, mid, L, R);

	else if(L > mid) Modify(rs[x], mid + 1, r, L, R);

	else Modify(ls[x], l, mid, L, mid), Modify(rs[x], mid + 1, r, mid + 1, R);

}

IL ll Query(RG int x, RG int ad, RG int l, RG int r, RG int L, RG int R){

	if(L <= l && R >= r) return sum[x] + 1LL * ad * (dis[pos[R]] - dis[fa[pos[L]]]);

	RG int mid = (l + r) >> 1;

	if(R <= mid) return Query(ls[x], ad + cov[x], l, mid, L, R);

	if(L > mid) return Query(rs[x], ad + cov[x], mid + 1, r, L, R);

	return Query(ls[x], ad + cov[x], l, mid, L, mid) + Query(rs[x], ad + cov[x], mid + 1, r, mid + 1, R);

}

IL ll Calc(RG int u, RG int i){

	RG ll ans = 0;

	while(top[u] ^ 1) ans += Query(rt[i], 0, 1, n, dfn[top[u]], dfn[u]), u = fa[top[u]];

	ans += Query(rt[i], 0, 1, n, 1, dfn[u]);

	return ans;

}

int main(RG int argc, RG char* argv[]){

    n = Read(); Q = Read(); A = Read();

	for(RG int i = 1; i <= n; ++i) o[i] = Read(), id[i] = i, fst[i] = -1;

	sort(id + 1, id + n + 1, Cmp); sort(o + 1, o + n + 1);

	for(RG int i = 1, aa, b, c; i < n; ++i) aa = Read(), b = Read(), c = Read(), Add(aa, b, c), Add(b, aa, c);

	Dfs1(1); cnt = 0; Dfs2(1, 1); cnt = 0; Build(rt[0], 1, n);

	for(RG int i = 1; i <= n; ++i){

		RG int u = id[i]; rt[i] = rt[i - 1];

		sdis[i] = sdis[i - 1] + dis[u];

		while(top[u] ^ 1) Modify(rt[i], 1, n, dfn[top[u]], dfn[u]), u = fa[top[u]];

		Modify(rt[i], 1, n, 1, dfn[u]);

	}

	for(RG ll L, R, u, a, b, ans = 0; Q; --Q){

		u = Read(); a = Read(); b = Read();

		L = min((a + ans) % A, (b + ans) % A);

		R = max((a + ans) % A, (b + ans) % A);

		L = lower_bound(o + 1, o + n + 1, L) - o;

		R = lower_bound(o + 1, o + n + 1, R + 1) - o - 1;

		ans = 1LL * (R - L + 1) * dis[u] + sdis[R] - sdis[L - 1] - 2LL * (Calc(u, R) - Calc(u, L - 1));

		printf("%lld\n", ans);

	}

    return 0;

}