# include <bits/stdc++.h>
using namespace std;
const int mod=1e9+7;
typedef long long ll;
typedef vector<int> Poly;
int n,sz[100],f[100][100];
ll fac[100],inv[100];
Poly P[100],P1[100];
vector<int> G[100];
template<typename T>
void upd(int &a,T b){a=(a+b)%mod;}
ll C(int n,int m){return fac[n]*inv[n-m]%mod*inv[m]%mod;}
Poly operator+(const Poly &a,const Poly &b)
{
	Poly ans(n+1);
	for(int i=0;i<=n;i++) ans[i]=(a[i]+b[i])%mod;
	return ans;
}
Poly operator*(Poly a,int b)
{
	for(int &i:a) i=(ll)i*b%mod;
	return a;
}
Poly operator*(const Poly &a,const Poly &b)
{
	Poly ans(n+1);
	for(int i=0;i<=n;i++) ans[i]=(ll)a[i]*b[i]%mod;
	return ans;
}
Poly Eval(const Poly &a)
{
	int sz=a.size();
	Poly ans(n+1);
	for(int i=0;i<=n;i++)
	{
		int v=0;
		for(int j=sz-1;j>=0;j--) v=((ll)v*i+a[j])%mod;
		ans[i]=v;
	}
	return ans;
}
Poly Lagrange(const Poly &a)
{
	Poly ans(n+1),mul(n+2);
	mul[0]=1;
	for(int i=0;i<=n;i++)
	{
		for(int j=i+1;j>=1;j--) mul[j]=((ll)(mod-i)*mul[j]+mul[j-1])%mod;
		mul[0]=(ll)(mod-i)*mul[0]%mod;
	}
	for(int i=0;i<=n;i++)
	{
		Poly F(n+1);
		for(int j=n,v=mul[j+1];j>=0;j--) F[j]=v,v=(mul[j]+(ll)v*i)%mod;
		ll w=((n-i)&1?mod-inv[n-i]:inv[n-i])*inv[i]%mod*a[i]%mod;
		for(int j=0;j<=n;j++) ans[j]=(ans[j]+w*F[j])%mod;
	}
	return ans;
}
int calc(Poly a)
{
	a=Lagrange(a);
	int ans=0;
	for(int i=0;i<=n;i++) upd(ans,a[i]*fac[i]);
	return ans;
}
void dfs(int u,int fa)
{
	sz[u]=1;
	for(int v:G[u])
		if(v!=fa) dfs(v,u);
	vector<int> vec;vec.push_back(u);
	for(int v:G[u])
		if(v!=fa) vec.push_back(v);
	int tot=vec.size();
	if(tot==1) return f[u][0]=0,f[u][1]=1,void();
	static Poly pre[100][100],suf[100][100],G[100];
	static int pres[100],sufs[100];
	pres[0]=sz[vec[0]];sufs[tot]=0;
	for(int i=0;i<=sz[vec[0]];i++) pre[0][i]=P[i]*C(sz[vec[0]],i);
	suf[tot][0]=Poly(n+1,1);
	for(int i=1;i<tot;i++)
	{
		pres[i]=pres[i-1]+sz[vec[i]];
		for(int j=0;j<=pres[i];j++) pre[i][j]=Poly(n+1);
		for(int j=0;j<=pres[i-1];j++)
			for(int k=0;k<=sz[vec[i]];k++)
				pre[i][j+k]=pre[i][j+k]+pre[i-1][j]*P[k]*C(sz[vec[i]],k);
	}
	for(int i=tot-1;i>=0;i--)
	{
		sufs[i]=sufs[i+1]+sz[vec[i]];
		for(int j=0;j<=sufs[i];j++) suf[i][j]=Poly(n+1);
		for(int j=0;j<=sufs[i+1];j++)
			for(int k=0;k<=sz[vec[i]];k++)
				suf[i][j+k]=suf[i][j+k]+suf[i+1][j]*P[k]*C(sz[vec[i]],k);
	}
	sz[u]=sufs[0];
	for(int i=1;i<=sz[u];i++) G[i]=Poly(n+1);
	for(int i=1;i<tot;i++)
	{
		int v=vec[i];
		static Poly F[100];
		for(int j=0;j<=sz[u]-sz[v];j++) F[j]=Poly(n+1);
		for(int j=0;j<=pres[i-1];j++)
			for(int k=0;k<=sufs[i+1];k++)
				F[j+k]=F[j+k]+pre[i-1][j]*suf[i+1][k];
		for(int j=0;j<=sz[u]-sz[v];j++)
			for(int k=1;k<=sz[v];k++)
				for(int l=1;l<=k;l++)
					G[j+k]=G[j+k]+F[j]*P1[k-l]*(C(sz[v]-l,k-l)*f[v][l]%mod);
	}
	for(int i=1;i<=sz[u];i++) G[i]=G[i]+suf[1][i-1];
	for(int i=1;i<=sz[u];i++) f[u][i]=calc(G[i]);
}
int main()
{
	int u,v;
	cin>>n;
	fac[0]=fac[1]=inv[0]=inv[1]=1;
	for(int i=2;i<=n;i++) fac[i]=fac[i-1]*i%mod;
	for(int i=2;i<=n;i++) inv[i]=(mod-mod/i)*inv[mod%i]%mod;
	for(int i=2;i<=n;i++) inv[i]=inv[i-1]*inv[i]%mod;
	for(int i=1;i<n;i++)
	{
		scanf("%d%d",&u,&v);
		G[u].push_back(v);G[v].push_back(u);
	}
	P1[0]=P[0]=Poly(n+1,1);
	for(int i=1;i<=n;i++)
	{
		P[i].resize(i+1);P1[i].resize(i+1);
		for(int j=1;j<=i;j++) P[i][j]=C(i-1,j-1)*fac[i]%mod*inv[j]%mod;
		for(int j=i-1;j>=0;j-=2) P[i][j]=mod-P[i][j];
		for(int j=0;j<i;j++) P1[i][j]=(P[i][j]+(ll)(mod-j-1)*P[i][j+1])%mod;
		P1[i][i]=P[i][i];
		P[i]=Eval(P[i]);P1[i]=Eval(P1[i]);
	}
	dfs(1,0);
	cout<<f[1][n]<<endl;
	return 0;
}