import sys
from collections import deque

mod = 10**9+7

input = sys.stdin.readline

N = int(input())
edge = [[] for i in range(N)]
for i in range(N-1):
    u,v,w = map(int,input().split())
    edge[u-1].append((v-1,w))
    edge[v-1].append((u-1,w))

deq = deque([0])
parent = [-1]*N
parent[0] = (0,0)
res = []
while deq:
    v = deq.popleft()
    res.append(v)
    for nv,w in edge[v]:
        if parent[nv]==-1:
            parent[nv] = (v,w)
            deq.append(nv)

res = res[::-1]
size = [1 for v in range(N)]
dist_sum = [0 for v in range(N)]
dist_sum_square = [0 for v in range(N)]
for v in res:
    for nv,w in edge[v]:
        if nv!=parent[v][0]:
            size[v] += size[nv]
            dist_sum[v] += dist_sum[nv] + w * size[nv]
            dist_sum[v] %= mod
            dist_sum_square[v] += dist_sum_square[nv] + 2 * w * dist_sum[nv] + size[nv] * w * w
            dist_sum_square[v] %= mod

#4 9 1 3
ans = [0 for v in range(N)]
ans[0] = dist_sum_square[0]
dist_sum_ans = [0 for v in range(N)]
dist_sum_ans[0] = dist_sum[0]
res = res[::-1]
for v in res:
    if v==0:
        continue
    pv,w = parent[v]
    tmp = ans[pv] - ( dist_sum_square[v] + 2 * w * dist_sum[v] + size[v] * w * w )
    tmp_dist_sum = dist_sum_ans[pv] - (w * size[v] + dist_sum[v])
    ans[v] = dist_sum_square[v] + tmp + 2 * w * (tmp_dist_sum) + (N - size[v]) * w * w
    dist_sum_ans[v] = dist_sum[v] + w * (N - size[v]) + tmp_dist_sum
    ans[v] %= mod
    dist_sum_ans[v] %= mod

S = 0
for v in range(N):
    S += ans[v]
    S %= mod

S *= pow(2,mod-2,mod)
print(S%mod)