|
|
A130054
|
|
Inverse Moebius transform of A023900.
|
|
5
|
|
|
1, 0, -1, -1, -3, 0, -5, -2, -3, 0, -9, 1, -11, 0, 3, -3, -15, 0, -17, 3, 5, 0, -21, 2, -7, 0, -5, 5, -27, 0, -29, -4, 9, 0, 15, 3, -35, 0, 11, 6, -39, 0, -41, 9, 9, 0, -45, 3, -11, 0, 15, 11, -51, 0, 27, 10, 17, 0, -57, -3, -59, 0, 15, -5, 33, 0, -65, 15, 21
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,5
|
|
COMMENTS
|
|
|
LINKS
|
|
|
FORMULA
|
A126988 * A130054 = d(n), A000005: (1, 2, 2, 3, 2, 4, 2, 4, 3, 4, ...).
Multiplicative with a(p^e) = 1 - (p-1) * e for prime p and e >= 0.
Dirichlet g.f.: (zeta(s))^2 / zeta(s-1).
Dirichlet inverse of A007431. (End)
|
|
MAPLE
|
with(numtheory): seq(add(d*mobius(d)*tau(n/d), d in divisors(n)), n=1..60); # Ridouane Oudra, Nov 17 2019
|
|
MATHEMATICA
|
b[n_] := Sum[d MoebiusMu[d], {d, Divisors[n]}];
a[n_] := Sum[b[n/d], {d, Divisors[n]}];
f[p_, e_] := 1-(p-1)*e; a[1] = 1; a[n_] := Times @@ f @@@ FactorInteger[n]; Array[a, 100] (* Amiram Eldar, Nov 23 2020 *)
|
|
PROG
|
b(n)={sumdivmult(n, d, d*moebius(d))}
(Magma) [&+[d*MoebiusMu(d)*NumberOfDivisors(n div d):d in Divisors(n)]:n in [1..70]]; // Marius A. Burtea, Nov 17 2019
|
|
CROSSREFS
|
|
|
KEYWORD
|
sign,mult,easy
|
|
AUTHOR
|
|
|
EXTENSIONS
|
Name changed and terms a(11) and beyond from Andrew Howroyd, Aug 03 2018
|
|
STATUS
|
approved
|
|
|
|