#include <bits/stdc++.h> 
using namespace std;
typedef long long ll;
typedef pair<int,int> Pint;
typedef pair<ll, ll> P;
//typedef pair<int, pair<int, int>> P;
//typedef tuple<int,int,int> T;
ll INFL = 1000000000000000010;//10^18 = 2^60
int INF = 2147483600;//10^9
ll MOD  = 1000000007;
vector<int> dy = {0,0,1,-1};
vector<int> dx = {1,-1,0,0};

int N;
double memo[110][110][110];
int X, Y, Z;


double dfs(int x, int y, int z){
    if(memo[x][y][z] != -1) return memo[x][y][z];
    memo[x][y][z] = 0;
    double exp_num = 1.0 * N / (N - z);//自己ループする回数
    if(x != N) memo[x][y][z] += (dfs(x+1, y, z) + exp_num) * (N - x) / (N - z);//０枚
    if(y != x) memo[x][y][z] += (dfs(x, y+1, z) + exp_num) * (x - y) / (N - z);//１枚
    if(z != y) memo[x][y][z] += (dfs(x, y, z+1) + exp_num) * (y - z) / (N - z);//２枚
    return memo[x][y][z];
}

int main(void){
    cin >> N;
    for(int i = 0; i < N; i++){
        int A; cin >> A;
        if(A >= 1) X++;
        if(A >= 2) Y++;
        if(A >= 3) Z++;
    }
    //X　1枚以上持っている枚数
    //Y　2枚以上持っている枚数
    //Z　3枚以上持っている枚数
    for(int i = 0; i < 110; i++){
        for(int j = 0; j < 110; j++){
            for(int k = 0; k < 110; k++){
                memo[i][j][k] = -1;
            }
        }
    }
    cout << fixed << setprecision(10) << dfs(X, Y, Z) << endl;
}