#include<bits/stdc++.h>
using namespace std;
using ll = long long;

int n;
double dp[110][110][110];
int cnt[4];

double rec(int x,int y,int z)
{//x一枚以上あるカードの種類
//y二枚以上あるカードの種類
//z三枚以上あるカードの種類
//それからコンプリートするために必要な回数の期待値

	if(dp[x][y][z]!=-1)return dp[x][y][z];
	dp[x][y][z] = 0;
	//z/nででないとき、でるまでの回数の期待値
	double exp_num = 1.0 * n/(n-z);
	//n-x　ゼロ枚の数
	if(x!=n)dp[x][y][z] += (rec(x+1,y,z)+exp_num)*(n-x)/(n-z);
	//x-y一枚のすう
	if(y!=x)dp[x][y][z] += (rec(x,y+1,z)+exp_num)*(x-y)/(n-z);
	//y−z 二枚のすう
	if(z!=y)dp[x][y][z] += (rec(x,y,z+1)+exp_num)*(y-z)/(n-z);
	return dp[x][y][z];
}


int main()
{
	cin >> n;
	for(int i = 0;i<n+1;i++)for(int j = 0;j<n+1;++j)for(int k = 0;k <n+1;k++)dp[i][j][k] = -1;
	dp[n][n][n] = 0;
	for(int i = 0;i<n;i++)
	{
		int x;
		cin >> x;
		cnt[min(x,3)]++;
	}
	cnt[2] += cnt[3];
	cnt[1] += cnt[2];
	printf("%.12lf\n",rec(cnt[1],cnt[2],cnt[3]));
	return 0;
}