#include #include #include using namespace std; using i32 = int32_t; using u32 = uint32_t; using i64 = int64_t; using u64 = uint64_t; #define rep(i,n) for(int i=0; i<(n); i++) struct heavy_light_decomposition { private: int N; vector P; vector PP; vector PD; vector D; public: vector rangeL; vector rangeR; heavy_light_decomposition(const vector>& E = { {} }) { N = E.size(); P.assign(N, -1); vector I = { 0 }; rep(i, I.size()) { int p = I[i]; for (int e : E[p]) if (P[p] != e) { I.push_back(e); P[e] = p; } } vector Z(N, 1); for (int i = N - 1; i >= 1; i--) Z[P[I[i]]] += Z[I[i]]; PP.resize(N); for (int i = 0; i < N; i++) PP[i] = i; vector nx(N, -1); for (int p : I) if (p != 0) { if (nx[P[p]] == -1) nx[P[p]] = p; if (Z[nx[P[p]]] < Z[p]) nx[P[p]] = p; } rep(i, N) if (nx[i] != -1) PP[nx[i]] = i; for (int p : I) if (p != 0) PP[p] = PP[PP[p]]; PD.assign(N, N); PD[0] = 0; for (int p : I) if (p != 0) PD[p] = min(PD[PP[p]], PD[P[p]] + 1); D.assign(N, 0); for (int p : I) D[p] = (p == 0) ? 0 : (D[P[p]] + 1); rangeL.assign(N, -1); rangeR.assign(N, -1); vector dfs; int ir = 0; dfs.push_back(0); while (dfs.size()) { int p = dfs.back(); dfs.pop_back(); if (p < 0) { rangeR[-1 - p] = ir; continue; } dfs.push_back(-1 - p); for (int e : E[p]) if (P[p] != e) if (e != nx[p]) dfs.push_back(e); if (nx[p] != -1) dfs.push_back(nx[p]); rangeL[p] = ir++; } } int lca(int u, int v) { if (PD[u] < PD[v]) swap(u, v); while (PD[u] > PD[v]) u = P[PP[u]]; while (PP[u] != PP[v]) { u = P[PP[u]]; v = P[PP[v]]; } return (D[u] > D[v]) ? v : u; } vector> getpath(int r, int c) { vector> res; while (PD[r] < PD[c]) { res.push_back({ rangeL[PP[c]], rangeL[c] + 1 }); c = P[PP[c]]; } if (PP[r] != PP[c]) return {}; if (D[r] > D[c]) return {}; res.push_back({ rangeL[r], rangeL[c] + 1 }); reverse(res.begin(), res.end()); return move(res); } }; class CoordinateCompress { using Elem = i64; static const Elem negInf = -1001001001001001001; vector> G; vector res; vector mRealval; Elem mMaxcoord; bool ok = true; void calc() { if (ok) return; sort(G.begin(), G.end()); res.resize(G.size()); mRealval.clear(); Elem x = negInf; int p = -1; rep(i, G.size()) { if (x != G[i].first) { x = G[i].first; mRealval.push_back(x); p++; } res[G[i].second] = p; } mMaxcoord = p; ok = true; } public: int push(Elem x) { ok = false; G.push_back({ x,(int)G.size() }); return G.back().second; } int operator[](int i) { calc(); return res[i]; } Elem realval(int x) { calc(); return mRealval[x]; } Elem maxcoord() { calc(); return mMaxcoord; } }; struct Edge{ int to; u32 h; }; int N; vector> E; vector> hldE; heavy_light_decomposition hld; vector P; vector Z; CoordinateCompress CCH; vector H; vector> black_stack; vector prev_black; vector block_size; vector root_size; i64 ans; void dfs(int p = 0, int pre = -1){ if(pre == -1){ black_stack.assign(N,{0}); prev_black.assign(N,0); block_size.assign(N,0); root_size.assign(N,N); ans = 0; } int h = H[p]; block_size[p] = Z[p]; prev_black[p] = black_stack[h].back(); if(prev_black[p] != 0) block_size[prev_black[p]] -= Z[p]; else root_size[h] -= Z[p]; for(int e : hldE[p]) if(e != pre){ black_stack[H[p]].push_back(e); dfs(e,p); black_stack[H[p]].pop_back(); } } int main(){ cin >> N; E.resize(N); hldE.resize(N); rep(i,N-1){ int u,v; u32 c; cin >> u >> v >> c; u--; v--; E[u].push_back({ v,c }); E[v].push_back({ u,c }); hldE[u].push_back(v); hldE[v].push_back(u); } hld = heavy_light_decomposition(hldE); vector I = {0}; P.assign(N,-1); H.assign(N,0); rep(i,I.size()){ int p = I[i]; for(auto e : E[p]){ if(P[p] == e.to) continue; I.push_back(e.to); P[e.to] = p; H[e.to] = H[p] ^ e.h; } } Z.assign(N,1); for(int i=N-1; i>=1; i--) Z[P[I[i]]] += Z[I[i]]; rep(i,N) H[i] = CCH.push(H[i]); rep(i,N) H[i] = CCH[H[i]]; dfs(); //rep(i,N) cout << prev_black[i] << " "; cout << endl; //rep(i,N) cout << block_size[i] << " "; cout << endl; i64 ans = 0; for(int p=1; p