#include<bits/stdc++.h>
using namespace std;
#define N 2505
#define M 55
#define int long long 
const int p=1e9+7;
int power(int a,int b){
    int ans=1;
    while(b){
        if(b&1)ans=ans*a%p;
        a=a*a%p;b>>=1;
    }
    return ans;
}
int f[M][N<<1],n,m,a[N],b[N],jc[N],inv[N],pw[N];
int C(int n,int m){
    if(n<m)return 0;
    return jc[n]*inv[m]%p*inv[n-m]%p;
}
signed main(){
    ios::sync_with_stdio(false);cin.tie(0);cout.tie(0);
    jc[0]=1;
    for(int i=1;i<=1500;++i)jc[i]=jc[i-1]*i%p;
    inv[1500]=power(jc[1500],p-2);pw[0]=1;
    for(int i=1;i<=1500;++i)pw[i]=pw[i-1]*((p+1)>>1)%p;
    for(int i=1500;i;--i)inv[i-1]=inv[i]*i%p;
    cin>>n;
    for(int i=1;i<=n;++i)cin>>a[i],b[a[i]]++;
    for(int i=50;i>=0;--i)b[i]+=b[i+1];
    for(int i=0;i<n;++i){
        for(int j=0,x=1;j<=i;++j,x=-x){
            f[0][i]+=x*C(n,j)%p*power(i-j+1,n)%p;
        }
        // cout<<f[0][i]<<" ";
        f[0][i]=(f[0][i]*inv[n]%p+p)%p;
    }
    for(int i=n-1;i;--i)f[0][i]=f[0][i]-f[0][i-1];
    for(int i=1;i<=50;++i){
        for(int j=0;j<=2000;++j)if(f[i-1][j])for(int k=0;k<=b[i];++k){
            (f[i][j+k>>1]+=f[i-1][j]*C(b[i],k)%p*pw[b[i]]%p)%=p;
        }
    }
    int ans=1;
    for(int i=1;i<=n;++i){
        int w=f[a[i]][0]%p;
        for(int j=a[i]+1;j<=50;++j)w=w*pw[b[j]]%p;
        ans=ans+p-w;
    }
    ans=(ans%p+p)%p;
    cout<<ans<<"\n";
}