#include <bits/stdc++.h>
using namespace std;

typedef long long ll;
typedef pair<int, int> pii;
typedef pair<ll, ll> pll;

const int MAXN = 1000;
const int MAXM = 50;

const int MOD = 1e9+7;
struct mint
{
    int x;
    mint() : x(0) {}
    mint(int _x) : x(_x) {}
    mint operator + (int ot) const { return x+ot >=MOD ? x+ot-MOD : x+ot; }
    mint operator - (int ot) const { return x<ot ? x+MOD-ot : x-ot; }
    mint operator - () const { return x ? MOD-x : 0; }
    mint operator * (int ot) const { return 1ll*ot*x%MOD; }
    mint operator += (int ot) { return *this = *this + ot; }
    mint operator -= (int ot) { return *this = *this - ot; }
    mint operator *= (int ot) { return *this = *this * ot; }
    operator int() const { return x; }
};

mint mpow(mint x, int a)
{
    mint ret=1;
    while(a)
    {
        if(a&1) ret*=x;
        x=x*x; a>>=1; 
    }
    return ret;
}
mint inv(mint x) { return mpow(x, MOD-2); }


int N, A[MAXM+10], B[MAXM+10];
mint fac[MAXN+10], ifac[MAXN+10], pows[MAXN+10], itp[MAXN+10];
mint dp[MAXM+10][MAXN+10];
mint ncr(int n, int r) { return fac[n]*ifac[r]*ifac[n-r]; }

int main()
{
    scanf("%d", &N);
    for(int i=1; i<=N; i++)
    {
        int t;
        scanf("%d", &t);
        A[t]++;
    }
    for(int i=MAXM; i>=0; i--) B[i]=B[i+1]+A[i+1];

    fac[0]=1;
    for(int i=1; i<=N; i++) fac[i]=fac[i-1]*i;
    ifac[N]=inv(fac[N]);
    for(int i=N; i>=1; i--) ifac[i-1]=ifac[i]*i;
    for(int i=1; i<=N; i++) pows[i]=mpow(i, N);
    itp[0]=1; itp[1]=inv(2);
    for(int i=2; i<=N; i++) itp[i]=itp[i-1]*itp[1];

    for(int k=0; k<=N; k++)
    {
        for(int i=0; i<=k; i++)
        {
            if(i%2==0) dp[0][k]+=pows[k-i]*ifac[i]*ifac[N-i];
            else dp[0][k]-=pows[k-i]*ifac[i]*ifac[N-i];
        }
    }
    for(int k=N; k>=1; k--) dp[0][k]-=dp[0][k-1];
    for(int k=1; k<=N; k++) dp[0][k-1]=dp[0][k];
    dp[0][N]=0;

    for(int i=1; i<=MAXM+1; i++)
    {
        for(int j=0; j<=N; j++)
        {
            for(int k=0; k<=B[i-1]; k++)
            {
                mint t;
                if(0<=2*j-k && 2*j-k<=N) t+=dp[i-1][2*j-k];
                if(0<=2*j+1-k && 2*j+1-k<=N) t+=dp[i-1][2*j+1-k];
                dp[i][j]+=t*ncr(B[i-1], k)*itp[B[i-1]];
            }
            // printf("%d %d : %d\n", i, j, dp[i][j].x);
        }
    }

    mint ans=1;
    for(int i=0; i<=MAXM; i++)
    {
        int sum=0;
        for(int j=i; j<=MAXM; j++) sum+=B[j];
        ans-=(mpow(itp[1], sum)*dp[i][0])*A[i];
    }
    printf("%d\n", ans.x);
}