Problem: Problem

A dolphin who lives in a certain aquarium will be rewarded if he jumps and goes through the $ N $ ring.

* Dolphin jumps from coordinates $ (0,0) $ and lands at $ (T,0) $.
* The trajectory of the jump is a parabola.
* The $ i $ th ring is determined to have passed through when the jump trajectory intersects the line segment connecting $ (X_i, L_i) $ and $ (X_i, H_i) $.
* $ 1 $ jump requires as much physical strength as the initial velocity.



Dolphins can jump as many times as they want. Let the gravitational acceleration vector be $ (0, -1) $ and find the minimum total physical strength required for the dolphin to pass through all the rings. However, suppose friction and air resistance are negligibly small.

Constraints

The input satisfies the following conditions.

* All inputs are integers.
* $ 1 \ le X_i <T \ le 10 ^ 6 $
* $ 1 \ le N \ le 10 ^ 5 $
* $ 1 \ le L_i <H_i \ le 10 ^ 6 $

Input

The input is given in the following format.


$ T $ $ N $
$ X_1 $ $ L_1 $ $ H_1 $
$ \ vdots $
$ X_N $ $ L_N $ $ H_N $


First, $ T $ and $ N $ are given to the $ 1 $ line. Then the $ N $ line is given the position of the $ i $ th ring, $ X_i $, $ L_i $, and $ H_i $.

Output

Output the answer on one line. If the absolute error or relative error is $ 10 ^ {-9} $ or less, the answer is judged to be correct.

Examples

Input

100 5
50 1 5
50 5 10
50 20 30
50 40 60
50 61 1000000


Output

48.6090201099


Input

64 15
38 133177 927361
48 177920 668766
12 680425 790550
43 6853 384115
17 214954 723798
62 63843 153825
28 399349 482937
2 336136 367001
33 138008 733496
6 203462 911631
58 321974 527734
17 696940 781678
55 265874 507640
41 56037 880001
34 279422 528651


Output

6087.909851326286

Solution (in CPP):

Solution:
#include<bits/stdc++.h>
using namespace std;

class fraction {
    long long numer; // bunshi
    long long denom; // bunbo
public:
    fraction(long long numer, long long denom) : numer(numer), denom(denom) {}
    int compare_to(const fraction& f) const {
        assert(denom > 0 && f.denom > 0);
        long long res = numer * f.denom - denom * f.numer;
        if (res > 0) return 1;
        else if (res < 0) return -1;
        else return 0;
    }
    long double to_float() const {
        return numer / (long double)denom;
    }
};

using P = pair<fraction, fraction>;


int main() {
    long long t;
    int n;
    cin >> t >> n;
    vector<P> upper, lower, middle;
    for (int i = 0; i < n; ++i) {
        long long x, a, b;
        cin >> x >> a >> b;
        fraction ta(a, x * (t - x));
        fraction tb(b, x * (t - x));
        if (a * t > x * (t - x)) upper.push_back({ta, tb});
        else if (b * t < x * (t - x)) lower.push_back({ta, tb});
        else middle.push_back({ta, tb});
    }
    sort(upper.begin(), upper.end(), [](const P& a, const P& b) { return a.first.compare_to(b.first) > 0; });
    sort(lower.begin(), lower.end(), [](const P& a, const P& b) { return a.second.compare_to(b.second) < 0; });

    long double ans = 0;
    fraction last1(1e6 + 10, 1), last2(-(1e6 + 10), 1);
    for (auto& i : upper) {
        if (last1.compare_to(i.second) > 0) {
            ans += sqrt(i.first.to_float() * t * t / 2 + 1 / (i.first.to_float() * 2));
            last1 = i.first;
        }
    }
    for (auto& i : lower) {
        if (last2.compare_to(i.first) < 0) {
            ans += sqrt(i.second.to_float() * t * t / 2 + 1 / (i.second.to_float() * 2));
            last2 = i.second;
        }
    }
    bool ad = false;
    for (auto& i : middle) {
        if ((i.second.compare_to(last1) < 0 || i.first.compare_to(last1) > 0) && (i.second.compare_to(last2) < 0 || i.first.compare_to(last2) > 0)) {
            ad = true;
            break;
        }
    }
    if (ad) ans += sqrt(t);
    cout << setprecision(24) << ans << '\n';
    return 0;
}