login

Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 61st year, we have over 378,000 sequences, and we’ve reached 11,000 citations (which often say “discovered thanks to the OEIS”).

A362763
Array read by antidiagonals: T(n,k) is the number of nonisomorphic k-sets of permutations of an n-set.
4
1, 1, 1, 0, 1, 1, 0, 0, 2, 1, 0, 0, 1, 3, 1, 0, 0, 0, 5, 5, 1, 0, 0, 0, 6, 23, 7, 1, 0, 0, 0, 5, 116, 89, 11, 1, 0, 0, 0, 3, 521, 2494, 484, 15, 1, 0, 0, 0, 1, 1931, 69366, 87984, 2904, 22, 1, 0, 0, 0, 0, 5906, 1592714, 15456557, 4250015, 22002, 30, 1
OFFSET
0,9
COMMENTS
Isomorphism is up to permutation of the elements of the n-set. Each permutation can be considered to be a set of disjoint directed cycles whose vertices cover the n-set. Permuting the elements of the n-set permutes each of the permutations in the k-set.
LINKS
Andrew Howroyd, Table of n, a(n) for n = 0..1325 (first 51 antidiagonals).
FORMULA
T(n,k) = 0 for k > n!.
T(n,k) = T(n, n!-k).
EXAMPLE
Array begins:
====================================================================
n/k| 0 1 2 3 4 5 6 ...
---+----------------------------------------------------------------
0 | 1 1 0 0 0 0 0 ...
1 | 1 1 0 0 0 0 0 ...
2 | 1 2 1 0 0 0 0 ...
3 | 1 3 5 6 5 3 1 ...
4 | 1 5 23 116 521 1931 5906 ...
5 | 1 7 89 2494 69366 1592714 30461471 ...
6 | 1 11 484 87984 15456557 2209040882 263190866673 ...
7 | 1 15 2904 4250015 5329123475 5366409944453 4503264576070573 ...
...
PROG
(PARI)
B(n, k) = {n!*k^n}
K(v)=my(S=Set(v)); prod(i=1, #S, my(k=S[i], c=#select(t->t==k, v)); B(c, k))
R(v, m)=concat(vector(#v, i, my(t=v[i], g=gcd(t, m)); vector(g, i, t/g)))
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}
T(n, k) = {if(n==0, k<=1, my(s=0); forpart(q=n, s+=permcount(q) * polcoef(exp(sum(m=1, k, K(R(q, m))*(x^m-x^(2*m))/m, O(x*x^k))), k)); s/n!)}
CROSSREFS
Columns k=0..3 are A000012, A000041, A362764, A362765.
Row sums are A362766.
Cf. A362644.
Sequence in context: A331385 A026729 A109466 * A259095 A326676 A076833
KEYWORD
nonn,tabl
AUTHOR
Andrew Howroyd, May 03 2023
STATUS
approved