//@formatter:off #include #define overload4(_1,_2,_3,_4,name,...) name #define rep1(i,n) for (ll i = 0; i < ll(n); ++i) #define rep2(i,s,n) for (ll i = ll(s); i < ll(n); ++i) #define rep3(i,s,n,d) for(ll i = ll(s); i < ll(n); i+=d) #define rep(...) overload4(__VA_ARGS__,rep3,rep2,rep1)(__VA_ARGS__) #define rrep(i,n) for (ll i = ll(n)-1; i >= 0; i--) #define all(a) a.begin(),a.end() #define rall(a) a.rbegin(),a.rend() #define pb push_back #define eb emplace_back #define vi vector #define vvi vector> #define vl vector #define vvl vector> #define vd vector #define vvd vector> #define vs vector #define vc vector #define vvc vector> #define vb vector #define vvb vector> #define vp vector

#define vvp vector> #ifdef __LOCAL #define debug(...) { cout << #__VA_ARGS__; cout << ": "; print(__VA_ARGS__); } #else #define debug(...) void(0) #endif #define INT(...) int __VA_ARGS__;scan(__VA_ARGS__) #define LL(...) ll __VA_ARGS__;scan(__VA_ARGS__) #define STR(...) string __VA_ARGS__;scan(__VA_ARGS__) #define CHR(...) char __VA_ARGS__;scan(__VA_ARGS__) #define DBL(...) double __VA_ARGS__;scan(__VA_ARGS__) #define LD(...) ld __VA_ARGS__;scan(__VA_ARGS__) using namespace std; using ll = long long; using P = pair; using LP = pair; template istream& operator>>(istream &is,pair &p) { return is >> p.first >> p.second; } template ostream& operator<<(ostream &os,const pair &p) { return os<<'{'< istream& operator>>(istream &is,vector &v) { for(T &t:v){is>>t;} return is; } template ostream& operator<<(ostream &os,const vector &v) { os<<'[';rep(i,v.size())os< void fin(T a) { cout << a << '\n'; exit(0); } template bool chmin(T& a,T b) {if(a > b){a = b; return true;} return false;} template bool chmax(T& a,T b) {if(a < b){a = b; return true;} return false;} void scan(){} template void scan(Head& head, Tail&... tail){ cin >> head; scan(tail...); } template void print(const T& t){ cout << t << '\n'; } template void print(const Head& head, const Tail&... tail){ cout<= mod) x -= mod; return *this; } constexpr mint &operator++() { return *this += mint(1); } constexpr mint &operator-=(const mint &a) { if ((x += mod - a.x) >= mod) x -= mod; return *this; } constexpr mint &operator--() { return *this -= mint(1); } constexpr mint &operator*=(const mint &a) { (x *= a.x) %= mod; return *this; } constexpr mint operator+(const mint &a) const { mint res(*this); return res += a; } constexpr mint operator-(const mint &a) const { mint res(*this); return res -= a; } constexpr mint operator*(const mint &a) const { mint res(*this); return res *= a; } constexpr mint pow(ll t) const { mint res = mint(1), a(*this); while (t > 0) { if (t & 1) res *= a; t >>= 1; a *= a; } return res; } // for prime mod constexpr mint inv() const { return pow(mod - 2); } constexpr mint &operator/=(const mint &a) { return *this *= a.inv(); } constexpr mint operator/(const mint &a) const { mint res(*this); return res /= a; } }; ostream &operator<<(ostream &os, const mint &a) { return os << a.x; } bool operator==(const mint &a, const mint &b) { return a.x == b.x; } bool operator!=(const mint &a, const mint &b) { return a.x != b.x; } bool operator==(const mint &a, const int &b) { return a.x == b; } bool operator!=(const mint &a, const int &b) { return a.x != b; } template class rerooting { int n; vvp tree; T identity; MERGE merge; ADDROOT addRoot; ADDEDGE addEdge; vector> dp; vector ans; T dfs(int v = 0, int p = -1) { T sum = identity; dp[v].resize(tree[v].size()); rep(i, tree[v].size()) { auto[u, w] = tree[v][i]; if (u == p) continue; dp[v][i] = dfs(u, v); sum = merge(sum, addEdge(dp[v][i], w)); } return addRoot(sum); } void dfs2(T dpP, int v = 0, int p = -1) { int sz = tree[v].size(); rep(i, sz) if (tree[v][i].first == p) dp[v][i] = dpP; vector sumL(sz + 1, identity), sumR(sz + 1, identity); rep(i, sz) sumL[i + 1] = merge(sumL[i], addEdge(dp[v][i], tree[v][i].second)); rrep(i, sz) sumR[i] = merge(sumR[i + 1], addEdge(dp[v][i], tree[v][i].second)); ans[v] = addRoot(sumL[sz]); rep(i, sz) { auto[u, w] = tree[v][i]; if (u == p) continue; T t = merge(sumL[i], sumR[i + 1]); dfs2(addRoot(t), u, v); } } void init() { dfs(); dfs2(identity); } public: rerooting(int n, vvp tree, T identity, MERGE merge, ADDROOT addRoot, ADDEDGE addEdge) : n(n), tree(tree), identity(identity), merge(merge), addRoot(addRoot), addEdge(addEdge), dp(n), ans(n) { init(); }; T get_ans(int i) { return ans[i]; } }; struct state { mint sq_sum, sum; int cnt; state(mint sq_sum = 0, mint sum = 0, int cnt = 0) : sq_sum(sq_sum), sum(sum), cnt(cnt) {} }; int main() { ios::sync_with_stdio(false); cin.tie(nullptr); cout.tie(nullptr); INT(n); vvp G(n); rep(_, n - 1) { INT(u, v, w); u--; v--; G[u].eb(v, w); G[v].eb(u, w); } auto merge = [](const state &a, const state &b) { return state(a.sq_sum + b.sq_sum, a.sum + b.sum, a.cnt + b.cnt); }; auto addRoot = [](const state &a) { return state(a.sq_sum, a.sum, a.cnt + 1); }; auto addEdge = [](const state &a, int weight) { return state(a.sq_sum + a.sum * weight * 2 + mint(weight) * weight * a.cnt, a.sum + mint(weight) * a.cnt, a.cnt); }; rerooting rt(n, G, state(0, 0, 0), merge, addRoot, addEdge); mint ans = 0; rep(i, n) ans += rt.get_ans(i).sq_sum; print(ans / 2); }