%I #68 Dec 04 2021 12:25:07
%S 1,2,-2,-28,-74,92,1324,3656,-4826,-70228,-197372,267896,3921724,
%T 11126936,-15347432,-225505648,-643622906,897078476,13214495764,
%U 37869162392,-53170602284,-784672445368,-2255295815192,3183829452272,47051201187676,135537088268792,-192142210448216
%N Expansion of 1/sqrt(1 - 4*x + 16*x^2).
%H T. D. Noe, <a href="/A012000/b012000.txt">Table of n, a(n) for n = 0..200</a>
%H 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 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.
%H Tony D. Noe, <a href="http://www.cs.uwaterloo.ca/journals/JIS/VOL9/Noe/noe35.html">On the Divisibility of Generalized Central Trinomial Coefficients</a>, Journal of Integer Sequences, Vol. 9 (2006), Article 06.2.7.
%F G.f.: 1/sqrt(1 - 4*x + 16*x^2). - _Vladeta Jovovic_, May 13 2003
%F Scaled Legendre polynomials evaluated at 1/2: 2^(2n)P(n, 1/2). - _Michael Somos_, Dec 03 2001
%F a(n) = (-1)^n*Sum_{k=0..n} C(n,k)^2*(-3)^k. - _Benoit Cloitre_, Oct 25 2003
%F a(n) = Sum_{k=0..floor(n/2)} C(n,k)*C(2(n-k),n)(-4)^k. - _Paul Barry_, Sep 08 2004
%F D-finite with recurrence: n*a(n) + 2*(1-2*n)*a(n-1) + 16*(n-1)*a(n-2) = 0. - _R. J. Mathar_, Nov 14 2011
%F G.f.: G(0), where G(k) = 1 + x*(1-x)*(4*k+1)/( 4*k+2 - x*(1-x)*(4*k+2)*(4*k+3)/(x*(1-x)*(4*k+3) + 4*(k+1)/G(k+1) )); (continued fraction). - _Sergei N. Gladkovskii_, Jul 06 2013
%F a(n) ~ 2^(2*n+1)*sin(Pi/12*(4*n+5))/(3^(1/4)*sqrt(Pi*n)). - _Vaclav Kotesovec_, Jul 30 2013
%F a(n) = [x^n] ((1 - x)*(1 + 3*x))^n. - _Peter Bala_, Aug 16 2015
%F a(n) = 4^n*hypergeom([-n,1+n], [1], 1/4). - _Peter Luschny_, May 09 2016
%F From _Peter Bala_, Nov 28 2021: (Start)
%F a(n) = 4^(2*n+1)*Sum_{k >= n} (-1)^(n+k)*C(k,n)^2*(1/3)^(k+1).
%F a(n) = (-1)^n*hypergeom([-n,-n], [1], -3).
%F a(n) = (4/3)*(16/3)^n*hypergeom([n+1,n+1], [1], -1/3). (End)
%e G.f. = 1 + 2*x - 2*x^2 - 28*x^3 - 74*x^4 + 92*x^5 + 1324*x^6 + 3656*x^7 + ...
%p a := n -> 4^n*hypergeom([-n,1+n],[1],1/4);
%p seq(simplify(a(n)),n=0..26); # _Peter Luschny_, May 09 2016
%t Table[ 2^(2n) LegendreP[ n, 1/2 ], {n, 12} ]
%o (PARI) {a(n) = 2^(2*n) * subst( pollegendre(n), x, 1/2)} /* _Michael Somos_, Dec 03 2001 */
%o (PARI) a(n) = polcoeff(((1 - x)*(1 + 3*x))^n, n); \\ _Michel Marcus_, Aug 16 2015
%Y Cf. A116091, A098332, A126869.
%K sign
%O 0,2
%A _Wouter Meeussen_
%E More terms from _Vladeta Jovovic_, May 13 2003