OFFSET
1,1
LINKS
Vincenzo Librandi, Table of n, a(n) for n = 1..1000
MATHEMATICA
ok [p_]:=Reduce[Mod[x^4 - 17, p]== 0, x, Integers]=!= False; Select[Prime[Range[180]], ok] (* Vincenzo Librandi, Sep 12 2012 *)
PROG
(Magma) [p: p in PrimesUpTo(800) | exists{x: x in ResidueClassRing(p) | x^4 eq 17}]; // Vincenzo Librandi, Sep 12 2012
(PARI) /* with workaround for bug with ispower( Mod(17, n), 4) "division by zero" */
select( n->(n==2) || (ispower( Mod(17, n), 4)), primes(1000) )
/* Joerg Arndt, Sep 12 2012 */
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
STATUS
approved