This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.
%I A116091 #47 Sep 08 2022 08:45:24 %S A116091 1,-2,-2,28,-74,-92,1324,-3656,-4826,70228,-197372,-267896,3921724, %T A116091 -11126936,-15347432,225505648,-643622906,-897078476,13214495764, %U A116091 -37869162392,-53170602284,784672445368,-2255295815192,-3183829452272,47051201187676,-135537088268792,-192142210448216 %N A116091 Expansion of 1/sqrt(1+4*x+16*x^2). %C A116091 Fourth binomial transform is expansion of 1/sqrt(1-4*x+16*x^2), A012000. %H A116091 G. C. Greubel, <a href="/A116091/b116091.txt">Table of n, a(n) for n = 0..1000</a> %H A116091 Hacène Belbachir and Abdelghani Mehdaoui, <a href="https://doi.org/10.2989/16073606.2020.1729269">Recurrence relation associated with the sums of square binomial coefficients</a>, Quaestiones Mathematicae (2021) Vol. 44, Issue 5, 615-624. %H A116091 Hacène Belbachir, Abdelghani Mehdaoui, and László Szalay, <a href="https://cs.uwaterloo.ca/journals/JIS/VOL22/Szalay/szalay42.html">Diagonal Sums in the Pascal Pyramid, II: Applications</a>, J. Int. Seq., Vol. 22 (2019), Article 19.3.5. %F A116091 E.g.f.: exp(-2*x)*Bessel_I(0, 2*sqrt(-3)*x). %F A116091 a(n) = Sum_{k=0..n} C(n,k)^2*(-3)^k. %F A116091 O.g.f.: P(-1/2,4*x) with the o.g.f. P(x,z):=1/sqrt(1-2*x*z+z^2) for the Legendre polynomials. _Wolfdieter Lang_, Mar 10 2011. %F A116091 G.f. A(x) = 1/(2*T(0)-4*x-1) where T(k)= 1 + 3*x/(1 - x/T(k+1)); (continued fraction, 2-step ). - _Sergei N. Gladkovskii_, Aug 23 2012 %F A116091 D-finite with recurrence: a(n+2) = -(16*(n+1)*a(n))/(n+2) - (2*(2*n+3)*a(n+1))/(n+2) with a(0)=1, a(1)=-2. - _Alexander R. Povolotsky_, Aug 23 2012 %F A116091 a(n) = (-4)^n*hypergeom([-n, 1+n], [1], 1/4). - _Peter Luschny_, May 09 2016 %F A116091 From _Peter Bala_, Nov 30 2021: (Start) %F A116091 a(n) = (-4)^n^P(n,1/2), where P(n,x) is the n-th Legendre polynomial. %F A116091 a(n) = (4/3)*(16^n)*Sum_{k >= n} C(k,n)^2*(-1/3)^k. %F A116091 a(n) = (-3)^n*hypergeom([-n, -n], [1], -1/3). %F A116091 a(n) = (4/3)*(-16/3)^n*hypergeom([n+1, n+1], [1], -1/3). %F A116091 a(n) = [x^n] ((1 + x)*(3 - x))^n. (End) %p A116091 a := n -> (-4)^n*hypergeom([-n,1+n],[1],1/4); %p A116091 seq(simplify(a(n)), n=0..26); # _Peter Luschny_, May 09 2016 %t A116091 Table[4^n*LegendreP[n,-1/2],{n,0,30}] (* _Vaclav Kotesovec_, Jul 23 2013 *) %t A116091 CoefficientList[Series[1/Sqrt[1+4x+16x^2],{x,0,30}],x] (* _Harvey P. Dale_, Jun 08 2015 *) %o A116091 (PARI) Vec(1/sqrt(1+4*x+16*x^2+O(x^30))) \\ _M. F. Hasler_, Aug 25 2012 %o A116091 (Magma) R<x>:=PowerSeriesRing(Rationals(), 30); Coefficients(R!( 1/Sqrt(1+4*x+16*x^2) )); // _G. C. Greubel_, May 09 2019 %o A116091 (Sage) (1/sqrt(1+4*x+16*x^2)).series(x, 30).coefficients(x, sparse=False) # _G. C. Greubel_, May 09 2019 %Y A116091 Cf. A012000, A098332, A126869. %K A116091 sign,easy %O A116091 0,2 %A A116091 _Paul Barry_, Feb 04 2006