#include<cstdio>
const int N=303;
int b[4],n,a;double f[N][N][N];
signed main(){
	scanf("%d",&n);
	for(int i=1;i<=n;i++)scanf("%d",&a),b[a]++;
	f[0][0][0]=0;
	for(int k=0;k<=b[3];k++)
	  for(int j=0;j+k<=b[2]+b[3];j++)
		for(int i=0;i+j+k<=b[1]+b[2]+b[3];i++)if(i+j+k){
		  f[i][j][k]=1.0*n/(i+j+k);
		  if(i)f[i][j][k]+=1.0*i/(i+j+k)*f[i-1][j][k];
		  if(j)f[i][j][k]+=1.0*j/(i+j+k)*f[i+1][j-1][k];
		  if(k)f[i][j][k]+=1.0*k/(i+j+k)*f[i][j+1][k-1];
		}
	printf("%.15lf",f[b[1]][b[2]][b[3]]);
	return 0;
}