login
Number of labeled acyclic digraphs with n nodes containing exactly n-1 points of in-degree zero.
22

%I #76 Feb 01 2024 12:10:19

%S 0,2,9,28,75,186,441,1016,2295,5110,11253,24564,53235,114674,245745,

%T 524272,1114095,2359278,4980717,10485740,22020075,46137322,96468969,

%U 201326568,419430375,872415206,1811939301,3758096356,7784628195,16106127330,33285996513

%N Number of labeled acyclic digraphs with n nodes containing exactly n-1 points of in-degree zero.

%C Convolution of 2^n+1 (A000051) and 2^n-1 (A000225). - _Graeme McRae_, Jun 07 2006

%C Let Q be a binary relation on the power set P(A) of a set A having n = |A| elements such that for all nonempty elements x,y of P(A), xRy if x is a proper subset of y and there are no z in P(A) such that x is a proper subset of z and z is a proper subset of y. Then a(n) = |Q|. - _Ross La Haye_, Feb 20 2008, Oct 21 2008

%C Also: convolution of A006589 with A000012 (i.e., partial sums of A006589). - _R. J. Mathar_, Jan 25 2009

%C The La Haye binary relation Q is more clearly stated as x is nonempty and y has one more element than x. If x is a k-set than the number of such pairs is binomial( n, k) * (n-k). - _Michael Somos_, Mar 29 2012

%C Select one of the n nodes of the digraph and select a nonempty subset of the rest to connect to the selected node. This can be done in n * (2^(n-1) - 1) ways. - _Michael Somos_, Mar 29 2012

%C Column 1 of A198204. - _Peter Bala_, Aug 01 2012

%C a(n) is the number of ternary sequences of length n that contain one 0 and at least one 1. For example, a(3)=9 since the sequences are the 3 permutations of 011 and the 6 permutations of 012. - _Enrique Navarrete_, Apr 05 2021

%C a(n) is also the number of multiplications required to compute the permanent of general n X n matrices using canonical trellis method (see Theorem 5, p. 10 in Kiah et al.). - _Stefano Spezia_, Nov 02 2021

%D F. Harary and E. M. Palmer, Graphical Enumeration, Academic Press, NY, 1973, p. 19, (1.6.4).

%D Gerta Rucker and Christoph Rucker, "Walk counts, Labyrinthicity and complexity of acyclic and cyclic graphs and molecules", J. Chem. Inf. Comput. Sci., 40 (2000), 99-106. See Table 1 on page 101. [From _Parthasarathy Nambi_, Sep 26 2008]

%H Vincenzo Librandi, <a href="/A058877/b058877.txt">Table of n, a(n) for n = 1..2001</a>

%H Jia Huang and Erkko Lehtonen, <a href="https://arxiv.org/abs/2401.15786">Associative-commutative spectra for some varieties of groupoids</a>, arXiv:2401.15786 [math.CO], 2024. See p. 12.

%H Han Mao Kiah, Alexander Vardy and Hanwen Yao, <a href="https://arxiv.org/abs/2107.07377">Computing Permanents on a Trellis</a>, arXiv:2107.07377 [cs.IT], 2021.

%H Szakács, Tamás <a href="https://doi.org/10.1515/cm-2017-0011">Convolution of second order linear recursive sequences. II.</a> Commun. Math. 25, No. 2, 137-148 (2017), remark 3

%H <a href="/index/Rec#order_04">Index entries for linear recurrences with constant coefficients</a>, signature (6,-13,12,-4).

%F a(n+1) = (n+1)*2^n - n - 1 = Sum_{j=0..n} (n+j)*2^(n-j-1) = A048493(n)-1 = Column sum of A062111. - _Henry Bottomley_, May 30 2001

%F From _R. J. Mathar_, Jan 25 2009: (Start)

%F G.f.: x^2*(2-3*x)/((1-2*x)*(1-x))^2.

%F a(n) = 6*a(n-1) - 13*a(n-2) + 12*a(n-3) - 4*a(n-4). (End)

%F a(n) = Sum_{k=1..n-1} binomial(n, k) * (n-k). - _Michael Somos_, Mar 29 2012

%F E.g.f: x*exp(x)*(exp(x)-1). - _Enrique Navarrete_, Apr 05 2021

%e G.f. = 2*x^2 + 9*x^3 + 28*x^4 + 75*x^5 + 186*x^6 + 441*x^7 + 1016*x^8 + ...

%p [seq (stirling2(n,2)*n,n=1..29)]; # _Zerinvary Lajos_, Dec 06 2006

%p a:=n->sum(k*binomial(n,k), k=2..n): seq(a(n), n=1..29); # _Zerinvary Lajos_, May 08 2007

%t Table[(n+1)*2^n-n-1, {n, 0, 200}] (* _Vladimir Joseph Stephan Orlovsky_, Jun 30 2011 *)

%t a[ n_] := Sum[ Binomial[ n, k] (n - k), {k, n-1}]; (* _Michael Somos_, Mar 29 2012 *)

%o (Sage) [stirling_number2(i,2)*i for i in range(1,26)] # _Zerinvary Lajos_, Jun 27 2008

%o (Sage) [(n+1)*gaussian_binomial(n,1,2) for n in range(0,29)] # _Zerinvary Lajos_, May 31 2009

%o (Magma) [(n+1)*2^n-n-1: n in [0..30]]; // _Vincenzo Librandi_, Sep 26 2011

%o (PARI) {a(n) = if( n<1, 0, n * (2^(n-1) - 1))} /* _Michael Somos_, Mar 29 2012 */

%Y Second column of A058876. Cf. A003025, A003026.

%Y Column k=1 of A133399.

%Y Cf. A198204.

%K nonn,easy

%O 1,2

%A _N. J. A. Sloane_, Jan 07 2001

%E More terms from _Vladeta Jovovic_, Apr 10 2001