login
A057840
Numbers k such that k | 8^k + 6^k + 4^k + 2^k + 1.
1
1, 3, 7, 9, 27, 49, 81, 133, 243, 267, 343, 729, 2187, 2401, 5999, 6561, 14063, 14337, 16807, 17253, 19683, 22329, 33323, 45619, 59049, 75573, 117649, 144531, 177147, 348519, 383913, 531441, 745339, 823543, 911853, 1594323, 2384883, 3610523, 3814857, 4782969, 5764801
OFFSET
1,2
MATHEMATICA
Select[ Range[ 10^6 ], Mod[ PowerMod[ 8, #, # ] + PowerMod[ 6, #, # ] + PowerMod[ 4, #, # ] + PowerMod[ 2, #, # ] + 1, # ] == 0 & ]
PROG
(PARI) is(k)=sum(i=1, 4, Mod(2*i, k)^k)==-1; \\ Jinyuan Wang, Jun 18 2026
(Python)
from itertools import count, islice
def A057840_gen(startvalue=1): # generator of terms >= startvalue
for k in count(max(startvalue, 1)):
a, b = pow(2, k, k), pow(3, k, k)
c = ((a+1)*a+b+1)%k
c = c*a%k
if c == k-1:
yield k
A057840_list = list(islice(A057840_gen(), 30)) # Chai Wah Wu, Jun 19 2026
CROSSREFS
KEYWORD
nonn,changed
AUTHOR
Robert G. Wilson v, Nov 09 2000
EXTENSIONS
a(36)-a(41) from Jinyuan Wang, Jun 18 2026
STATUS
approved