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!)
A332466 a(n) = n! * Sum_{d|n} mu(d) / d!. 1
1, 1, 5, 12, 119, 241, 5039, 20160, 302400, 1784161, 39916799, 160332480, 6227020799, 43571848321, 1078831353601, 10461394944000, 355687428095999, 2143016754278400, 121645100408831999, 1196177491129420800, 42565648051390464001, 562000335730215782401 (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!).

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

MAPLE

with(numtheory):

a:= n-> n! * add(mobius(d)/d!, d=divisors(n)):

seq(a(n), n=1..23);  # Alois P. Heinz, Feb 13 2020

MATHEMATICA

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

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

PROG

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

CROSSREFS

Cf. A008683, A057625, A068107, A099740, A132958, A332467.

Sequence in context: A267271 A009754 A096314 * A323565 A195538 A330218

Adjacent sequences:  A332463 A332464 A332465 * A332467 A332468 A332469

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 May 28 12:02 EDT 2020. Contains 334681 sequences. (Running on oeis4.)