OFFSET
1,1
COMMENTS
Union of 2, 3, A068231 (primes congruent to 11 modulo 12), primes p == 1 (mod 4) such that 3^((p-1)/4) == 1 (mod p). - Jianing Song, Jun 22 2025
LINKS
Vincenzo Librandi, Table of n, a(n) for n = 1..1000
MATHEMATICA
ok [p_]:=Reduce[Mod[x^4- 3, 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 3}]; // Vincenzo Librandi, Sep 11 2012
(PARI) isA040101(p) = isprime(p) && (p==2 || p==3 || p%12==11 || (p%4==1 && Mod(3, p)^((p-1)/4) == 1)) \\ Jianing Song, Jun 22 2025
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
STATUS
approved
