|
|
A007427
|
|
Moebius transform applied twice to sequence 1,0,0,0,....
(Formerly M0198)
|
|
66
|
|
|
1, -2, -2, 1, -2, 4, -2, 0, 1, 4, -2, -2, -2, 4, 4, 0, -2, -2, -2, -2, 4, 4, -2, 0, 1, 4, 0, -2, -2, -8, -2, 0, 4, 4, 4, 1, -2, 4, 4, 0, -2, -8, -2, -2, -2, 4, -2, 0, 1, -2, 4, -2, -2, 0, 4, 0, 4, 4, -2, 4, -2, 4, -2, 0, 4, -8, -2, -2, 4, -8, -2, 0, -2, 4, -2, -2, 4, -8, -2, 0, 0
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,2
|
|
COMMENTS
|
|a(n)| is the number of ways to write n as a product of 2 squarefree numbers (i.e., number of ways to write n = x*y with 1 <= x <= n, 1 <= y <= n, x and y squarefree). - Benoit Cloitre, Jan 01 2003
|
|
REFERENCES
|
Tom M. Apostol, Introduction to Analytic Number Theory, Springer-Verlag, 1976, page 30.
N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).
|
|
LINKS
|
|
|
FORMULA
|
Dirichlet g.f.: 1/zeta(s)^2.
Multiplicative function with a(p^e) = binomial(2, e)*(-1)^e for p prime and e >= 0.
Dirichlet g.f. of |a(n)|: zeta(s)^2/zeta(2s)^2 (conjectured). - Ralf Stephan, Jul 05 2013. The conjecture is correct because 1+Sum_{e>=1} binomial(2,e)/p^(e*s) = (p^s+1)^2/p^2s, whose product over p is zeta(s)^2/zeta(2s)^2. - Michael Shamos
Sum_{k=1..n} abs(a(k)) ~ (n/zeta(2)^2) * (log(n) + 2*gamma - 1 - 4*zeta'(2)/zeta(2)), where gamma is Euler's constant (A001620). - Amiram Eldar, Dec 24 2023
|
|
EXAMPLE
|
G.f. = x - 2*x^2 - 2*x^3 + x^4 - 2*x^5 + 4*x^6 - 2*x^7 + x^9 + 4*x^10 + ...
We have a(3^1) = C(2, 1)*(-1)^1 = -2, a(3^2) = C(2, 2)*(-1)^2 = 1, and a(3^m) = C(2, m)*(-1)^m = 0 for m >= 3. - Petros Hadjicostas, Jun 07 2019
|
|
MAPLE
|
möbius := proc(a) local b, i, mo: b := NULL:
mo := (m, n) -> `if`(irem(m, n) = 0, numtheory:-mobius(m/n), 0);
for i to nops(a) do b := b, add(mo(i, j)*a[j], j=1..i) od: [b] end:
|
|
MATHEMATICA
|
f[n_] := Plus @@ Times @@@ (MoebiusMu[{#, n/#}] & /@ Divisors@n); Array[f, 105] (* Robert G. Wilson v *)
a[n_] := DivisorSum[n, MoebiusMu[#]*MoebiusMu[n/#]&]; Array[a, 80] (* Jean-François Alcover, Dec 01 2015 *)
|
|
PROG
|
(PARI) {a(n) = if( n<1, 0, direuler(p=2, n, (1 - X)^2)[n])}; /* Michael Somos, Nov 15 2002 */
(PARI) {a(n) = if(n<1, 0, sumdiv(n, d, moebius(d) * moebius(n/d)))}; /* Michael Somos, Nov 15 2002 */
(PARI) a(n)=if(n>1, my(f=factor(n)[, 2], s=sum(i=1, #f, f[i]==1)); if(vecmax(f)>2, 0, (-1)^s<<s), 1) \\ Charles R Greathouse IV, Jun 28 2011
(Haskell)
a007427 n = sum $ zipWith (*) mds $ reverse mds where
mds = a225817_row n
|
|
CROSSREFS
|
|
|
KEYWORD
|
sign,easy,nice,mult
|
|
AUTHOR
|
|
|
EXTENSIONS
|
Added a proof of Stephan's conjecture about the Dirichlet g.f. of |a(n)|.
|
|
STATUS
|
approved
|
|
|
|