OFFSET
1,1
COMMENTS
LINKS
MATHEMATICA
ok[p_]:= Reduce[Mod[x^44 - 2, p] == 0, x, Integers] =!= False; Select[Prime[Range[300]], ok] (* Vincenzo Librandi, Sep 14 2012 *)
PROG
(Magma) [p: p in PrimesUpTo(1100) | exists(t){x : x in ResidueClassRing(p) | x^44 eq 2}]; // Vincenzo Librandi, Sep 14 2012
(PARI)
N=10^4; default(primelimit, N);
ok(p, r, k)={ return ( Mod(r, p)^((p-1)/gcd(k, p-1)) == 1 ); }
forprime(p=2, N, if (ok(p, 2, 44), print1(p, ", ")));
/* Joerg Arndt, Sep 21 2012 */
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
STATUS
approved