login

Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 61st year, we have over 378,000 sequences, and we’ve reached 11,000 citations (which often say “discovered thanks to the OEIS”).

A051406
a(n) = (3^n+1) * (3^(n+1)+1) / 8.
2
1, 5, 35, 287, 2501, 22265, 199655, 1794707, 16145801, 145292525, 1307573675, 11767985927, 105911341901, 953200482785, 8578799562095, 77209181709947, 694882592342801, 6253943201945045, 56285488430084915, 506569394708502767, 4559124548889740501
OFFSET
0,2
FORMULA
a(n) = 13*a(n-1) - 39*a(n-2) + 27*a(n-3). G.f.: -(1-8*x+9*x^2)/((x-1)*(3*x-1)*(9*x-1)). [R. J. Mathar, Jul 25 2009]
PROG
(PARI) Vec((1 - 8*x + 9*x^2) / ((1 - x)*(1 - 3*x)*(1 - 9*x)) + O(x^30)) \\ Colin Barker, Jul 16 2017
CROSSREFS
Sequence in context: A369157 A002294 A369128 * A000356 A370844 A379103
KEYWORD
nonn,easy
AUTHOR
STATUS
approved