#include #include using namespace std; using namespace atcoder; const long nPrime = 1000000007; //const long nPrime = 998244353; typedef long long ll; struct UnionFind { vector size, parents; UnionFind() {} UnionFind(int n) { size.resize(n, 1); parents.resize(n, 0); for (int i = 0; i < n; i++) { MakeTree(i); } } void MakeTree(int x) { parents[x] = x; } bool IsSame(int x, int y) { return FindRoot(x) == FindRoot(y); } void Unite(int x, int y) { x = FindRoot(x); y = FindRoot(y); if(x==y){ return; } if (size[x] > size[y]) { parents[y] = x; size[x] += size[y]; size[y] = -1; } else { parents[x] = y; size[y] += size[x]; size[x] = -1; } } int FindRoot(int x) { if (x != parents[x]) parents[x] = FindRoot(parents[x]); return parents[x]; } int GetSize(int x){ return size[FindRoot(x)]; } }; const long long nMod = nPrime; class ModInt { long long x; public: ModInt(long long x=0) : x((x%nMod+nMod)%nMod) {} ModInt operator-() const { return ModInt(-x); } ModInt& operator+=(const ModInt& a) { if ((x += a.x) >= nMod) x -= nMod; return *this; } ModInt& operator-=(const ModInt& a) { if ((x += nMod-a.x) >= nMod) x -= nMod; return *this; } ModInt& operator*=(const ModInt& a) { (x *= a.x) %= nMod; return *this; } ModInt operator+(const ModInt& a) const { ModInt res(*this); return res+=a; } ModInt operator-(const ModInt& a) const { ModInt res(*this); return res-=a; } ModInt operator*(const ModInt& a) const { ModInt res(*this); return res*=a; } ModInt Pow(long long t) const { if (!t) return 1; ModInt a = Pow(t>>1); a *= a; if (t&1) a *= *this; return a; } // for prime nMod ModInt Inv() const { return Pow(nMod-2); } ModInt& operator/=(const ModInt& a) { return (*this) *= a.Inv(); } ModInt operator/(const ModInt& a) const { ModInt res(*this); return res/=a; } friend ostream& operator<<(ostream& os, const ModInt& m){ os << m.x; return os; } }; void MakeTable(vector &viInv, vector &viFctr, vector &viFctrInv, bool bUseFctr = false, const long nTableSize = 1000000){ viInv.resize(nTableSize+1, 1); for(long i = 1; i < nTableSize+1; i++){ viInv[i] /= i; } if(!bUseFctr){ return; } viFctr.resize(nTableSize+1, 1); viFctrInv.resize(nTableSize+1, 1); for(long i = 1; i < nTableSize+1; i++){ viFctr[i] = viFctr[i-1] * i; } viFctrInv[nTableSize] /= viFctr[nTableSize]; for(long i = nTableSize; i >= 1; i--){ viFctrInv[i-1] = viFctrInv[i]*i; } } int main() { long n; cin >> n; vector,long>> viEdge(n-1); for(long i = 0; i < n-1; i++){ long a,b,c; cin >> a >> b >> c; a--,b--; viEdge[i] = make_pair(make_pair(a,b),c); } ModInt nAns = 0; vector vzGroup(31,UnionFind(n)); for(long i = 0; i <= 30; i++){ for(long j = 0; j < n-1; j++){ long a = viEdge[j].first.first; long b = viEdge[j].first.second; long c = viEdge[j].second; if((c&(1<