cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

A331508 Array read by antidiagonals: A(n,k) is the number of nonisomorphic T_0 n-regular set multipartitions (multisets of sets) on a k-set.

This page as a plain text file.
%I A331508 #13 Jan 16 2024 17:33:08
%S A331508 1,1,1,0,1,1,0,1,1,1,0,1,2,1,1,0,1,5,3,1,1,0,1,11,12,4,1,1,0,1,26,66,
%T A331508 25,5,1,1,0,1,68,445,278,44,6,1,1,0,1,177,4279,5532,966,73,7,1,1,0,1,
%U A331508 497,53340,200589,53535,2957,112,8,1,1,0,1,1476,846254,11662671,7043925,431805,8149,166,9,1,1
%N A331508 Array read by antidiagonals: A(n,k) is the number of nonisomorphic T_0 n-regular set multipartitions (multisets of sets) on a k-set.
%C A331508 An n-regular set multipartition is a finite multiset of nonempty sets in which each element appears in n blocks.
%C A331508 A set multipartition is T_0 if for every two distinct elements there exists a block containing one but not the other element.
%C A331508 A(n,k) is the number of nonequivalent binary matrices with k distinct columns and any number of nonzero rows with n ones in every column up to permutation of rows and columns.
%C A331508 A(n,k) is the number of non-isomorphic set-systems with k parts each of size n.
%H A331508 Andrew Howroyd, <a href="/A331508/b331508.txt">Table of n, a(n) for n = 0..152</a> (first 17 antidiagonals)
%F A331508 A306019(n) = Sum_{d|n} A(n/d, d).
%e A331508 Array begins:
%e A331508 ===============================================
%e A331508 n\k | 0 1 2  3    4      5       6        7
%e A331508 ----+------------------------------------------
%e A331508   0 | 1 1 0  0    0      0       0        0 ...
%e A331508   1 | 1 1 1  1    1      1       1        1 ...
%e A331508   2 | 1 1 2  5   11     26      68      177 ...
%e A331508   3 | 1 1 3 12   66    445    4279    53340 ...
%e A331508   4 | 1 1 4 25  278   5532  200589 11662671 ...
%e A331508   5 | 1 1 5 44  966  53535 7043925 ...
%e A331508   6 | 1 1 6 73 2957 431805 ...
%e A331508   ...
%e A331508 The A(2,3) = 5 matrices are:
%e A331508   [1 0 0]  [1 1 0]  [1 1 1]  [1 1 0]  [1 1 0]
%e A331508   [1 0 0]  [1 0 0]  [1 0 0]  [1 0 1]  [1 0 1]
%e A331508   [0 1 0]  [0 1 0]  [0 1 0]  [0 1 0]  [0 1 1]
%e A331508   [0 1 0]  [0 0 1]  [0 0 1]  [0 0 1]
%e A331508   [0 0 1]  [0 0 1]
%e A331508   [0 0 1]
%o A331508 (PARI)
%o A331508 WeighT(v)={Vec(exp(x*Ser(dirmul(v, vector(#v, n, (-1)^(n-1)/n))))-1, -#v)}
%o A331508 permcount(v) = {my(m=1, s=0, k=0, t); for(i=1, #v, t=v[i]; k=if(i>1&&t==v[i-1], k+1, 1); m*=t*k; s+=t); s!/m}
%o A331508 K(q, t, k)={WeighT(Vec(sum(j=1, #q, my(g=gcd(t, q[j])); g*x^(q[j]/g)) + O(x*x^k), -k))[k]}
%o A331508 T(n,k)={my(m=n*k, s=0); if(m==0, k<=1, forpart(q=m, my(g=sum(t=1, k, K(q, t, n)*x^t/t) + O(x*x^k)); s+=permcount(q)*polcoef(exp(g - subst(g,x,x^2)), k)); s/m!)}
%o A331508 { for(n=0, 6, for(k=0, 5, print1(T(n, k), ", ")); print) } \\ _Andrew Howroyd_, Jan 16 2024
%Y A331508 Rows n=1..6 are A000012, A000664, A058790, A058791, A058792, A058793.
%Y A331508 Cf. A283877, A306019, A330942, A331126, A331277, A331461, A331509, A331510.
%K A331508 nonn,tabl
%O A331508 0,13
%A A331508 _Andrew Howroyd_, Jan 18 2020