OFFSET
1,2
LINKS
Natalia L. Skirrow, counting endofunctions on [n] with j trees and k leaves: interpreting a Lagrange inversion result, Mathematics StackExchange (2026)
FORMULA
E.g.f.: 1/(1 - A(x,y)) where A(x,y) is the e.g.f. for A055302.
T(n,k) = (n!/k!) * Sum_{j=1..k} (n-1)!/(n-j)! * j * Stirling2(n-j,n-k). - Natalia L. Skirrow, Feb 13 2026
Sum_{k=1..n} k * T(n,k) = n * A089901(n-1). - Natalia L. Skirrow, Feb 13 2026
EXAMPLE
T(3,3) = 6 because we have: (1,2,3),(2,1,3),(3,2,1),(1,3,2),(2,3,1),(3,1,2). In these 6 functions represented as a word there are 3 (all) elements with zero nonrecurrent elements mapped to them.
Triangle begins:
nk| 1 2 3 4 5 6
--+-----------------------------
1| 1
2| 2 2
3| 6 15 6
4| 24 108 100 24
5|120 840 1340 705 120
6|720 7200 17400 15150 5466 720
...
MATHEMATICA
nn=6; Map[Select[#, #>0&]&, Drop[Range[0, nn]!CoefficientList[Series[1/(1- (-x+x y-ProductLog[-Exp[x (-1+y)] x])), {x, 0, nn}], {x, y}], 1]]//Grid
PROG
(Python)
from math import factorial as fact
from sympy.functions.combinatorial.numbers import stirling
A231536=lambda n, k: fact(n)//fact(k)*sum(fact(n-1)//fact(n-j)*j*stirling(n-j, n-k) for j in range(1, k+1)) # Natalia L. Skirrow, Feb 13 2026
CROSSREFS
KEYWORD
nonn,tabl
AUTHOR
Geoffrey Critzer, Nov 10 2013
STATUS
approved
