#include using namespace std; struct iofast_t { iofast_t() { ios::sync_with_stdio(false); cin.tie(nullptr); } } iofast; struct uns_t {} uns; template auto vec(Element init, Head arg, Args ...args) { if constexpr (sizeof...(Args) == 0) return vector(arg, init); else return vector(arg, vec(init, args...)); } template auto vec(uns_t, Head arg, Args ...args) { return vec(Element(), arg, args...); } template auto distance(const Container &c, decltype(begin(c)) iter) { return distance(begin(c), iter); } template ::value_type>> auto isort(RIter first, RIter last, Compare comp = Compare()) { vector i(distance(first, last)); iota(begin(i), end(i), 0); sort(begin(i), end(i), [&](auto x, auto y) { return comp(*(first + x), *(first + y)); }); return i; } template typename, typename = void_t<>> struct detect : false_type {}; template typename Check> struct detect>> : true_type {}; template typename Check> constexpr inline bool detect_v = detect::value; template using has_member_sort = decltype(declval().sort()); template > auto sorted(Container c, Compare comp = Compare()) { if constexpr (detect_v) { c.sort(comp); return c; } else { sort(begin(c), end(c), comp); return c; } } template > auto uniqued(Container c, Compare comp = Compare()) { c.erase(unique(begin(c), end(c), comp), end(c)); return c; } template > T &chmin(T &l, T r, Compare &&f = less()) { return l = min(l, r, f); } template > T &chmax(T &l, T r, Compare &&f = less()) { return l = max(l, r, f); } template constexpr auto fix(F &&f) noexcept { return [f = std::tuple(std::forward(f))](auto &&...args) mutable { return std::get<0>(f)(fix(std::get<0>(f)), std::forward(args)...); }; } #include #include using mint = atcoder::modint1000000007; int main() { int n, m, x; cin >> n >> m >> x; atcoder::dsu dsu(n); auto g = vec>(uns, n, 0); for (int i = 0; i < m; ++i) { int x, y, z; cin >> x >> y >> z; --x; --y; if (!dsu.same(x, y)) { g[x].push_back({ y, z }); g[y].push_back({ x, z }); dsu.merge(x, y); } } auto sub = vec(uns, n); fix([&](auto f, int v, int p = -1) -> int { sub[v] = 1; for (auto [u, _] : g[v]) { if (u == p) { continue; } sub[v] += f(u, v); } return sub[v]; })(0); mint ans = 0; fix([&](auto f, int v, int p = -1, int z = -1) -> void { if (0 <= p) { ans += mint(x).pow(z) * sub[v] * (n - sub[v]); } for (auto [u, w] : g[v]) { if (u == p) { continue; } f(u, v, w); } })(0); cout << ans.val() << endl; }