Problem: Malek is a rich man. He also is very generous. That's why he decided to split his money between poor people. A charity institute knows n poor people numbered from 1 to n. The institute gave Malek q recommendations. A recommendation is a segment of people like [l, r] which means the institute recommended that Malek gives one dollar to every person whose number is in this segment.

However this charity has very odd rules about the recommendations. Because of those rules the recommendations are given in such a way that for every two recommendation [a, b] and [c, d] one of the following conditions holds: 

  * The two segments are completely disjoint. More formally either a ≤ b < c ≤ d or c ≤ d < a ≤ b
  * One of the two segments are inside another. More formally either a ≤ c ≤ d ≤ b or c ≤ a ≤ b ≤ d. 



The goodness of a charity is the value of maximum money a person has after Malek finishes giving his money. The institute knows for each recommendation what is the probability that Malek will accept it. They want to know the expected value of goodness of this charity. So they asked you for help.

You have been given the list of recommendations and for each recommendation the probability of it being accepted by Malek. You have also been given how much money each person initially has. You must find the expected value of goodness.

Input

In the first line two space-separated integers n, q (1 ≤ n ≤ 105, 1 ≤ q ≤ 5000) are given.

In the second line n space-separated integers a1, a2, ..., an (0 ≤ ai ≤ 109) are given meaning that person number i initially has ai dollars. 

Each of the next q lines contains three space-separated numbers li, ri, pi (1 ≤ li ≤ ri ≤ n, 0 ≤ p ≤ 1) where li and ri are two integers describing the segment of recommendation and pi is a real number given with exactly three digits after decimal point which is equal to probability of Malek accepting this recommendation.

Note that a segment may appear several times in recommendations.

Output

Output the sought value. Your answer will be considered correct if its absolute or relative error is less than 10 - 6.

Examples

Input

5 2
1 7 2 4 3
1 3 0.500
2 2 0.500


Output

8.000000000


Input

5 2
281 280 279 278 282
1 4 1.000
1 4 0.000


Output

282.000000000


Input

3 5
1 2 3
1 3 0.500
2 2 0.250
1 2 0.800
1 1 0.120
2 2 0.900


Output

4.465000000

Solution (in CPP):

Solution:
#include <bits/stdc++.h>
using namespace std;
template <class T>
inline void rd(T &x) {
  char ch;
  x = 0;
  bool fl = false;
  while (!isdigit(ch = getchar())) (ch == '-') && (fl = true);
  for (x = (ch ^ '0'); isdigit(ch = getchar()); x = x * 10 + (ch ^ '0'))
    ;
  (fl == true) && (x = -x);
}
template <class T>
inline void output(T x) {
  if (x / 10) output(x / 10);
  putchar(x % 10 + '0');
}
template <class T>
inline void ot(T x) {
  if (x < 0) putchar('-'), x = -x;
  output(x);
  putchar(' ');
}
template <class T>
inline void prt(T a[], int st, int nd) {
  for (register int i = st; i <= nd; ++i) ot(a[i]);
  putchar('\n');
}
namespace Modulo {
const int mod = 998244353;
inline int ad(int x, int y) { return x + y >= mod ? x + y - mod : x + y; }
inline int sub(int x, int y) { return ad(x, mod - y); }
inline int mul(int x, int y) { return (long long)x * y % mod; }
inline void inc(int &x, int y) { x = ad(x, y); }
inline void inc2(int &x, int y) { x = mul(x, y); }
inline int qm(int x, int y = mod - 2) {
  int ret = 1;
  while (y) {
    if (y & 1) ret = mul(x, ret);
    x = mul(x, x);
    y >>= 1;
  }
  return ret;
}
template <class... Args>
inline int ad(const int a, const int b, const Args &...args) {
  return ad(ad(a, b), args...);
}
template <class... Args>
inline int mul(const int a, const int b, const Args &...args) {
  return mul(mul(a, b), args...);
}
}  // namespace Modulo
namespace Miracle {
const int N = 1e5 + 5000 + 6;
int n, m;
int a[N];
struct node {
  int nxt, to;
} e[2 * N];
int hd[N], cnt;
void add(int x, int y) {
  e[++cnt].nxt = hd[x];
  e[cnt].to = y;
  hd[x] = cnt;
}
vector<pair<int, double> > f[5005], tmp;
int del[N], dc;
int id[N];
int cur;
int nc() {
  int r = dc ? del[dc--] : ++cur;
  f[r].clear();
  return r;
}
void dele(int r) { del[++dc] = r; }
vector<double> c[N];
int tot;
void merge(vector<pair<int, double> > &f, vector<pair<int, double> > &g) {
  if (f.empty())
    f.swap(g);
  else {
    tmp.clear();
    int dw = max(f.front().first, g.front().first);
    int up = max(f.back().first, g.back().first);
    for (register int i = dw; i <= up; ++i) {
      tmp.push_back(make_pair(i, 0.0));
    }
    for (const auto &x : f) {
      for (const auto &y : g) {
        tmp[max(x.first, y.first) - dw].second += x.second * y.second;
      }
    }
    f.swap(tmp);
  }
}
void dfs(int x) {
  id[x] = nc();
  for (register int i = hd[x]; i; i = e[i].nxt) {
    int y = e[i].to;
    dfs(y);
    merge(f[id[x]], f[id[y]]);
    dele(id[y]);
  }
  if (!hd[x]) {
    f[id[x]].push_back(make_pair(a[x], 1.00));
  }
  for (const auto &p : c[x]) {
    f[id[x]].push_back(make_pair(f[id[x]].back().first + 1, 0));
    for (register int i = f[id[x]].size() - 1; i > 0; --i) {
      f[id[x]][i].second =
          f[id[x]][i].second * (1 - p) + f[id[x]][i - 1].second * p;
    }
    f[id[x]][0].second *= (1 - p);
  }
}
struct qs {
  int l, r;
  double p;
  bool friend operator<(qs a, qs b) {
    if (a.r - a.l == b.r - b.l) return a.l < b.l;
    return a.r - a.l < b.r - b.l;
  }
} q[N];
int fa[N], R[N];
int fin(int x) { return fa[x] == x ? x : fa[x] = fin(fa[x]); }
int main() {
  rd(n);
  rd(m);
  for (register int i = 1; i <= n; ++i) rd(a[i]);
  for (register int i = 1; i <= n; ++i) {
    fa[i] = 0;
  }
  for (register int i = 1; i <= m; ++i) {
    rd(q[i].l);
    rd(q[i].r);
    scanf("%lf", &q[i].p);
  }
  ++m;
  q[m].l = 1;
  q[m].r = n;
  q[m].p = 0;
  sort(q + 1, q + m + 1);
  tot = n;
  for (register int i = 1; i <= m; ++i) {
    if (q[i].l == q[i].r) {
      c[q[i].l].push_back(q[i].p);
    } else {
      if (i == 1 || (q[i].l != q[i - 1].l) || (q[i].r != q[i - 1].r)) {
        ++tot;
        int j = q[i].l;
        while (j <= q[i].r) {
          if (!fa[j]) {
            fa[j] = tot;
            R[j] = q[i].r;
            add(tot, j);
            ++j;
          } else {
            int nxt = R[j] + 1;
            add(tot, fa[j]);
            fa[j] = tot;
            R[j] = q[i].r;
            j = nxt;
          }
        }
        c[tot].push_back(q[i].p);
      } else {
        c[tot].push_back(q[i].p);
      }
    }
  }
  dfs(tot);
  double ans = 0;
  for (const auto &y : f[id[tot]]) {
    ans += y.first * y.second;
  }
  printf("%.10lf", ans);
  return 0;
}
}  // namespace Miracle
signed main() {
  Miracle::main();
  return 0;
}