#include<bits/stdc++.h>
#define ff first
#define ss second
#define pb push_back
typedef long long ll;
using namespace std;
typedef pair<int,int> pii;

const int mod=998244353;
inline int add(int x,int y){int ret=x+y;if(ret>=mod)ret-=mod;return ret;}
inline int sub(int x,int y){int ret=x-y;if(ret<0)ret+=mod;return ret;}
inline int mul(int x,int y){return ((ll)x*y)%mod;}
inline int step(int base,int pw){int ret=1;while(pw){if(pw&1)ret=mul(ret,base);base=mul(base,base);pw>>=1;}return ret;}
inline int invv(int x){return step(x,mod-2);}


const int maxn=1e5+10;

vector<int>divs[maxn],primes[maxn];
int n,a[maxn],sum;
double p[maxn],e[maxn];

void prek_divs(){

    for(int i=1;i<maxn;i++){
        for(int j=1;i*j<maxn;j++)divs[i*j].pb(i);

        if(i!=1 && divs[i].size()==2){
            for(int j=1;i*j<maxn;j++)primes[i*j].pb(i);
        }

    }

}

int main(){

    ///freopen("test.txt","r",stdin);

    prek_divs();

    scanf("%d",&n);
    for(int i=1;i<=n;i++){
        scanf("%d",&a[i]);
        sum+=a[i];
    }
    for(int i=1;i<=n;i++){
        p[i]=(double)a[i]/sum;
    }

    for(int i=n;i>=1;i--){

        unordered_map<int,int>mapa;
        mapa.max_load_factor(0.25);
        for(int j=0;j<primes[i].size();j++){
            int pom=1;
            int p=primes[i][j];
            while(i%(pom*p)==0)pom*=p;
            mapa[p]=pom;
        }

        double q=0;
        double suma=0;
        for(int j=0;j<divs[i].size();j++){
            int id=divs[i][j];
            q+=p[id];
        }

        for(int j=2;i*j<=n;j++){

            int g=__gcd(i,j);
            int d=1;
            for(int k=0;k<primes[g].size();k++){
                int id=primes[g][k];
                d*=mapa[id];
            }

            int pomd=d*j;
            d=i/d;

            double ps=0;
            for(int k=0;k<divs[d].size();k++){
                int id=divs[d][k];
                ps+=p[id*pomd];
            }

            suma+=ps*(1+e[i*j]);
        }

        e[i]=(suma+q)/(1-q);
    }

    cout<<fixed<<setprecision(11)<<e[1]<<endl;

    return 0;
}