/* cerberus97 - Hanit Banga */

#include <iostream>
#include <iomanip>
#include <cassert>
#include <cmath>
#include <cstdio>
#include <cstring>
#include <cstdlib>
#include <map>
#include <set>
#include <queue>
#include <stack>
#include <vector>
#include <algorithm>
#include <numeric>

using namespace std;

#define pb push_back
#define fast_cin() ios_base::sync_with_stdio(false); cin.tie(NULL)

typedef long long ll;
typedef long double ld;
typedef pair <int, int> pii;
typedef pair <ll, ll> pll;

const int N = 110, X = 510, mod = 1e9 + 7;

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

int mul(ll x, int y) {
	return (x * y) % mod;
}

int l[N], dp[N][N][X];

int main() {
	fast_cin();
	int n, x;
	cin >> n >> x;
	for (int i = 1; i <= n; ++i) {
		cin >> l[i];
	}
	sort(l + 1, l + 1 + n, greater<int>());
	dp[0][0][0] = 1;
	for (int i = 0; i < n; ++i) {
		for (int k = 0; k * l[i] <= x and k <= i; ++k) {
			for (int sum = 0; sum <= x; ++sum) {
				int cur = dp[i][k][sum];
				if (!cur) {
					continue;
				}
				if (sum + l[i + 1] <= x) {
					add(dp[i + 1][k + 1][sum + l[i + 1]], cur);
				}
				add(dp[i + 1][k][sum], mul(cur, sum - k * (l[i + 1] - 1)));
				for (int y = 0; y < l[i + 1] - 1 and sum + y <= x; ++y) {
					add(dp[i + 1][k][sum + y + 1], mul(cur, 2 * k));
					add(dp[i + 1][k - 1][sum + y], mul(cur, mul(l[i + 1] - y - 1, k * (k - 1))));
				}
			}
		}
	}
	int ans = 0, fact = 1;
	for (int k = 1; k <= n; ++k) {
		fact = mul(fact, k);
		add(ans, mul(fact, dp[n][k][x]));
	}
	cout << ans << endl;
}