Reminder: The OEIS is hiring a new managing editor, and the application deadline is January 26.
%I #15 Mar 07 2022 07:55:18
%S 1,1,2,6,24,840,720,5040,40320,59814720,3628800,83701537920,479001600,
%T 26980643289600,2642646473026560,1307674368000,20922789888000,
%U 41837259585747225600,6402373705728000,598354114828973074790400,18160977780223038067507200
%N Number of endofunctions on [n] such that the LCM of their cycle lengths equals n.
%H Alois P. Heinz, <a href="/A290961/b290961.txt">Table of n, a(n) for n = 1..389</a>
%F a(n) = A222029(n,n).
%p b:= proc(n, m) option remember; `if`(n=0, x^m, add((j-1)!*
%p b(n-j, ilcm(m, j))*binomial(n-1, j-1), j=1..n))
%p end:
%p a:= n-> add(coeff(b(j, 1), x, n)*n^(n-j)*binomial(n-1, j-1), j=0..n):
%p seq(a(n), n=1..25);
%t b[n_, m_] := b[n, m] = If[n == 0, x^m, Sum[(j - 1)!*
%t b[n - j, LCM[m, j]]*Binomial[n - 1, j - 1], {j, 1, n}]];
%t a[n_] := Sum[Coefficient[b[j, 1], x, n]*n^(n-j)*Binomial[n-1, j-1], {j, 0, n}];
%t Table[a[n], {n, 1, 25}] (* _Jean-François Alcover_, Mar 07 2022, after _Alois P. Heinz_ *)
%Y Main diagonal of A222029.
%Y Cf. A074351 (the same for permutations).
%K nonn
%O 1,3
%A _Alois P. Heinz_, Aug 15 2017