OFFSET
1,2
LINKS
Michael S. Branicky, Table of n, a(n) for n = 1..10000
EXAMPLE
1158 is a term because its digital sum is 1+1+5+8 = 15 and "15" is a substring of 1158 and of 1158^3 = 1552836312.
PROG
(Python)
def sd(n): return sum(map(int, str(n)))
def ok(n): s = str(sd(n)); return s in str(n) and s in str(n**3)
print([k for k in range(1, 5411) if ok(k)]) # Michael S. Branicky, Jan 31 2022
CROSSREFS
KEYWORD
nonn,base
AUTHOR
Claudio Meller, Jun 24 2009
EXTENSIONS
a(47) and beyond from Michael S. Branicky, Jan 31 2022
STATUS
approved