#include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include using namespace std; using Int = long long; template ostream &operator<<(ostream &os, const pair &a) { return os << "(" << a.first << ", " << a.second << ")"; }; template ostream &operator<<(ostream &os, const vector &as) { const int sz = as.size(); os << "["; for (int i = 0; i < sz; ++i) { if (i >= 256) { os << ", ..."; break; } if (i > 0) { os << ", "; } os << as[i]; } return os << "]"; } template void pv(T a, T b) { for (T i = a; i != b; ++i) cerr << *i << " "; cerr << endl; } template bool chmin(T &t, const T &f) { if (t > f) { t = f; return true; } return false; } template bool chmax(T &t, const T &f) { if (t < f) { t = f; return true; } return false; } #define COLOR(s) ("\x1b[" s "m") struct Hld { int n, rt; // needs to be tree // vertex lists // modified in build(rt) (parent removed, heavy child first) vector> graph; vector sz, par, dep; int zeit; vector dis, fin, sid; // head vertex (minimum depth) in heavy path vector head; Hld() : n(0), rt(-1), zeit(0) {} explicit Hld(int n_) : n(n_), rt(-1), graph(n), zeit(0) {} void ae(int u, int v) { assert(0 <= u); assert(u < n); assert(0 <= v); assert(v < n); graph[u].push_back(v); graph[v].push_back(u); } void dfsSz(int u) { sz[u] = 1; for (const int v : graph[u]) { auto it = std::find(graph[v].begin(), graph[v].end(), u); if (it != graph[v].end()) graph[v].erase(it); par[v] = u; dep[v] = dep[u] + 1; dfsSz(v); sz[u] += sz[v]; } } void dfsHld(int u) { dis[u] = zeit++; const int deg = graph[u].size(); if (deg > 0) { int vm = graph[u][0]; int jm = 0; for (int j = 1; j < deg; ++j) { const int v = graph[u][j]; if (sz[vm] < sz[v]) { vm = v; jm = j; } } swap(graph[u][0], graph[u][jm]); head[vm] = head[u]; dfsHld(vm); for (int j = 1; j < deg; ++j) { const int v = graph[u][j]; head[v] = v; dfsHld(v); } } fin[u] = zeit; } void build(int rt_) { assert(0 <= rt_); assert(rt_ < n); rt = rt_; sz.assign(n, 0); par.assign(n, -1); dep.assign(n, -1); dep[rt] = 0; dfsSz(rt); zeit = 0; dis.assign(n, -1); fin.assign(n, -1); head.assign(n, -1); head[rt] = rt; dfsHld(rt); assert(zeit == n); sid.assign(n, -1); for (int u = 0; u < n; ++u) sid[dis[u]] = u; } friend ostream &operator<<(ostream &os, const Hld &hld) { const int maxDep = *max_element(hld.dep.begin(), hld.dep.end()); vector ss(2 * maxDep + 1); int pos = 0, maxPos = 0; for (int j = 0; j < hld.n; ++j) { const int u = hld.sid[j]; const int d = hld.dep[u]; if (hld.head[u] == u) { if (j != 0) { pos = maxPos + 1; ss[2 * d - 1].resize(pos, '-'); ss[2 * d - 1] += '+'; } } else { ss[2 * d - 1].resize(pos, ' '); ss[2 * d - 1] += '|'; } ss[2 * d].resize(pos, ' '); ss[2 * d] += std::to_string(u); if (maxPos < static_cast(ss[2 * d].size())) { maxPos = ss[2 * d].size(); } } for (int d = 0; d <= 2 * maxDep; ++d) os << ss[d] << '\n'; return os; } bool contains(int u, int v) const { return (dis[u] <= dis[v] && dis[v] < fin[u]); } int lca(int u, int v) const { assert(0 <= u); assert(u < n); assert(0 <= v); assert(v < n); for (; head[u] != head[v]; ) (dis[u] > dis[v]) ? (u = par[head[u]]) : (v = par[head[v]]); return (dis[u] > dis[v]) ? v : u; } int jumpUp(int u, int d) const { assert(0 <= u); assert(u < n); assert(d >= 0); if (dep[u] < d) return -1; const int tar = dep[u] - d; for (u = head[u]; ; u = head[par[u]]) { if (dep[u] <= tar) return sid[dis[u] + (tar - dep[u])]; } } int jump(int u, int v, int d) const { assert(0 <= u); assert(u < n); assert(0 <= v); assert(v < n); assert(d >= 0); const int l = lca(u, v); const int du = dep[u] - dep[l], dv = dep[v] - dep[l]; if (d <= du) { return jumpUp(u, d); } else if (d <= du + dv) { return jumpUp(v, du + dv - d); } else { return -1; } } // [u, v) or [u, v] template void doPathUp(int u, int v, bool inclusive, F f) const { assert(0 <= u); assert(u < n); assert(0 <= v); assert(v < n); assert(contains(v, u)); for (; head[u] != head[v]; u = par[head[u]]) f(dis[head[u]], dis[u] + 1); if (inclusive) { f(dis[v], dis[u] + 1); } else { if (v != u) f(dis[v] + 1, dis[u] + 1); } } // not path order, include lca(u, v) or not template void doPath(int u, int v, bool inclusive, F f) const { assert(0 <= u); assert(u < n); assert(0 <= v); assert(v < n); const int l = lca(u, v); doPathUp(u, l, false, f); doPathUp(v, l, inclusive, f); } // find deepest true for pred: [u, root] -> bool, increasing // -1 if !pred(rt) template int findUp(int u, Pred pred) const { assert(0 <= u); assert(u < n); for (; ~u; ) { const int h = head[u]; if (pred(h)) { int lo = dis[h], hi = dis[u] + 1; for (; lo + 1 < hi; ) { const int mid = (lo + hi) / 2; (pred(sid[mid]) ? lo : hi) = mid; } return sid[lo]; } u = par[h]; } return -1; } // (vs, ps): compressed tree // vs: DFS order (sorted by dis) // vs[ps[x]]: the parent of vs[x] // ids[vs[x]] = x, not set for non-tree vertex vector ids; pair, vector> compress(vector us) { // O(n) first time ids.resize(n, -1); std::sort(us.begin(), us.end(), [&](int u, int v) -> bool { return (dis[u] < dis[v]); }); us.erase(std::unique(us.begin(), us.end()), us.end()); int usLen = us.size(); assert(usLen >= 1); for (int x = 1; x < usLen; ++x) us.push_back(lca(us[x - 1], us[x])); std::sort(us.begin(), us.end(), [&](int u, int v) -> bool { return (dis[u] < dis[v]); }); us.erase(std::unique(us.begin(), us.end()), us.end()); usLen = us.size(); for (int x = 0; x < usLen; ++x) ids[us[x]] = x; vector ps(usLen, -1); for (int x = 1; x < usLen; ++x) ps[x] = ids[lca(us[x - 1], us[x])]; return make_pair(us, ps); } }; //////////////////////////////////////////////////////////////////////////////// struct Set { // max{ceil(log_64(n)), 1} int log64N, n; vector a[6]; explicit Set(int n_ = 0) : n(n_) { assert(n >= 0); int m = n ? n : 1; for (int d = 0; ; ++d) { m = (m + 63) >> 6; a[d].assign(m, 0); if (m == 1) { log64N = d + 1; break; } } } bool empty() const { return !a[log64N - 1][0]; } bool contains(int x) const { return (a[0][x >> 6] >> (x & 63)) & 1; } void insert(int x) { for (int d = 0; d < log64N; ++d) { const int q = x >> 6, r = x & 63; a[d][q] |= 1ULL << r; x = q; } } void erase(int x) { for (int d = 0; d < log64N; ++d) { const int q = x >> 6, r = x & 63; if ((a[d][q] &= ~(1ULL << r))) break; x = q; } } // max s.t. <= x (or -1) int prev(int x) const { if (x > n - 1) x = n - 1; for (int d = 0; d <= log64N; ++d) { if (x < 0) break; const int q = x >> 6, r = x & 63; const unsigned long long lower = a[d][q] << (63 - r); if (lower) { x -= __builtin_clzll(lower); for (int e = d; --e >= 0; ) x = x << 6 | (63 - __builtin_clzll(a[e][x])); return x; } x = q - 1; } return -1; } // min s.t. >= x (or n) int next(int x) const { if (x < 0) x = 0; for (int d = 0; d < log64N; ++d) { const int q = x >> 6, r = x & 63; if (static_cast(q) >= a[d].size()) break; const unsigned long long upper = a[d][q] >> r; if (upper) { x += __builtin_ctzll(upper); for (int e = d; --e >= 0; ) x = x << 6 | __builtin_ctzll(a[e][x]); return x; } x = q + 1; } return n; } }; template void bAdd(vector &bit, int pos, const T &val) { const int bitN = bit.size(); for (int x = pos; x < bitN; x |= x + 1) bit[x] += val; } template T bSum(const vector &bit, int pos) { T ret = 0; for (int x = pos; x > 0; x &= x - 1) ret += bit[x - 1]; return ret; } template T bSum(const vector &bit, int pos0, int pos1) { return bSum(bit, pos1) - bSum(bit, pos0); } int N; vector A, B; vector C; Hld hld; Set on; vector bit; int dist(int l, int r) { const int u = hld.sid[l]; const int v = hld.sid[r]; return hld.dep[u] + hld.dep[v] - 2 * hld.dep[hld.lca(u, v)]; } void add(int l, int r, int sig) { if (0 <= l && r < N) bAdd(bit, l, sig * dist(l, r)); } void add(int u) { const int m = hld.dis[u]; const int l = on.prev(m); const int r = on.next(m); on.insert(m); add(l, r, -1); add(l, m, +1); add(m, r, +1); } void rem(int u) { const int m = hld.dis[u]; on.erase(m); const int l = on.prev(m); const int r = on.next(m); add(l, m, -1); add(m, r, -1); add(l, r, +1); } int calc(int L, int R) { const int l = on.next(L); if (l >= R) return 0; const int r = on.prev(R - 1); return (bSum(bit, l, r) + dist(l, r)) / 2 + 1; } int calc(int L0, int R0, int L1, int R1) { const int l0 = on.next(L0); if (l0 >= R0) return calc(L1, R1); const int l1 = on.next(L1); if (l1 >= R1) return calc(L0, R0); const int r0 = on.prev(R0 - 1); const int r1 = on.prev(R1 - 1); return (bSum(bit, l0, r0) + dist(r0, l1) + bSum(bit, l1, r1) + dist(l0, r1)) / 2 + 1; } int main() { for (; ~scanf("%d", &N); ) { A.resize(N - 1); B.resize(N - 1); for (int i = 0; i < N - 1; ++i) { scanf("%d%d", &A[i], &B[i]); --A[i]; --B[i]; } C.resize(N); for (int u = 0; u < N; ++u) scanf("%d", &C[u]); hld = Hld(N); for (int i = 0; i < N - 1; ++i) hld.ae(A[i], B[i]); hld.build(0); on = Set(N); bit.assign(N, 0); for (int u = 0; u < N; ++u) if (C[u]) add(u); int Q; scanf("%d", &Q); for (int q = 0; q < Q; ++q) { int O; scanf("%d", &O); if (O == 1) { int V; scanf("%d", &V); --V; (C[V] ^= 1) ? add(V) : rem(V); } else if (O == 2) { // X-Y-* int X, Y; scanf("%d%d", &X, &Y); --X; --Y; int ans = 0; if (X == Y) { ans = calc(0, N); } else { X = hld.jump(Y, X, 1); if (X == hld.par[Y]) { ans = calc(hld.dis[Y], hld.fin[Y]); } else { ans = calc(0, hld.dis[X], hld.fin[X], N); } } printf("%d\n", ans); } else { assert(false); } } } return 0; }