A351708 a(n) = Sum_{p|n, p prime} (n/p)!.
0, 1, 1, 2, 1, 8, 1, 24, 6, 122, 1, 744, 1, 5042, 126, 40320, 1, 363600, 1, 3628824, 5046, 39916802, 1, 479041920, 120, 6227020802, 362880, 87178291224, 1, 1307677997520, 1, 20922789888000, 39916806, 355687428096002, 5160, 6402374184729600, 1
Offset: 1
Examples
a(6) = 8; a(6) = Sum_{p|6} (6/p)! = (6/2)! + (6/3)! = 3*2 + 2*1 = 8.
Links
- Robert Israel, Table of n, a(n) for n = 1..897
Programs
-
Maple
f:= proc(n) local p; add((n/p)!, p = numtheory:-factorset(n)) end proc: map(f, [$1..40]); # Robert Israel, Dec 05 2022
-
Mathematica
Table[DivisorSum[n, Factorial[n/#] &, PrimeQ], {n, 36}] (* Michael De Vlieger, Dec 06 2022 *)
Formula
a(A000040(n)) = 1.