%I #29 Jun 20 2022 08:46:34
%S 1,1,2,3,8,16,38,80,180,384,840,1792,3856,8192,17440,36864,77888,
%T 163840,344192,720896,1507584,3145728,6554112,13631488,28312576,
%U 58720256,121636864,251658240,520097792,1073741824,2214600704,4563402752,9395257344,19327352832
%N Numbers of isomers of unbranched a-4-catapolypentagons - see Brunvoll reference for precise definition.
%C Numbers of polypentagons with two connected internal vertices.
%H Colin Barker, <a href="/A121133/b121133.txt">Table of n, a(n) for n = 1..1000</a>
%H J. Brunvoll, S. J. Cyvin and B. N. Cyvin, <a href="http://dx.doi.org/10.1016/0166-1280(95)04463-9">Isomer enumeration of polygonal systems...</a>, J. Molec. Struct. (Theochem), 364 (1996), 1-13, Table 10.
%H S. J. Cyvin et al., <a href="http://dx.doi.org/10.1021/ci00013a027">Theory of polypentagons</a>, J. Chem. Inf. Comput. Sci., 33 (1993), 466-474.
%H <a href="/index/Rec#order_04">Index entries for linear recurrences with constant coefficients</a>, signature (4,-2,-8,8).
%F G.f.: x*(2*x^6+2*x^5-5*x^3+3*x-1)/((2*x-1)^2*(2*x^2-1)). - _Colin Barker_, Nov 30 2012
%F From _Colin Barker_, Oct 28 2016: (Start)
%F a(n) = 2^(n-5)*(n+2) for n>3 and even.
%F a(n) = 2^(n-5)*(n+2)+2^((n-5)/2) for n>3 and odd.
%F a(n) = 4*a(n-1)-2*a(n-2)-8*a(n-3)+8*a(n-4) for n>7.
%F (End)
%p H := proc(r,alpha,q) local rhalf,alphahalf ; rhalf := floor(r/2) ; alphahalf := floor(alpha/2) ; (binomial(rhalf-1,alphahalf-1)*(q-3)+binomial(rhalf-1,alphahalf))*(q-3)^(rhalf-alphahalf-1) ; end:
%p J := proc(r,alpha,q) (binomial(r-2,alpha-2)*(q-3)^2+2*binomial(r-2,alpha-1)*(q-3)+binomial(r-2,alpha))*(q-3)^(r-alpha-2) ; end: Ifunc := proc(r,alpha,q) J(r,alpha,q)/4+binomial(2,r-alpha)/4+ (1+(-1)^(r+alpha)+(1+(-1)^alpha)*(1-(-1)^r)/2)*H(r,alpha,q)/4 ; end:
%p A121133 := proc(n) if n = 1 then 1; else Ifunc(n,1,5) ; fi ; end: for n from 1 to 80 do printf("%d,",A121133(n)) ; od: # _R. J. Mathar_, Aug 07 2008
%t Rest@ CoefficientList[Series[x (2 x^6 + 2 x^5 - 5 x^3 + 3 x - 1)/((2 x - 1)^2*(2 x^2 - 1)), {x, 0, 34}], x] (* _Michael De Vlieger_, Oct 28 2016 *)
%t LinearRecurrence[{4,-2,-8,8},{1,1,2,3,8,16,38},40] (* _Harvey P. Dale_, Feb 10 2019 *)
%o (PARI) Vec(x*(2*x^6+2*x^5-5*x^3+3*x-1)/((2*x-1)^2*(2*x^2-1)) + O(x^50)) \\ _Colin Barker_, Oct 28 2016
%K nonn,easy
%O 1,3
%A _N. J. A. Sloane_, Aug 13 2006, Sep 24 2006
%E Edited and extended by _R. J. Mathar_, Aug 07 2008