|
%I #36 Jun 07 2023 21:15:00
%S 1,8,32,96,256,640,1536,3584,8192,18432,40960,90112,196608,425984,
%T 917504,1966080,4194304,8912896,18874368,39845888,83886080,176160768,
%U 369098752,771751936,1610612736,3355443200,6979321856,14495514624,30064771072,62277025792
%N Expansion of (1 + 2*x)^2/(1 - 2*x)^2.
%H Vincenzo Librandi, <a href="/A241204/b241204.txt">Table of n, a(n) for n = 0..1000</a>
%H <a href="/index/Rec#order_02">Index entries for linear recurrences with constant coefficients</a>, signature (4,-4).
%F a(n) = 2^(2+n)*n for n>0. - _Colin Barker_, Apr 23 2014
%F a(n) = 4*a(n-1)-4*a(n-2) for n>2. - _Colin Barker_, Apr 23 2014
%F From _Amiram Eldar_, Jan 13 2021: (Start)
%F Sum_{n>=1} 1/a(n) = log(2)/4.
%F Sum_{n>=1} (-1)^(n+1)/a(n) = log(3/2)/4. (End)
%F E.g.f.: 1 + 8*x*exp(x). - _G. C. Greubel_, Jun 07 2023
%p A241204:= n->`if`(n=0, 1, 2^(n+2)*n); seq(A241204(n), n=0..20); # _Wesley Ivan Hurt_, Apr 22 2014
%t Table[2^(n+2)*n + Boole[n==0], {n,0,40}] (* _G. C. Greubel_, Jun 07 2023 *)
%o (Magma) R<x>:=PowerSeriesRing(Integers(), 41); Coefficients(R!((1+2*x)^2/(1-2*x)^2));
%o (PARI) Vec((2*x+1)^2/(2*x-1)^2 + O(x^100)) \\ _Colin Barker_, Apr 22 2014
%o (Sage)
%o def A241204(i):
%o if i==0: return 1
%o else: return 2^(2+i)*i;
%o [A241204(n) for n in (0..30)] # _Bruno Berselli_, Apr 23 2014
%Y Subsequence of A008574.
%K nonn,easy
%O 0,2
%A _Ilya Lopatin_ and _Juri-Stepan Gerasimov_, Apr 17 2014
|
