# Greetings from The On-Line Encyclopedia of Integer Sequences! http://oeis.org/

Search: id:a057271
Showing 1-1 of 1

%I A057271 #16 Jan 16 2022 17:46:36
%S A057271 1,0,2,1,0,0,6,20,15,6,1,0,0,0,24,234,672,908,792,495,220,66,12,1,0,0,
%T A057271 0,0,120,2544,16880,55000,111225,161660,183006,167660,125945,77520,
%U A057271 38760,15504,4845,1140,190,20,1
%N A057271 Triangle T(n,k) of number of digraphs with a source and a sink on n labeled nodes and k arcs, k=0,1,..,n*(n-1).
%D A057271 V. Jovovic, G. Kilibarda, Enumeration of labeled initially-finally connected digraphs, Scientific review, Serbian Scientific Society, 19-20 (1996), p. 245.
%H A057271 Andrew Howroyd, <a href="/A057271/b057271.txt">Table of n, a(n) for n = 1..2680</a> (rows 1..20)
%H A057271 V. Jovovic and G. Kilibarda, <a href="http://dx.doi.org/10.1016/S0012-365X(00)00112-6">Enumeration of labeled quasi-initially connected digraphs</a>, Discrete Math., 224 (2000), 151-163.
%H A057271 R. W. Robinson, <a href="http://cobweb.cs.uga.edu/~rwr/publications/components.pdf">Counting digraphs with restrictions on the strong components</a>, Combinatorics and Graph Theory '95 (T.-H. Ku, ed.), World Scientific, Singapore (1995), 343-354.
%e A057271 Triangle starts:
%e A057271 [1] 1;
%e A057271 [2] 0,2,1;
%e A057271 [3] 0,0,6,20,15,6,1;
%e A057271 [4] 0,0,0,24,234,672,908,792,495,220,66,12,1;
%e A057271   ...
%e A057271 The number of digraphs with a source and a sink on 3 labeled nodes is 48 = 6+20+15+6+1.
%o A057271 (PARI) \\ Following Eqn 20 in the Robinson reference.
%o A057271 Z(p,f)={my(n=serprec(p,x)); serconvol(p, sum(k=0, n-1, x^k*f(k), O(x^n)))}
%o A057271 G(e,p)={Z(p, k->1/e^(k*(k-1)/2))}
%o A057271 U(e,p)={Z(p, k->e^(k*(k-1)/2))}
%o A057271 DigraphEgf(n,e)={sum(k=0, n, e^(k*(k-1))*x^k/k!, O(x*x^n) )}
%o A057271 StrongD(n,e=2)={-log(U(e, 1/G(e, DigraphEgf(n, e))))}
%o A057271 InitFinally(n, e=2)={my(S=StrongD(n, e)); Vec(serlaplace( S - S^2 + exp(S) * U(e, G(e, S*exp(-S))^2*G(e, DigraphEgf(n,e))) ))}
%o A057271 row(n)={Vecrev(InitFinally(n, 1+'y)[n]) }
%o A057271 { for(n=1, 5, print(row(n))) } \\ _Andrew Howroyd_, Jan 16 2022
%Y A057271 Row sums give A049524.
%Y A057271 The unlabeled version is A057278.
%Y A057271 Cf. A057272, A057273, A057274, A057275, A062735, A350791.
%K A057271 nonn,tabf
%O A057271 1,3
%A A057271 _Vladeta Jovovic_, Goran Kilibarda, Sep 14 2000

# Content is available under The OEIS End-User License Agreement: http://oeis.org/LICENSE