OFFSET
1,1
COMMENTS
See A246057 for the van der Poorten et al. reference and a comment.
153 is the Armstrong number A005188(10). [Typo corrected by Jeremy Tan, Feb 25 2023]
LINKS
Colin Barker, Table of n, a(n) for n = 1..333
Index entries for linear recurrences with constant coefficients, signature (1111,-112110,1111000,-1000000).
FORMULA
From Colin Barker, Feb 08 2017: (Start)
G.f.: 9*x*(17 - 550*x + 33500*x^2) / ((1 - x)*(1 - 10*x)*(1 - 100*x)*(1 - 1000*x)).
a(n) = (-2 + 2^(1+n)*5^n - 100^n + 1000^n) / 6.
a(n) = 1111*a(n-1) - 112110*a(n-2) + 1111000*a(n-3) - 1000000*a(n-4) for n>4. (End)
EXAMPLE
1^3 + 5^3 + 3^3 = 153, 16^3 + 50^3 + 33^3 = 165033, 166^3 + 500^3 + 333^3 = 166500333, ...
MATHEMATICA
Table[FromDigits@ Join[ReplacePart[ConstantArray[6, n], 1 -> 1], ReplacePart[ConstantArray[0, n], 1 -> 5], ConstantArray[3, n]], {n, 12}] (* Michael De Vlieger, Feb 08 2017 *)
PROG
(PARI) Vec(9*x*(17 - 550*x + 33500*x^2) / ((1 - x)*(1 - 10*x)*(1 - 100*x)*(1 - 1000*x)) + O(x^15)) \\ Colin Barker, Feb 08 2017
(PARI) a(n) = (((10^n - 4)/6)^3) + ((10^n/2)^3) + (((10^n - 1)/3)^3) \\ Jean-Jacques Vaudroz, Aug 11 2024
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Wolfdieter Lang, Feb 07 2017
STATUS
approved