#include <bits/stdc++.h>

using namespace std;
typedef long long ll;
#define int ll

#define REP(i,n) for(int i=0;i<n;++i)
#define ALL(name) name.begin(),name.end()
#define SORT(name) sort(name.begin(), name.end())
#define ZERO(p) memset(p, 0, sizeof(p))
#define MINUS(p) memset(p, -1, sizeof(p))
#if 1
#  define DBG(fmt, ...) printf(fmt, ##__VA_ARGS__)
#else
#  define DBG(fmt, ...)
#endif

const ll LLINF = (1LL<<60);
const int INF = (1LL<<30);
const double DINF = std::numeric_limits<double>::infinity();
const int MOD = 1000000007;
#define MAX_N 100010

class NumberOfCase {
    public:
        // 初期化 O(n log m) m は素数
        // n : 求めたい組み合わせの最大値, m : 剰余(素数)
        NumberOfCase(ll n, ll m) : m_Size(n), m_Mod(m) {
            m_Fact.resize(m_Size + 1);
            m_FInv.resize(m_Size + 1);
            m_Inv.resize(m_Size + 1);
            // 0! ～ m_Size! までの階乗を保管する配列の作成 (mod m_Mod をとってある)
            m_Fact[0] = 1;
            for(int i = 1; i <= m_Size; ++i) { m_Fact[i] = (m_Fact[i-1] * i) % m_Mod; }
            // m_Size! の逆元 1/m_Size! を二分累乗法で求める
            // フェルマーの小定理より a^-1 = a^(p-2) mod p -> (m_Size!)^-1 = (m_Size!)^(p-2) mod p
            m_FInv[m_Size] = Pow(m_Fact[m_Size], m_Mod-2);
            // 以下 1 / (N-1)! = N / N! であることを利用
            for(int i = m_Size - 1; i >= 0; --i) { m_FInv[i] = (m_FInv[i + 1] * (i + 1)) % m_Mod; }
            // 逆元の作成（割り算のとき利用）
            m_Inv[1] = 1;
            for(int i = 2; i <= m_Size; ++i) { m_Inv[i] = m_Inv[m_Mod % i] * (m_Mod - m_Mod / i) % m_Mod; }
        }
        // x の y 乗を求める O(log MOD)
        ll Pow(ll x, ll y) {
            if (y == 0) return 1;
            ll res = Pow(x, y / 2);
            return res * res % m_Mod * (y & 1 ? x : 1) % m_Mod;
        }
        // 階乗 n! の取得 O(1)
        ll Factorial(ll n) { return m_Fact[n]; }
        // n! の逆元の取得 O(1)
        ll FInverse(ll n) { return m_FInv[n]; }
        // 逆元の取得 O(1) or O(log MOD)
        // m_Size を超える値も可能。ただし O(log M)
        ll Inverse(ll n) {
            if(n < m_Size) { return m_Inv[n]; }
            else { return Pow(n, m_Mod - 2); }
        }
        // 剰余を考慮した割り算 n / m の結果を返す
        ll Div(ll n, ll m) { return (n * m_Inv[m]) % m_Mod; }
        // 順列 nPr の取得 O(1)
        ll Permutation(ll n, ll r) { return (m_Fact[n] * m_FInv[n - r]) % m_Mod; }
        // 組み合わせ nCr の取得 O(1)
        ll Combination(ll n, ll r) {
            ll ret = (m_Fact[n] * m_FInv[n - r]) % m_Mod;
            ret = (ret * m_FInv[r]) % m_Mod;
            return ret;
        }
        // 重複組み合わせ nHr の取得 O(1)
        ll Homogen(ll n, ll r) { return Combination(r+n-1, r); }
    private:
        const ll    m_Size;     // 導出する範囲 (10^5 くらいまでにして欲しい)
        const ll    m_Mod;      // 剰余
        vector<ll>  m_Fact;     // 階乗
        vector<ll>  m_FInv;     // 階乗の逆元
        vector<ll>  m_Inv;      // 逆元
};

ll N, K;
vector<ll> graph[MAX_N];
NumberOfCase noc(MAX_N, MOD);

// 頂点 i まで決定済みで、それ以降の子を見た時のパターン数
ll dfs(ll i, ll parent) {
    ll ret = 0;
    if(i == 0) {
        // 根の時は親がいない
        ll child_cnt = graph[i].size();
        if(K - 1 < child_cnt) { return 0; }
        ret = noc.Permutation(K - 1, child_cnt);
    }
    else {
        ll child_cnt = graph[i].size() - 1;
        if(K - 2 < child_cnt) { return 0; }
        ret = noc.Permutation(K - 2, child_cnt);
    }
    for(auto& child : graph[i]) {
        if(child == parent) { continue; }
        ret *= dfs(child, i);
        ret %= MOD;
    }
    return ret;
}

signed main() {
    cin.tie(0); cout.tie(0);
    ios::sync_with_stdio(false);

    cin >> N >> K;
    REP(i, N-1) {
        ll a, b;
        cin >> a >> b;
        a--; b--;
        graph[a].push_back(b);
        graph[b].push_back(a);
    }

    ll ans = (K * dfs(0, -1)) % MOD;
    printf("%lld\n", ans);

    return 0;
}


Measured memory_bytes: 12800