Problem: Farmer John is hosting a tennis tournament with his n cows. Each cow has a skill level si, and no two cows having the same skill level. Every cow plays every other cow exactly once in the tournament, and each cow beats every cow with skill level lower than its own.

However, Farmer John thinks the tournament will be demoralizing for the weakest cows who lose most or all of their matches, so he wants to flip some of the results. In particular, at k different instances, he will take two integers ai, bi (ai < bi) and flip all the results between cows with skill level between ai and bi inclusive. That is, for any pair x, y <image> he will change the result of the match on the final scoreboard (so if x won the match, the scoreboard will now display that y won the match, and vice versa). It is possible that Farmer John will change the result of a match multiple times. It is not guaranteed that ai and bi are equal to some cow's skill level.

Farmer John wants to determine how balanced he made the tournament results look. In particular, he wants to count the number of triples of cows (p, q, r) for which the final leaderboard shows that cow p beats cow q, cow q beats cow r, and cow r beats cow p. Help him determine this number.

Note that two triples are considered different if they do not contain the same set of cows (i.e. if there is a cow in one triple that is not in the other).

Input

On the first line are two space-separated integers, n and k (3 ≤ n ≤ 105; 0 ≤ k ≤ 105). On the next line are n space-separated distinct integers, s1, s2, ..., sn (1 ≤ si ≤ 109), denoting the skill levels of the cows. On the next k lines are two space separated integers, ai and bi (1 ≤ ai < bi ≤ 109) representing the changes Farmer John made to the scoreboard in the order he makes it.

Output

A single integer, containing the number of triples of cows (p, q, r) for which the final leaderboard shows that cow p beats cow q, cow q beats cow r, and cow r beats cow p.

Please do not use the %lld specifier to read or write 64-bit integers in С++. It is preferred to use the cin, cout streams or the %I64d specifier.

Examples

Input

3 2
1 2 3
1 2
2 3


Output

1


Input

5 3
5 9 4 1 7
1 7
2 8
3 9


Output

3

Note

In the first sample, cow 3 > cow 1, cow 3 > cow 2, and cow 2 > cow 1. However, the results between cows 1 and 2 and cows 2 and 3 are flipped, so now FJ's results show that cow 1 > cow 2, cow 2 > cow 3, and cow 3 > cow 1, so cows 1, 2, and 3 form a balanced triple. 

Solution (in CPP):

Solution:
#include <bits/stdc++.h>
using namespace ::std;
const long long maxn = 1e5 + 500;
const long long mod = 1e9 + 7;
const long long inf = 1e9 + 500;
long long a[maxn];
vector<long long> vec[maxn];
vector<long long> sar[maxn];
long long out[maxn];
long long seg[maxn * 4];
bool lazy[maxn * 4];
void shift(long long a, long long tol) {
  if (lazy[a]) {
    long long v = tol / 2;
    lazy[2 * a] = !lazy[2 * a];
    lazy[2 * a + 1] = !lazy[2 * a + 1];
    seg[2 * a] = v - seg[2 * a];
    seg[2 * a + 1] = tol - v - seg[2 * a + 1];
    lazy[a] = 0;
  }
}
void update(long long L, long long R, long long node, long long l,
            long long r) {
  if (l <= L && R <= r) {
    seg[node] = R - L - seg[node];
    lazy[node] = !lazy[node];
    return;
  }
  if (R <= l || r <= L) {
    return;
  }
  shift(node, R - L);
  long long mid = (L + R) / 2;
  update(L, mid, 2 * node, l, r);
  update(mid, R, 2 * node + 1, l, r);
  seg[node] = seg[2 * node] + seg[node * 2 + 1];
}
long long finds(long long L, long long R, long long node, long long l,
                long long r) {
  if (l <= L && R <= r) {
    return seg[node];
  }
  if (R <= l || r <= L) {
    return 0;
  }
  shift(node, R - L);
  long long mid = (L + R) / 2;
  return finds(L, mid, 2 * node, l, r) + finds(mid, R, 2 * node + 1, l, r);
}
long long doaz(long long a) { return a * (a - 1) / 2; }
int main() {
  long long n, m;
  cin >> n >> m;
  for (long long i = 0; i < n; i++) {
    cin >> a[i];
  }
  sort(a, a + n);
  for (long long i = 0; i < m; i++) {
    long long c, b;
    cin >> b >> c;
    b = lower_bound(a, a + n, b) - a;
    c = upper_bound(a, a + n, c) - a - 1;
    if (b != c) {
      vec[b].push_back(c);
      sar[c].push_back(b);
    }
  }
  for (long long i = 0; i < n; i++) {
    for (long long j = 0; j < vec[i].size(); j++) {
      long long v = vec[i][j];
      update(0, n, 1, i, v + 1);
    }
    out[i] = (n - i - 1) - finds(0, n, 1, i + 1, n);
  }
  fill(lazy, lazy + maxn * 4, 0);
  fill(seg, seg + maxn * 4, 0);
  long long ans = 0;
  for (long long i = n - 1; i >= 0; i--) {
    for (long long j = 0; j < sar[i].size(); j++) {
      long long v = sar[i][j];
      update(0, n, 1, v, i);
    }
    out[i] += finds(0, n, 1, 0, i);
    ans += doaz(out[i]);
  }
  cout << n * (n - 1) * (n - 2) / 6 - ans;
}