OFFSET
1,1
COMMENTS
LINKS
MAPLE
select(p -> isprime(p) and nops([msolve(x^63-2, p)])>0, [2, seq(2*i+1, i=1..1000)]); # Robert Israel, Nov 03 2014
MATHEMATICA
ok[p_]:= Reduce[Mod[x^63 - 2, p] == 0, x, Integers] =!= False; Select[Prime[Range[150]], ok] (* Vincenzo Librandi, Sep 15 2012 *)
PROG
(Magma) [p: p in PrimesUpTo(600) | exists(t){x : x in ResidueClassRing(p) | x^63 eq 2}]; // Vincenzo Librandi, Sep 15 2012
(PARI) N=10^4;
ok(p, r, k)={ return ( (p==r) || (Mod(r, p)^((p-1)/gcd(k, p-1))==1) ); }
forprime(p=2, N, if (ok(p, 2, 63), print1(p, ", ")));
/* Joerg Arndt, Sep 21 2012 */
(Python)
from itertools import islice
from sympy import nextprime, is_nthpow_residue
def A049595_gen(startvalue=2): # generator of terms >= startvalue
p = max(startvalue-1, 1)
while (p:=nextprime(p)):
if is_nthpow_residue(2, 63, p):
yield p
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
STATUS
approved