OFFSET
1,1
COMMENTS
Primes p == 3 (mod 4) are precisely the rational primes in the ring of Gaussian integers.
5 is a quadratic residue of integers modulo p for p being a term of this sequence. (See A385165).
LINKS
Jianing Song, Table of n, a(n) for n = 1..10000
EXAMPLE
5479 is a term since (2+-i)^125081 == 1 (mod 5479), and 125081 is divisible by neither 2 nor 3.
PROG
(PARI) ord(p) = my(d = divisors((p+1)*znorder(Mod(5, p)))); for(i=1, #d, if(Mod([2, -1; 1, 2], p)^d[i] == 1, return(d[i]))) \\ for a prime p == 3 (mod 4), returns ord(2+-i, p)
isA385188(p) = isprime(p) && p%4==3 && ord(p)%2 && ord(p)%3
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Jianing Song, Jun 20 2025
STATUS
approved
