#include <bits/stdc++.h>

namespace RealComp {
    
    const double EPS = std::numeric_limits<double>::epsilon();
    
    struct isEqual {
        bool operator()(double a, double b) const {
            return std::abs(a-b) < 2 * EPS * std::max(std::abs(a), std::abs(b));
        }
    };
    
    struct isLess {
        bool operator()(double a, double b) const {
            return !isEqual()(a,b) && a < b;
        }
    };
}

using namespace std;

set<double, RealComp::isLess> st;

int main() {
    double eps = 3e-16;
    double b = 1.0, a = b - eps, c = b + eps;
    cout <<  std::fixed << std::setprecision(20) << a << " " << b << " " << c << "\n"; // prints 0.99999999999999966693 1.00000000000000000000 1.00000000000000022204
    cout << RealComp::isEqual()(a,b) << " " << RealComp::isEqual()(a,c) << " " <<RealComp::isEqual()(b,c) << "\n"; // prints 1 0 1
    
    st.insert(a);
    st.insert(b);
    st.insert(c);
    cout << st.size() << "\n"; // prints 2
    return 0;
}