#include<bits/stdc++.h>
using namespace std;
using ll=long long;
const int N=1e3;
const int mod=1e9+7;
int add(int x,int y){ return (x+=y)>=mod?x-mod:x; }
int sub(int x,int y){ return (x-=y)<0?x+mod:x; }
int ksm(ll x,int tp,int s=1){
    for(;tp;x=x*x%mod,tp>>=1) if(tp&1) s=x*s%mod;
    return s;
}
int C[N+5][N+5];
int ikp[N+5];
int n;
int a[N+5];
int f[55][N+5];
int g[N*2+5];
int ans;
void get0(int *f){
    static int h[N+5];
    int iv=1;
    for(int i=1;i<=n;i++) iv=1ll*iv*i%mod;
    iv=ksm(iv,mod-2);
    for(int i=1;i<=n;i++)
        h[i]=ksm(i,n,iv);
    for(int i=1;i<=n;i++){
        for(int j=0;j<=i;j++)
            if(j&1) f[i]=(f[i]-1ll*C[n][j]*h[i-j])%mod;
            else f[i]=(f[i]+1ll*C[n][j]*h[i-j])%mod;
        if(f[i]<0)
            f[i]+=mod;
    }
    for(int i=0;i<n;i++)
        f[i]=sub(f[i+1],f[i]);
    f[n]=0;
}
int main(){
    scanf("%d",&n);
    for(int i=1;i<=n;i++)
        scanf("%d",&a[i]);
    ikp[0]=1;
    for(int i=1;i<=n;i++)
        ikp[i]=ksm(2,mod-2,ikp[i-1]);
    C[0][0]=1;
    for(int i=1;i<=n;i++){
        C[i][0]=C[i][i]=1;
        for(int j=1;j<i;j++)
            C[i][j]=add(C[i-1][j],C[i-1][j-1]);
    }
    get0(f[0]);
    for(int i=0;i<=50;i++){
        int m=0;
        for(int j=1;j<=n;j++){
            if(a[j]==i)
                ans=add(ans,f[i][0]);
            if(a[j]>i)
                m++;
        }
        ans=1ll*ans*ikp[m]%mod;
        memset(g,0,sizeof g);
        for(int j=0;j<=n;j++)
            for(int k=0;k<=m;k++)
                g[j+k]=(g[j+k]+1ll*f[i][j]*C[m][k])%mod;
        for(int j=0;j<n*2;j++)
            f[i+1][j/2]=(f[i+1][j/2]+1ll*g[j]*ikp[m])%mod;
        // printf("%d\n",ans);
    }
    printf("%d\n",sub(1,ans));
    return 0;
}