|
|
A271102
|
|
a(n) is multiplicative with a(p^e) = -1 if e=2, a(p^e) = 0 if e=1 or e>2.
|
|
12
|
|
|
1, 0, 0, -1, 0, 0, 0, 0, -1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, -1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, -1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1
|
|
LINKS
|
|
|
FORMULA
|
Dirichlet g.f.: 1 / zeta(2*s).
|
|
EXAMPLE
|
G.f. = x - x^4 - x^9 - x^25 + x^36 - x^49 + x^100 - x^121 - x^169 + ...
|
|
MATHEMATICA
|
Table[Times @@ Apply[Times, FactorInteger[n] /. {p_, e_} /; p > 1 :> If[e == 2, -1, 0]], {n, 105}] (* Michael De Vlieger, Jul 29 2017 *)
Table[DivisorSum[n, Abs[MoebiusMu[#]]*MoebiusMu[n/#] &], {n, 1, 100}] (* Vaclav Kotesovec, Apr 08 2019 *)
|
|
PROG
|
(PARI) {a(n) = if( n<1, 0, direuler( p=2, n, 1 - X^2 )[n])};
(PARI) {a(n) = my(A, p, e); if( n<1, 0, A = factor(n); prod(k=1, matsize(A)[1], [p, e] = A[k, ]; -(e==2)))};
|
|
CROSSREFS
|
|
|
KEYWORD
|
sign,mult
|
|
AUTHOR
|
|
|
EXTENSIONS
|
|
|
STATUS
|
approved
|
|
|
|