// 
#include <bits/stdc++.h>
#define int long long
#define ll long long
#define db long double
#define wdnmd const int mod=::mod;
#define eps 1e-12
#define lb(i) ((i)&(-(i)))
#define wr(x,ch) write(x),putchar(ch)
using namespace std;
typedef pair<int,int> pai;
#define x first
#define y second
#define gh() getchar()
inline long long read() {
	char ch=gh();
	long long x=0;
	char t=0;
	while(ch<'0'||ch>'9')   t|=ch=='-',ch=gh();
	while(ch>='0'&&ch<='9') x=x*10+(ch^48),ch=gh();
	return t?-x:x;
}
void write(ll x) {
    if(x<0) putchar('-'),x=-x;
    if(x>=10) write(x/10);
    putchar(x%10^48);
    return ;
}
const int NN=5e5+5,N=3005,M=55,mod=1e9+7;
int fac[NN],ifac[NN],pw[NN],ipw[NN];
int cnt[M],n,a[N],dp[M][N],ans;
int qpow(int x,int y=mod-2) {
	int ans=1;
	while(y) {
		if(y&1) ans=ans*x%mod;
		x=x*x%mod,y>>=1;
	}
	return ans;
}
int C(int x,int y) {
	if(x<y||y<0) return 0;
	return fac[x]*ifac[y]%mod*ifac[x-y]%mod;
}
void init(int n) {
	fac[0]=pw[0]=1;
	for(int i=1;i<=n;i++) fac[i]=fac[i-1]*i%mod,pw[i]=pw[i-1]*2%mod;
	ifac[n]=qpow(fac[n]),ipw[n]=qpow(pw[n]);
	for(int i=n-1;i>=0;i--) ifac[i]=ifac[i+1]*(i+1)%mod,ipw[i]=ipw[i+1]*2%mod;
	return ;
}
signed main() {
	init(5e5);
	n=read();
	for(int i=1;i<=n;i++) {
		a[i]=read();
		for(int j=0;j<a[i];j++) cnt[j]++;
	}
	for(int i=1;i<=n;i++) 
		for(int j=0;j<=i;j++) 
			dp[0][i]=(dp[0][i]+(j&1?mod-1:1)*qpow(i-j,n)%mod*ifac[n]%mod*C(n,j)%mod)%mod;
	for(int i=n+1;i<=3000;i++) dp[0][i]=1;
//	for(int i=1;i<=n;i++) wr(dp[0][i],' ');puts("");
	for(int i=1;i<=50;i++) 
		for(int j=1;j<=3000;j++) 
			for(int k=0;k<=cnt[i-1]&&k<=j*2;k++) 
				dp[i][j]=(dp[i][j]+(j*2-k>=3000?1:dp[i-1][j*2-k])*ipw[cnt[i-1]]%mod*C(cnt[i-1],k)%mod)%mod;
//    for(int i=1;i<=50;i++){
//        for(int j=1;j<=3000;j++){
//            for(int u=0;u<=cnt[i-1];u++){
//                int k1=ipw[cnt[i-1]]*C(cnt[i-1],u)%mod;
//                if(j*2-u<=3000&&j*2-u>=0)k1=k1*dp[i-1][j*2-u]%mod;
//                if(j*2-u<0)k1=0;
//                dp[i][j]=(dp[i][j]+k1)%mod;
//            }
//        }
//    }
	ans=1;
	for(int i=1;i<=n;i++) {
		int tmp=dp[a[i]][1];
		for(int j=a[i];j<=50;j++) tmp=tmp*ipw[cnt[j]]%mod;
		ans=(ans-tmp+mod)%mod;
	}
	wr(ans,'\n');
    return 0;
}