#include<bits/stdc++.h>
using namespace std;
const int M=4005;
const int P=998244353,G=3,G_i=332748118;
int n,tot;
int hd[M],to[M<<1],nx[M<<1];
void add(int u,int v) {
    nx[++tot]=hd[u];
    to[tot]=v;
    hd[u]=tot;
}
int dp[2][M][M];
int Sz[M];
int A[M<<1],B[M<<1];
int C[M][M];
int F[M<<1];
int Mx=1,cnt,now,sw=1;
void Add(int &x,int y) {
    x+=y;(x>=P)&&(x-=P);
}
int fpow(int a,int cnt) {
	int res=1;
	for(int i=cnt; i; i>>=1,a=1LL*a*a%P)if(i&1)res=1LL*res*a%P;
	return res;
}
int pos[M<<1];
void NTT(int *a,int f) {
	for(int i=0; i<Mx; ++i)if(i<pos[i])swap(a[pos[i]],a[i]);
	for(int len=1; len<Mx; len<<=1) {
		now=fpow(f?G:G_i,(P-1)/(len<<1));
		for(int L=0; L<Mx; sw=1,L+=len<<1) {
			for(int R=L; R<L+len; ++R,sw=1LL*sw*now%P) {
				int xx=a[R],yy=1LL*a[R+len]*sw%P;
				a[R]=xx+yy;(a[R]>=P)&&(a[R]-=P);
				a[R+len]=xx-yy;(a[R+len]<0)&&(a[R+len]+=P);
			}
		}
	}
}
int Sum[2][M<<1];
void Solve(int now,int Fa) {
    dp[1][now][1]=1;Sz[now]=1;
    for(int i=hd[now]; i; i=nx[i]) {
        int To=to[i];
        if(To==Fa)continue;
        Solve(To,now);
//        cout<<"To:"<<To<<" "<<endl;
        int Minn=min(Sz[now],Sz[To]);
        memset(Sum,0,sizeof Sum);
        for(int k=0; k<=Minn; ++k) {
            Mx=1;cnt=0;
            while(Mx<=Sz[now]+Sz[To]-k*2)Mx<<=1,cnt++;
	        for(int j=0; j<Mx; ++j)pos[j]=(pos[j>>1]>>1)|((j&1)<<cnt-1);

            for(int j=0; j<Mx; ++j)A[j]=B[j]=0;
            for(int j=k; j<=Sz[now]; ++j)if(j)A[j-k]=1LL*dp[1][now][j]*C[j-1][k]%P;
            for(int j=k; j<=Sz[To]; ++j)B[j-k]=1LL*(dp[1][To][j]+dp[0][To][j])*C[j][k]%P*F[k]%P;
            NTT(A,1),NTT(B,1);
	        for(int j=0; j<Mx; ++j)A[j]=1LL*A[j]*B[j]%P;
	        NTT(A,0);
            int inv_Mx=fpow(Mx,P-2);
	        for(int j=0; j<=Sz[now]+Sz[To]-k*2; ++j)Sum[1][j]=(Sum[1][j]+1LL*A[j]*inv_Mx)%P;
	        
	        
            for(int j=0; j<Mx; ++j)A[j]=0;
            for(int j=k; j<=Sz[now]; ++j)A[j-k]=1LL*dp[0][now][j]*C[j][k]%P;
            NTT(A,1);
            for(int j=0; j<Mx; ++j)A[j]=1LL*A[j]*B[j]%P;
            NTT(A,0);
//            if(To==3&&now==1)cout<<k<<"A[0]:"<<A[1]<<" "<<Sum[0][1]<<" "<<Sum[0][1]<<endl;
            for(int j=0; j<=Sz[now]+Sz[To]-k*2; ++j)Sum[0][j]=(Sum[0][j]+1LL*A[j]*inv_Mx)%P;
//            if(To==3&&now==1)cout<<k<<"A[0]:"<<A[1]<<" "<<Sum[0][1]<<" "<<Sum[0][1]<<endl;
            
//            if(To==2&&k==1)cout<<Sum[0][0]<<" WWW"<<" "<<dp[0][1][1]<<" "<<dp[1][2][1]<<" "<<C[1][1]<<endl;
			for(int j=0; j<Mx; ++j)A[j]=B[j]=0;
            for(int j=k+1; j<=Sz[now]; ++j)A[j-k-1]=1LL*dp[1][now][j]*C[j-1][k]%P;
            for(int j=k+1; j<=Sz[To]; ++j) {
				B[j-k-1]=(1LL*dp[1][To][j]*(j-1)%P*C[j-1][k]%P*F[k]+1LL*dp[0][To][j]%P*C[j][k+1]%P*F[k+1])%P;
//				if(now==2&&k==0&&j==2)cout<<"test:"<<1LL*dp[0][To][j]*j%P*C[j][k+1]%P*F[k+1]<<" "<<dp[0][To][j]<<" "<<B[1]<<endl;
			}
//			if(now==2)cout<<"A[0]:"<<A[0]<<" B[1]"<<B[1]<<" k:"<<k<<endl;
            NTT(A,1),NTT(B,1);
	        for(int j=0; j<Mx; ++j)A[j]=1LL*A[j]*B[j]%P;
	        NTT(A,0);
//            if(To==3&&now==1)cout<<k<<"A[0]:"<<A[1]<<" "<<Sum[0][1]<<" "<<Sum[0][1]<<endl;
	        for(int j=0; j<Mx; ++j)Sum[0][j]=(Sum[0][j]+1LL*A[j]*inv_Mx)%P;
//			if(now==2)cout<<"k:"<<k<<" "<<Sum[0][1]<<endl;
//            if(To==3&&now==1)cout<<k<<"A[0]:"<<A[1]<<" "<<Sum[0][1]<<" "<<Sum[0][1]<<endl;

        }
        for(int j=0; j<=Sz[now]+Sz[To]; ++j) dp[0][now][j]=Sum[0][j],dp[1][now][j]=Sum[1][j];
//        if(now==1)cout<<dp[0][1][1]<<" ASDS"<<endl;
//        if(To==3)cout<<"SADS"<<dp[0][now][1]<<" "<<endl;
		Sz[now]+=Sz[To];
    }
//    cout<<now<<"  "<<endl;
//    for(int i=0; i<=Sz[now]; ++i)cout<<dp[1][now][i]<<" "<<" "<<dp[0][now][i]<<endl;
}
int main() {
//	freopen("a.in","r",stdin);
    scanf("%d",&n);
    for(int i=1; i<n+n; ++i) {
        int u,v;
        scanf("%d%d",&u,&v);
        add(u,v);add(v,u);
    }
    F[0]=1;
	for(int i=1; i<=n+n; ++i)F[i]=1LL*F[i-1]*i%P;
    for(int i=0; i<=n+n; ++i) {
        C[i][0]=1;
        for(int j=1; j<=i; ++j) {
            C[i][j]=C[i-1][j];
            Add(C[i][j],C[i-1][j-1]);
        }
    }
    Solve(1,0);
    printf("%d",dp[0][1][0]);
    return 0;
}