#include<bits/stdc++.h>
#define ll long long
#define int long long
using namespace std;

const int mod = 1e9+7;
const int maxn = 105;
inline void add(int &a,const int &b){ a+=b;if(a>=mod)a-=mod; }
inline void dec(int &a,const int &b){ a-=b;if(a<0)a+=mod; }

int pw(int x,int k)
{
	int re=1;
	for(;k;k>>=1,x=(ll)x*x%mod) if(k&1)
		re=(ll)re*x%mod;
	return re;
}

int n,m;
int a[maxn];
vector<int>v1;
int f[maxn][maxn*maxn/2];
vector< vector<int> >G;
vector<int>Gcov;
int ansf[2][maxn*maxn/2];

signed main()
{
	ios_base::sync_with_stdio(false);
	//freopen("tmp.in","r",stdin);
	
	int Tcase; cin>>Tcase;
	while(Tcase--)
	{
		cin>>n>>m;
		
		{
			v1.clear();
			vector< vector<int> >_; G.swap(_);
			Gcov.clear();
		}
		v1.push_back(0);
		for(int i=1;i<=n;i++) 
		{
			cin>>a[i];
			if(a[i]==1) v1.push_back(i);
		}
		v1.push_back(n+1);
		
		for(int i=0,sz=(int)v1.size();i+1<sz;i++)
		{
			int x=v1[i],y=v1[i+1];
			int len= y-1-x;
			if(len>0)
			{
				int flag=0;
				for(int j=x+1;j<y;j++) if(a[j]==2) flag=1;
				
				int u= len*(len+1)/2;
		//		for(int j=x;j<=y;j++) for(int c=0;c<=u;c++)
		//			f[j][c]=0;
				f[x][0]=1;
				
				for(int l=x;l<y;l++) for(int c=0;c<=u;c++) if(f[l][c])
				{
					int &temp= f[l][c];
					for(int r=l+1;r<=y;r++) if(a[r]==2 || r==y)
					{
						int rl= r-l-1;
						int rc= rl*(rl+1)/2;
						if(c+rc<=u) dec(f[r][c+rc],temp);
					}
					temp=0;
				}
				
				vector<int>g(u+1);
				for(int c=0;c<=u;c++) if(f[y][c])
				{
					int &temp=f[y][c];
					dec(g[c],f[y][c]);
					temp=0;
				}
				G.emplace_back(g);
				Gcov.emplace_back(flag);
				
		/*		{
					cerr<<"[ "<<x<<" , "<<y<<" ]:";
					for(int c=0;c<=u;c++) cerr<<g[c]<<" \n"[c==u];
				}*/
			}
		}
		
		for(int c=0;c<=n*(n+1)/2;c++) ansf[0][c]=ansf[1][c]=0;
		int now=0;
		ansf[now][0]=1; int ansn=0;
		for(int gi=0,sz=(int)G.size();gi<sz;gi++)
		{
			auto &g=G[gi];
			int flag=Gcov[gi];
			
			int gn=g.size();
			now=!now;
			for(int ic=0;ic<=ansn;ic++) if(ansf[!now][ic])
			{
				int &temp=ansf[!now][ic];
		//		if(flag) dec(ansf[now][ic],temp);
				for(int jc=0;jc<gn;jc++) if(g[jc])
				{
					add(ansf[now][ic+jc],(ll)temp*g[jc]%mod);
				}
				temp=0;
			}
			ansn+= gn-1;
		}
		
/*		{
			cerr<<"ansf: ";
			for(int c=0;c<=n*(n+1)/2;c++)
				cerr<<ansf[now][c]<<" \n"[c==n*(n+1)/2];
		}*/
		
		int ans=0;
		for(int c=0;c<=n*(n+1)/2;c++)
			add(ans,(ll)ansf[now][c]*pw(c,m)%mod);
		cout<<ans<<'\n';
	}
	
	return 0;
}