|
%I #36 Nov 16 2017 17:18:38
%S 1,1,10,159,3496,98345,3373056,136535455,6371523712,336784920849,
%T 19888195110400,1297716672601151,92721494240225280,
%U 7199830049013964921,603715489091812335616,54366622743565012989375,5233114241479255004839936,536180296483497244155041825
%N Number of pairs of functions (f, g) on a set of n elements into itself satisfying f(g(x)) = f(x).
%H Alois P. Heinz, <a href="/A239761/b239761.txt">Table of n, a(n) for n = 0..345</a>
%F a(n) ~ 5^(-1/4) * ((1+sqrt(5))/2)^(3*n-1/2) * n^n / exp(2*n/(1+sqrt(5))). - _Vaclav Kotesovec_, Aug 07 2014
%F a(n) = Sum_{k = 1..n} A060281(n,k) n^k. - _David Einstein_, Oct 31 2016
%F a(n) = n! * [x^n] 1/(1 + LambertW(-x))^n. - _Ilya Gutkovskiy_, Oct 03 2017
%p b:= proc(n, i) option remember; `if`(n=0, 1, `if`(i<1, 0,
%p add(b(n-j, i-1)*binomial(n, j)*j^j, j=0..n)))
%p end:
%p a:= n-> b(n$2):
%p seq(a(n), n=0..20); # _Alois P. Heinz_, Jul 17 2014
%t f4[n_] := Sum[n^k Sum[Binomial[n - 1, j]*n^(n - 1 - j)*StirlingS1[j + 1, k] *(-1)^(j + k + 1), {j, 0, n - 1}], {k, 1, n}] (* _David Einstein_, Oct 31 2016 *)
%Y Cf. A181162, A239769, A239773.
%Y Column k=1 of A245910.
%Y Cf. A001622, A295188.
%K nonn
%O 0,3
%A _Chad Brewbaker_, Mar 26 2014
%E a(6)-a(7) from _Giovanni Resta_, Mar 28 2014
%E a(8)-a(17) from _Alois P. Heinz_, Jul 17 2014
|
