cp's OEIS Frontend

A060052 Triangle read by rows: T(n,k) gives number of r-bicoverings of an n-set with k blocks, n >= 2, k = 3..n+floor(n/2).

%I A060052 #22 Jan 30 2020 18:14:09
%S A060052 1,1,4,0,15,25,3,0,30,222,226,40,0,30,1230,3670,2706,535,15,0,0,5040,
%T A060052 39900,69450,40405,8141,420,0,0,15120,345240,1254960,1498035,722275,
%U A060052 142877,9730,105,0,0,30240,2492280,18587520,40701780,36450820,15031204,2871240,226828,5040
%N A060052 Triangle read by rows: T(n,k) gives number of r-bicoverings of an n-set with k blocks, n >= 2, k = 3..n+floor(n/2).
%C A060052 A bicovering is r-bicovering if intersection of every two blocks contains at most one element.
%D A060052 I. P. Goulden and D. M. Jackson, Combinatorial Enumeration, John Wiley and Sons, N.Y., 1983.
%H A060052 Andrew Howroyd, <a href="/A060052/b060052.txt">Table of n, a(n) for n = 2..1802</a> (rows n=2..50)
%H A060052 <a href="/A060052/a060052.pdf">Table</a>
%F A060052 E.g.f.: A(x, y) = exp(-x-1/2*x^2*y)*Sum_{i>=0} (1+y)^binomial(i, 2)*x^i/i!.
%F A060052 T(n, k) = (n!/k!) * A276640(k, n). - _David Pasino_, Sep 22 2016
%F A060052 T(n,k) = 0 for n > binomial(k,2). - _Andrew Howroyd_, Jan 30 2020
%e A060052 Triangle starts:
%e A060052 [1],
%e A060052 [1, 4],
%e A060052 [0, 15, 25, 3],
%e A060052 [0, 30, 222, 226, 40],
%e A060052 [0, 30, 1230, 3670, 2706, 535, 15],
%e A060052 [0, 0, 5040, 39900, 69450, 40405, 8141, 420],
%e A060052 [0, 0, 15120, 345240, 1254960, 1498035, 722275, 142877, 9730, 105],
%e A060052 [0, 0, 30240, 2492280, 18587520, 40701780, 36450820, 15031204, 2871240, 226828, 5040],
%e A060052 ...
%o A060052 (PARI) \\ returns k-th column as vector.
%o A060052 C(k)=if(k<3, [], Vecrev(serlaplace(polcoef(exp(-x-1/2*x^2*y + O(x*x^k))*sum(i=0, 3*k\2, (1+y)^binomial(i, 2)*x^i/i!), k))/y)) \\ _Andrew Howroyd_, Jan 30 2020
%o A060052 (PARI)
%o A060052 T(n)={my(m=(3*n\2), y='y + O('y^(n+1))); my(g=exp(-x-1/2*x^2*y + O(x*x^m))*sum(k=0, m, (1+y)^binomial(k, 2)*x^k/k!)); Mat([Col(serlaplace(p), -n) | p<-Vec(g)[2..m+1]])}
%o A060052 { my(A=T(8)); for(n=2, matsize(A)[1], print(A[n, 3..3*n\2])) } \\ _Andrew Howroyd_, Jan 30 2020
%Y A060052 Row sums are A060053.
%Y A060052 Column sums are A060051.
%Y A060052 Cf. A059443, A059530, A060092, A060487, A060492, A276640, A331039.
%K A060052 easy,nonn,tabf
%O A060052 2,3
%A A060052 _Vladeta Jovovic_, Feb 15 2001
%E A060052 Zeros inserted into data by _Andrew Howroyd_, Jan 30 2020