OFFSET
1,2
COMMENTS
All terms are multiples of 3.
n is in the sequence iff 10*n is. Are there infinitely many terms not divisible by 10? - Robert Israel, Nov 20 2015
LINKS
Robert Israel, Table of n, a(n) for n = 1..320
EXAMPLE
126 is in the sequence because 126^3 = 2000376 and 3*126 = 378 have the same digit sum: 18.
MAPLE
select(n -> convert(convert(n^3, base, 10), `+`)=convert(convert(3*n, base, 10), `+`), 3*[$0..10^5]); # Robert Israel, Nov 20 2015
MATHEMATICA
Select[Range[0, 50000], Total[IntegerDigits[3 #]] == Total[IntegerDigits[#^3]] &]
PROG
(Magma) [n: n in [0..50000] | &+Intseq(3*n) eq &+Intseq(n^3)];
(PARI) for(n=0, 1e5, if(sumdigits(n^3)==sumdigits(3*n), print1(n, ", "))) \\ Altug Alkan, Nov 20 2015
CROSSREFS
KEYWORD
nonn,base,easy
AUTHOR
Vincenzo Librandi, Nov 18 2015
STATUS
approved