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;
const int N = 1e5 + 10, LOG = 20, Q = 5010;
pair<pair<int, int>, long double> s[N];
int rmq[N][LOG];
vector<int> adj[N];
vector<pair<int, long double> > pi[Q];
int n, q, a[N];
bool vis[N];
void build() {
  for (int j = 0; j < LOG; j++) {
    for (int i = 0; i < n; i++) {
      if (j == 0) {
        rmq[i][0] = a[i];
        continue;
      }
      if (i + (1 << j) <= n)
        rmq[i][j] = max(rmq[i][j - 1], rmq[i + (1 << (j - 1))][j - 1]);
    }
  }
}
int get(int l, int r) {
  if (l == r) return 0;
  int dis = r - l, ret = -1;
  ;
  for (int i = 0; i < LOG; i++) {
    if (dis & (1 << i)) {
      ret = max(ret, rmq[l][i]);
      l += (1 << i);
    }
  }
  return ret;
}
vector<pair<int, long double> > merge1(vector<pair<int, long double> >& a,
                                       vector<pair<int, long double> >& b) {
  if (a.size() == 0) return b;
  if (b.size() == 0) return a;
  long double ps1[Q], ps2[Q];
  vector<pair<int, long double> > ret;
  int sz1 = a.size(), sz2 = b.size();
  ps1[0] = ps2[0] = 0;
  for (int i = 0; i < sz1; i++) {
    ps1[i + 1] = ps1[i] + a[i].second;
  }
  for (int i = 0; i < sz2; i++) {
    ps2[i + 1] = ps2[i] + b[i].second;
  }
  for (int i = 0; i < sz1; i++) {
    int value = a[i].first;
    int ind =
        lower_bound(b.begin(), b.end(), pair<int, long double>(value, -1)) -
        b.begin();
    long double p = a[i].second, ans = 0;
    ans += p * ps2[ind];
    if (ind < sz2 && b[ind].first == value) {
      ans += p * b[ind].second;
      ans += b[ind].second * ps1[i];
    }
    ret.push_back({value, ans});
  }
  for (int i = 0; i < sz2; i++) {
    int value = b[i].first;
    int ind =
        lower_bound(a.begin(), a.end(), pair<int, long double>(value, -1)) -
        a.begin();
    long double p = b[i].second, ans = 0;
    if (ind < sz1 && value == a[ind].first) continue;
    ans += p * ps1[ind];
    ret.push_back({value, ans});
  }
  sort(ret.begin(), ret.end());
  return ret;
}
vector<pair<int, long double> > merge2(vector<pair<int, long double> >& a,
                                       long double p) {
  int sz = a.size();
  int mnval = a[0].first, mxval = a[sz - 1].first + 1;
  vector<pair<int, long double> > ret;
  for (int i = mnval; i <= mxval; i++) {
    int ind = lower_bound(a.begin(), a.end(), pair<int, long double>(i, -1)) -
              a.begin();
    long double ans = 0;
    if (ind < sz && a[ind].first == i)
      ans += (long double)(1 - p) * a[ind].second;
    ind--;
    if (ind >= 0 && a[ind].first == i - 1)
      ans += (long double)p * a[ind].second;
    if (ans > 0) ret.push_back({i, ans});
  }
  return ret;
}
void dfs(int v, int p = -1) {
  int ls = s[v].first.first, maxi = -1;
  vis[v] = true;
  for (auto u : adj[v]) {
    if (u == p) continue;
    if (vis[u]) continue;
    dfs(u, v);
    int l = s[u].first.first, r = -s[u].first.second;
    if (ls < l) {
      maxi = max(maxi, get(ls, l));
      ls = r + 1;
    }
    pi[v] = merge1(pi[v], pi[u]);
  }
  if (ls <= -s[v].first.second)
    maxi = max(maxi, get(ls, -s[v].first.second + 1));
  if (maxi != -1) {
    vector<pair<int, long double> > tmp;
    tmp.push_back({maxi, (long double)1});
    pi[v] = merge1(pi[v], tmp);
  }
  pi[v] = merge2(pi[v], s[v].second);
  return;
}
int main() {
  ios_base::sync_with_stdio(0);
  cin.tie(0);
  cin >> n >> q;
  for (int i = 0; i < n; i++) cin >> a[i];
  build();
  for (int i = 0; i < q; i++) {
    cin >> s[i].first.first >> s[i].first.second >> s[i].second;
    s[i].first.first--;
    s[i].first.second--;
    s[i].first.second *= (-1);
  }
  s[q].first.first = 0;
  s[q].first.second = 1 - n;
  s[q].second = 0;
  sort(s, s + q + 1);
  for (int i = 0; i <= q; i++) {
    for (int j = 0; j <= q; j++) {
      if (i == j) continue;
      int l1 = s[i].first.first, r1 = -s[i].first.second, l2 = s[j].first.first,
          r2 = -s[j].first.second;
      if (l1 < l2 && r2 <= r1)
        adj[i].push_back(j);
      else if (l1 <= l2 && r2 < r1)
        adj[i].push_back(j);
      else if (l1 == l2 && r1 == r2 && i < j)
        adj[i].push_back(j);
    }
  }
  dfs(0);
  long double ans = 0;
  for (auto X : pi[0]) {
    long long value = X.first;
    long double p = X.second;
    ans += (long double)p * value;
  }
  cout << fixed << setprecision(12) << ans << '\n';
}