#include<bits/stdc++.h>
using namespace std;
int ct[4],mod;
int mul(int x,int y){return 1ll*x*y%mod;}
void add(int &x,const int y){if((x+=y)>=mod)x-=mod;}
void sub(int &x,const int y){if((x+=(mod-y))>=mod)x-=mod;}
int qpow(int a,int b){
	int c=1;
	while(b){
		if(b&1)c=mul(c,a);
		a=mul(a,a),b>>=1;
	}return c;
}
int inv(int a){return qpow(a,mod-2);}
struct Complex{int a,b;};
Complex operator*(Complex a,Complex b){
	Complex c;
	c.a=(1ll*a.a*b.a+mod-1ll*a.b*b.b%mod)%mod;
	c.b=(1ll*a.a*b.b+1ll*a.b*b.a)%mod;
	return c;
}
Complex operator+(Complex a,Complex b){
	add(a.a,b.a),add(a.b,b.b);
	return a;
}
Complex operator-(Complex a,Complex b){
	sub(a.a,b.a),sub(a.b,b.b);
	return a;
}
Complex qpow(Complex a,int b){
	Complex c;
	c.a=1,c.b=0;
	while(b){
		if(b&1)c=c*a;
		a=a*a;b>>=1;
	}return c;
}
Complex ans,t[3][4];
int vis[3],col[3];
void dfs(int o,int p,Complex prod){
	if(p==3){
		int nn=ct[0]+ct[1]+ct[2]+ct[3];
		if(!o){
			if(vis[1])prod=prod*qpow(Complex{inv(2),inv(2)},nn);
			if(vis[2])prod=prod*qpow(Complex{inv(2),mod-inv(2)},nn);
		}else{
			if((ct[1]+ct[3])&1)prod=prod*Complex{mod-1,0};
			if(vis[1])prod=prod*qpow(Complex{inv(2),mod-inv(2)},nn);
			if(vis[2])prod=prod*qpow(Complex{inv(2),inv(2)},nn);	
		}
		ans=ans+prod;
		return;
	}
	vis[p]=0,dfs(o,p+1,prod);
	vis[p]=1;
	for(int i=0;i<4;++i)
	  if(ct[i]){
	  	--ct[i];
	  	dfs(o,p+1,prod*Complex{ct[i]+1,0}*t[p][i]);
	  	++ct[i];
	  }
	vis[p]=0;
}
void solve(){
	int ni;cin>>ni;
	ans.a=ans.b=0;int n=0;
	for(int i=0;i<4;++i)scanf("%d",&ct[i]),n+=ct[i];
	vis[1]=vis[2]=0;
	dfs(0,1,Complex{1,0});
	dfs(1,1,Complex{1,0});
	Complex now;now.a=inv(4),now.b=0;
	while(n--)ans=ans*now;
    printf("%d\n",ans.a);
}
int main(){
	int T;cin>>T>>mod;
	t[1][0]=t[2][0]=Complex{1,0};
	t[1][2]=t[2][2]=Complex{mod-1,0};
	t[1][1]=t[2][3]=Complex{0,mod-1};
	t[1][3]=t[2][1]=Complex{0,1};
	while(T--)solve();
	return 0;
}