%I #27 Nov 14 2024 02:10:58
%S 0,1,1,6,1,180,1,5040,20160,378000,1,59875200,1,6235669440,
%T 47221574400,1307674368000,1,533531142144000,1,126713646259200000,
%U 1219830034655232000,51090956251003468800,1,38778025108327464960000,25852016738884976640000
%N a(n) = Sum_{p|n, p prime} (n-1)!/(p-1)!.
%H Robert Israel, <a href="/A352012/b352012.txt">Table of n, a(n) for n = 1..449</a>
%F E.g.f.: -Sum_{p prime} log(1-x^p)/p!.
%F a(n) = 1 if and only if n is prime.
%p f:= proc(n) local p;
%p add( (n-1)!/(p-1)!, p = numtheory:-factorset(n))
%p end proc:
%p map(f, [$1..30]): # _Robert Israel_, Nov 14 2024
%t a[1] = 0; a[n_] := (n - 1)! * Plus @@ (1/(FactorInteger[n][[;; , 1]] - 1)!); Array[a, 25] (* _Amiram Eldar_, Mar 01 2022 *)
%o (PARI) a(n) = sumdiv(n, d, isprime(d)*(n-1)!/(d-1)!);
%o (PARI) my(N=40, x='x+O('x^N)); concat(0, Vec(serlaplace(-sum(k=1, N, isprime(k)*log(1-x^k)/k!))))
%o (PARI) a(n) = my(f=factor(n)); sum(k=1, #f~, (n-1)!/(f[k,1]-1)!); \\ _Michel Marcus_, Mar 01 2022
%Y Cf. A087906, A352004, A352014.
%K nonn
%O 1,4
%A _Seiichi Manyama_, Feb 28 2022