#include <cassert>
#include <cmath>
#include <utility>

template <typename Domain>
auto midpoint(Domain a, Domain b) { return (a + b) / 2; }

template <typename Domain>
auto distance(Domain a, Domain b) { using std::abs; return abs(a - b); }

template <typename Function, typename Domain, typename Distance>
std::pair<Domain, Domain> 
root_bisection(Function f, Domain a, Domain b, Distance epsilon, int steps = 50)
{
  assert(f(a) * f(b) < 0);
  auto ya = f(a);

  while (distance(a, b) > epsilon && steps-- > 0)
  {
    auto m = midpoint(a, b);
    auto ym = f(m);
    if (ya * ym < 0) { b = m; }
    else { a = m; ya = ym; }
  }
  return {a, b};
}