#include<bits/stdc++.h>
using namespace std;
#define FOR(i,s,t) for(int i=(s),_t=(t); i<=_t; ++i)
#define DOR(i,s,t) for(int i=(s),_t=(t); i>=_t; --i)
typedef long long ll;
const int N=2005;
const int Mod=1e9+7;
ll Fast(ll x,int b) {
    ll t=1;
    for(; b; b>>=1,x=x*x%Mod) if(b&1) t=t*x%Mod;
    return t;
}
ll f[N][N];
ll g[N];
int main() {
    ios::sync_with_stdio(false);
    cin.tie(nullptr);
    int n;
    cin >> n;
    ll inv2=Fast(2,Mod-2);
    f[1][1]=1;
    FOR(i,2,n) FOR(j,1,i) {
        f[i][j]=f[i-1][j-1]*(j*(j-1)/2)%Mod*inv2+f[i-1][j]*((i-j)*(i-j-1)/2)+f[i-1][j]*(j*(i-j))%Mod*inv2;
        f[i][j]%=Mod;
        f[i][j]=f[i][j]*Fast(i*(i-1)/2,Mod-2)%Mod;
    }
    g[0]=1;
    FOR(i,1,n) {
        int x;
        cin >> x;
        DOR(j,i-1,0) {
            g[j+1]=(g[j+1]+g[j]*x)%Mod;
        }
    }
    ll res=0;
    FOR(i,1,n) {
        res=(res+f[n][i]*g[i])%Mod;
        // cerr << i << ' ' << f[n][i] << endl;
    }
    cout << res << '\n';
    return 0;
}