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

int n,x,a[3];
double dp[305][305][305];

int main()
{
	scanf("%d",&n);
	for(int i=1;i<=n;i++) scanf("%d",&x),a[x-1]++;
	for(int k=0;k<=a[2];k++)
	{
		for(int j=0;j<=a[2]+a[1]-k;j++)
		{
			for(int i=0;i<=n-j-k;i++)
			{
				if(i||j||k) dp[i][j][k]=1.0*n/(i+j+k);
				if(i) dp[i][j][k]+=1.0*i/(i+j+k)*dp[i-1][j][k];
				if(j) dp[i][j][k]+=1.0*j/(i+j+k)*dp[i+1][j-1][k];
				if(k) dp[i][j][k]+=1.0*k/(i+j+k)*dp[i][j+1][k-1];
			}
		}
	}
	printf("%.12lf\n",dp[a[0]][a[1]][a[2]]);
}