login
The OEIS Foundation is supported by donations from users of the OEIS and by a grant from the Simons Foundation.

 

Logo


Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A332467 a(n) = n! * Sum_{d|n} mu(d) / (d!)^(n/d). 1
1, 1, 5, 18, 119, 611, 5039, 37800, 361200, 3515149, 39916799, 471148524, 6227020799, 86497207369, 1307505443245, 20841060240000, 355687428095999, 6389731861649136, 121645100408831999, 2430526115576719732, 51090759661943327041, 1123451899297246814569 (list; graph; refs; listen; history; text; internal format)
OFFSET

1,3

LINKS

Table of n, a(n) for n=1..22.

FORMULA

E.g.f.: Sum_{k>=1} Sum_{j>=1} mu(j) * x^(k*j) / (j!)^k).

E.g.f.: Sum_{k>=1} mu(k) * x^k / (k! - x^k).

a(n) ~ n!. - Vaclav Kotesovec, Feb 16 2020

MATHEMATICA

Table[n! DivisorSum[n, MoebiusMu[#]/(#!)^(n/#) &], {n, 1, 22}]

nmax = 22; CoefficientList[Series[Sum[MoebiusMu[k] x^k/(k! - x^k), {k, 1, nmax}], {x, 0, nmax}], x] Range[0, nmax]! // Rest

PROG

(PARI) a(n)={sumdiv(n, d, moebius(d)*n!/(d!)^(n/d))} \\ Andrew Howroyd, Feb 13 2020

(MAGMA) [Factorial(n)*&+[MoebiusMu(d) /(Factorial(d))^(n div d):d in Divisors(n)]:n in [1..22]]; // Marius A. Burtea, Feb 13 2020

CROSSREFS

Cf. A008683, A061095, A327587, A332466.

Sequence in context: A009442 A332784 A302435 * A036688 A009348 A009365

Adjacent sequences:  A332464 A332465 A332466 * A332468 A332469 A332470

KEYWORD

nonn

AUTHOR

Ilya Gutkovskiy, Feb 13 2020

STATUS

approved

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recent
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified July 8 04:16 EDT 2020. Contains 335504 sequences. (Running on oeis4.)