%I #6 Aug 21 2014 04:35:32
%S 5040,408240,22680000,1106859600,51732172800,2408384618640,
%T 113960430904320,5541379593750000,278592031202284800,
%U 14529619059476320800,787422201081850414080,44373594768472437642720,2600326096882824360960000,158404803877320370312773600
%N Number of endofunctions on [n] where the largest cycle length equals 8.
%H Alois P. Heinz, <a href="/A246218/b246218.txt">Table of n, a(n) for n = 8..200</a>
%F a(n) ~ (8*exp(761/280) - 7*exp(363/140)) * n^(n-1). - _Vaclav Kotesovec_, Aug 21 2014
%p with(combinat):
%p b:= proc(n, i) option remember; `if`(n=0, 1, `if`(i<1, 0,
%p add((i-1)!^j*multinomial(n, n-i*j, i$j)/j!*
%p b(n-i*j, i-1), j=0..n/i)))
%p end:
%p A:= (n, k)-> add(binomial(n-1, j-1)*n^(n-j)*b(j, min(j, k)), j=0..n):
%p a:= n-> A(n, 8) -A(n, 7):
%p seq(a(n), n=8..25);
%Y Column k=8 of A241981.
%K nonn
%O 8,1
%A _Alois P. Heinz_, Aug 19 2014
|