OFFSET
0,2
COMMENTS
For n>0, numbers such that sqrt(a(n)) has the continued fraction {k;[1,1,1,2k]}, where the part in [] is repeated and k is of the form 3m+2 (A016789). - Bruno Berselli, May 30 2013
For n >= 1, the continued fraction expansion of sqrt(4*a(n)) is [6n-1; {3, 3n-1, 3, 12n-2}]. - Magus K. Chu, Sep 18 2022
LINKS
G. C. Greubel, Table of n, a(n) for n = 0..5000
Index entries for linear recurrences with constant coefficients, signature (3,-3,1).
FORMULA
a(n) = a(n-1) + 18*n - 11 with n>0, a(0)=0. - Vincenzo Librandi, Nov 22 2010
a(0)=0, a(1)=7, a(2)=32, a(n) = 3*a(n-1) - 3*a(n-2) + a(n-3). - Harvey P. Dale, Jul 07 2012
From G. C. Greubel, Mar 11 2022: (Start)
G.f.: x*(7 - 11*x)/(1-x)^3.
E.g.f.: x*(7 + 9*x)*exp(x). (End)
Sum_{n>=1} 1/a(n) = -(psi(7/9)+gamma)/2 = (A354640-A001620)/2 = 0.22000753... - R. J. Mathar, Apr 22 2024
MATHEMATICA
Table[n(9n-2), {n, 0, 50}] (* or *) LinearRecurrence[{3, -3, 1}, {0, 7, 32}, 50] (* Harvey P. Dale, Jul 07 2012 *)
PROG
(PARI) a(n)=n*(9*n-2) \\ Charles R Greathouse IV, Jun 17 2017
(Magma) [n*(9*n-2): n in [0..60]]; // G. C. Greubel, Mar 11 2022
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
STATUS
approved