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

const int maxn=1e5+10;
const int mod=1e9+7;

int a[60];
int jc[maxn];
int ans;

void dfs(int rem,int n,int tmp)
{
	if(rem>=50||!n)
	{
		ans=(ans+1ll*tmp*jc[n])%mod;
		return;
	}
	for(int i=1;i<=rem;i++)
		if(a[i])
		{
			a[i]--;
			dfs(rem+i,n-1,1ll*tmp*(a[i]+1)%mod);
			a[i]++;
		}
}

int main()
{
	int n,rem;
	scanf("%d%d",&n,&rem);
	jc[0]=1;
	for(int i=1;i<=n;i++) jc[i]=1ll*jc[i-1]*i%mod;
	for(int i=1;i<=n;i++)
	{
		int x;
		scanf("%d",&x);
		a[x]++;
	}
	dfs(rem,n-a[0],1);
	for(int i=1;i<=a[0];i++) ans=1ll*ans*(n-i+1)%mod;
	printf("%d\n",ans);
	return 0;
}