%time a= int(1938)^int(99484)
def legendre(a,p): for i in range(p): if a % p == i^2 % p: is_sqr_modp=3 else: is_sqr_modp=2 return is_sqr_modp
legendre(2,7)
b=range ︠ab674799-b8f0-4e1a-88c0-bfb94c086a44︠