#include <stdio.h>
#include <string.h>
#include <assert.h>

#include <vector>
#include <algorithm>

#define eb emplace_back

#define LOG(FMT...) fprintf(stderr, FMT)

using namespace std;

typedef long long LL;
constexpr int N = 105, mod = (int)1e9 + 7;

constexpr int reduce(int x) {
	return x >= mod ? x - mod : x;
}

constexpr int neg(int x) {
	return x ? mod - x : 0;
}

void Add(int &x, int y) {
	if ((x += y) >= mod) {
		x -= mod;
	}
}

constexpr int inv(int x, int k = mod - 2) {
	int r = 1;
	while (k) {
		if (k & 1) {
			r = x * (LL)r % mod;
		}
		x = x * (LL)x % mod;
		k >>= 1;
	}
	return r;
}

struct Pair {
	int x, y;
	Pair() {}
	Pair(int _x, int _y) : x(_x), y(_y) {}
	bool operator < (const Pair &t) const {
		return y < t.y;
	}
};
struct Poly {
	vector<Pair> v;
	Poly() : v{Pair(1, 0)} {}
	Poly(int k) : v{Pair(1, k * (k + 1) / 2)} {}
	Poly(vector<Pair> &t) : v(t) {}
	Pair &operator[] (int k) {
		return v[k];
	}
	Pair operator[] (int k) const {
		return v[k];
	}
	int size() const {
		return (int)v.size();
	}
	Poly operator * (const Poly &t) const {
		vector<Pair> ret;
		for (int i = 0; i < v.size(); ++i) {
			for (int j = 0; j < t.size(); ++j) {
				ret.eb(v[i].x * (LL)t[j].x % mod, reduce(v[i].y + t[j].y));
			}
		}
		return ret;
	}
	Poly operator - (const Poly &t) const {
		vector<Pair> ret(v.size() + t.size());
		for (int i = 0; i < v.size(); ++i) {
			ret[i] = v[i];
		}
		for (int i = 0; i < t.size(); ++i) {
			auto &[x, y]  = ret[i + v.size()];
			x = neg(t[i].x), y = t[i].y;
		}
		sort(ret.begin(), ret.end());
		int k = 0;
		for (int i = 0, j; (j = i) < ret.size(); i = j) {
			int coef = 0;
			while (j < ret.size() && ret[i].y == ret[j].y) {
				Add(coef, ret[i].x);
				++j;
			}
			ret[k++] = Pair(coef, ret[i].y);
		}
		ret.erase(ret.begin() + k, ret.end());
		return ret;
	}
	int val(int k) const {
		int ret = 0;
		for (const auto &[x, y] : v) {
			ret = (ret + x * (LL)inv(y, k)) % mod;
		}
		return ret;
	}
};

int a[N], n, m;
bool was[N][N];
Poly f[N][N];

Poly dp(int l, int r) {
	if (l == r + 1) {
		return Poly();
	}
	if (was[l][r]) {
		return f[l][r];
	}
	was[l][r] = true;
	Poly &ret = f[l][r];
	ret = Poly(r - l + 1);
	for (int k = l; k <= r; ++k) {
		if (a[k] == 1) {
			LOG("ERR\n");
			while (true);
		}
		if (a[k] == 2) {
			ret = ret - dp(l, k - 1) * Poly(r - k);
		}
	}
	return ret;
}

void solve() {
	scanf("%d%d", &n, &m);
	for (int i = 1; i <= n; ++i) {
		scanf("%d", a + i);
	}
	Poly ret;
	for (int i = 1, j; (j = i) <= n; i = j + 1) {
		while (j <= n && a[j] != 1) {
			j += 1;
		}
		if (a[i] != 1) {
			ret = ret * dp(i, j - 1);
		}
	}
	int res = ret.val(m);
	// LOG("%d\n", ret.size());
	// for (auto [x, y] : ret.v) LOG("%d %d\n", x, y);
	/* for (int i = 1; i <= m; ++i) {
		res = res * (LL)i % mod;
	} */
	printf("%d\n", res);
}

int main() {
	int task;
	for (scanf("%d", &task); task--; ) {
		solve();
	}
	return 0;
}