def solve(N, A):
    '''
    dp[i][j][k] 0枚のカードがi種類、1枚のカードがj種類、2枚のカードがk種類のときに、
    　　　　　　　　　　全ての種類のカードで3枚以上入手するまで引く回数の期待値とする。
    dp[i][j][k] =  1 + dp[i][j][k]*(N-i-j-k)/N
                     + dp[i][j][k-1]*k/N
                     + dp[i][j-1][k+1]*j/N
                     + dp[i-1][j+1][k]*i/N
    ->
    dp[i][j][k] = (N + dp[i][j][k-1]*k + dp[i][j-1][k+1]*j + dp[i-1][j+1][k]*i)/(i+j+k)
    '''
    global memo
    n0 = A.count(0)
    n1 = A.count(1)
    n2 = A.count(2)
    ns = n0 + n1 + n2 + 1
    memo = [[[-1]*ns for j in range(ns)] for k in range(ns)]
    memo[0][0][0] = 0
    return dp(n0, n1, n2, n0+n1+n2, N)

memo = []
def dp(n0, n1, n2, n_sum, N):
    global memo
    if n0 + n1 + n2 > n_sum or n0 < 0 or n1 < 0 or n2 < 0:
        return 0
    v = memo[n0][n1][n2]
    if v != -1:
        return v
    v = (N + dp(n0, n1, n2 - 1, n_sum, N) * n2
           + dp(n0, n1 - 1, n2 + 1, n_sum, N) * n1
           + dp(n0 - 1, n1 + 1, n2, n_sum, N) * n0)/(n0 + n1 + n2)
    memo[n0][n1][n2] = v
    return v
    
N = int(input())
As = list(map(int, input().split()))
print(solve(N, As))