OFFSET
1,2
EXAMPLE
For 311, (3+1+1) + (3*1*1) = 5 + 3 = 8 = 2^3. Hence, 311 is a member of the sequence.
MATHEMATICA
sdpdQ[n_]:=Module[{idn=IntegerDigits[n]}, IntegerQ[Surd[Total[ idn]+ Times@@idn, 3]]]; Select[Range[0, 1500], sdpdQ] (* Harvey P. Dale, Apr 05 2016 *)
PROG
(Python)
def DS(n):
s = 0
for i in str(n):
s += int(i)
return s
def DP(n):
p = 1
for i in str(n):
p *= int(i)
return p
for n in range(2000):
k = 0
while k**3 <= DS(n)+DP(n):
if k**3 == DS(n)+DP(n):
print(n, end=', ')
break
else:
k += 1
## Edited by Derek Orr, Apr 10 2015
(PARI) for(n=0, 2000, d=digits(n); p=prod(i=1, #d, d[i]); s=sumdigits(n); if(ispower(s+p, 3), print1(n, ", "))) \\ Derek Orr, Apr 10 2015
CROSSREFS
KEYWORD
nonn,base,easy
AUTHOR
Derek Orr, Sep 19 2013
EXTENSIONS
More terms from Derek Orr, Apr 10 2015
STATUS
approved
