%I M4458 N1889 #102 Jul 07 2024 12:10:53
%S 1,1,7,127,4369,243649,20036983,2280356863,343141433761,
%T 65967241200001,15773461423793767,4591227123230945407,
%U 1598351733247609852849,655782249799531714375489,313160404864973852338669783,172201668512657346455126457343,108026349476762041127839800617281
%N a(n) = Sum_{k=0..n-1} binomial(2*n,2*k)*a(k)*a(n-k-1).
%C Number of increasing rooted triangular cacti of 2n+1 nodes. (In an increasing rooted graph, nodes are numbered and numbers increase as you move away from the root.)
%C a(n) is (2n)!/2^n times the n-th coefficient in the series for inverf(2x/sqrt(Pi)). - _Paul Barry_, Apr 12 2010
%C Number of ordered bilabeled increasing trees with 2n labels. - _Markus Kuba_, Nov 17 2014
%C Limit_{n->oo} (a(n)/(n!)^2)^(1/n) = 8/Pi. - _Vaclav Kotesovec_, Nov 19 2014
%C From _David Callan_, Jul 21 2017: (Start)
%C Conjectures:
%C a(n) is the Hafnian of the triangular array (u(i,j))_{1 <= i < j <= 2n} with u(i,j)=i. The Hafnian is the same as the Pfaffian except without the alternating signs just as the permanent of a matrix is the determinant without the signs.
%C a(n) is the total weight of Dyck n-paths with weight defined as follows. Given a Dyck path, for each upstep, record its position in the path and the height of its upper endpoint; then multiply together all of these positions and heights. For example, the Dyck 4-path P = UUDUUDDD has upsteps in positions 1,2,4,5 ending at heights 1,2,2,3 respectively, and hence weight(P) = 480. (In fact the positions determine the heights because, for the k-th upstep, position + height = 2k.) (End)
%D F. Bergeron, G. Labelle and P. Leroux, Combinatorial Species and Tree-Like Structures, Camb. 1998, cf. Chapter 5.
%D N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence).
%D N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).
%H T. D. Noe, <a href="/A002067/b002067.txt">Table of n, a(n) for n = 0..50</a>
%H L. Carlitz, <a href="http://projecteuclid.org/euclid.pjm/1103035736">The inverse of the error function</a>, Pacific J. Math., 13 (1963), 459-470. [See Eq. 1.3 and Section 6.]
%H D. Dominici, <a href="http://arxiv.org/abs/math/0607230">Asymptotic analysis of the derivatives of the inverse error function</a>, arXiv:math/0607230 [math.CA], 2006-2007.
%H A. J. E. M. Janssen, <a href="https://www.eurandom.tue.nl/wp-content/uploads/2023/03/mrt2023_final.pdf">Analysis of a constrained initial value for an ODE arising in the study of a power-flow model</a>, Eindhoven Univ. Tech. (Netherlands, 2023).
%H Markus Kuba and Alois Panholzer, <a href="http://arxiv.org/abs/1411.4587">Combinatorial families of multilabelled increasing trees and hook-length formulas</a>, arXiv:1411.4587 [math.CO], 2014.
%H Wikipedia, <a href="http://en.wikipedia.org/wiki/Error_function">Error Function</a>
%H <a href="/index/Ca#cacti">Index entries for sequences related to cacti</a>
%F a(n) = b(2n+1), where e.g.f. of b satisfies B'(x)=exp(B(x)^2/2).
%F a(n) = (2n)! * A092676(n) / (2^n*A092677(n)). - _Paul Barry_, Apr 12 2010
%F a(n) = 1/2^n * A026944(n+1). Let D denote the operator g(x) -> (1/sqrt(2))*d/dx(exp(x^2)*g(x)). Then a(n) = D^(2*n)(1) evaluated at x = 0. - _Peter Bala_, Sep 08 2011
%F E.g.f. B(x)=Sum_{n>=1} a(n-1)*x^(2*n)/(2*n)! satisfies differential equation B''(x) - B(x)*B''(x) - 1 = 0, B'(0)=1/2. - _Vladimir Kruchinin_, Aug 12 2019
%F E.g.f. satisfies: A(x) = exp( Integral A(x)*B(x) dx ), where A(x) = Sum_{n>=0} a(n)*x^(2*n)/(2*n)! and B(x) = Sum_{n>=0} a(n)*x^(2*n+1)/(2*n+1)!, and the constant of integration is zero. - _Paul D. Hanna_, Jun 02 2015 [formula revised by _Paul D. Hanna_, Jul 06 2024 following a suggestion from _Petros Hadjicostas_]
%e E.g.f.: A(x) = 1 + x^2/2! + 7*x^4/4! + 127*x^6/6! + 4369*x^8/8! + ...
%p a:=proc(n) option remember; if n <= 0 then RETURN(1); else RETURN( add( binomial(2*n,2*k)*a(k)*a(n-k-1), k=0..n-1 ) ); fi; end;
%t max = 16; se = Series[ InverseErf[ 2*x/Sqrt[Pi] ], {x, 0, 2*max+1} ]; a[n_] := (2n+1)!/2^n*Coefficient[ se, x, 2*n+1]; Table[ a[n], {n, 0, max} ] (* _Jean-François Alcover_, Mar 07 2012, after _Paul Barry_ *)
%o (PARI) /* E.g.f. A(x) = exp( Integral A(x) * Integral A(x) dx dx ): */
%o {a(n) = local(A=1+x); for(i=1,n, A = exp( intformal( A * intformal( A + x*O(x^n)) ) ) ); n!*polcoeff(A,n)}
%o for(n=0,20,print1(a(2*n),", ")) \\ _Paul D. Hanna_, Jun 02 2015
%o (PARI) /* By definition: */
%o {a(n) = if(n==0,1,sum(k=0,n-1, binomial(2*n,2*k)*a(k)*a(n-k-1)))}
%o for(n=0,20,print1(a(n),", ")) \\ _Paul D. Hanna_, Jun 02 2015
%Y The sequence of fractions A092676/A132467 is closely related.
%Y Periods: A122149, A122159.
%K nonn,eigen,easy,nice
%O 0,3
%A _N. J. A. Sloane_
%E Alternate description, formula and comment from _Christian G. Bower_
%E New definition and more terms from _Vladeta Jovovic_, Oct 22 2005