Problem: You are given a permutation of n numbers p1, p2, ..., pn. We perform k operations of the following type: choose uniformly at random two indices l and r (l ≤ r) and reverse the order of the elements pl, pl + 1, ..., pr. Your task is to find the expected value of the number of inversions in the resulting permutation.

Input

The first line of input contains two integers n and k (1 ≤ n ≤ 100, 1 ≤ k ≤ 109). The next line contains n integers p1, p2, ..., pn — the given permutation. All pi are different and in range from 1 to n.

The problem consists of three subproblems. The subproblems have different constraints on the input. You will get some score for the correct submission of the subproblem. The description of the subproblems follows.

  * In subproblem G1 (3 points), the constraints 1 ≤ n ≤ 6, 1 ≤ k ≤ 4 will hold. 
  * In subproblem G2 (5 points), the constraints 1 ≤ n ≤ 30, 1 ≤ k ≤ 200 will hold. 
  * In subproblem G3 (16 points), the constraints 1 ≤ n ≤ 100, 1 ≤ k ≤ 109 will hold. 

Output

Output the answer with absolute or relative error no more than 1e - 9.

Examples

Input

3 1
1 2 3


Output

0.833333333333333


Input

3 4
1 3 2


Output

1.458333333333334

Note

Consider the first sample test. We will randomly pick an interval of the permutation (1, 2, 3) (which has no inversions) and reverse the order of its elements. With probability <image>, the interval will consist of a single element and the permutation will not be altered. With probability <image> we will inverse the first two elements' order and obtain the permutation (2, 1, 3) which has one inversion. With the same probability we might pick the interval consisting of the last two elements which will lead to the permutation (1, 3, 2) with one inversion. Finally, with probability <image> the randomly picked interval will contain all elements, leading to the permutation (3, 2, 1) with 3 inversions. Hence, the expected number of inversions is equal to <image>.

Solution (in CPP):

Solution:
#include <bits/stdc++.h>
const int INF = 1e9 + 7;
int ps[101];
double weight[101][101];
double px0[101][101];
double px1[101][101];
double py0[101][101];
double py1[101][101];
double pd0[101][101];
double pd1[101][101];
double tmp[101][101];
struct Segment {
  int x0, x1;
  int m, b;
};
const std::vector<Segment>& lower_envelope(int bn1, int b0, int b1, int x0,
                                           int x1) {
  static std::vector<Segment> res;
  res.clear();
  int t0 = b0 - b1, t1 = bn1 - b0;
  if (t0 < t1) {
    res.push_back({x0, t0 - 1, 1, b1});
    res.push_back({t0, t1 - 1, 0, b0});
    res.push_back({t1, x1, -1, bn1});
  } else {
    int t2 = (bn1 - b1 + 1) / 2;
    res.push_back({x0, t2 - 1, 1, b1});
    res.push_back({t2, x1, -1, bn1});
  }
  for (auto& it : res) {
    it.x0 = std::max(it.x0, x0);
    it.x1 = std::min(it.x1, x1);
  }
  return res;
}
int main() {
  int N, K;
  scanf("%d %d", &N, &K);
  for (int i = 0; i < N; i++) {
    scanf("%d", &ps[i]);
  }
  for (int i = 0; i < N; i++) {
    for (int j = 0; j < N; j++) {
      if (ps[i] < ps[j]) {
        weight[i][j] = 1;
      } else {
        weight[j][i] = 1;
      }
    }
  }
  long long start = clock();
  long long timer0 = 0, timer1 = 0, timer2 = 0;
  while (K-- && (clock() - start < 1.9 * CLOCKS_PER_SEC)) {
    timer0 -= clock();
    for (int x = 0; x < N; x++) {
      for (int y = 0; y < N; y++) {
        if (x == y) continue;
        auto transfer = [x, y, N](int x2, int y2, double v) {
          tmp[x2][y2] += weight[x][y] * v;
        };
        if (x < y) {
          transfer(
              x, y,
              ((x + 1) * x + (y - x) * (y - x - 1) + (N - y) * (N - y - 1)) /
                  2.0);
        } else {
          transfer(
              x, y,
              ((y + 1) * y + (x - y) * (x - y - 1) + (N - x) * (N - x - 1)) /
                  2.0);
        }
        if (x < y) {
          {
            auto& segs = lower_envelope(y, std::min(x + 1, y - x), 1, 0, y - 1);
            for (auto seg : segs) {
              if (seg.x0 > seg.x1) continue;
              px0[seg.x0][y] += seg.b * weight[x][y];
              px1[seg.x0][y] += seg.m * weight[x][y];
              px0[seg.x1 + 1][y] -= seg.b * weight[x][y];
              px1[seg.x1 + 1][y] -= seg.m * weight[x][y];
            }
          }
          {
            auto& segs =
                lower_envelope(N, std::min(y - x, N - y), -x, x + 1, N - 1);
            for (auto seg : segs) {
              if (seg.x0 > seg.x1) continue;
              py0[x][seg.x0] += seg.b * weight[x][y];
              py1[x][seg.x0] += seg.m * weight[x][y];
              py0[x][seg.x1 + 1] -= seg.b * weight[x][y];
              py1[x][seg.x1 + 1] -= seg.m * weight[x][y];
            }
          }
        } else {
          {
            auto& segs =
                lower_envelope(N, std::min(x - y, N - x), -y, y + 1, N - 1);
            for (auto seg : segs) {
              if (seg.x0 > seg.x1) continue;
              px0[seg.x0][y] += seg.b * weight[x][y];
              px1[seg.x0][y] += seg.m * weight[x][y];
              px0[seg.x1 + 1][y] -= seg.b * weight[x][y];
              px1[seg.x1 + 1][y] -= seg.m * weight[x][y];
            }
          }
          {
            auto& segs = lower_envelope(x, std::min(y + 1, x - y), 1, 0, x - 1);
            for (auto seg : segs) {
              if (seg.x0 > seg.x1) continue;
              py0[x][seg.x0] += seg.b * weight[x][y];
              py1[x][seg.x0] += seg.m * weight[x][y];
              py0[x][seg.x1 + 1] -= seg.b * weight[x][y];
              py1[x][seg.x1 + 1] -= seg.m * weight[x][y];
            }
          }
        }
        {
          auto& segs = lower_envelope(
              std::min(N, N - x + y), std::min({x + 1, y + 1, N - x, N - y}),
              std::min(1, 1 + x - y), std::max(0, y - x),
              std::min(N - 1, y - x + N - 1));
          for (auto seg : segs) {
            if (seg.x0 > seg.x1) continue;
            pd0[seg.x0][seg.x0 + x - y] += seg.b * weight[x][y];
            pd1[seg.x0][seg.x0 + x - y] += seg.m * weight[x][y];
            pd0[seg.x1 + 1][seg.x1 + 1 + x - y] -= seg.b * weight[x][y];
            pd1[seg.x1 + 1][seg.x1 + 1 + x - y] -= seg.m * weight[x][y];
          }
        }
      }
    }
    timer0 += clock();
    timer1 -= clock();
    for (int x = 0; x < N; x++) {
      for (int y = 0; y < N; y++) {
        if (x > 0) px0[x][y] += px0[x - 1][y];
        if (x > 0) px1[x][y] += px1[x - 1][y];
        if (y > 0) py0[x][y] += py0[x][y - 1];
        if (y > 0) py1[x][y] += py1[x][y - 1];
        if (x > 0 && y > 0) pd0[x][y] += pd0[x - 1][y - 1];
        if (x > 0 && y > 0) pd1[x][y] += pd1[x - 1][y - 1];
      }
    }
    for (int x = 0; x < N; x++) {
      for (int y = 0; y < N; y++) {
        tmp[x][y] += px0[x][y];
        tmp[x][y] += x * px1[x][y];
        tmp[x][y] += py0[x][y];
        tmp[x][y] += y * py1[x][y];
        tmp[x][y] += pd0[x][y];
        tmp[x][y] += x * pd1[x][y];
      }
    }
    for (int x = 0; x < N; x++) {
      for (int y = 0; y < N; y++) {
        tmp[x][y] /= N * (N + 1) / 2.0;
      }
    }
    timer1 += clock();
    timer2 -= clock();
    memcpy(weight, tmp, sizeof(tmp));
    memset(tmp, 0, sizeof(tmp));
    memset(px0, 0, sizeof(px0));
    memset(px1, 0, sizeof(px1));
    memset(py0, 0, sizeof(py0));
    memset(py1, 0, sizeof(py1));
    memset(pd0, 0, sizeof(pd0));
    memset(pd1, 0, sizeof(pd1));
    timer2 += clock();
  }
  double ans = 0;
  for (int i = 0; i < N; i++) {
    for (int j = 0; j < i; j++) {
      ans += weight[i][j];
    }
  }
  printf("%.10lf\n", ans);
}