OFFSET
1,3
COMMENTS
4 is the only zeroless number < 10^7 that is a member of this sequence and A229185 (Numbers incremented by their digit product produce a cube).
EXAMPLE
167 - 1*6*7 = 125 = 5^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
More terms and prepended a(1) = 0 from Derek Orr, Mar 12 2015
STATUS
approved