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”).

A037480
Base-3 digits are, in order, the first n terms of the periodic sequence with initial period 1,2.
3
1, 5, 16, 50, 151, 455, 1366, 4100, 12301, 36905, 110716, 332150, 996451, 2989355, 8968066, 26904200, 80712601, 242137805, 726413416, 2179240250, 6537720751, 19613162255, 58839486766, 176518460300, 529555380901
OFFSET
1,2
FORMULA
a(n) = (5*3^n + (-1)^n - 6)/8 . - Paul D. Hanna, Sep 23 2007
a(n) = 3*a(n-1)+a(n-2)-3*a(n-3). G.f.: x*(2*x+1) / ((x-1)*(x+1)*(3*x-1)). - Colin Barker, Apr 30 2014
MAPLE
a[0]:=0:a[1]:=1:for n from 2 to 50 do a[n]:=2*a[n-1]+3*a[n-2]+3 od: seq(a[n], n=1..33); # Zerinvary Lajos, Dec 14 2008
MATHEMATICA
CoefficientList[Series[(2 x + 1)/((x - 1) (x + 1) (3 x - 1)), {x, 0, 40}], x] (* Vincenzo Librandi, May 01 2014 *)
Table[FromDigits[PadRight[{}, n, {1, 2}], 3], {n, 30}] (* or *) LinearRecurrence[ {3, 1, -3}, {1, 5, 16}, 30] (* Harvey P. Dale, Dec 15 2019 *)
PROG
(PARI) Vec(x*(2*x+1)/((x-1)*(x+1)*(3*x-1)) + O(x^100)) \\ Colin Barker, Apr 30 2014
CROSSREFS
Sequence in context: A082001 A084356 A007806 * A027108 A304168 A317817
KEYWORD
nonn,base,easy
STATUS
approved