OFFSET
1,3
LINKS
Derek Orr, Table of n, a(n) for n = 1..461
EXAMPLE
512*(5+1+2) = 4096 = 16^3. Hence, 512 is a term of the sequence.
MATHEMATICA
Select[Range@ 1000000, IntegerQ@ Power[# Plus @@ IntegerDigits@ #, 1/3] == True &] (* Michael De Vlieger, Mar 23 2015 *)
PROG
(Python)
def DS(n):
return sum(int(i) for i in str(n))
def a(n):
k = 0
nDSn = n * DS(n)
while k <= n:
if k**3 == nDSn:
return True
if k**3 > nDSn:
return False
k += 1
[n for n in range(10**5) if a(n)]
# Simplified by Derek Orr, Mar 22 2015
(Sage)
n=100000 # change n for more terms
[x for x in [0..n] if floor((x*sum(Integer(x).digits(base=10)))^(1/3))==(x*sum(Integer(x).digits(base=10)))^(1/3)] # Tom Edgar, Sep 21 2013
(PARI) for(n=0, 10^6, if((n==0) || ispower(n*sumdigits(n), 3), print1(n, ", "))) \\ Derek Orr, Mar 22 2015
CROSSREFS
KEYWORD
nonn,base,easy
AUTHOR
Derek Orr, Sep 19 2013
STATUS
approved