#include<iostream>
#include<iomanip>
using namespace std;
double dp[101][101][101];
bool vis[101][101][101];
int N;
double f(int x,int y,int z)
{
	if(x<0||y<0||z<0||x+y+z>N)return 0;
	if(vis[x][y][z])return dp[x][y][z];
	if(!x&&!y&&!z)return 0;
	dp[x][y][z]=(x*f(x-1,y+1,z)+y*f(x,y-1,z+1)+z*f(x,y,z-1)+N)/(double)(x+y+z);
	vis[x][y][z]=true;
	return dp[x][y][z];
}
main()
{
	cin>>N;
	int I=0,J=0,K=0;
	for(int i=0;i<N;i++)
	{
		int A;cin>>A;
		if(A==0)I++;
		else if(A==1)J++;
		else if(A==2)K++;
	}
	cout<<fixed<<setprecision(16)<<f(I,J,K)<<endl;
}