#include<bits/stdc++.h>
using namespace std;
constexpr double sb[5][5]={{0.5000000000,0.6000000000,0.6428571429,0.6666666667,0.6818181818,},{1.0000000000,1.2000000000,1.2857142857,1.3333333333,1.3636363636,},{1.5000000000,1.8000000000,1.9285714286,2.0000000000,2.0454545455,},{2.0000000000,2.4000000000,2.5714285714,2.6666666667,2.7272727273,},{2.5000000000,3.0000000000,3.2142857143,3.3333333333,3.4090909091,},};
int main(){
    std::ios::sync_with_stdio(0);cin.tie(0);cout.tie(0);
    double ans=0;
    int n,cnt[5]={};cin>>n;
    for(int i=1,x;i<=n;i++)cin>>x,cnt[x-1]++;
    for(int i=0;i<5;i++) ans+=cnt[i]*3*(i+1)+(cnt[i])*(cnt[i]-1.0)*sb[i][i];
    for(int i=0;i<5;i++)for(int j=i+1;j<5;j++)ans+=cnt[i]*(sb[i][j]+sb[j][i])*cnt[j];
    cout<<fixed<<setprecision(10)<<ans<<'\n';
}