a binary search on floating point numbers / 二分探索 (浮動小数点数)
(utils/binary_search_float.hpp)

• View this file on GitHub
• Last update: 2020-05-23 03:04:39+09:00
• Include: #include "utils/binary_search_float.hpp"

Depends on

• utils/macros.hpp

Code

#pragma once
#include <cassert>
#include "utils/macros.hpp"

/**
 * @brief a binary search on floating point numbers / 二分探索 (浮動小数点数)
 */
template <int Iteration = 64, typename UnaryPredicate>
double binsearch_float(double l, double r, UnaryPredicate p) {
    assert (l <= r);
    REP (i, Iteration) {
        double m = (l + r) / 2;
        (p(m) ? r : l) = m;
    }
    return r;
}

#line 2 "utils/binary_search_float.hpp"
#include <cassert>
#line 2 "utils/macros.hpp"
#define REP(i, n) for (int i = 0; (i) < (int)(n); ++ (i))
#define REP3(i, m, n) for (int i = (m); (i) < (int)(n); ++ (i))
#define REP_R(i, n) for (int i = (int)(n) - 1; (i) >= 0; -- (i))
#define REP3R(i, m, n) for (int i = (int)(n) - 1; (i) >= (int)(m); -- (i))
#define ALL(x) std::begin(x), std::end(x)
#line 4 "utils/binary_search_float.hpp"

/**
 * @brief a binary search on floating point numbers / 二分探索 (浮動小数点数)
 */
template <int Iteration = 64, typename UnaryPredicate>
double binsearch_float(double l, double r, UnaryPredicate p) {
    assert (l <= r);
    REP (i, Iteration) {
        double m = (l + r) / 2;
        (p(m) ? r : l) = m;
    }
    return r;
}

Back to top page