OFFSET
1,2
COMMENTS
Whenever 15|a(n) (n = 32, 36, 40, 86, 100, 135, 143, 194, 197, 201), then a(n+1) = a(n) + 1. Zero also satisfies the definition (n = Sum_{i=1..k} d[i]^k where d[1..k] are the base-15 digits of n), but this sequence only considers positive terms. - M. F. Hasler, Nov 22 2019
LINKS
Joseph Myers, Table of n, a(n) for n = 1..202 (the full list of terms, from Winter)
Henk Koppelaar and Peyman Nasehpour, On Hardy's Apology Numbers, arXiv:2008.08187 [math.NT], 2020.
Eric Weisstein's World of Mathematics, Narcissistic Number
D. T. Winter, Table of Armstrong Numbers
MATHEMATICA
Select[Range[10^7], # == Total[IntegerDigits[#, 15]^IntegerLength[#, 15]] &] (* Michael De Vlieger, Nov 04 2020 *)
PROG
(PARI) select( is_A161952(n)={n==vecsum([d^#n|d<-n=digits(n, 15)])}, [1..10^5]) \\ M. F. Hasler, Nov 22 2019
CROSSREFS
KEYWORD
base,fini,full,nonn
AUTHOR
Joseph Myers, Jun 22 2009
EXTENSIONS
Terms sorted in increasing order by Pontus von Brömssen, Mar 03 2019
STATUS
approved