#include<bits/stdc++.h>
using namespace std;

typedef long long ll;
const int N=110,M=110*110,mod=1e9+7;
int f[4*N][M],ID[N][N],a[N],sum;

int mi(int x,int t){
    int d=1;
    while(t){
        if(t%2) d=(ll)d*x%mod;
        x=(ll)x*x%mod;t/=2;
    }
    return d;
}
int ni(int x) {return mi(x,mod-2);}

int dfs(int l,int r){
    if(ID[l][r]) return ID[l][r];
    int mx=(r-l+1)*(r-l+2)/2,now=++sum;
    for(int i=0;i<=mx;i++) f[now][i]=0;
    if(l>r){
        f[now][0]=1;
        return now;
    }
    bool pd=true;
    for(int i=l;i<=r;i++)
        if(a[i]==1){
            int lc=dfs(l,i-1),rc=dfs(i+1,r);
            int mx1=(i-l)*(i-l+1)/2,mx2=(r-i)*(r-i+1)/2;
            for(int j=0;j<=mx1;j++)
                for(int k=0;k<=mx2;k++)
                    f[now][j+k]=(f[now][j+k]+(ll)f[lc][j]*f[rc][k])%mod;
            pd=false;break;
        }
    if(!pd) return now;
    for(int i=l;i<=r;i++)
        if(a[i]==2){
            int son=dfs(i+1,r);
            int mx1=(i-l)*(i-l+1)/2,mx2=(r-i)*(r-i+1)/2;
            for(int j=0;j<=mx2;j++)
                f[now][mx1+j]=(f[now][mx1+j]-f[son][j]+mod)%mod;
        }
    f[now][mx]=(f[now][mx]+1)%mod;
    return now;
}

int main()
{

    int T;scanf("%d",&T);

    while(T--){
        int n,m;sum=0;
        scanf("%d%d",&n,&m);
        for(int i=1;i<=n;i++) scanf("%d",&a[i]);
        for(int i=0;i<=n+1;i++)
            for(int j=0;j<=n+1;j++) ID[i][j]=0;
        int rt=dfs(1,n),mx=n*(n+1)/2,Ans=0;
        for(int i=1;i<=mx;i++) Ans=(Ans+(ll)f[rt][i]*mi(i,m))%mod;
        printf("%d\n",Ans);
    }

    return 0;
}