#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
typedef unsigned long long ull;
typedef pair<int,int> pii;
#define debug(fmt,...) fprintf(stderr, fmt,##__VA_ARGS__)
#define debugln(fmt,...) fprintf(stderr,"[%d] " fmt,__LINE__,##__VA_ARGS__)
#define TIME (1e3*clock()/CLOCKS_PER_SEC)
bool memBeg;
template<typename T> void chkmin(T &x,T y) {x=min(x,y);}
template<typename T> void chkmax(T &x,T y) {x=max(x,y);}
constexpr int mod=1e9+7;
void inc(int &x,int y) {x+=y; x-=(x>=mod)*mod;}
void dec(int &x,int y) {x-=y; x+=(x<0)*mod;}
constexpr int add(int x,int y) {return (x+y>=mod)?x+y-mod:x+y;}
constexpr int sub(int x,int y) {return (x<y)?x-y+mod:x-y;}
constexpr int quick_pow(int x,ll times,int ret=1) {
    for(;times;times>>=1,x=1ll*x*x%mod) if(times&1) ret=1ll*ret*x%mod;
    return ret;
}
constexpr int maxn=805;
int n,head[maxn],to[maxn<<1],nxt[maxn<<1],tot;
void add_edge(int u,int v) {
    to[tot]=v; nxt[tot]=head[u]; head[u]=tot++;
    to[tot]=u; nxt[tot]=head[v]; head[v]=tot++;
}
vector<int> flr[maxn<<1][maxn];
void dfs0(int u,int pre,int anc,int dep) {
    flr[anc][dep].emplace_back(u);
    for(int i=head[u];~i;i=nxt[i]) {
        if((i^1)==pre) continue;
        dfs0(to[i],i,anc,dep+1);
    }
}
int wys[maxn<<1][maxn],nwys[maxn<<1][maxn],sav[maxn<<1][maxn],fac[maxn],ifac[maxn],C[maxn][maxn];
int binom(int x,int y) {
    if(x<0||y<0||x<y) return 0;
    return C[x][y];
}
bool memEn;
void fl() {
    freopen("tree.in","r",stdin);
    freopen("tree.out","w",stdout);
}
int main() {
    debug("%.24lf\n",fabs(&memEn-&memBeg)/1024.0/1024.0);
    // fl();
// #ifdef LOCAL
//     freopen("debug","w",stderr);
// #endif
    scanf("%d",&n);
    if(n==1) return puts("0"),0;
    memset(head,-1,sizeof(head));
    for(int i=1,u,v;i<n;i++) {
        scanf("%d%d",&u,&v);
        add_edge(u,v);
    }
    fac[0]=1;
    for(int i=1;i<=n;i++) fac[i]=1ll*fac[i-1]*i%mod;
    ifac[n]=quick_pow(fac[n],mod-2);
    for(int i=n-1;i>=0;i--) ifac[i]=1ll*ifac[i+1]*(i+1)%mod;
    for(int i=0;i<=n;i++) {
        C[i][0]=1;
        for(int j=1;j<=i;j++) {
            C[i][j]=add(C[i-1][j-1],C[i-1][j]);
        }
    }
    for(int i=0;i<tot;i++) dfs0(to[i],i,i,0);
    for(int i=0;i<tot;i++) wys[i][1]=1;
    int ret=0;
    for(int i=0;i<=(n-2)/2;i++) {
        memset(nwys,0,sizeof(nwys));
        memset(sav,0,sizeof(sav));
        for(int j=0;j<tot;j++) {
            int u=to[j^1],v=to[j];
            if(!flr[j^1][i+1].empty()&&!flr[j][i+1].empty()) {
                int b=(int)flr[j][i].size()-1;
                for(int a=1;a<=(int)flr[j^1][i].size();a++) {
                    if(!wys[j][a]) continue;
                    int A=(int)flr[j^1][i].size()-a;
                    for(int c=1;c<=(int)flr[j^1][i+1].size();c++) {
                        int lhs=1ll*fac[A]*fac[c]%mod*binom((int)flr[j^1][i+1].size(),c)%mod,rhs=fac[b];
                        int coef=1ll*lhs*rhs%mod*binom(A+c+b,b)%mod;
                        nwys[j][c]=(nwys[j][c]+1ll*wys[j][a]*coef%mod*(int)flr[j][i+1].size())%mod;
                    }
                    for(int d=1;d<=(int)flr[j][i+1].size();d++) {
                        int lhs=fac[A],rhs=1ll*fac[b]*fac[d]%mod*binom((int)flr[j][i+1].size(),d)%mod;
                        int coef=1ll*lhs*rhs%mod*binom(A+b+d,A)%mod;
                        nwys[j^1][d]=(nwys[j^1][d]+1ll*wys[j][a]*coef%mod*(int)flr[j^1][i+1].size())%mod;
                    }
                }
            }
            if(!flr[j][i+2].empty()) {
                int ext1=(int)flr[j][i].size()-1;
                static int hlp[maxn];
                memset(hlp,0,sizeof(int)*(1+flr[j^1][i].size()));
                int rhs=fac[ext1];
                for(int dt=0;dt<(int)flr[j^1][i].size();dt++) {
                    int coef=1ll*fac[dt]*rhs%mod;
                    if(dt) {
                        coef=1ll*coef*binom(dt+ext1-1,dt)%mod;
                    }
                    for(int k=1;k+dt<=(int)flr[j^1][i].size();k++) {
                        hlp[k+dt]=(hlp[k+dt]+1ll*wys[j][k]*coef%mod*binom((int)flr[j^1][i].size()-k,dt))%mod;
                    }
                }
                for(int e=head[v];~e;e=nxt[e]) {
                    if((e^1)==j) continue;
                    int off=flr[j][i+1].size()-flr[e][i].size();
                    for(int k=1;k<=(int)flr[j^1][i].size();k++) {
                        sav[e][k]=(sav[e][k]+1ll*hlp[k]*fac[flr[j^1][i].size()-k+off])%mod;
                    }
                }
            }
            if(!flr[j^1][i+2].empty()) {
                static int hlp[maxn];
                memset(hlp,0,sizeof(int)*(1+flr[j][i].size()));
                for(int a=1;a<=(int)flr[j^1][i].size();a++) {
                    if(!wys[j][a]) continue;
                    int ext1=flr[j^1][i].size()-a;
                    int lhs=fac[ext1];
                    for(int dt=0;dt<(int)flr[j][i].size();dt++) {
                        int coef=1ll*lhs*fac[dt]%mod;
                        if(dt) {
                            coef=1ll*coef*binom(dt+ext1-1,dt)%mod;
                        }
                        hlp[dt+1]=(hlp[dt+1]+1ll*wys[j][a]*coef%mod*binom((int)flr[j][i].size()-1,dt))%mod;
                    }
                }
                for(int e=head[u];~e;e=nxt[e]) {
                    if(e==j) continue;
                    int off=flr[j^1][i+1].size()-flr[e][i].size();
                    for(int k=1;k<=(int)flr[j][i].size();k++) {
                        sav[e][k]=(sav[e][k]+1ll*hlp[k]*fac[flr[j][i].size()-k+off])%mod;
                    }
                }
            }
            if(flr[j^1][i+1].empty()&&flr[j][i+1].empty()) {
                for(int a=1;a<=(int)flr[j^1][i].size();a++) {
                    if(wys[j][a]) {
                        int len=(int)flr[j^1][i].size()-a+(int)flr[j][i].size()-1;
                        ret=(ret+1ll*wys[j][a]*fac[len])%mod;
                    }
                }
            } else if(flr[j^1][i+1].empty()&&flr[j][i+2].empty()) {
                for(int a=1;a<=(int)flr[j^1][i].size();a++) {
                    if(wys[j][a]) {
                        int A=(int)flr[j^1][i].size()-a,b=(int)flr[j][i].size()-1,d=flr[j][i+1].size();
                        int coef=1ll*fac[A]*fac[b]%mod*fac[d]%mod*binom(A+b+d,A)%mod;
                        ret=(ret+1ll*wys[j][a]*coef)%mod;
                    }
                }
            } else if(flr[j^1][i+2].empty()&&flr[j][i+1].empty()) {
                for(int a=1;a<=(int)flr[j^1][i].size();a++) {
                    if(wys[j][a]) {
                        int A=(int)flr[j^1][i].size()-a,c=flr[j^1][i+1].size(),b=(int)flr[j][i].size()-1;
                        int coef=1ll*fac[A]*fac[b]%mod*fac[c]%mod*binom(A+b+c,b)%mod;
                        ret=(ret+1ll*wys[j][a]*coef)%mod;
                    }
                }
            }
            // debug("nwys[2][1] = %d\n",nwys[2][1]);
        }
        for(int j=0;j<tot;j++) {
            static int g[maxn],tmp[maxn];
            int len=flr[j^1][i+1].size();
            for(int dt=0;dt<len;dt++) {
                g[dt]=1ll*fac[(int)flr[j][i].size()+dt]*ifac[dt]%mod*flr[j][i+1].size()%mod;
            }
            for(int k=1;k<=len;k++) tmp[k]=0;
            for(int k=1;k<=len;k++) {
                for(int dt=0;k+dt<=len;dt++) {
                    tmp[k+dt]=(tmp[k+dt]+1ll*sav[j][k]*g[dt])%mod;
                }
            }
            for(int k=1;k<=len;k++) {
                nwys[j][k]=(nwys[j][k]+1ll*tmp[k]*ifac[len-k])%mod;
            }
        }
        for(int j=0;j<tot;j++) {
            memcpy(wys[j],nwys[j],sizeof(int)*(n+1));
        }
    }
    printf("%d\n",ret);
    return 0;
}
/*
g++ template.cpp -o template -std=c++14 -Wall -Wextra -Wshadow -O2 -fsanitize=address,undefined
*/