#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
#define REP(i, n)           for(int(i)=0;(i)<(n);++(i))
const int MOD = 1000000007;

int k[4];

double dfs(int N, int k0, int k1, int k2, int k3, int cnt){
    if(k0 == 0 && k1 == 0 && k2 == 0) return 0;
    if(cnt>100) return 0;

    double res = 1;
    if(k0 > 0) res += dfs(N, k0-1, k1+1, k2,   k3    , cnt+1) * k0 / N;
    if(k1 > 0) res += dfs(N, k0,   k1-1, k2+1, k3    , cnt+1) * k1 / N;
    if(k2 > 0) res += dfs(N, k0,   k1,   k2-1, k3 + 1, cnt+1) * k2 / N;
    if(k3 > 0) res += dfs(N, k0,   k1,   k2,   k3    , cnt+1) * k3 / N;
    return res;
}

int main(){
    int N; cin >> N;
    REP(i,N){
        int A; cin >> A;
        if(A < 3) k[A]++;
    }

    cout << dfs(N, k[0],k[1],k[2],N-k[0]-k[1]-k[2],0) << endl;
    return 0;
}