#include <bits/stdc++.h>

#define N 500010 

using namespace std;

const int mod = 1000000007 ;

char *p1, *p2, buf[100000];

typedef unsigned long long ll;

ll bs1 = 20021214 ;

ll bs2 = 20030315 ;

#define nc() (p1 == p2 && (p2 = (p1 = buf) + fread(buf, 1, 100000, stdin), p1 == p2) ? EOF : *p1 ++ )

int rd() {
    int x = 0;
    char c = nc();
    while (c < 48) {
        c = nc();
    }
    while (c > 47) {
        x = (((x << 2) + x) << 1) + (c ^ 48), c = nc();
    }
    return x;
}

int n, inv[N], to[N << 1], nxt[N << 1], head[N], tot, sz[N], w[N], root, rt1, rt2, flag, f[2][N], tmp[N];

ll hs[N];

inline void add(int x, int y) {
    to[ ++ tot] = y;
    nxt[tot] = head[x];
    head[x] = tot;
}

void getroot(int p, int fa) {
    sz[p] = 1;
    w[p] = 0;
    for (int i = head[p]; i; i = nxt[i]) {
        if(to[i] != fa) {
            getroot(to[i], p);
            sz[p] += sz[to[i]];
            w[p] = max(w[p], sz[to[i]]);
        }
    }
    w[p] = max(w[p], n - sz[p]);
    if (w[p] < w[root]) {
        root = p;
    }
}

int C(int n, int m) {
    int re = 1;
    for (int i = n - m + 1; i <= n; i ++ ) {
        re = (ll)re * i % mod;
    }
    for (int i = 1; i <= m; i ++ ) {
        re = (ll)re * inv[i] % mod;
    }
    return re;
}

inline bool cmp(int i, int j) {
    return hs[i] < hs[j];
}

void dfs(int p, int fa) {
    sz[p] = f[0][p] = f[1][p] = 1;
    for (int i = head[p]; i; i = nxt[i]) {
        if (to[i] != fa) {
            dfs(to[i], p);
            sz[p] += sz[to[i]];
        }
    }
    int len = 0;
    for (int i = head[p]; i; i = nxt[i]) {
        if (to[i] != fa) {
            tmp[ ++ len] = to[i];
        }
    }
    sort(tmp + 1, tmp + len + 1, cmp);
    for (int i = 1, j = 1; i <= len; i = j) {
        while (j <= len && hs[tmp[j]] == hs[tmp[i]]) {
            j ++ ;
        }
        f[0][p] = (ll)f[0][p] * C(f[0][tmp[i]] + f[1][tmp[i]] + j - i - 1, j - i) % mod;
        f[1][p] = (ll)f[1][p] * C(f[0][tmp[i]] + j - i - 1, j - i) % mod;
    }
    hs[p] = bs2 * len + sz[p];
    for (int i = 1; i <= len; i ++ ) {
        hs[p] = (hs[p] * bs1) ^ hs[tmp[i]];
    }
}

int main() {
    n = rd();
    inv[0] = inv[1] = 1;
    for (int i = 2; i <= n; i ++ ) {
        inv[i] = (ll)inv[mod % i] * (mod - mod / i) % mod;
    }
    for (int i = 1; i < n; i ++ ) {
        int x = rd(), y = rd();
        add(x, y), add(y, x);
    }
    w[0] = n;
    getroot(1, 0);
    getroot(root, 0);
    for (int i = head[root], last = 0; i; last = i, i = nxt[i]) {
        if (sz[to[i]] * 2 == n) {
            n ++ ;
            if (i == head[root]) {
                head[root] = nxt[i];
            }
            else {
                nxt[last] = nxt[i];
            }
            for (int j = head[to[i]], lst = 0; j; lst = j, j = nxt[j]) {
                if (to[j] == root) {
                    if (j == head[to[i]]) {
                        head[to[i]] = nxt[j];
                    }
                    else {
                        nxt[lst] = nxt[j];
                    }
                    break;
                }
            }
            add(n, root);
            add(root, n);
            add(n, to[i]);
            add(to[i], n);
            rt1 = root;
            rt2 = to[i];
            root = n;
            flag = 1;
            break;
        }
    }
    dfs(root, 0);
    if (!flag) {
        cout << (f[0][root] + f[1][root]) % mod << endl ;
    }
    else if (hs[rt1] == hs[rt2]) {
        cout << (f[0][root] - C(f[1][rt1] + 1, 2) + mod) % mod;
    }
    else {
        cout << ((ll)f[0][rt1] * f[0][rt2] + (ll)f[0][rt1] * f[1][rt2] + (ll)f[1][rt1] * f[0][rt2]) % mod << endl ;
    }
}