#include <bits/stdc++.h>
using namespace std;
#define ll long long
ll mod;
inline void md(ll &x){if(x>=mod) x-=mod;}
inline ll _md(ll x){return x>=mod?x-mod:x;}
inline ll qpow(ll x,ll y)
{
    ll res=1;while(y){if(y&1ll) res=res*x%mod; x=x*x%mod,y>>=1;} return res;
}
struct Com
{
    ll x,y;Com(ll _x=0,ll _y=0){x=_x,y=_y;}
    inline Com operator+(const Com &z)const{return Com(_md(x+z.x),_md(y+z.y));}
    inline Com operator-(const Com &z)const{return Com(_md(x+mod-z.x),_md(y+mod-z.y));}
    inline Com operator*(const Com &z)const{return Com((x*z.x+mod-y*z.y%mod)%mod,(x*z.y+y*z.x)%mod);}
    inline Com operator*(ll k)const{return Com(x*k%mod,y*k%mod);}
    inline Com conj()const{return Com(x,_md(mod-y));}
    inline Com operator/(const Com &z)const{return (*this)*z.conj()*qpow((z.x*z.x+z.y*z.y)%mod,mod-2);}
}rt[4];
inline Com qpow(Com x,ll y)
{
    Com res(1);while(y){if(y&1ll) res=res*x; x=x*x,y>>=1;} return res;
}
int ct[4],n,a[4][4],s[4];
inline void calc(Com &res,ll K)
{
    Com ans(1);
    for(int i=0;i<4;i++)
    {
        int sum=0;
        for(int j=0;j<4;j++)
        {
            if(a[j][i]<0) return;
            sum+=a[j][i];
        }   
        if(sum!=ct[i]) return;
    }
    for(int i=0;i<4;i++)
    {
        s[i]=0;int nws=0;
        for(int j=0;j<4;j++)
            nws+=a[i][j]*j,s[i]+=a[i][j];
        ans=ans/qpow(rt[i],nws);
        ans=ans*qpow(Com(1)+rt[i]*rt[i],1ll*s[i]*(s[i]-1)/2);
    }
    for(int i=0;i<4;i++)
        for(int j=i+1;j<4;j++)
            ans=ans*qpow(Com(1)+rt[i]*rt[j],1ll*s[i]*s[j]);
    res=res+ans*K;
}
inline void work()
{
    cin>>n>>ct[0]>>ct[1]>>ct[2]>>ct[3];
    Com ans(0,0);
    if(n<=2)
    {
        for(int i=0;i<2;i++)
            for(int j=0;j<2;j++)
                for(int p1=0;p1<4;p1++)
                    for(int p2=0;p2<4;p2++)
                    {
                        if(!i&&p1) continue;
                        if(!j&&p2) continue;
                        memset(a,0,sizeof(a));
                        ll K=1;
                        if(i&&j) K=(p1==p2?1ll*ct[p1]*(ct[p1]-1)%mod:1ll*ct[p1]*ct[p2]%mod);
                        else if(i) K=ct[p1];
                        else if(j) K=ct[p2];
                        calc(ans,K);
                    }
    }
    for(int i=0;i<2;i++)
        for(int j=0;j<2;j++)
            for(int k=0;k<2;k++)
                for(int p1=0;p1<4;p1++)
                    for(int p2=0;p2<4;p2++)
                    {
                        if(!j&&p1) continue;
                        if(!k&&p2) continue;
                        memset(a,0,sizeof(a));
                        for(int p=0;p<4;p++)
                            a[i?0:2][p]=ct[p];
                        if(j) a[1][p1]++,a[i?0:2][p1]--;
                        if(k) a[3][p2]++,a[i?0:2][p2]--;
                        ll K=1;
                        if(j&&k) K=(p1==p2?1ll*ct[p1]*(ct[p1]-1)%mod:1ll*ct[p1]*ct[p2]%mod);
                        else if(j) K=ct[p1];
                        else if(k) K=ct[p2];
                        calc(ans,K);
                    }
    ans=ans*qpow(qpow(2,2ll*n+1ll*n*(n-1)/2ll),mod-2);
    cout<<ans.x<<"\n";
}
int main()
{
    ios::sync_with_stdio(false);
    cin.tie(0),cout.tie(0);
    int T;cin>>T>>mod;
    rt[0]=Com(1,0),rt[1]=Com(0,1);
    rt[2]=Com(mod-1,0),rt[3]=Com(0,mod-1);
    while(T--) work();


    return 0;
}