%I #65 Oct 25 2022 06:32:14
%S 1,-3,25,-243,2601,-29403,344569,-4141875,50737129,-630663003,
%T 7930793025,-100681224075,1288236350025,-16592960274075,
%U 214939203248025,-2797935722568243,36578032462268649,-480000660000226875,6320012816203363489,-83462977778600141643,1105193229806740453201
%N G.f.: Sum_{n>=0} x^n*Sum_{k=0..n} (-1)^k * C(n,k)^2 * Sum_{j=0..k} C(k,j)^2 * x^j.
%C The g.f.s formed from a(2*n)^(1/2) and (-a(2*n+1)/3)^(1/2) are:
%C A245926: sqrt( (1-x + sqrt(1-14*x+x^2)) / (2*(1-14*x+x^2)) );
%C A245927: sqrt( (1-x - sqrt(1-14*x+x^2)) / (6*x*(1-14*x+x^2)) ).
%C Lim_{n->infinity} a(n+1)/a(n) = -(7 + 4*sqrt(3)).
%H Seiichi Manyama, <a href="/A245925/b245925.txt">Table of n, a(n) for n = 0..500</a>
%H A. Bostan, S. Boukraa, J.-M. Maillard, and J.-A. Weil, <a href="http://arxiv.org/abs/1507.03227">Diagonals of rational functions and selected differential Galois groups</a>, arXiv preprint arXiv:1507.03227 [math-ph], 2015.
%F G.f.: Sum_{n>=0} x^n / (1+x)^(2*n+1) * ( Sum_{k=0..n} C(n,k)^2*(-x)^k )^2.
%F G.f.: 1 / AGM(1-x^2, sqrt(1+14*x^2+x^4)), where AGM(x,y) = AGM((x+y)/2, sqrt(x*y)) is the arithmetic-geometric mean.
%F a(2*n) = A245926(n)^2.
%F a(2*n+1) = (-3)*A245927(n)^2.
%F a(n) = Sum_{k=0..n} Sum_{j=0..2*n-2*k} (-1)^(j+k) * C(2*n-k,j+k)^2 * C(j+k,k)^2.
%F D-finite with recurrence: n^2*(2*n-3)*a(n) = -(2*n-1)*(13*n^2 - 26*n + 10)*a(n-1) + (2*n-3)*(13*n^2 - 26*n + 10)*a(n-2) + (n-2)^2*(2*n-1)*a(n-3). - _Vaclav Kotesovec_, Aug 16 2014
%F a(n) ~ (-1)^n * (2+sqrt(3)) * (7+4*sqrt(3))^n / (4*Pi*n). - _Vaclav Kotesovec_, Aug 16 2014
%F a(n) = (-1)^n*Sum_{k=0..n} binomial(2*(n-k), n-k)*binomial(2*n-k, k)^2. - _Peter Luschny_, Aug 17 2014
%F a(n) = (-1)^n*binomial(2*n,n)*hyper4F3([-n,-n,-n,-n+1/2],[1,-2*n,-2*n], 4). - _Peter Luschny_, Aug 17 2014
%F a(n) = Sum_{k=0..n} (-1)^k * C(2*k, k)^2 * C(n+k, n-k). - _Paul D. Hanna_, Aug 17 2014
%F a(n) = (-1)^n*hypergeom([-n, -n, n + 1, n + 1], [1/2, 1, 1], 1/4). - _Peter Luschny_, Mar 14 2018
%F a(n) = Legendre_P(n, sqrt(-3))^2. - _Peter Bala_, Dec 22 2020
%F G.f.: Sum_{n >= 0} (-1)^n*binomial(2*n,n)^2*x^n/(1-x)^(2*n+1). - _Peter Bala_, Feb 07 2022
%F From _Peter Bala_, Apr 05 2022: (Start)
%F a(n) = Sum_{k = 0..n} (-1)^(n+k)*binomial(2*n-k,k)*binomial(2*n-2*k,n-k)^2.
%F a(n) = (-1)^n * Sum_{k = 0..n} binomial(2*k,k)*binomial(n+k,n-k)^2.
%F a(n) = (-1)^n*binomial(2*n,n)^2*hypergeom([-n,-n,-n,], [-2*n,-n+1/2], 1/4). (End)
%e G.f.: A(x) = 1 - 3*x^2 + 25*x^4 - 243*x^6 + 2601*x^8 - 29403*x^10 + ...
%e where the g.f. is given by the binomial series:
%e A(x) = 1 + x*(1 - (1+x)) + x^2*(1 - 2^2*(1+x) + (1+2^2*x+x^2))
%e + x^3*(1 - 3^2*(1+x) + 3^2*(1+2^2*x+x^2) - (1+3^2*x+3^2*x^2+x^3))
%e + x^4*(1 - 4^2*(1+x) + 6^2*(1+2^2*x+x^2) - 4^2*(1+3^2*x+3^2*x^2+x^3) + (1+4^2*x+6^2*x^2+4^2*x^3+x^4))
%e + x^5*(1 - 5^2*(1+x) + 10^2*(1+2^2*x+x^2) - 10^2*(1+3^2*x+3^2*x^2+x^3) + 5^2*(1+4^2*x+6^2*x^2+4^2*x^3+x^4) - (1+5^2*x+10^2*x^2+10^2*x^3+5^2*x^4+x^5))
%e + x^6*(1 - 6^2*(1+x) + 15^2*(1+2^2*x+x^2) - 20^2*(1+3^2*x+3^2*x^2+x^3) + 15^2*(1+4^2*x+6^2*x^2+4^2*x^3+x^4) - 6^2*(1+5^2*x+10^2*x^2+10^2*x^3+5^2*x^4+x^5) + (1+6^2*x+15^2*x^2+20^2*x^3+15^2*x^4+6^2*x^5+x^6)) + ...
%e in which the coefficients of odd powers of x vanish.
%e We can also express the g.f. by the binomial series identity:
%e A(x) = 1/(1+x) + x/(1+x)^3*(1-x)^2 + x^2/(1+x)^5*(1 - 2^2*x + x^2)^2
%e + x^3/(1+x)^7*(1 - 3^2*x + 3^2*x^2 - x^3)^2
%e + x^4/(1+x)^9*(1 - 4^2*x + 6^2*x^2 - 4^2*x^3 + x^4)^2
%e + x^5/(1+x)^11*(1 - 5^2*x + 10^2*x^2 - 10^2*x^3 + 5^2*x^4 - x^5)^2
%e + x^6/(1+x)^13*(1 - 6^2*x + 15^2*x^2 - 20^2*x^3 + 15^2*x^4 - 6^2*x^5 + x^6)^2 + ...
%p A245925 := n -> (-1)^n*add(binomial(2*(n-k), n-k)*binomial(2*n-k, k)^2, k=0..n); seq(A245925(n), n=0..20); # _Peter Luschny_, Aug 17 2014
%t Table[Sum[Sum[(-1)^(j+k) * Binomial[2*n - k, j + k]^2 * Binomial[j + k, k]^2, {j, 0, 2*n - 2*k}], {k, 0, n}], {n, 0, 20}] (* _Vaclav Kotesovec_, Aug 16 2014 after _Paul D. Hanna_ *)
%t a[n_] := (-1)^n*HypergeometricPFQ[{-n, -n, n + 1, n + 1}, {1/2, 1, 1}, 1/4];
%t Table[a[n], {n, 0, 20}] (* _Peter Luschny_, Mar 14 2018 *)
%o (PARI) /* By definition: */
%o {a(n)=polcoeff(sum(m=0, n, x^m*sum(k=0, m, (-1)^k*binomial(m, k)^2*sum(j=0, k, binomial(k, j)^2*x^j)+x*O(x^n))), n)}
%o for(n=0, 20, print1(a(2*n), ", "))
%o (PARI) /* From alternate g.f.: */
%o {a(n)=local(A=1);A=sum(m=0,n,x^m/(1+x)^(2*m+1)*sum(k=0,m,binomial(m,k)^2*(-x)^k)^2+x*O(x^n));polcoeff(A,n)}
%o for(n=0,20,print1(a(2*n),", "))
%o (PARI) /* From formula for a(n); printing only nonzero terms: */
%o {a(n)=sum(k=0, n\2, sum(j=0, n-2*k, (-1)^(j+k)*binomial(n-k, j+k)^2*binomial(j+k, k)^2))}
%o for(n=0, 20, print1(a(2*n), ", "))
%o (PARI) /* From formula for a(n) (nonzero terms): */
%o {a(n)=sum(k=0, n, sum(j=0, 2*n-2*k, (-1)^(j+k)*binomial(2*n-k,j+k)^2*binomial(j+k, k)^2))}
%o for(n=0, 20, print1(a(n), ", "))
%o (PARI) /* Formula for a(n), after _Peter Luschny_ and _Robert Israel_: */
%o {a(n) = (-1)^n * sum(k=0,n, binomial(2*k, k) * binomial(n+k, n-k)^2)}
%o for(n=0, 20, print1(a(n), ", "))
%o (PARI) /* Simpler formula for a(n): */
%o {a(n) = sum(k=0, n, (-1)^k * binomial(2*k, k)^2 * binomial(n+k, n-k) )}
%o for(n=0, 20, print1(a(n), ", "))
%o (PARI) /* Using AGM: */
%o {a(n)=polcoeff( 1 / agm(1-x, sqrt(1+14*x+x^2 +x*O(x^n))), n)}
%o for(n=0,20,print1(a(n),", "))
%o (Sage)
%o A245925 = lambda n: (-1)^n*sum(binomial(2*(n-k), n-k)*binomial(2*n-k, k)^2 for k in (0..n))
%o [A245925(n) for n in range(21)] # _Peter Luschny_, Aug 17 2014
%Y Cf. A245926, A245927, A245929, A243945, A249891.
%K sign,easy
%O 0,2
%A _Paul D. Hanna_, Aug 15 2014