|
%I #39 May 23 2019 08:56:31
%S 1,0,0,2,0,2,0,2,0,2,0,4,0,2,0,2,0,2,0,4,0,2,0,4,0,2,0,4,0,4,0,2,0,2,
%T 0,4,0,2,0,4,0,4,0,4,0,2,0,4,0,2,0,4,0,2,0,4,0,2,0,8,0,2,0,2,0,4,0,4,
%U 0,4,0,4,0,2,0,4,0,4,0,4,0,2,0,8,0,2,0
%N a(n) = Sum_{d|n} möbius(d)^(n/d).
%H Reinhard Zumkeller, <a href="/A264782/b264782.txt">Table of n, a(n) for n = 1..10000</a>
%F a(n) = Sum_{d|n} möbius(d)^(n/d).
%F For odd n, a(n)=0.
%F For n = 2 * p1^k1 * p2^k2 * ... * pr^kr, a(n) = 2^r.
%F For n = 2^m * p1^k1 * p2^k2 * ... * pr^kr, a(n) = 2^(r+1) if m > 1.
%F a(2n) = A034444(n) for n > 1.
%F From _Gevorg Hmayakyan_, Dec 31 2016: (Start)
%F If b(n) = Sum_{d|n} möbius(d)^d, then b(n) = (A209229(n)+1)*((-1)^n + 1)/2*a(2*n)/2, for n > 1.
%F Dirichlet g.f.: -1 + 2^(-s) + (2^(-s) Zeta[s]^2)/Zeta[2s]. (End)
%F Sum_{k=1..n} a(k) ~ 3*n / Pi^2 * (log(n) - 1 + 2*gamma - log(2) - 12*Zeta'(2)/Pi^2), where gamma is the Euler-Mascheroni constant A001620. - _Vaclav Kotesovec_, Feb 02 2019
%F G.f.: Sum_{k>=1} mu(k)*x^k/(1 - mu(k)*x^k). - _Ilya Gutkovskiy_, May 23 2019
%e a(1) = 1, a(p) = mu(1)^p + mu(p)^1 = 0.
%e a(p1*p2) = mu(1)^p1*p2 + mu(p1)^p2 + mu(p2)^p1 + mu(p1*p2) = 1+(-1)+(-1)+1 = 0.
%e a(2*p) = mu(1)^2*p + mu(2)^p + mu(p)^2 + mu(2*p) = 1+(-1)+1+1 = 2.
%t Table[Sum[MoebiusMu[d]^(n/d), {d, Divisors@ n}], {n, 87}] (* _Michael De Vlieger_, Nov 25 2015 *)
%o (PARI) a(n) = sumdiv(n, d, moebius(n/d)^d);
%o (Haskell)
%o a264782 n = sum $ zipWith (^) (map a008683 divs) (reverse divs)
%o where divs = a027750_row n
%o -- _Reinhard Zumkeller_, Dec 19 2015
%o (Perl) use ntheory ":all"; sub a264782 { my $n=shift; divisor_sum($n, sub { moebius($_[0]) ** ($n/$_[0]) }); } # _Dana Jacobsen_, Dec 29 2015
%Y Cf. A008683, A027750, A034444.
%K nonn,easy
%O 1,4
%A _Gevorg Hmayakyan_, Nov 24 2015
|
