OFFSET
0,8
COMMENTS
A cycle type set of [k] requires that n be at least k*(k+1)/2 hence in the array k ranges from 0 to floor((sqrt(1+8*n)-1)/2).
REFERENCES
P. Flajolet and R. Sedgewick, Analytic Combinatorics, Cambridge University Press, 2009.
LINKS
Alois P. Heinz, Rows n = 0..460, flattened
Marko Riedel, Maple code for sequence.
Marko Riedel, Derivation of recurrence and generating functions.
FORMULA
T(n,k) = n! * rho(n,k) where rho(n,k) = Sum_{p=1..floor(n/k)} rho(n-k*p,k-1)/p!/k^p with rho(n,0) = delta(n,0).
EXAMPLE
Triangle begins:
n/k 0 1 2 3 4
+-------------------------
0| 1
1| 0 1
2| 0 1
3| 0 1 3
4| 0 1 6
5| 0 1 25
6| 0 1 60 120
7| 0 1 231 420
8| 0 1 658 2800
9| 0 1 2619 2016
10| 0 1 8550 105840 151200
...
MAPLE
b:= proc(n, i) option remember; `if`(n=0, x^(i-1), add(b(n-i*j, i+1)*
combinat[multinomial](n, n-i*j, i$j)/j!*(i-1)!^j, j=1..n/i))
end:
T:= n-> (p-> seq(coeff(p, x, i), i=0..degree(p)))(b(n, 1)):
seq(T(n), n=0..15); # Alois P. Heinz, Feb 01 2026
MATHEMATICA
b[n_, i_] := b[n, i] = If[n == 0, x^(i-1), Sum[b[n-i*j, i+1]* Multinomial@@Join[{n-i*j}, Table[i, {j}]]/j!*(i-1)!^j, {j, 1, n/i}]];
T[n_] := CoefficientList[b[n, 1], x];
Table[T[n], {n, 0, 15}] // Flatten (* Jean-François Alcover, Feb 14 2026, after Alois P. Heinz *)
CROSSREFS
KEYWORD
nonn,tabf
AUTHOR
Marko Riedel, Jan 31 2026
STATUS
approved
