#include<bits/stdc++.h>
using namespace std;
const int maxn=4010;
const int mod=998244353;
int fir[maxn],to[2*maxn],Next[2*maxn],edge;
long long cnt[maxn],dp[maxn][maxn][2],mul[maxn];
void add(int x,int y){
	to[++edge]=y;Next[edge]=fir[x];fir[x]=edge;
}
void dfs(int now,int fa){
	dp[now][0][0]=cnt[now]=1;
	for(int i=fir[now];i;i=Next[i]){
		int z=to[i];
		if(z==fa) continue;
		dfs(z,now);
		for(int j=cnt[now]/2;j>=0;j--){
			for(int l=cnt[z]/2;l>=0;l--){
				if(l>0){
					(dp[now][j+l][0]+=dp[now][j][0]*(dp[z][l][0]+dp[z][l][1])%mod)%=mod;
					(dp[now][j+l][1]+=dp[now][j][1]*(dp[z][l][0]+dp[z][l][1])%mod)%=mod;
				}
				(dp[now][j+l+1][1]+=dp[now][j][0]*dp[z][l][0]%mod)%=mod;
			}
		}
		cnt[now]+=cnt[z];
	}
}
int main(){
	int n;
	cin>>n;
	for(int i=1;i<=2*n-1;i++){
		int u,v;cin>>u>>v;
		add(u,v);add(v,u);
	}
	mul[0]=1;
	dfs(1,0);
	for(int i=1;i<=n;i++) mul[i]=mul[i-1]*(2*i-1)%mod;
	long long ans=0;
	for(int i=0;i<=n;i++){
		if(i&1) (ans-=(dp[1][i][0]+dp[1][i][1])*mul[n-i]%mod)%=mod;
		else (ans+=(dp[1][i][0]+dp[1][i][1])*mul[n-i]%mod)%=mod;
	}
	cout<<ans<<endl;
	return 0;
}