OFFSET
1,1
LINKS
Vincenzo Librandi, Table of n, a(n) for n = 1..1000
MATHEMATICA
ok [p_]:=Reduce[Mod[x^4 - 10, p]== 0, x, Integers]=!= False; Select[Prime[Range[180]], ok] (* Vincenzo Librandi, Sep 12 2012 *)
PROG
(PARI) isA040119(p)={r=0; for(m=0, p-1, if(Mod(m, p)^4==Mod(10, p), r=1)); r}
(Magma) [p: p in PrimesUpTo(800) | exists(t){x : x in ResidueClassRing(p) | x^4 eq 10}]; // Vincenzo Librandi, Sep 12 2012
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
EXTENSIONS
Definition corrected and PARI program added by Michael B. Porter, Oct 13 2009
STATUS
approved