# Greetings from The On-Line Encyclopedia of Integer Sequences! http://oeis.org/

Search: id:a352058
Showing 1-1 of 1

%I A352058 #16 Mar 02 2022 12:02:39
%S A352058 0,1,1,6,1,120,1,840,10080,30240,1,3659040,1,17297280,3632428800,
%T A352058 259459200,1,1490852563200,1,845092863014400,3379030566912000,
%U A352058 28158588057600,1,2565331847811532800,1077167364120207360000,64764752532480000,5001134190558105600000
%N A352058 a(n) = Sum_{p|n, p prime} n!/(p!*(n/p)!).
%F A352058 E.g.f.: Sum_{p prime} (exp(x^p) - 1)/p!.
%F A352058 a(n) = 1 if and only if n is prime.
%t A352058 a[1] = 0; a[n_] := Plus @@ (n!/((p=FactorInteger[n][[;;,1]])!*(n/p)!)); Array[a, 30] (* _Amiram Eldar_, Mar 02 2022 *)
%o A352058 (PARI) a(n) = my(f=factor(n)); sum(k=1, #f~, n!/(f[k, 1]!*(n/f[k, 1])!));
%o A352058 (PARI) my(N=40, x='x+O('x^N)); concat(0, Vec(serlaplace(sum(k=1, N, isprime(k)*(exp(x^k)-1)/k!))))
%Y A352058 Cf. A121860, A352012, A352059, A352060.
%K A352058 nonn
%O A352058 1,4
%A A352058 _Seiichi Manyama_, Mar 02 2022

# Content is available under The OEIS End-User License Agreement: http://oeis.org/LICENSE