#include<bits/stdc++.h>
using namespace std;
#define ll long long
#define MP make_pair
mt19937 rnd(time(0));
const int MAXN=5005;
const int MOD=998244353;
bool Mbe;
int m,n,siz[MAXN];
vector<int> edg[MAXN];
ll f[MAXN][MAXN],g[MAXN][MAXN],h[MAXN][MAXN],tf[MAXN],tg[MAXN],th[MAXN];
ll p[MAXN],fac[MAXN],s;
// 向父亲延申（不算自己），不延展，延展但孤立点
ll q[MAXN][MAXN];// 把 i 个有区别元素分为 j 个有区别链的方案
bool Med;
void add(ll &x,ll y){x=(x+y)%MOD;}
void dfs(int u,int fa){
	siz[u]=1;
	h[u][0]=1;
	for(int v:edg[u]) if(v!=fa){
		dfs(v,u);
		memset(tf,0,sizeof(tf));memset(tg,0,sizeof(tg));memset(th,0,sizeof(th));
		for(int i=0;i<=siz[u];i++) for(int j=0;j<=siz[v];j++){
			add(tf[i+j],f[u][i]*g[v][j]%MOD+h[u][i]*f[v][j]%MOD);
			add(th[i+j],h[u][i]*g[v][j]%MOD);
			add(tg[i+j],g[u][i]*g[v][j]%MOD);
			add(tg[i+j+1],2*f[u][i]*f[v][j]%MOD);
		}
		swap(f[u],tf);swap(g[u],tg);swap(h[u],th);
		siz[u]+=siz[v];
	}
	for(int i=0;i<=siz[u];i++){
		add(g[u][i+1],2*f[u][i]+h[u][i]);
		add(f[u][i],h[u][i]);
	}
}
int main(){
	ios::sync_with_stdio(false);
	fac[0]=p[0]=1;
	for(int i=1;i<MAXN;i++) fac[i]=fac[i-1]*i%MOD;
	cin>>m;
	q[0][0]=1;
	for(int i=1;i<MAXN;i++) for(int j=1;j<=i;j++){
		q[i][j]=((i+j-1)*q[i-1][j]+q[i-1][j-1])%MOD;
	}
	p[0]=1;
	for(int i=1;i<=m;i++){
		cin>>n;
		for(int i=1;i<=n;i++) edg[i].clear();
		for(int i=1;i<=n;i++){
			memset(f[i],0,sizeof(f[i]));
			memset(g[i],0,sizeof(g[i]));
			memset(h[i],0,sizeof(h[i]));
		}
		for(int j=1;j<n;j++){
			int u,v;cin>>u>>v;
			edg[u].push_back(v);
			edg[v].push_back(u);
		}
		dfs(1,0);
		for(int i=1;i<=siz[1];i++) for(int j=1;j<i;j++) 
			add(g[1][j],((i-j)&1?MOD-1:1)*q[i][j]%MOD*g[1][i]%MOD);
		for(int i=s+n;i>=0;i--){
			p[i]=0;
			for(int j=1;j<=n&&j<=i;j++) add(p[i],g[1][j]*p[i-j]%MOD);
		}
		s+=n;
	}
	ll ans=0;
	for(int i=1;i<=s;i++) ans+=fac[i-1]*p[i]%MOD;
	cout<<ans%MOD<<'\n';
	return 0;
}