#include <bits/stdc++.h>

using namespace std;

typedef long long i64;
typedef unsigned long long u64;
typedef double lf;

const int mod = 1'000'000'000 + 7;

template<typename T>
T add(T x) {
    return x;
}
template<typename T, typename... Ts>
T add(T x, Ts... y) {
    T res = x + add(y...);
    if (res >= mod)
        res -= mod;
    return res;
}
template<typename T, typename... Ts>
T sub(T x, Ts... y) {
    return add(x, mod - add(y...));
}
template<typename T, typename... Ts>
void udd(T &x, Ts... y) {
    x = add(x, y...);
}
template<typename T>
T mul(T x) {
    return x;
}
template<typename T, typename... Ts>
T mul(T x, Ts... y) {
    return (x * 1ll * mul(y...)) % mod;
}
template<typename T, typename... Ts>
void uul(T &x, Ts... y) {
    x = mul(x, y...);
}

int bin(int a, i64 deg) {
    int r = 1;
    while (deg) {
        if (deg & 1)
            uul(r, a);
        deg >>= 1;
        uul(a, a);
    }
    return r;
}
int inv(int x) {
    assert(x);
    return bin(x, mod - 2);
}

const int maxb = 50;

void solve() {
    int n;
    cin >> n;
    vector<int> cnt(maxb + 1);
    vector<int> cnt_bit(maxb + 1);
    int maxv = 0;
    for (int i = 0; i < n; i++) {
        int x;
        cin >> x;
        cnt[x]++;
        maxv = max(maxv, x);
        for (int j = 0; j < x; ++j) {
            cnt_bit[j]++;
        }
    }

    int half = inv(2);
    
    vector<int> pows(maxb * maxb + 1);
    pows[0] = 1;
    for (int i = 1; i <= maxb * maxb; ++i) {
        pows[i] = mul(pows[i - 1], half);
    }
    
    vector< vector<int> > choose(maxb + 1, vector<int>(maxb + 1));
    choose[0][0] = 1;
    for (int i = 1; i <= maxb; ++i) {
        choose[i][0] = 1;
        for (int j = 1; j <= i; ++j) {
            choose[i][j] = add(choose[i - 1][j], choose[i - 1][j - 1]);
        }
    }
    vector<int> cdf(n + 1);
    int inv_fact = 1;
    for (int i = 0; i < n; ++i) {
        uul(inv_fact, inv(i + 1));
    }
    for (int i = 0; i <= n; ++i) {
        int sgn = 1;
        for (int k = 0; k <= i; ++k) {
            udd(cdf[i], mul(sgn, choose[n][k], bin(sub(i, k), n)));
            sgn = sub(0, sgn);
        }
        uul(cdf[i], inv_fact);
    }

    vector< vector<int> > dp(maxv + 1, vector<int>(2 * n + 1));
    for (int i = 0; i < n; ++i) {
        dp[0][i] = sub(cdf[i + 1], cdf[i]);
    }
    for (int i = 1; i <= maxv; ++i) {
        for (int j = 0; j <= 2 * n; ++j) {
            int &res = dp[i][j];
            for (int k = 0; k <= cnt_bit[i - 1]; ++k) {
                int v = 0;
                if (0 <= 2 * j - k && 2 * j - k <= 2 * n) {
                    udd(v, dp[i - 1][2 * j - k]);
                }
                if (0 <= 2 * j - k + 1 && 2 * j - k + 1 <= 2 * n) {
                    udd(v, dp[i - 1][2 * j - k + 1]);
                }
                uul(v, choose[cnt_bit[i - 1]][k], pows[cnt_bit[i - 1]]);
                udd(res, v);
            }
        }
    }
    vector<int> prob_ge(maxv + 1);
    for (int i = 0; i <= maxv; ++i) {
        int high = 0;
        for (int j = i; j <= maxv; ++j) {
            high += cnt_bit[j];
        }
        prob_ge[i] = sub(1, pows[high]);
        int v = 0;
        for (int j = 1; j <= 2 * n; ++j) {
            udd(v, dp[i][j]);
        }
        udd(prob_ge[i], mul(v, pows[high]));
    }
    int ans = 0;
    for (int i = 0; i <= maxv; ++i) {
        udd(ans, mul(cnt[i], sub(1, prob_ge[i])));
    }
    cout << sub(1, ans) << '\n';
}

signed main() {
#ifndef ONLINE_JUDGE
    freopen("input.txt", "r", stdin);
#endif
    ios_base::sync_with_stdio(false);
    int t = 1;
    // cin >> t;
    for (int i = 0; i < t; i++) {
        solve();
    }
}