OFFSET
1,2
LINKS
Michael S. Branicky, Table of n, a(n) for n = 1..10000
EXAMPLE
1145 is a term because its digital sum is 1+1+4+5 = 11 and "11" is a substring of 1145, of 1145^2 = 1311025 and of 1145^3 = 1501123625.
MATHEMATICA
sddQ[n_]:=With[{d=IntegerDigits[Total[IntegerDigits[n]]]}, SequenceCount[IntegerDigits[n], d]>0&&SequenceCount[IntegerDigits[n^2], d]>0&&SequenceCount[IntegerDigits[n^3], d]>0]; Select[Range[11000], sddQ] (* Harvey P. Dale, Jun 16 2025 *)
PROG
(Python)
def sd(n): return sum(map(int, str(n)))
def ok(n): s = str(sd(n)); return all(s in str(t) for t in [n, n**2, n**3])
print([k for k in range(1, 10140) if ok(k)]) # Michael S. Branicky, Jan 31 2022
CROSSREFS
KEYWORD
nonn,base
AUTHOR
Claudio Meller, Jun 24 2009
EXTENSIONS
a(21) and beyond from Michael S. Branicky, Jan 31 2022
STATUS
approved
