#include<bits/stdc++.h>
#define ll long long
using namespace std;
const int N=2e5+10;
int n,m;
int vis[N],imp[N];
int dep[N],num[N],id;
struct Ready{
	struct tree{
		int v,nxt;
	}e[N];
	int hd[N],tot;
	void add(int u,int v){e[++tot]={v,hd[u]},hd[u]=tot;}
	int fa[N],top[N],son[N],sz[N];
	void clr(){for(int i=1;i<=n;i++)hd[i]=0;tot=id=0;}
	void dfs(int u){
		sz[u]=1,son[u]=0;
		for(int i=hd[u];i;i=e[i].nxt){
			int v=e[i].v;
			if(v==fa[u])continue;
			fa[v]=u,dep[v]=dep[u]+1;
			dfs(v);
			sz[u]+=sz[v];
			if(sz[son[u]]<sz[v])son[u]=v;
		}
	}
	void dfs2(int u,int h){
		num[u]=++id,top[u]=h;
		if(son[u])dfs2(son[u],h);
		for(int i=hd[u];i;i=e[i].nxt){
			int v=e[i].v;
			if(v==fa[u]||v==son[u])continue;
			dfs2(v,v);
		}
	}
	int lca(int u,int v){
		while(top[u]!=top[v]){
			if(num[u]<num[v])swap(u,v);
			u=fa[top[u]];
		}
		return num[u]<num[v]?u:v;
	}
}tr;
int x[N];
bool cmp(int a,int b){return num[a]<num[b];}
int st[N],tp;
int Fa[N],isa[N],isb[N],isc[N];
int stk[N];
int cnt;
int nd[N],Cnt;
void add(int u,int v,int k){
	if(vis[v]!=k)vis[v]=k,nd[++Cnt]=v;
	if(vis[u]!=k)vis[u]=k,nd[++Cnt]=u;
	if(isa[v]==k)isa[u]=k;
	if(isb[v]==k)isb[u]=k;
	Fa[v]=u;
	stk[++cnt]=v;
}
int T;
signed main(){
	scanf("%d",&T);
	for(;T--;){
		scanf("%d%d",&n,&m);
		tr.clr();
		for(int i=1;i<=n;i++)imp[i]=vis[i]=isa[i]=isb[i]=isc[i]=Fa[i]=stk[i]=st[i]=dep[i]=num[i]=0;
		for(int fa,i=2;i<=n;i++){
			scanf("%d",&fa);
			tr.add(fa,i);
		}
		dep[1]=1;
		tr.dfs(1),tr.dfs2(1,1);
		for(int l=1,k1,k2,k3;l<=m;l++){
			scanf("%d%d%d",&k1,&k2,&k3);
			int k=k1+k2+k3;
			for(int i=1,t;i<=k1;i++)scanf("%d",&t),x[i]=t,imp[t]=isa[t]=l;
			for(int i=k1+1,t;i<=k1+k2;i++)scanf("%d",&t),x[i]=t,imp[t]=isb[t]=l;
			for(int i=k1+k2+1,t;i<=k1+k2+k3;i++)scanf("%d",&t),x[i]=t,imp[t]=isc[t]=l;
			sort(x+1,x+k+1,cmp);
			k=unique(x+1,x+k+1)-x-1;
			st[tp=1]=1,cnt=0,Cnt=0;
			for(int i=1;i<=k;i++){
				if(x[i]==1)continue;
				int y=tr.lca(st[tp],x[i]);
				if(y^st[tp]){
					while(tp>1&&num[y]<num[st[tp-1]])add(st[tp-1],st[tp],l),--tp;
					add(y,st[tp],l);
					if(num[y]>num[st[tp-1]])st[tp]=y;
					else --tp;
				}
				st[++tp]=x[i];
			}
			for(int i=tp-1;i>=1;i--)add(st[i],st[i+1],l);
			Fa[1]=0;
			for(;cnt;--cnt){
				int v=stk[cnt],u=Fa[v];
				if(isc[u]==l)isc[v]=l;
			}
//			for(int i=1;i<=Cnt;i++)printf("Node:%d %d %d %d %d\n",x[i],isa[x[i]],isb[x[i]],isc[x[i]],Fa[x[i]]);
			for(int i=1;i<=Cnt;i++)if(isa[nd[i]]==l&&isb[nd[i]]==l&&isc[Fa[nd[i]]]==l)cnt+=dep[nd[i]]-dep[Fa[nd[i]]]-1;
			for(int i=1;i<=Cnt;i++)if(isa[nd[i]]==l&&isb[nd[i]]==l&&isc[nd[i]]==l)cnt++;
			printf("%d\n",cnt);
		}
	}
}