// 
#include<bits/stdc++.h>
#define ll long long
#define fi first
#define se second
#define pb push_back
#define db double
#define inf 0x3f3f3f3f3f3f3f3f
using namespace std;
const int mod=1e9+7;
const int N=5005;
int n,sz[N];
ll f[N][N],fac[N],inv[N];
vector<int>G[N];
ll fpow(ll x,ll y=mod-2){
    ll z(1);
    for(;y;y>>=1){
        if(y&1)z=z*x%mod;
        x=x*x%mod;
    }return z;
}
void init(int n){
    fac[0]=1;for(int i=1;i<=n;i++)fac[i]=fac[i-1]*i%mod;
    inv[n]=fpow(fac[n]);for(int i=n;i>=1;i--)inv[i-1]=inv[i]*i%mod;
}
ll binom(int x,int y){
    if(x<0||y<0||x<y)return 0;
    return fac[x]*inv[y]%mod*inv[x-y]%mod;
}
void dfs(int u,int topf){
    f[u][1]=1,sz[u]=1;
    for(auto v:G[u]){
        if(v==topf)continue;
        dfs(v,u);
        for(int i=sz[u];i>=1;i--){
            for(int j=sz[v];j>=1;j--){
                f[u][i+j]=(f[u][i+j]+f[u][i]*f[v][j]%mod*binom(i+j-1,j))%mod;
            }
            f[u][i]=f[u][i]*f[v][0]%mod;
        }
        sz[u]+=sz[v];
    }
    for(int i=sz[u]-1;i>=0;i--)f[u][i]=(f[u][i]+f[u][i+1])%mod;
}
int main(){
    //freopen("data.in","r",stdin);
    ios::sync_with_stdio(false);
    cin.tie(0),cout.tie(0);
    cin>>n,init(n);
    for(int i=1;i<n;i++){
        int x,y;cin>>x>>y;
        G[x].pb(y),G[y].pb(x);
    }
    dfs(1,0);
    cout<<(f[1][0]+mod)%mod;
}