A116423 Binomial transform of A006053.
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
Keywords
Examples
a(5) = 26 = 1*0 + 1*4 + 4*1 + 4*3 + 6*1 = 4 + 4 + 12 + 6 = 26.
Links
- Jinyuan Wang, Table of n, a(n) for n = 1..1000
- Andrew B. Hudson, Illustration of the first 7 terms as a spiral tiling of similar trapezoids.
- Peter Steinbach, Golden Fields: A Case for the Heptagon, Mathematics Magazine, Vol. 70, No. 1, Feb. 1997.
- Index entries for linear recurrences with constant coefficients, signature (4,-3,-1)
Crossrefs
Cf. A006053.
Programs
-
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
-
Mathematica
LinearRecurrence[{4, -3, -1}, {0, 1, 3}, 40] (* Vincenzo Librandi, Jul 11 2019 *)
-
PARI
concat(0, Vec(x^2*(1-x)/(1-4*x+3*x^2+x^3) + O(x^50))) \\ Michel Marcus, Sep 13 2014
Formula
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)
Extensions
More terms from R. J. Mathar, Apr 02 2008
More terms from Michel Marcus, Sep 13 2014
Comments