let shift_right_big_int x y = (* idem division rounding to -oo *) div_big_int x (power_int_positive_big_int 2 y)