OFFSET
0,2
COMMENTS
LINKS
Colin Barker, Table of n, a(n) for n = 0..500
Editors, L'Intermédiaire des Mathématiciens, Query 4500: The equation x(x+1)/2 = y*(y+1)/3, L'Intermédiaire des Mathématiciens, 22 (1915), 255-260 (I).
Editors, L'Intermédiaire des Mathématiciens, Query 4500: The equation x(x+1)/2 = y*(y+1)/3, L'Intermédiaire des Mathématiciens, 22 (1915), 255-260 (II).
Editors, L'Intermédiaire des Mathématiciens, Query 4500: The equation x(x+1)/2 = y*(y+1)/3, L'Intermédiaire des Mathématiciens, 22 (1915), 255-260 (III).
Editors, L'Intermédiaire des Mathématiciens, Query 4500: The equation x(x+1)/2 = y*(y+1)/3, L'Intermédiaire des Mathématiciens, 22 (1915), 255-260 (IV).
Index entries for linear recurrences with constant coefficients, signature (99,-99,1).
FORMULA
G.f. 10*x / ((1-x)*(x^2-98*x+1)). - R. J. Mathar, Dec 20 2011
a(n) = 99*a(n-1)-99*a(n-2)+a(n-3) for n>2. - Colin Barker, Mar 02 2016
a(n) = (-10+(5-2*sqrt(6))*(49+20*sqrt(6))^(-n)+(5+2*sqrt(6))*(49+20*sqrt(6))^n)/96. - Colin Barker, Mar 07 2016
EXAMPLE
3*0 = 2*0.
3*10 = 2*15.
3*990 = 2*1485.
3*97020 = 2*145530.
MATHEMATICA
LinearRecurrence[{99, -99, 1}, {0, 10, 990}, 20] (* Harvey P. Dale, Feb 25 2018 *)
PROG
(PARI) concat(0, Vec(10*x/((1-x)*(1-98*x+x^2)) + O(x^40))) \\ Colin Barker, Mar 02 2016
(Magma) m:=30; R<x>:=PowerSeriesRing(Integers(), m); [0] cat Coefficients(R!(10*x/((1-x)*(x^2-98*x+1)))); // G. C. Greubel, Jul 15 2018
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Charlie Marion, Dec 15 2011
STATUS
approved