OFFSET
1,2
COMMENTS
Also, positive integers k that divide A014931(k).
LINKS
Harvey P. Dale, Table of n, a(n) for n = 1..2000
MATHEMATICA
Join[{1}, Select[Range[1000], PowerMod[16, #, #]==1&]] (* Harvey P. Dale, Jun 12 2024 *)
PROG
(Python)
A014957_list = [n for n in range(1, 10**6) if n == 1 or pow(16, n, n) == 1] # Chai Wah Wu, Mar 25 2021
CROSSREFS
KEYWORD
nonn
AUTHOR
EXTENSIONS
Edited by Max Alekseyev, Sep 10 2011
STATUS
approved