|
%I #14 Apr 28 2018 18:14:01
%S 1,1,1,1,2,6,16,36,75,164,401,1046,2718,6878,17200,43486,112202,
%T 293540,770535,2019891,5296670,13942944,36902130,98097968,261456388,
%U 697970447,1866383507,5001333169,13432923544,36154294520,97475330092,263188299372
%N G.f. satisfies: A(x) = Sum_{n>=0} x^(n^2) * A(x)^(n^2).
%C Apparently: Number of Dyck n-paths with each ascent length being a square number. [_David Scambler_, May 09 2012]
%F G.f. satisfies: A(x) = B(x*A(x)) where B(x) = A(x/B(x)) = Sum_{n>=0} x^(n^2),
%F where a(n) = [x^n] B(x)^(n+1)/(n+1) for n>=0.
%F G.f.: A(x) = (1/x)*Series_Reversion( x / Sum_{n>=0} x^(n^2) ).
%F From _Paul D. Hanna_, Apr 24 2010: (Start)
%F SPECIAL VALUES:
%F . at x = 2*exp(-Pi)/(1+Pi^(1/4)/gamma(3/4)) = 0.04142369369176926261...
%F . A(x) = B(exp(-Pi)) = (1+Pi^(1/4)/gamma(3/4))/2 = 1.043217405606654...
%F RADIUS OF CONVERGENCE r:
%F . at r = 0.3529672118496605771445592553666318566205464502456806...,
%F . A(r) = 1.9530374869760035836323161721583051467541841357702661...,
%F where r and A(r) are given by:
%F . r = z/B(z) and
%F . A(r) = B(z) = Sum_{n>=0} z^(n^2)
%F such that z is the real root nearest the origin that satisfies:
%F . B(z) - z*B'(z) = 0, which has solution:
%F . z = 0.689358196415787767209694723600383373645983284157633311584643...
%F Here, B(z) = Sum_{n>=0} z^(n^2), the partial Jacobi theta_3 function.
%F (End)
%e G.f.: A(x) = 1 + x + x^2 + x^3 + 2*x^4 + 6*x^5 + 16*x^6 + 36*x^7 +...
%e A(x)^4 = 1 + 4*x + 10*x^2 + 20*x^3 + 39*x^4 + 88*x^5 + 228*x^6 +...
%e A(x)^9 = 1 + 9*x + 45*x^2 + 165*x^3 + 504*x^4 + 1404*x^5 +...
%e A(x)^16 = 1 + 16*x + 136*x^2 + 816*x^3 + 3892*x^4 + 15824*x^5 +...
%e where
%e A(x) = 1 + x*A(x) + x^4*A(x)^4 + x^9*A(x)^9 + x^16*A(x)^16 +...
%e A(x) = (1/x)*Series_Reversion(x/(1 + x + x^4 + x^9 + x^16 +...)).
%t f[x_, y_, d_] := f[x, y, d] = If[x < 0 || y < x, 0, If[x == 0 && y == 0, 1, f[x-1, y, 0] + f[x, y - If[d == 0, 1, Sqrt[d]*2 + 1], If[d == 0, 1, Sqrt[d]*2 + 1 + d]]]]; Table[f[n, n, 0], {n, 0, 31}] (* _David Scambler_, May 09 2012 *)
%o (PARI) {a(n)=local(A=1+x+x*O(x^n));for(i=1,n,A=sum(m=0,sqrtint(n),x^(m^2)*A^(m^2)));polcoeff(A,n)}
%o (PARI) seq(n)={Vec(serreverse(x/(1 + sum(i=1, sqrtint(n), x^(i^2))) + O(x*x^n)))} \\ _Andrew Howroyd_, Apr 28 2018
%Y Cf. A157133, A157134, A157135.
%K nonn
%O 0,5
%A _Paul D. Hanna_, Feb 24 2009
| 