OFFSET
1,1
COMMENTS
Union of 2, 5, A122869 (primes congruent to 11 or 19 modulo 20), and primes p == 1 (mod 4) such that 5^((p-1)/4) == 1 (mod p). - Jianing Song, Jun 20 2025
LINKS
Vincenzo Librandi, Table of n, a(n) for n = 1..1000
MATHEMATICA
ok [p_]:=Reduce[Mod[x^4- 5, p] == 0, x, Integers] =!= False; Select[Prime[Range[200]], ok] (* Vincenzo Librandi, Sep 11 2012 *)
PROG
(Magma) [p: p in PrimesUpTo(800) | exists(t){x : x in ResidueClassRing(p) | x^4 eq 5}]; // Vincenzo Librandi, Sep 11 2012
(PARI) isA040105(p) = isprime(p) && (p==2 || p==5 || p%20==11 || p%20==19 || (p%4==1 && Mod(5, p)^((p-1)/4) == 1)) \\ Jianing Song, Jun 20 2025
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
STATUS
approved
