faster ``is_quadratic_residue``;
in ``sqrt_mod`` changed default value of ``limit`` to 1;
changed a test and added another one