#include<bits/stdc++.h> using namespace std; const int N=5e3+10,mod=998244353,iv2=(mod+1)/2; void Add(int &a,int b){a+=b;if(a>=mod) a-=mod;}; int n; vector<int> E[N]; int f[N][N][3]; int tmp[N][3],res[N]; int fac[N],siz[N]; void dfs(int u,int fa){ f[u][1][0]=1; int cur=1; for(auto v:E[u]){ if(v==fa) continue; dfs(v,u); memset(tmp,0,sizeof(tmp)); for(int i=1;i<=cur;i++){ for(int j=1;j<=siz[v];j++){ int sumj=(1ll*f[v][j][0]+1ll*f[v][j][1]+1ll*f[v][j][2])%mod; Add(tmp[i+j][2],1ll*f[u][i][2]*sumj%mod); Add(tmp[i+j][1],1ll*f[u][i][1]*sumj%mod); Add(tmp[i+j][0],1ll*f[u][i][0]*sumj%mod); Add(tmp[i+j-1][2],1ll*f[u][i][1]*f[v][j][1]%mod*iv2%mod); Add(tmp[i+j-1][2],1ll*f[u][i][1]*f[v][j][0]%mod); Add(tmp[i+j-1][1],1ll*f[u][i][0]*f[v][j][1]%mod); Add(tmp[i+j-1][1],2ll*f[u][i][0]*f[v][j][0]%mod); } }cur+=siz[v]; for(int i=1;i<=cur;i++){ f[u][i][0]=tmp[i][0]; f[u][i][1]=tmp[i][1]; f[u][i][2]=tmp[i][2]; } }siz[u]=cur; } int main(){ cin>>n; for(int i=1;i<n;i++){ int u,v;cin>>u>>v; E[u].push_back(v); E[v].push_back(u); } fac[0]=1;for(int i=1;i<=n;i++) fac[i]=1ll*fac[i-1]*i%mod; dfs(1,0); for(int i=1;i<=n;i++){ Add(res[i],f[1][i][0]); Add(res[i],f[1][i][1]); Add(res[i],f[1][i][2]); }int ans=0; for(int i=0;i<n;i++){ int w=(i&1)?(mod-1):1; Add(ans,1ll*w*res[n-i]%mod*fac[n-i]%mod); }cout<<ans<<endl; return 0; }