#include<bits/stdc++.h>
using namespace std;
const long long N=1e9+7;
int n,a[100010],cn[51];
long long ans[100010],t1,an;
void dfs(int p){
	if (p>=a[n]){
		int t=0;
		for (int i=1;i<=50;i++)t+=cn[i];
		an=(an+t1*ans[t]%N)%N;
		return ;
	}
	for (int i=1;i<=p;i++)
		if (cn[i]){
			long long tt=t1;t1=1ll*t1*cn[i]%N;cn[i]--;
			dfs(p+i);
			t1=tt;cn[i]++;
		}
}
int main(){
	cin>>n;
	for (int i=0;i<=n;i++)cin>>a[i];
//	memset(cn,0,sizeof(cn));
	sort(a+1,a+n+1);
	for (int i=1;i<=n;i++)cn[a[i]]++;
	ans[1]=1;
	for (int i=2;i<=n;i++){
		ans[i]=ans[i-1]*i%N;//cout<<ans[i]<<"\n";
	}
	t1=1;an=0;dfs(a[0]);
	for (int i=n-cn[0]+1;i<=n;i++)an=1ll*an*i%N;
	cout<<an;
	return 0;
}