|
|
A307420
|
|
Dirichlet g.f.: zeta(2*s) * zeta(3*s) / zeta(s).
|
|
3
|
|
|
1, -1, -1, 1, -1, 1, -1, 0, 1, 1, -1, -1, -1, 1, 1, 0, -1, -1, -1, -1, 1, 1, -1, 0, 1, 1, 0, -1, -1, -1, -1, 0, 1, 1, 1, 1, -1, 1, 1, 0, -1, -1, -1, -1, -1, 1, -1, 0, 1, -1, 1, -1, -1, 0, 1, 0, 1, 1, -1, 1, -1, 1, -1, 1, 1, -1, -1, -1, 1, -1, -1, 0, -1, 1, -1
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1
|
|
COMMENTS
|
|
|
LINKS
|
|
|
FORMULA
|
Multiplicative with a(p^e) = 1 if e == 0 or 2 (mod 6), -1 if e == 1 (mod 6), and 0 otherwise. - Amiram Eldar, Dec 25 2022
|
|
MATHEMATICA
|
Table[DivisorSum[n, Boole[IntegerQ[#^(1/3)]] * LiouvilleLambda[n/#]&], {n, 1, 100}]
Table[DivisorSum[n, Mod[DivisorSigma[0, #], 3, -1] * Boole[IntegerQ[Sqrt[n/#]]] &], {n, 1, 100}]
f[p_, e_] := Switch[Mod[e, 6], 0, 1, 1, -1, 2, 1, 3, 0, 4, 0, 5, 0]; a[1] = 1; a[n_] := Times @@ f @@@ FactorInteger[n]; Array[a, 100] (* Amiram Eldar, Dec 25 2022 *)
|
|
PROG
|
(PARI) for(n=1, 100, print1(direuler(p=2, n, 1/((1+X)*(1-X^3)))[n], ", ")) \\ Vaclav Kotesovec, Jun 14 2020
|
|
CROSSREFS
|
|
|
KEYWORD
|
sign,mult
|
|
AUTHOR
|
|
|
STATUS
|
approved
|
|
|
|