Problem: The Little Elephant loves permutations of integers from 1 to n very much. But most of all he loves sorting them. To sort a permutation, the Little Elephant repeatedly swaps some elements. As a result, he must receive a permutation 1, 2, 3, ..., n.

This time the Little Elephant has permutation p1, p2, ..., pn. Its sorting program needs to make exactly m moves, during the i-th move it swaps elements that are at that moment located at the ai-th and the bi-th positions. But the Little Elephant's sorting program happened to break down and now on every step it can equiprobably either do nothing or swap the required elements.

Now the Little Elephant doesn't even hope that the program will sort the permutation, but he still wonders: if he runs the program and gets some permutation, how much will the result of sorting resemble the sorted one? For that help the Little Elephant find the mathematical expectation of the number of permutation inversions after all moves of the program are completed.

We'll call a pair of integers i, j (1 ≤ i < j ≤ n) an inversion in permutatuon p1, p2, ..., pn, if the following inequality holds: pi > pj.

Input

The first line contains two integers n and m (1 ≤ n, m ≤ 1000, n > 1) — the permutation size and the number of moves. The second line contains n distinct integers, not exceeding n — the initial permutation. Next m lines each contain two integers: the i-th line contains integers ai and bi (1 ≤ ai, bi ≤ n, ai ≠ bi) — the positions of elements that were changed during the i-th move.

Output

In the only line print a single real number — the answer to the problem. The answer will be considered correct if its relative or absolute error does not exceed 10 - 6.

Examples

Input

2 1
1 2
1 2


Output

0.500000000


Input

4 3
1 3 2 4
1 2
2 3
1 4


Output

3.000000000

Solution (in CPP):

Solution:
#include <bits/stdc++.h>
using namespace std;
vector<bool> H[1004][1004];
vector<double> dp[1004][1004];
int i, j, k, n, m, x, y, z, a[(int)(1e5 + 10)];
pair<int, int> B[(int)(1e5 + 10)];
vector<pair<int, int> > C[(int)(1e5 + 10)];
double f(int x, int y, int t) {
  int p1 =
      lower_bound(C[x].begin(), C[x].end(), make_pair(t, 0)) - C[x].begin();
  int p2 =
      lower_bound(C[y].begin(), C[y].end(), make_pair(t, 0)) - C[y].begin();
  double &r = dp[x][y][p1 + p2];
  if (H[x][y][p1 + p2]) return r;
  H[x][y][p1 + p2] = 1;
  int l2 = C[y].size();
  int l1 = C[x].size();
  if (p1 == l1 and p2 == l2) return r = (x > y);
  if ((p2 == l2 and p1 != l1) or
      (p1 != l1 and C[x][p1].first < C[y][p2].first)) {
    r = f(C[x][p1].second, y, C[x][p1].first + 1) / 2;
  } else if ((p1 == l1 and p2 != l2) or
             (p2 != l2 and C[y][p2].first < C[x][p1].first))
    r = f(x, C[y][p2].second, C[y][p2].first + 1) / 2;
  else if (p1 != l1 and p2 != l2 and C[x][p1].first == C[y][p2].first) {
    r = f(y, x, C[x][p1].first + 1) / 2;
  }
  int k = m;
  if (p1 != l1) k = C[x][p1].first;
  if (p2 != l2) k = min(k, C[y][p2].first);
  return r += f(x, y, k + 1) / 2;
}
int main() {
  cin >> n >> m;
  for (i = 1; i <= n; i++) scanf("%d", a + i);
  for (i = 1; i <= m; i++) {
    scanf("%d %d", &B[i].first, &B[i].second);
    C[B[i].first].push_back(make_pair(i, B[i].second));
    C[B[i].second].push_back(make_pair(i, B[i].first));
  }
  double ans = 0;
  for (i = 1; i <= n; i++)
    for (j = 1; j <= n; j++) {
      H[i][j].resize(C[i].size() + C[j].size() + 2);
      dp[i][j].resize(C[i].size() + C[j].size() + 2);
    }
  for (i = 1; i <= n; i++)
    for (j = i + 1; j <= n; j++) {
      if (a[i] < a[j])
        ans += f(i, j, 1);
      else
        ans += 1 - f(i, j, 1);
    }
  cout << setprecision(10) << fixed << ans << '\n';
}