// {{{
#include <algorithm>
#include <array>
#include <bitset>
#include <cassert>
#include <cmath>
#include <cstdio>
#include <iostream>
#include <iterator>
#include <list>
#include <map>
#include <numeric>
#include <queue>
#include <random>
#include <set>
#include <sstream>
#include <stack>
#include <string>
#include <tuple>
#include <utility>
#include <vector>
// }}}

using namespace std;

using ll = long long;
using ull = unsigned long long;
using ld = long double;

static constexpr int mod = (int)1e9 + 7;
static constexpr int inf = 100100100;
static constexpr ll linf = 1e18;
static constexpr double eps = 1e-9;
static constexpr double pi = 3.14159265359;

#define rep(i, n) for (ll i = 0; i < n; ++i)
#define rrep(i, n) for (ll i = n; i >= 0; --i)
#define all(c) begin(c), end(c)
#define rall(c) rbegin(c), rend(c)
#define pb push_back
#define ist insert
#define fst first
#define snd second


// {{{
template <typename T>
struct modint {
    static_assert(std::is_integral<T>::value, "T must be integral");
    modint(T value = 0, T mod = (T)1e9 + 7) : value(value % mod), mod(mod) {}
    modint operator +(const modint& rhs) const {
        return modint(*this) += rhs;
    }
    modint operator -(const modint& rhs) const {
        return modint(*this) -= rhs;
    }
    modint operator *(const modint& rhs) const {
        return modint(*this) *= rhs;
    }
    modint operator /(const modint& rhs) const {
        return modint(*this) /= rhs;
    }
    modint& operator +=(const modint& rhs) {
        value += rhs.value;
        if (value >= mod) {
            value -= mod;
        }
        return *this;
    }
    modint& operator -=(const modint& rhs) {
        if (value < rhs.value) {
            value += mod;
        }
        value -= rhs.value;
        return *this;
    }
    modint& operator *=(const modint& rhs) {
        value = value * rhs.value % mod;
        return *this;
    }
    modint& operator /=(const modint& rhs) {
        modint tmp = rhs;
        T exp = mod - 2;
        while (exp) {
            if (exp % 2) {
                *this *= tmp;
            }
            tmp *= tmp;
            exp /= 2;
        }
        return *this;
    }
    bool operator ==(const modint& rhs) const {
        return value == rhs.value;
    }
    bool operator !=(const modint& rhs) const {
        return !(*this == rhs);
    }
    bool operator >(const modint& rhs) const {
        return value > rhs.value;
    }
    bool operator <(const modint& rhs) const {
        return value < rhs.value;
    }
    bool operator >=(const modint& rhs) const {
        return !(*this < rhs);
    }
    bool operator <=(const modint& rhs) const {
        return !(*this > rhs);
    }
    friend ostream& operator <<(ostream &out, const modint& m) {
        out << m.value;
        return out;
    }
    friend istream& operator >>(istream &in,  modint& m) {
        in >> m.value;
        m.value %= m.mod;
        return in;
    }
    T value;
    T mod;
};
// }}}

modint<ll> dfs(ll K, const vector<vector<ll>>& G, ll cur, ll par) {
    ll k;
    if (par == -1) {
        k = K - 1;
    } else {
        k = K - 2;
    }
    if (K < G[cur].size()) {
        return 0;
    }
    modint<ll> c = 1;
    for (auto e : G[cur]) {
        if (e == par) continue;
        c *= k;
        --k;
    }
    for (auto e : G[cur]) {
        if (e == par) continue;
        c *= dfs(K, G, e, cur);
    }
    return c;
}

int main() {
    // cin.tie(0);
    // ios_base::sync_with_stdio(false);
    ll N, K;
    cin >> N >> K;

    vector<vector<ll>> G(N);
    rep (i, N - 1) {
        ll a, b;
        cin >> a >> b;
        --a; --b;
        G[a].pb(b);
        G[b].pb(a);
    }
    modint<ll> ans = dfs(K, G, 0, -1) * K;
    cout << ans << endl;
}


Measured memory_bytes: 11648