%I #34 Jan 17 2022 14:02:43
%S 1,3,1,54,9,1,3834,243,18,1,1027080,20790,675,30,1,1067308488,6364170,
%T 67635,1485,45,1,4390480193904,7543111716,23031540,171045,2835,63,1,
%U 72022346388181584,35217115838604,30469951764,63580545,370440,4914,84,1
%N Triangular array read by rows. T(n,k) is the number of digraphs with n labeled nodes having exactly k undirected (or weak) components, n >= 1, 1 <= k <= n.
%C The Bell transform of A003027(n+1). For the definition of the Bell transform see A264428. - _Peter Luschny_, Jan 18 2016
%H Alois P. Heinz, <a href="/A189898/b189898.txt">Rows n = 1..32, flattened</a>
%F E.g.f. for column k: log(A(x))^k/k! where A(x) is the e.g.f. for A053763.
%e 1
%e 3 1
%e 54 9 1
%e 3834 243 18 1
%e 1027080 20790 675 30 1
%p T:= (n, k)-> coeff(series(log(add(2^(i^2-i) *x^i/i!, i=0..n))^k /k!,
%p x, n+1), x, n) *n!:
%p seq(seq(T(n, k), k=1..n), n=1..10); # _Alois P. Heinz_, May 01 2011
%t a= Sum[4^Binomial[n,2]x^n/n!,{n,0,10}];
%t Transpose[Map[Drop[#, 1] &,Table[Range[0, 10]! CoefficientList[Series[Log[a]^n/n!, {x, 0, 10}], x], {n, 1, 10}]]] // Grid
%o (Sage) # uses[bell_matrix from A264428, A003027]
%o # Adds a column 1,0,0,0, ... at the left side of the triangle.
%o bell_matrix(lambda n: A003027(n+1), 10) # _Peter Luschny_, Jan 18 2016
%Y Column 1 = A003027, row sums = A053763, lower diagonal = A045943.
%K nonn,tabl
%O 1,2
%A _Geoffrey Critzer_, May 01 2011
%E Name clarified by _Andrew Howroyd_, Jan 11 2022