OFFSET
1,2
LINKS
Carole Dubois, Table of n, a(n) for n = 1..226 (terms 1..100 from Michael S. Branicky)
Carole Dubois, Pin plot of A133509
FORMULA
PROG
(Python)
def ok(n):
d, lim = 1, 1
while lim < n*9*d: d, lim = d+1, lim*10
return not any(sum(map(int, str(k**n))) == k for k in range(2, lim+1))
for k in range(195):
if ok(k): print(k, end=", ") # Michael S. Branicky, Jul 06 2022
CROSSREFS
KEYWORD
nonn,base
AUTHOR
Farideh Firoozbakht, Dec 04 2007
EXTENSIONS
Description improved by T. D. Noe, Nov 26 2008
Extension by T. D. Noe, Nov 26 2008
Edited by Charles R Greathouse IV, Aug 02 2010
a(1) = 0 and a(46) and beyond from Michael S. Branicky, Jul 06 2022
STATUS
approved