OFFSET
0,3
LINKS
G. C. Greubel, Table of n, a(n) for n = 0..1000
Marco Abrate, Stefano Barbero, Umberto Cerruti, and Nadir Murru, Polynomial sequences on quadratic curves, Integers, Vol. 15, 2015, #A38.
S. Barbero, U. Cerruti, and N. Murru, On polynomial solutions of the Diophantine equation (x + y - 1)^2 = wxy, Rendiconti Sem. Mat. Univ. Pol. Torino (2020) Vol. 78, No. 1, 5-12.
Index entries for linear recurrences with constant coefficients, signature (10,-10,1).
FORMULA
a(n) = 2*(T(n, 9/2)-1)/7 with twice Chebyshev's polynomials of the first kind evaluated at x=9/2: 2*T(n, 9/2) = A056918(n) = ((9 + sqrt(77))^n + (9 - sqrt(77))^n)/2^n.
a(n) = 9*a(n-1) - a(n-2) + 2, n >= 2, a(0)=0, a(1)=1.
a(n) = 10*a(n-1) - 10*a(n-2) + a(n-3), n >= 3, a(0)=0, a(1)=1, a(2)=11.
G.f.: x*(1+x)/((1-x)*(1-9*x+x^2)) = x*(1+x)/(1-10*x+10*x^2-x^3) (from the Stephan link, see A092184).
MATHEMATICA
LinearRecurrence[{10, -10, 1}, {0, 1, 11}, 30] (* Harvey P. Dale, Jan 27 2012 *)
PROG
(PARI) my(x='x+O('x^30)); concat([0], Vec(x*(1+x)/((1-x)*(1-9*x+x^2)))) \\ G. C. Greubel, May 24 2019
(Magma) I:=[0, 1, 11]; [n le 3 select I[n] else 10*Self(n-1)-10*Self(n-2) + Self(n-3): n in [1..30]]; // G. C. Greubel, May 24 2019
(Sage) (x*(1+x)/((1-x)*(1-9*x+x^2))).series(x, 30).coefficients(x, sparse=False) # G. C. Greubel, May 24 2019
(GAP) a:=[0, 1, 11];; for n in [4..30] do a[n]:=10*a[n-1]-10*a[n-2]+ a[n-3]; od; a; # G. C. Greubel, May 24 2019
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Wolfdieter Lang, Oct 18 2004
STATUS
approved