OFFSET
0,2
LINKS
Vincenzo Librandi, Table of n, a(n) for n = 0..10000
Index entries for linear recurrences with constant coefficients, signature (3,-3,1).
FORMULA
G.f.: (1+118*x+81*x^2)/(1-x)^3. - Bruno Berselli, Jul 30 2011
a(n) = a(n-1) + 40*(5*n-2), n > 0; a(0)=1. - Miquel Cerda, Oct 30 2016
a(n) = A017281(n)^2. - Michel Marcus, Oct 30 2016
E.g.f.: (1 +120*x +100*x^2)*exp(x). - G. C. Greubel, Dec 24 2022
MATHEMATICA
(* Programs from Michael De Vlieger, Mar 30 2017 *)
Table[(10 n+1)^2, {n, 0, 35}]
FoldList[#1 + 200 #2 - 80 &, 1, Range@ 35]
CoefficientList[Series[(1+118x+81x^2)/(1-x)^3, {x, 0, 35}], x] (* End *)
LinearRecurrence[{3, -3, 1}, {1, 121, 441}, 40] (* Harvey P. Dale, Sep 21 2017 *)
PROG
(Magma) [(10*n+1)^2: n in [0..35]]; // Vincenzo Librandi, Jul 30 2011
(PARI) for(n=0, 35, print1((10*n+1)^2", ")); \\ Bruno Berselli, Jul 30 2011
(SageMath) [(10*n+1)^2 for n in range(51)] # G. C. Greubel, Dec 24 2022
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
EXTENSIONS
More terms from Bruno Berselli, Jul 30 2011
STATUS
approved