#include <iostream>
#include <vector>
using ull = unsigned long long;
constexpr ull M = 1000000007;
int main()
{
    std::size_t N;
    std::cin >> N;
    std::vector<ull> L(N + 1, 1), R(N + 1, 1);
    for (std::size_t i = 1; i <= N; i++) { L[i] = L[i - 1] * (ull)(i == 1 ? 1 : i - 1) % M; }
    for (std::size_t i = N; i >= 1; i--) { R[i - 1] = R[i] * (ull)(i) % M; }
    std::vector<std::vector<std::size_t>> g(N);
    for (std::size_t i = 0; i < N - 1; i++) {
        std::size_t u, v;
        std::cin >> u >> v, u--, v--;
        g[u].push_back(v), g[v].push_back(u);
    }
    std::vector<std::size_t> dnum(N + 1, 0);
    auto dfs = [&](auto&& self, const std::size_t s, const std::size_t p, const std::size_t d) -> void {
        dnum[d]++;
        for (const std::size_t to : g[s]) {
            if (to == p) { continue; }
            self(self, to, s, d + 1);
        }
    };
    dfs(dfs, 0, N, 1);
    ull ans = 0;
    for (std::size_t d = 1; d <= N; d++) { (ans += (ull)dnum[d] * (L[d] * R[d] % M) % M) %= M; }
    std::cout << ans << std::endl;
    return 0;
}