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;
double* lp = new double[1000000];
pair<int, int>* lpvp = new pair<int, int>[1000000];
vector<double> ir;
vector<pair<pair<int, int>, pair<double, vector<int> > > > plp;
double tt(int l, int r, int ind) {
  lpvp[ind] = make_pair(l, r);
  if (l == r) {
    lp[ind] = ir[l];
    return lp[ind];
  }
  int mid = (l + r) / 2;
  lp[ind] = max(tt(l, mid, ind * 2), tt(mid + 1, r, ind * 2 + 1));
  return lp[ind];
}
double mv(int l, int r, int ind) {
  if (lpvp[ind].first > r || lpvp[ind].second < l) {
    return -1;
  }
  if (lpvp[ind].first >= l && lpvp[ind].second <= r) {
    return lp[ind];
  }
  return max(mv(l, r, ind * 2), mv(l, r, ind * 2 + 1));
}
vector<double> tml(int ind) {
  vector<double> vv;
  int l = plp[ind].first.first;
  int r = plp[ind].first.second;
  double t = plp[ind].second.first;
  vector<int> ch = plp[ind].second.second;
  vv.push_back(1 - t);
  vv.push_back(t);
  if (ch.size() == 0) {
    return vv;
  }
  vector<double> ncvd;
  int hn;
  vector<double> chal;
  int cl;
  int cr;
  int cm;
  vector<double> cvd;
  int crm = -1;
  for (int i = 0; i < ch.size(); i++) {
    ncvd.clear();
    hn = -1;
    chal.clear();
    chal = tml(ch[i]);
    cl = plp[ch[i]].first.first;
    cr = plp[ch[i]].first.second;
    cm = mv(cl, cr, 1);
    if (i == 0) {
      cvd = chal;
      crm = cm;
      continue;
    }
    int fild[cvd.size()][chal.size()];
    double prob[cvd.size()][chal.size()];
    for (int a = 0; a < cvd.size(); a++) {
      for (int b = 0; b < chal.size(); b++) {
        int cma = crm + a - max(crm, cm);
        int cmb = cm + b - max(crm, cm);
        int fab = max(cma, cmb);
        fild[a][b] = fab;
        prob[a][b] = cvd[a] * chal[b];
        if (fild[a][b] > hn) {
          hn++;
          ncvd.push_back(prob[a][b]);
        } else {
          ncvd[fild[a][b]] += prob[a][b];
        }
      }
    }
    cvd = ncvd;
    crm = max(crm, cm);
  }
  ncvd.clear();
  hn = -1;
  chal.clear();
  chal = vv;
  cl = l;
  cr = r;
  cm = mv(cl, cr, 1);
  int fild[cvd.size()][chal.size()];
  double prob[cvd.size()][chal.size()];
  for (int a = 0; a < cvd.size(); a++) {
    for (int b = 0; b < chal.size(); b++) {
      int cma = crm + a - max(crm, cm);
      int cmb = cm + b - max(crm, cm);
      int fab = 0;
      if (cma >= 0) {
        fab += cma;
      }
      if (cmb >= 0) {
        fab += cmb;
      }
      fild[a][b] = fab;
      prob[a][b] = cvd[a] * chal[b];
      if (fild[a][b] > hn) {
        hn++;
        ncvd.push_back(prob[a][b]);
      } else {
        ncvd[fild[a][b]] += prob[a][b];
      }
    }
  }
  cvd = ncvd;
  crm = max(crm, cm);
  return cvd;
}
int main() {
  int n, q;
  cin >> n >> q;
  for (int i = 0; i < n; i++) {
    double d;
    cin >> d;
    ir.push_back(d);
  }
  tt(0, n - 1, 1);
  vector<pair<pair<int, int>, double> > sl;
  sl.push_back(make_pair(make_pair(0, (n - 1) * -1), 0));
  for (int i = 0; i < q; i++) {
    int vas, par;
    double ton;
    cin >> vas >> par >> ton;
    vas--;
    par--;
    sl.push_back(make_pair(make_pair(vas, -1 * par), ton));
  }
  sort(sl.begin(), sl.end());
  int ci = 0;
  for (int i = 0; i <= q; i++) {
    int l = sl[i].first.first;
    int r = sl[i].first.second * -1;
    double t = sl[i].second;
    vector<int> vi;
    if (ci == 0) {
      plp.push_back(make_pair(make_pair(l, r), make_pair(t, vi)));
      ci++;
      continue;
    }
    for (int j = ci - 1; j >= 0; j--) {
      int el = plp[j].first.first;
      int er = plp[j].first.second;
      if (el <= l && er >= r) {
        plp[j].second.second.push_back(ci);
        plp.push_back(make_pair(make_pair(l, r), make_pair(t, vi)));
        ci++;
        break;
      }
    }
  }
  vector<double> vv = tml(0);
  double vas = 0;
  for (int i = 0; i < vv.size(); i++) {
    vas += vv[i] * i;
  }
  vas += mv(0, n - 1, 1);
  cout.precision(15);
  cout << vas << endl;
  cin >> n;
}