#include<bits/stdc++.h>
using namespace std;
typedef long long LL;
typedef pair<int, int> pi;
typedef double db;
#define V vector
#define pb push_back
#define fi first
#define se second
#define all1(x) (x).begin() + 1, (x).end()
const int INF = 2e9 + 10;

int main() {
    int n; cin >> n;
    V<int> w(n + 1);
    for(int i = 1; i <= n; i++) cin >> w[i];
    int sum = accumulate(all1(w), 0);
    V<V<pi>> dis(n + 1);

    for(int i = 1; i <= n; i++) {
        int now = i;
        for(int j = 2; j * j <= now; j++) {
            int pe = 1;
            while(now % j == 0) {
                pe *= j;
                now /= j;
            }
            if(pe != 1) dis[i].pb({j, pe});
        }
        if(now != 1) dis[i].pb({now, now});
    }

    map<pi, int> mp;
    V<int> tmp(n + 1);
    for(int x = 1; x <= n; x++) {
        for(int y = x; y <= n; y += x) {
            tmp[y] = 0;
        }
        for(int y = x; y <= n; y += x) {
            for(int z = y; z <= n; z += y) tmp[z] += w[y];
        }
        for(int y = x; y <= n; y += x) {
            mp[{x, y}] = tmp[y];
        }
    }

    auto get_lb = [&](int i, int j) -> int {
        int l = 0, r = 0;
        int ans = 1;
        while(r != dis[j].size()) {
            if(l == dis[i].size() || dis[i][l] != dis[j][r]) {
                if(l != dis[i].size() && dis[i][l].fi == dis[j][r].fi) l++;
                ans *= dis[j][r].se;
                r++;
            } else {
                l++, r++;
            }
        }
        return ans;
    };

    V<db> dp(n + 1);
    for(int i = n; i >= 1; i--) {
        db pself = (db)mp[{get_lb(i, i), i}] / sum;
        db tmp = 1;
        for(int j = 2 * i; j <= n; j += i) {
            tmp += (db)mp[{get_lb(i, j), j}] / sum * dp[j];
        }
        dp[i] = tmp * (1 / (1 - pself));
    }

    cout << fixed << setprecision(20) << dp[1] - 1 << endl;
}