A014957 Positive integers k that divide 16^k - 1.

1, 3, 5, 9, 15, 21, 25, 27, 39, 45, 55, 63, 75, 81, 105, 117, 125, 135, 147, 155, 165, 171, 189, 195, 205, 225, 243, 273, 275, 315, 333, 351, 375, 405, 441, 465, 495, 507, 513, 525, 567, 585, 605, 609, 615, 625, 657, 675, 729, 735, 775, 819, 825, 855, 903

Offset: 1

Author

Olivier Gérard

Keywords

• nonn

Comments

Also, positive integers k that divide A014931(k).

Links

• Harvey P. Dale, Table of n, a(n) for n = 1..2000

Crossrefs

Cf. A067945, A014945, A067946, A014946, A067947, A014949, A068382, A014950, A068383, A014951, A116621, A014956, A177805, A177807, A128358, A128360

Cf. A128396, A177916

Programs

Mathematica

Join[{1},Select[Range[1000],PowerMod[16,#,#]==1&]] (* Harvey P. Dale, Jun 12 2024 *)

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

Extensions

Edited by Max Alekseyev, Sep 10 2011