|
|
A292562
|
|
Odd numbers n such that 2^d(n) == 2^n (mod n) where d = A000005.
|
|
0
|
|
|
1, 9, 19313, 226703, 345023, 673663, 929783, 1697609, 1790159, 1878073, 2391991, 3441743, 8026103, 9520441, 9770153, 10751841, 30652223, 35482433, 40642103, 41839231, 50358247, 57233047, 65046479, 69601193, 85246441, 110944089, 117368839, 119611073, 122527433
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,2
|
|
COMMENTS
|
Square terms are 1, 3^2, 3^2*1093^2, 3^2*3511^2 and there are no other squares k^2 with k < 3*10^9. Note that 1093 and 3511 are Wieferich primes (A001220). - Giovanni Resta, Sep 19 2017
|
|
LINKS
|
|
|
EXAMPLE
|
19313 = 7*31*89 is a term because d(19313) = 2^3 and 7*31*89 divides 2^(7*31*89) - 2^8.
|
|
MATHEMATICA
|
Select[Range[1, 10^7, 2], PowerMod[2, DivisorSigma[0, #], #] == PowerMod[2, #, #] &] (* Michael De Vlieger, Sep 19 2017 *)
|
|
PROG
|
(PARI) isok(n) = n%2==1 && Mod(2, n)^n==2^numdiv(n);
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn
|
|
AUTHOR
|
|
|
STATUS
|
approved
|
|
|
|