``sqrt_mod(a, p)`` returns a root ``x <= p // 2``