A354022 a(n) = n! * Sum_{d|n} mu(n/d) / d!.

1, -1, -5, -11, -119, 241, -5039, -1679, -60479, 1784161, -39916799, 218877121, -6227020799, 43571848321, 1078831353601, -518918399, -355687428095999, 1058152455360001, -121645100408831999, 1115079416638387201, 42565648051390464001, 562000335730215782401

Offset: 1

Author

Ilya Gutkovskiy, May 14 2022

Keywords

• sign

Links

• Winston de Greef, Table of n, a(n) for n = 1..448

Crossrefs

Cf. A008683, A057625, A062794, A068107, A099740, A132958, A327243, A332466.

Programs

Mathematica

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

PARI

a(n)=n! * sumdiv(n, d, moebius(n/d) / d!) \\ Winston de Greef, Sep 19 2023

Formula

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

Sum_{n>=1} a(n) * x^n / (n! * (1 - x^n)) = exp(x) - 1.