%I #22 Sep 08 2022 08:45:27
%S 1,0,12,8,108,160,960,2144,9040,24832,89664,270976,916416,2885120,
%T 9500160,30412288,99084544,319299584,1035979776,3347073024,
%U 10842246144,35064422400,113514577920,367253348352,1188632055808,3846143410176,12447083347968,40278203727872
%N Expansion of g.f.: 1/( (1+2*x)*(1-2*x-4*x^2)*(1-2*x^2)^2 ).
%C Based on characteristic polynomial of a square-within-a-square bonding graph.
%H G. C. Greubel, <a href="/A121961/b121961.txt">Table of n, a(n) for n = 0..1000</a>
%H <a href="/index/Rec#order_07">Index entries for linear recurrences with constant coefficients</a>, signature (0,12,8,-36,-32,32,32).
%F From _G. C. Greubel_, Oct 05 2019: (Start)
%F a(2*n) = 4^(n+1)*(Fibonacci(2*n-2) + 1) - 2^n*(3*n-1).
%F a(2*n+1) = 2^(2*n+3)*(Fibonacci(2*n-1) - 1) + 2^(n+2)*n. (End)
%p seq(coeff(series(1/((1-8*x^2-8*x^3)*(1-2*x^2)^2), x, n+1), x, n), n = 0 .. 40); # _G. C. Greubel_, Oct 05 2019
%t M = {{0,1,0,1,1,0,0,1}, {1,0,1,0,1,1,0,0}, {0,1,0,1,0,1,1,0}, {1,0,1,0,0, 0,1,1}, {1,1,0,0,0,0,0,0}, {0,1,1,0,0,0,0,0}, {0,0,1,1,0,0,0,0}, {1,0,0, 1,0,0,0,0}}; f[x_] = Det[M - x*IdentityMatrix[8]]; Table[ SeriesCoefficient[ Series[x/(x^10*f[1/x]), {x, 0, 30}], n], {n, 0, 30}]
%t LinearRecurrence[{0,12,8,-36,-32,32,32},{1,0,12,8,108,160,960},40] (* _Harvey P. Dale_, May 28 2017 *)
%t Table[If[EvenQ[n], 4^(n/2 +1)*(Fibonacci[n-2] +1) - 2^(n/2)*(3*n/2 -1), 2^(n+2)*(Fibonacci[n-2] -1) + 2^((n+1)/2)*(n-1)], {n, 0, 40}] (* _G. C. Greubel_, Oct 05 2019 *)
%o (PARI) my(x='x+O('x^40)); Vec(1/((1-8*x^2-8*x^3)*(1-2*x^2)^2)) \\ _G. C. Greubel_, Oct 05 2019
%o (Magma) R<x>:=PowerSeriesRing(Integers(), 40); Coefficients(R!( 1/((1-8*x^2-8*x^3)*(1-2*x^2)^2) )); // _G. C. Greubel_, Oct 05 2019
%o (Sage)
%o def A121961_list(prec):
%o P.<x> = PowerSeriesRing(ZZ, prec)
%o return P(1/((1-8*x^2-8*x^3)*(1-2*x^2)^2)).list()
%o A121961_list(30) # _G. C. Greubel_, Oct 05 2019
%o (GAP) a:=[1,0,12,8,108,160,960];; for n in [8..40] do a[n]:=12*a[n-2] +8*a[n-3]-36*a[n-4]-32*a[n-5]+32*a[n-6]+32*a[n-7]; od; a; # _G. C. Greubel_, Oct 05 2019
%Y Cf. A000045.
%K nonn,easy
%O 0,3
%A _Roger L. Bagula_, Sep 02 2006
%E Edited by _N. J. A. Sloane_, Feb 01 2007
|