#include <bits/stdc++.h>
using namespace std;
#define For(Ti,Ta,Tb) for(auto Ti=(Ta);Ti<=(Tb);++Ti)
#define Dec(Ti,Ta,Tb) for(auto Ti=(Ta);Ti>=(Tb);--Ti)
#define debug(...) fprintf(stderr,__VA_ARGS__)
#define range(Tx) begin(Tx),end(Tx)
using ll=long long;
const int N=805,Mod=1e9+7;
int n,cnt[N*2][N];
struct Edge{
	int u,v;
}edge[N*2];
vector<pair<int,int>> e[N];
ll fac[N],C[N][N],f[2][N*2][N],g[2][N*2][N];
ll binom(int n,int m){
	return n<m||m<0?0:C[n][m];
}
void dfs(int u,int fa,int d,int (&cnt)[N]){
	++cnt[d];
	for(auto [v,i]:e[u]){
		if(v!=fa){
			dfs(v,u,d+1,cnt);
		}
	}
}
void get_diam(int u,int fa,int d,int& x,int& xd){
	if(d>xd){
		x=u;
		xd=d;
	}
	for(auto [v,i]:e[u]){
		if(v!=fa){
			get_diam(v,u,d+1,x,xd);
		}
	}
}
ll combine(int i,int j){
	return i>0?fac[i+j]-j*fac[i+j-1]%Mod+Mod:0;
}
int main(){
	cin.tie(nullptr)->sync_with_stdio(false);
	// assert(freopen("tree.in","r",stdin));
	// assert(freopen("tree.out","w",stdout));
	cin>>n;
	if(n==1){
		cout<<"1\n";
		return 0;
	}
	For(i,1,n-1){
		int u,v;
		cin>>u>>v;
		edge[i]={u,v};
		edge[i+n-1]={v,u};
		e[u].emplace_back(v,i);
		e[v].emplace_back(u,i+n-1);
	}
	For(i,1,(n-1)*2){
		dfs(edge[i].v,edge[i].u,0,cnt[i]);
	}
	fac[0]=1;
	For(i,1,n){
		fac[i]=fac[i-1]*i%Mod;
	}
	For(i,0,n){
		For(j,0,i){
			C[i][j]=j?(C[i-1][j]+C[i-1][j-1])%Mod:1;
		}
	}
	int len=-1;
	{
		int x=0,xd=-1,y=0;
		get_diam(1,0,0,x,xd);
		get_diam(x,0,0,y,len);
	}
	int rt=0;
	if(len%2==0){
		For(i,1,(n-1)*2){
			int i_=(i>n-1?i-(n-1):i+(n-1));
			if(cnt[i][len/2]&&cnt[i_][len/2-1]){
				rt=i;
			}
		}
	}else{
		For(i,1,(n-1)*2){
			int i_=(i>n-1?i-(n-1):i+(n-1));
			if(cnt[i][len/2]&&cnt[i_][len/2]){
				rt=i;
			}
		}
	}
	For(j,1,(n-1)*2){
		f[0][j][0]=1;
		g[0][j][0]=2;
	}
	For(i,1,len/2){
		For(j,1,(n-1)*2){
			if(!cnt[j][i]) continue;
			auto [u,v]=edge[j];
			int j_=(j>n-1?j-(n-1):j+(n-1));
			fill(f[i&1][j],f[i&1][j]+cnt[j][i],0);
			for(auto [x,k]:e[v]){
				if(x==u||!cnt[k][i-1]) continue;
				int s=cnt[j_][i-1]+cnt[j][i]-cnt[k][i-1];
				For(l,0,cnt[j][i]-1){
					For(l_,0,min(l,cnt[k][i-1]-1)){
						For(y,0,cnt[k][i-1]-1-l_){
							(f[i&1][j][l]+=g[~i&1][k][l_+y]*binom(l_+y,l_)%Mod*binom(cnt[j][i]-cnt[k][i-1],l-l_)%Mod*combine(s-(l-l_),y))%=Mod;
						}
					}
				}
			}
			For(l,0,cnt[j][i]-1){
				For(l_,0,cnt[j_][i-1]-1){
					(f[i&1][j][l]+=g[~i&1][j_][l_]*binom(cnt[j][i],l)%Mod*combine(cnt[j][i]-l,l_))%=Mod;
				}
			}
		}
		For(j,1,(n-1)*2){
			if(!cnt[j][i]) continue;
			fill(g[i&1][j],g[i&1][j]+cnt[j][i],0);
			int j_=(j>n-1?j-(n-1):j+(n-1));
			For(l,0,cnt[j][i]-1){
				For(l_,0,cnt[j][i]-1-l){
					(g[i&1][j][l]+=f[i&1][j][l+l_]*binom(l+l_,l_)%Mod*combine(cnt[j_][i],l_))%=Mod;
				}
			}
			For(l,0,cnt[j][i]-1){
				For(l_,0,cnt[j_][i]-1){
					(g[i&1][j][l]+=f[i&1][j_][l_]*binom(cnt[j][i],l)%Mod*combine(cnt[j][i]-l,l_))%=Mod;
				}
			}
		}
	}
	if(len%2==0){
		cout<<f[(len/2)&1][rt][0]<<'\n';
	}else{
		cout<<g[(len/2)&1][rt][0]<<'\n';
	}
	return 0;
}