login

Reminder: The OEIS is hiring a new managing editor, and the application deadline is January 26.

A363552
Möbius function of rank 4: a(n) = lambda(n) = A008836(n) if n is biquadratefree (A046100) and 0 otherwise.
3
1, -1, -1, 1, -1, 1, -1, -1, 1, 1, -1, -1, -1, 1, 1, 0, -1, -1, -1, -1, 1, 1, -1, 1, 1, 1, -1, -1, -1, -1, -1, 0, 1, 1, 1, 1, -1, 1, 1, 1, -1, -1, -1, -1, -1, 1, -1, 0, 1, -1, 1, -1, -1, 1, 1, 1, 1, 1, -1, 1, -1, 1, -1, 0, 1, -1, -1, -1, 1, -1, -1, -1, -1, 1, -1
OFFSET
1
COMMENTS
Dirichlet inverse of A353519.
LINKS
Masato Kobayashi, Möbius functions of higher rank and Dirichlet series, arXiv:2108.01822 [math.NT], 2021.
FORMULA
a(n) = A008836(n) * A307430(n).
abs(a(n)) = A307430(n).
Multiplicative with a(p^e) = (-1)^e if e <= 3, and 0 otherwise.
Dirichlet g.f.: Product_{p prime} (1 - 1/p^s + 1/p^(2*s) - 1/p^(3*s)) = zeta(2*s)/(zeta(s)*zeta(4*s)).
MATHEMATICA
f[p_, e_] := If[e < 4, (-1)^e, 0]; a[1] = 1; a[n_] := Times @@ f @@@ FactorInteger[n]; Array[a, 100]
PROG
(PARI) f(e) = if(e < 4, (-1)^e, 0);
a(n) = vecprod(apply(f, factor(n)[, 2]));
CROSSREFS
Other generalizations of the Möbius function: A053864, A053865, A053981, A189021, A189022, A189023.
Sequence in context: A167851 A340369 A307445 * A307430 A053865 A189022
KEYWORD
sign,mult,easy
AUTHOR
Amiram Eldar, Jun 10 2023
STATUS
approved