OFFSET
1,2
LINKS
Charles R Greathouse IV, Table of n, a(n) for n = 1..10000 (first 100 terms from Bajpai)
EXAMPLE
19683 is in the sequence because 19683 divided by its digital sum (1+9+6+8+3 = 27) gives 729 which is also a cube: 729 = 9^3.
46656 is in the sequence because 46656 divided by its digital sum (4+6+6+5+6 = 27) gives 1728 which is also a cube: 1728 = 12^3.
MAPLE
with(StringTools):KD := proc() local a, b, d, e; a:=n^3; b:=add( i, i = convert((a), base, 10))(a); d:=a/b; e:=evalf(d^(1/3)); if e=floor(e) then RETURN (a); fi; end: seq(KD(), n=1..200);
PROG
(PARI)
digsum(n) = d=eval(Vec(Str(n))); sum(i=1, #d, d[i])
s=[]; for(n=1, 200, d=digsum(n^3); if(n^3%d==0 && ispower(n^3\d, 3), s=concat(s, n^3))); s \\ Colin Barker, Jan 22 2014
CROSSREFS
KEYWORD
base,nonn
AUTHOR
K. D. Bajpai, Jan 19 2014
STATUS
approved