OFFSET
0,9
COMMENTS
Two permutations x,y on [n] commute if x*y = y*x.
LINKS
Andrew Howroyd, Table of n, a(n) for n = 0..1325 (first 51 antidiagonals).
Tad White, Counting Free Abelian Actions, arXiv preprint arXiv:1304.2830 [math.CO], 2013.
FORMULA
T(n,k) = n!*A362826(n,k) for k > 0.
EXAMPLE
Array begins:
========================================================
n/k| 0 1 2 3 4 5 6 ...
---+----------------------------------------------------
0 | 1 1 1 1 1 1 1 ...
1 | 1 1 1 1 1 1 1 ...
2 | 1 2 4 8 16 32 64 ...
3 | 1 6 18 48 126 336 918 ...
4 | 1 24 120 504 2016 7944 31200 ...
5 | 1 120 840 4680 24720 130440 699840 ...
6 | 1 720 7920 66240 516240 3968640 30672720 ...
7 | 1 5040 75600 856800 9122400 97030080 1050336000 ...
...
PROG
(PARI)
EulerT(v)={Vec(exp(x*Ser(dirmul(v, vector(#v, n, 1/n))))-1, -#v)}
M(n, m=n)={my(v=vector(m+1), u=vector(n, n, n==1), f=vector(n, n, n!)); v[1]=vectorv(n+1, i, 1); for(j=1, #v-1, my(t=EulerT(u)); v[j+1]=vectorv(n+1, i, i--; if(i, f[i]*t[i], 1)); u=dirmul(u, vector(n, n, n^(j-1)))); Mat(v)}
{ my(A=M(7)); for(n=1, #A, print(A[n, ])) }
CROSSREFS
KEYWORD
nonn,tabl
AUTHOR
Andrew Howroyd, May 08 2023
STATUS
approved