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

A116423
Binomial transform of A006053.
6
0, 1, 3, 9, 26, 74, 209, 588, 1651, 4631, 12983, 36388, 101972, 285741, 800660, 2243445, 6286059, 17613241, 49351342, 138279586, 387451077, 1085614208, 3041824015, 8523002359, 23880923183, 66912861640, 187485674652, 525323190505, 1471922876424, 4124236259529
OFFSET
1,3
COMMENTS
a(n)/a(n-1) tends to 2.801... = 1 + 2*cos(Pi/7).
A(n) := a(n+1)*(-1)^(n+1) appears in the following formula for the nonpositive powers of rho*sigma, where rho:=2*cos(Pi/7) and sigma:=sin(3*Pi/7)/sin(Pi/7) = rho^2-1 are the ratios of the smaller and larger diagonal length to the side length in a regular 7-gon (heptagon). See the Steinbach reference where the basis <1,rho,sigma> is used in an extension of the rational field. (rho*sigma)^(-n) = C(n) + B(n)*rho + A(n)*sigma, n >= 0, with C(n)= A085810(n)*(-1)^n, and B(n)= A181880(n-2)*(-1)^n. For the nonnegative powers see A120757(n), |A122600(n-1)| and A181879(n), respectively. See also a comment under A052547.
This sequence is constructible as a spiral tiling of similar trapezoids, as follows: start with an isosceles trapezoid with side lengths 3,1,4,1. Each new trapezoid is rotated and scaled so one leg fills all unoccupied space on the short base of the previous trapezoid. a(n) is given by the length of the n-th trapezoid's legs. This process is identical to the recursion relation added by R. J. Mathar in the Formula section. See the Links section for an illustration. - Andrew B. Hudson, Jun 19 2019
FORMULA
Binomial transform of A006053 starting with A006053(1): (0, 1, 1, 3, 4, 9, 14, ...).
From R. J. Mathar, Apr 02 2008: (Start)
O.g.f.: x^2(1-x)/(1 - 4x + 3x^2 + x^3).
a(n) = 4*a(n-1) - 3*a(n-2) - a(n-3). (End)
EXAMPLE
a(5) = 26 = 1*0 + 1*4 + 4*1 + 4*3 + 6*1 = 4 + 4 + 12 + 6 = 26.
MATHEMATICA
LinearRecurrence[{4, -3, -1}, {0, 1, 3}, 40] (* Vincenzo Librandi, Jul 11 2019 *)
PROG
(PARI) concat(0, Vec(x^2*(1-x)/(1-4*x+3*x^2+x^3) + O(x^50))) \\ Michel Marcus, Sep 13 2014
(Magma) I:=[0, 1, 3]; [n le 3 select I[n] else 4*Self(n-1)-3*Self(n-2)-Self(n-3): n in [1..30]]; // Vincenzo Librandi, Jul 11 2019
CROSSREFS
Cf. A006053.
Sequence in context: A258911 A268093 A127911 * A077845 A291000 A276068
KEYWORD
nonn
AUTHOR
Gary W. Adamson, Feb 14 2006
EXTENSIONS
More terms from R. J. Mathar, Apr 02 2008
More terms from Michel Marcus, Sep 13 2014
STATUS
approved