OFFSET
1,1
LINKS
Vincenzo Librandi, Table of n, a(n) for n = 1..200
Index entries for linear recurrences with constant coefficients, signature (33711,-33711,1).
FORMULA
a(n+2) = 33710*a(n+1)-a(n)+724722.
G.f.: 43*x*(27*x^2-16855*x-26) / ((x-1)*(x^2-33710*x+1)). - Colin Barker, Oct 21 2014
EXAMPLE
a(1)=1118 because the first relation is : (1118+43)^3-1118^3=12943^2.
MATHEMATICA
CoefficientList[Series[43 (27 x^2 - 16855 x - 26)/((x - 1) (x^2 - 33710 x + 1)), {x, 0, 20}], x] (* Vincenzo Librandi, Oct 21 2014 *)
LinearRecurrence[{33711, -33711, 1}, {1118, 38413663, 1294925303334}, 20] (* Harvey P. Dale, Feb 22 2017 *)
PROG
(PARI) Vec(43*x*(27*x^2-16855*x-26)/((x-1)*(x^2-33710*x+1)) + O(x^20)) \\ Colin Barker, Oct 21 2014
(Magma) I:=[1118, 38413663, 1294925303334]; [n le 3 select I[n] else 33711*Self(n-1)-33711*Self(n-2)+Self(n-3): n in [1..10]]; // Vincenzo Librandi, Oct 21 2014
CROSSREFS
KEYWORD
easy,nonn
AUTHOR
Richard Choulet, Oct 08 2008
EXTENSIONS
Editing from Colin Barker, Oct 21 2014
STATUS
approved
