#include <iostream>
#include <vector>
#include <algorithm>
#include <iomanip>
using namespace std;

#define RREP(i,s,e) for (i = s; i >= e; i--)
#define rrep(i,n) RREP(i,n,0)
#define REP(i,s,e) for (i = s; i < e; i++)
#define rep(i,n) REP(i,0,n)
#define INF 1e8

typedef long long ll;

double dp[101][101][101];
double n;

double e(int a, int b, int c) {
    if (dp[a][b][c] > 0)
        return dp[a][b][c];
    if (a+b+c > 0) {
        dp[a][b][c] = n;
        if (a > 0) dp[a][b][c] += e(a-1,b+1,c)*a;
        if (b > 0) dp[a][b][c] += e(a,b-1,c+1)*b;
        if (c > 0) dp[a][b][c] += e(a,b,c-1)*c;
        dp[a][b][c] /= a+b+c;
    }
    return dp[a][b][c];
}

int main() {
    int i, a, b, c;
    cin >> n;
    a = b = c = 0;
    rep (i,n) {
        int x;
        cin >> x;
        if (x == 0)
            a++;
        else if (x == 1)
            b++;
        else if (x == 2)
            c++;
    }
    cout << fixed << setprecision(7);
    cout << e(a,b,c) << endl;
    return 0;
}