%I #15 Jul 20 2023 10:39:05
%S 1,0,1,1,6,21,110,904,4312,74400,731412,5600761,128196024,792051157,
%T 18696610816,264267572121,7136433698464,57948743342529,
%U 2228312959187256,22463157401776612,681974906329502904,15395459281239915282,463374873030990445252,6091833036158810701465
%N Number of endofunctions on [n] such that the number of elements that are mapped to i is either 0 or a prime divisor of i.
%H Alois P. Heinz, <a href="/A364328/b364328.txt">Table of n, a(n) for n = 0..444</a>
%e a(0) = 1: ().
%e a(2) = 1: (22).
%e a(3) = 1: (333).
%e a(4) = 6: (4422), (4242), (4224), (2442), (2424), (2244).
%e a(5) = 21: (55555), (44333), (43433), (43343), (43334), (34433), (34343), (34334), (33443), (33434), (33344), (33322), (33232), (33223), (32332), (32323), (32233), (23332), (23323), (23233), (22333).
%p b:= proc(n, i) option remember; `if`(n=0, 1, `if`(i<1, 0, b(n, i-1)+add(
%p `if`(d>n, 0, b(n-d, i-1)*binomial(n, d)), d=numtheory[factorset](i))))
%p end:
%p a:= n-> b(n$2):
%p seq(a(n), n=0..23);
%Y Cf. A000040, A178682, A334370, A364327, A364344.
%K nonn
%O 0,5
%A _Alois P. Heinz_, Jul 18 2023