aleaxit / gmpy

General Multi-Precision arithmetic for Python 2.6+/3+ (GMP, MPIR, MPFR, MPC)
https://gmpy2.readthedocs.io/en/latest/
GNU Lesser General Public License v3.0
510 stars 86 forks source link

Description of `is_euler_prp` #427

Closed haru-44 closed 1 year ago

haru-44 commented 1 year ago

I’m sorry to bother you again.

https://github.com/aleaxit/gmpy/blob/2beaae28d2274f32e11ebd1512de4755821045bb/src/gmpy_mpz_prp.c#L134-L135

This description is different from the immediately preceding comment. I think the correct description should be a**((n-1)/2) == (a/n) (mod n) [(a/n) is the Jacobi symbol].

haru-44 commented 12 months ago

thanks