OFFSET
1,2
EXAMPLE
681 + 6*8*1 = 729 = 9^3.
MATHEMATICA
Select[Range[0, 100000], IntegerQ[(# + Times @@ IntegerDigits[#])^(1/3)] &] (* T. D. Noe, Sep 16 2013 *)
PROG
(Python)
def DP(n):
p = 1
for i in str(n):
p *= int(i)
return p
for n in range(10**4):
k = 0
P = n + DP(n)
while P >= k**3:
if P == k**3:
print(n, end=', ')
break
k += 1
# Simplified by Derek Orr, Mar 12 2015
(PARI) for(n=0, 10^5, d=digits(n); P=n+prod(i=1, #d, d[i]); if(ispower(P, 3), print1(n, ", "))) \\ Derek Orr, Mar 12 2015
CROSSREFS
KEYWORD
nonn,easy,base
AUTHOR
Derek Orr, Sep 15 2013
EXTENSIONS
Prepended a(1) = 0 from Derek Orr, Mar 12 2015
STATUS
approved