#include <bits/stdc++.h>
using namespace std;
const long long MAXN = 100005;
long long hjj[4 * MAXN];
long long hj[4 * MAXN];
long long sum(long long *t, long long v, long long tl, long long tr,
              long long l, long long r) {
  if (l > r) return 0LL;
  if (l == tl && r == tr) return (long long)t[v];
  long long tm = (tl + tr) / 2;
  return sum(t, v * 2, tl, tm, l, min(r, tm)) +
         sum(t, v * 2 + 1, tm + 1, tr, max(l, tm + 1), r);
}
void upd(long long *t, long long v, long long tl, long long tr, long long pos,
         long long x) {
  if (tl == tr) {
    t[v] = x;
    return;
  }
  long long tm = (tl + tr) / 2;
  if (pos <= tm)
    upd(t, v * 2, tl, tm, pos, x);
  else
    upd(t, v * 2 + 1, tm + 1, tr, pos, x);
  t[v] = t[v * 2] + t[v * 2 + 1];
}
int main() {
  long long n;
  long long inv_total = 0;
  long long inv = 0;
  long long d = 0;
  cin >> n;
  for (long long i = 1; i <= n; i++) {
    long long p;
    cin >> p;
    p--;
    d += sum(hj, 1, 0, n - 1, p + 1, n - 1);
    inv_total += d;
    inv += sum(hjj, 1, 0, n - 1, p + 1, n - 1);
    upd(hj, 1, 0, n - 1, p, i);
    upd(hjj, 1, 0, n - 1, p, 1);
  }
  double randi = 0;
  double rand_total = 0;
  for (long long i = 2; i <= n; i++) {
    randi += (double)(i - 1) / 2;
    rand_total += randi * (double)(n - i + 1);
  }
  double ans = inv + (rand_total - (double)inv_total) /
                         (double)((long long)n * (long long)(n + 1) / 2);
  cout << setprecision(16) << fixed << ans;
  return 0;
}
