// 
#include<bits/stdc++.h>
#define sgn(x) ((x)%2==1?(mod-1):1)
using namespace std;
const int N=3005,mod=1000000007,lim=1500,B=50;
inline void add(int&x,int y) {
	x=(x+y)%mod;
}
inline int qkpow(int x,int y) {
	int res=1;
	while (y) {
		if (y&1) res=1ll*res*x%mod;
		x=1ll*x*x%mod,y>>=1;
	}
	return res;
}
int n,ans=1,a[N],keg[N],f[B+5][N];
int fac[N],ifac[N],pw[N];
inline void init(int n) {
	fac[0]=1;
	for (int i=1;i<=n;i++) fac[i]=1ll*fac[i-1]*i%mod;
	ifac[n]=qkpow(fac[n],mod-2);
	for (int i=n-1;i>=0;i--) ifac[i]=1ll*ifac[i+1]*(i+1)%mod;
	pw[0]=1;
	int iv=qkpow(2,mod-2);
	for (int i=1;i<=n;i++) pw[i]=1ll*pw[i-1]*iv%mod;
}
inline int C(int n,int m) {
	return 1ll*fac[n]*ifac[m]%mod*ifac[n-m]%mod;
}
int main() {
	scanf("%d",&n);
	for (int i=1;i<=n;i++) {
		scanf("%d",&a[i]);
		for (int j=0;j<a[i];j++) keg[j]++;
	}
	init(n); 
	for (int i=0;i<=n;i++) {
		for (int j=0;j<=i;j++) {
			add(f[0][i],1ll*sgn(j)*C(n,j)%mod*qkpow(i-j,n)%mod);
		}
		f[0][i]=1ll*f[0][i]*ifac[n]%mod;
	}
	for (int i=n;i;i--) add(f[0][i],mod-f[0][i-1]);
	for (int i=0;i<=lim;i++) f[0][i]=f[0][i+1];
	for (int i=0;i<=B;i++) {
		for (int j=0;j<=lim;j++) {
			for (int k=0;k<=keg[i];k++) {
				add(f[i+1][(j+k)>>1],1ll*f[i][j]*C(keg[i],k)%mod*pw[keg[i]]%mod);
			}
		}
	}
	for (int i=1;i<=n;i++) {
		int res=f[a[i]][0];
		for (int j=a[i];j<=B;j++) {
			res=1ll*res*pw[keg[j]]%mod;
		}
		add(ans,mod-res);
	}
	printf("%d",ans);
	return 0;
}

/*

3
0 2 0

*/