#include using namespace std; const int KMAX = 10; const int INF = 1000000; const long long MOD = 998244353; struct affine{ long long a, b; affine(){ a = 1; b = 0; } affine(int a, int b): a(a), b(b){ } }; affine composite(affine A, affine B){ return affine(A.a * B.a % MOD, (A.b * B.a + B.b) % MOD); } int value(affine A, int x){ return (A.a * x + A.b) % MOD; } template struct dual_segment_tree{ int N; vector ST; function f; T E; dual_segment_tree(int n, function f, T E): f(f), E(E){ N = 1; while (N < n){ N *= 2; } ST = vector(N * 2 - 1, E); } void push(int i){ if (i < N - 1){ ST[i * 2 + 1] = f(ST[i * 2 + 1], ST[i]); ST[i * 2 + 2] = f(ST[i * 2 + 2], ST[i]); ST[i] = E; } } T operator [](int k){ int v = 0; for (int i = N / 2; i >= 1; i >>= 1){ push(v); if ((k & i) == 0){ v = v * 2 + 1; } else { v = v * 2 + 2; } } return ST[v]; } void range_apply(int L, int R, T x, int i, int l, int r){ if (r <= L || R <= l){ } else if (L <= l && r <= R){ ST[i] = f(ST[i], x); } else { push(i); int m = (l + r) / 2; range_apply(L, R, x, i * 2 + 1, l, m); range_apply(L, R, x, i * 2 + 2, m, r); } } void range_apply(int L, int R, T x){ range_apply(L, R, x, 0, 0, N); } }; int main(){ ios_base::sync_with_stdio(false); cin.tie(nullptr); int N, Q; cin >> N >> Q; vector> E(N); for (int i = 0; i < N - 1; i++){ int A, B; cin >> A >> B; A--; B--; E[A].push_back(B); E[B].push_back(A); } vector X(N); for (int i = 0; i < N; i++){ cin >> X[i]; } int N2 = N + KMAX; E.resize(N2); E[0].push_back(N); E[N].push_back(0); for (int i = 0; i < KMAX - 1; i++){ E[N + i].push_back(N + i + 1); E[N + i + 1].push_back(N + i); } vector p(N2, -1); vector> c(N2); vector d(N2, 0); vector b; queue q; q.push(N2 - 1); while (!q.empty()){ int v = q.front(); q.pop(); b.push_back(v); for (int w : E[v]){ if (w != p[v]){ p[w] = v; c[v].push_back(w); d[w] = d[v] + 1; q.push(w); } } } reverse(b.begin(), b.end()); vector pk(N2); for (int i = 0; i < N; i++){ if (d[i] >= KMAX){ pk[i] = i; for (int j = 0; j < KMAX; j++){ pk[i] = p[pk[i]]; } } } vector>> ck(KMAX + 1, vector>(N2)); for (int i = 0; i < N2; i++){ ck[0][i].push_back(i); } for (int i = 0; i < KMAX; i++){ for (int j = 0; j < N2; j++){ for (int k : c[j]){ for (int l : ck[i][k]){ ck[i + 1][j].push_back(l); } } } } vector sz(N2, 1); for (int v : b){ for (int &w : c[v]){ sz[v] += sz[w]; if (sz[w] > sz[c[v][0]]){ swap(w, c[v][0]); } } } reverse(b.begin(), b.end()); vector id(N2, 0); for (int v : b){ for (int w : c[v]){ if (w == c[v][0]){ id[w] = id[v] + 1; } else { id[w] = 0; } } } vector> hc(KMAX + 1, vector(N2, -1)); for (int i = 0; i < N2; i++){ hc[0][i] = i; for (int j = 0; j < KMAX; j++){ if (hc[j][i] != -1){ if (!c[hc[j][i]].empty()){ hc[j + 1][i] = c[hc[j][i]][0]; } } } } vector in(N); int t = 0; vector in2(N), out2(N); vector hld(N); int t2 = 0; auto dfs1 = [&](auto dfs1, int v) -> void { if (v < N){ in2[v] = t; hld[v] = t2; t2++; } if (id[v] >= KMAX && v < N){ in[v] = t; t++; } for (int w : c[v]){ dfs1(dfs1, w); } if (v < N){ out2[v] = t; } }; dfs1(dfs1, N2 - 1); vector in3(N), out3(N); auto dfs2 = [&](auto dfs2, int v) -> void { for (int w : ck[KMAX][v]){ if (id[w] < KMAX){ in[w] = t; t++; } } if (v < N){ in3[v] = t; } for (int w : c[v]){ dfs2(dfs2, w); } if (v < N){ out3[v] = t; } }; dfs2(dfs2, N2 - 1); vector next(N, 0); for (int v : b){ if (v < N){ for (int w : c[v]){ if (in[w] == in[v] + 1){ next[w] = next[v]; } else { next[w] = w; } } } } vector> mn(KMAX + 1, vector(N2, INF)); vector> mx(KMAX + 1, vector(N2, -1)); for (int i = 0; i < N; i++){ int v = i; for (int j = 0; j <= KMAX; j++){ if (id[i] < KMAX){ mn[j][v] = min(mn[j][v], in[i]); mx[j][v] = max(mx[j][v], in[i]); } if (v == N2 - 1){ break; } v = p[v]; } } dual_segment_tree ST(N, composite, affine()); for (int i = 0; i < Q; i++){ int T; cin >> T; if (T == 1){ int V; cin >> V; V--; cout << value(ST[in[V]], X[V]) << '\n'; } if (T == 2){ int V, K; long long C, D; cin >> V >> K >> C >> D; V--; for (int j = 0; j <= K; j++){ for (int k = K - j; k >= max(K - j - 1, 0); k--){ if (mn[k][V] != INF){ ST.range_apply(mn[k][V], mx[k][V] + 1, affine(C, D)); } if (hc[k][V] != -1 && hc[k][V] < N){ if (id[hc[k][V]] >= KMAX){ ST.range_apply(in[hc[k][V]], in[hc[k][V]] + 1, affine(C, D)); } } } V = p[V]; } } if (T == 3){ int V; long long C, D; cin >> V >> C >> D; V--; ST.range_apply(in2[V], out2[V], affine(C, D)); for (int j = 0; j <= KMAX; j++){ if (mn[j][V] != INF){ ST.range_apply(mn[j][V], mx[j][V] + 1, affine(C, D)); } } ST.range_apply(in3[V], out3[V], affine(C, D)); } if (T == 4){ int U, V; long long C, D; cin >> U >> V >> C >> D; U--; V--; while (true){ if (hld[U] > hld[V]){ swap(U, V); } if (next[U] == next[V]){ ST.range_apply(in[U], in[V] + 1, affine(C, D)); break; } ST.range_apply(in[next[V]], in[V] + 1, affine(C, D)); V = p[next[V]]; } } } }