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.
%I A219585 #33 Sep 16 2019 04:29:24 %S A219585 1,1,1,1,1,1,1,2,1,1,1,5,5,2,1,1,15,40,17,2,1,1,52,457,364,46,3,1,1, %T A219585 203,6995,14595,2897,123,4,1,1,877,136771,937776,407287,21369,323,5,1, %U A219585 1,4140,3299218,88507276,107652681,10200931,148257,809,6,1 %N A219585 Number A(n,k) of k-partite partitions of {n}^k into distinct k-tuples; square array A(n,k), n>=0, k>=0, read by antidiagonals. %C A219585 A(n,k) is the number of factorizations of m^n into distinct factors where m is a product of k distinct primes. A(2,2) = 5: (2*3)^2 = 36 has 5 factorizations into distinct factors: 36, 3*12, 4*9, 2*18, 2*3*6. %H A219585 Andrew Howroyd, <a href="/A219585/b219585.txt">Table of n, a(n) for n = 0..209</a> %F A219585 A(n,k) = [(Product_{j=1..k} x_j)^n] 1/2 * Product_{i_1,...,i_k>=0} (1+Product_{j=1..k} x_j^i_j). %e A219585 A(1,3) = 5: [(1,1,1)], [(1,1,0),(0,0,1)], [(1,0,1),(0,1,0)], [(1,0,0),(0,1,0),(0,0,1)], [(0,1,1),(1,0,0)]. %e A219585 A(3,2) = 17: [(3,3)], [(3,0),(0,3)], [(3,2),(0,1)], [(2,3),(1,0)], [(3,1),(0,2)], [(2,2),(1,1)], [(1,3),(2,0)], [(2,1),(1,2)], [(2,1),(1,1),(0,1)], [(3,0),(0,2),(0,1)], [(2,2),(1,0),(0,1)], [(2,1),(0,2),(1,0)], [(1,2),(2,0),(0,1)], [(1,2),(1,1),(1,0)], [(0,3),(2,0),(1,0)], [(2,0),(1,1),(0,2)], [(2,0),(0,2),(1,0),(0,1)]. %e A219585 Square array A(n,k) begins: %e A219585 1, 1, 1, 1, 1, 1, 1, ... %e A219585 1, 1, 2, 5, 15, 52, 203, ... %e A219585 1, 1, 5, 40, 457, 6995, 136771, ... %e A219585 1, 2, 17, 364, 14595, 937776, 88507276, ... %e A219585 1, 2, 46, 2897, 407287, 107652681, ... %e A219585 1, 3, 123, 21369, 10200931, 10781201973, ... %e A219585 1, 4, 323, 148257, 233051939, ... %e A219585 1, 5, 809, 970246, 4909342744, ... %t A219585 f[n_, k_] := f[n, k] = 1/2 Product[Sum[O[x[j]]^(n+1), {j, 1, k}]+1+ Product[x[j]^i[j], {j, 1, k}], Evaluate[Sequence @@ Table[{i[j], 0, n}, {j, 1, k}]]]; %t A219585 a[0, _] = a[_, 0] = 1; a[n_, k_] := SeriesCoefficient[f[n, k], Sequence @@ Table[{x[j], 0, n}, {j, 1, k}]]; %t A219585 Table[Print[a[n-k, k]]; a[n-k, k], {n, 0, 9}, {k, n, 0, -1}] // Flatten (* _Jean-François Alcover_, Dec 11 2013, updated Sep 16 2019 *) %o A219585 (PARI) %o A219585 EulerT(v)={Vec(exp(x*Ser(dirmul(v,vector(#v,n,1/n))))-1, -#v)} %o A219585 D(p, n, k)={my(v=vector(n)); for(i=1, #p, v[p[i]]++); EulerT(v)[n]^k/prod(i=1, #v, i^v[i]*v[i]!)} %o A219585 T(n, k)={my(m=n*k+1, q=Vec(exp(intformal(O(x^m) - x^n/(1-x)))/(1+x))); if(n==0, 1, (-1)^m*sum(j=0, m, my(s=0); forpart(p=j, s+=(-1)^#p*D(p, n, k), [1, n]); s*q[#q-j])/2)} \\ _Andrew Howroyd_, Dec 16 2018 %Y A219585 Columns k=0..5 give: A000012, A000009, A219554, A219560, A219561, A219565. %Y A219585 Rows n=0..3 give: A000012, A000110, A094574, A319591. %Y A219585 Cf. A188445, A219727 (partitions of {n}^k into k-tuples), A318286. %K A219585 nonn,tabl %O A219585 0,8 %A A219585 _Alois P. Heinz_, Nov 23 2012