#include<bits/stdc++.h>
#define rep(i, n) for(int i = 0; i < (n); ++i)
using namespace std;
typedef long long ll;
typedef long double ld;
typedef unsigned long long ull;
typedef vector<int> vi;
typedef pair<int, int> pii;
typedef pair<ll, int> pli;
typedef pair<ll, ll> pll;
typedef array<int, 3> ai3;
const int inf = 0x3f3f3f3f;
const ll INF = 0x3f3f3f3f3f3f3f3fll;
const int Mod = 1e9 + 7;
inline int sign(int a){ return (a&1) ? (Mod-1) : 1; }
inline void uadd(int &a, int b){ a += b-Mod; a += (a>>31) & Mod; }
inline void usub(int &a, int b){ a -= b, a += (a>>31) & Mod; }
inline void umul(int &a, int b){ a = (int)(1ll * a * b % Mod); }
inline int add(int a, int b){ a += b-Mod; a += (a>>31) & Mod; return a; }
inline int sub(int a, int b){ a -= b, a += (a>>31) & Mod; return a; }
inline int mul(int a, int b){ a = (int)(1ll * a * b % Mod); return a; }
int qpow(int b, ll p){ int ret = 1; while(p){ if(p&1) umul(ret, b); umul(b, b), p >>= 1; } return ret; }
const int fN = 10010;
int fact[fN], invfact[fN], pw2[fN], invpw2[fN], inv[fN];
void initfact(int n){
	pw2[0] = 1; for(int i = 1; i <= n; ++i) pw2[i] = mul(pw2[i-1], 2);
	invpw2[0] = 1; for(int i = 1; i <= n; ++i) invpw2[i] = mul(invpw2[i-1], (Mod+1) / 2);
	fact[0] = 1; for(int i = 1; i <= n; ++i) fact[i] = mul(fact[i-1], i);
	invfact[n] = qpow(fact[n], Mod-2); for(int i = n; i > 0; --i) invfact[i-1] = mul(invfact[i], i);
	for(int i = 1; i <= n; ++i) inv[i] = mul(invfact[i], fact[i-1]);
}
int binom(int n, int m){ return (m < 0 || m > n) ? 0 : mul(fact[n], mul(invfact[m], invfact[n-m])); }
int qbinom[2077][2077];
const double pi = acos(-1);
template<typename T> inline void chmax(T &_a, T _b){ (_a<_b) ? (_a=_b) : 0; }
template<typename T> inline void chmin(T &_a, T _b){ (_b<_a) ? (_a=_b) : 0; }
mt19937_64 rng(58);
inline int myrand(int l, int r){ return (int)(rng() % (r-l+1)) + l; }

int n, a[1010];
int cnt[55], suf[55], prepw[55];

int pr[1010];
int dp[55][2020];

int cdp = 0, cop = 0;

int main(){
	//freopen("j.in", "r", stdin);
	ios::sync_with_stdio(false);
	cin.tie(0);
	
	cin >> n;

	initfact(10000);
	rep(i, n + 2) rep(j, n + 2) qbinom[i][j] = binom(i, j);
	for(int i = 1; i <= n; ++i){
		for(int j = 0; j <= i; ++j){
			uadd(pr[i], mul(binom(n, j), mul(sign(j), mul(qpow(i - j, n), invfact[n]))));
		}
	}
	//for(int j = 1; j <= n; ++j) cout << pr[j] << " ";
	//cout << "\n";

	rep(i, n) cin >> a[i], ++cnt[a[i]];
	for(int i = 50; i >= 0; --i) suf[i] = cnt[i] + suf[i + 1];
	rep(i, 51) --suf[i];

	prepw[0] = 1;
	for(int i = 1; i <= 50; ++i) prepw[i] = mul(prepw[i-1], pw2[suf[i] + 1]);
	for(int i = 0; i <= 50; ++i) for(int d = 0; d <= 2*n; ++d){
		//cerr << i << " " << d << endl;
		++cdp;
		int &ret = dp[i][d];
		if(!i){
			ret = (d >= n ? 1 : pr[d + 1]);
		} else {
			ret = 0;
			for(int x = 0; x <= min(2*d + 1, suf[i]); ++x){
				++cop;
				int coef = qbinom[suf[i]][x];
				int sum = add((2*d - x >= 0) ? (2*d - x > 2*n ? prepw[i-1] : dp[i - 1][2*d - x]) : 0, (2*d + 1 - x >= 0) ? (2*d + 1 - x > 2*n ? prepw[i-1] : dp[i - 1][2*d + 1 - x]) : 0);
				uadd(ret, mul(sum, coef));
			}
		}
	}

	int ans = 0;
	rep(i, 51){
		int ret = dp[i][0];
		//cout << i << ": " << ret << endl;
		usub(ans, mul(ret, cnt[i]));
		//cout << ans << endl;
	}
	for(int i = 1; i <= 50; ++i)
		umul(ans, invpw2[suf[i] + 1]);
	uadd(ans, 1);

	//for(int i = 0; i <= 3; ++i){
	//	for(int j = 0; j <= 2*n; ++j)
	//		cout << i << " " << j << ": " << dp[i][j] << "\n";
	//}

	cout << ans << "\n";

	cerr << 1. * clock() / CLOCKS_PER_SEC << endl;
	cerr << cdp << " " << cop << endl;

	return 0;
}