OFFSET
1,1
LINKS
Paolo Xausa, Table of n, a(n) for n = 1..250
Jamal Agbanwa, Generating Infinitely Many Taxicab Numbers, Figshare, (preprint) 2025.
Index entries for linear recurrences with constant coefficients, signature (2716,-37830,2716,-1).
FORMULA
a(n) = x(n-1)^3 + y(n-1)^3 = u(n-1)^3 + w(n-1)^3 where:
x(n) = (-54 + 9*((15 - 7*sqrt(3))*(7 - 4*sqrt(3))^n + (15 + 7*sqrt(3))*(7 + 4*sqrt(3))^n))/4 + 27,
y(n) = (-162 + 9*((7 - 5*sqrt(3))*(7 - 4*sqrt(3))^n + (7 + 5*sqrt(3))*(7 + 4*sqrt(3))^n))/4,
u(n) = (-54 + 9*((15 - 7*sqrt(3))*(7 - 4*sqrt(3))^n + (15 + 7*sqrt(3))*(7 + 4*sqrt(3))^n))/4, and
w(n) = (-162 + 9*((7 - 5*sqrt(3))*(7 - 4*sqrt(3))^n + (7 + 5*sqrt(3))*(7 + 4*sqrt(3))^n))/4 + 81.
G.f.: 530712*x*(1387 - 38121*x + 2715*x^2 - x^3)/((1 - 2702*x + x^2)*(1 - 14*x + x^2)). - Andrew Howroyd, Feb 25 2026
EXAMPLE
For n = 8, a(8) = x(7)^3 + y(7)^3 = ((-54 + 9*((15 - 7*sqrt(3))*(7 - 4*sqrt(3))^7 + (15 + 7*sqrt(3))*(7 + 4*sqrt(3))^7))/4 + 27)^3 + ((-162 + 9((7 - 5*sqrt(3))*(7 - 4*sqrt(3))^7 + (7 + 5*sqrt(3))*(7 + 4*sqrt(3))^7))/4)^3 = 285011938432925772595065675864.
MATHEMATICA
LinearRecurrence[{2716, -37830, 2716, -1}, {530712, 736097544, 1979009657352, 5347145100161592}, 10] (* Paolo Xausa, Mar 12 2026 *)
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Jamal Agbanwa, Feb 25 2026
STATUS
approved
