Function _integer_sqrt_of_frac_rto

Lib/statistics.py:1705–1709

Square root of n/m, rounded to the nearest integer using round-to-odd.

(n: int, m: int)

def _integer_sqrt_of_frac_rto(n: int, m: int) -> int:
    """Square root of n/m, rounded to the nearest integer using round-to-odd."""
    # Reference: https://www.lri.fr/~melquion/doc/05-imacs17_1-expose.pdf
    a = math.isqrt(n // m)
    return a | (a*a*m != n)


# For 53 bit precision floats, the bit width used in