login
The OEIS is supported by the many generous donors to the OEIS Foundation.

 

Logo
Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A074722 a(n) = Sum_{d divides n} phi(n/d)*(-1)^bigomega(d). 4

%I #33 Aug 28 2023 08:25:52

%S 1,0,1,2,3,0,5,2,5,0,9,2,11,0,3,6,15,0,17,6,5,0,21,2,17,0,13,10,27,0,

%T 29,10,9,0,15,10,35,0,11,6,39,0,41,18,15,0,45,6,37,0,15,22,51,0,27,10,

%U 17,0,57,6,59,0,25,22,33,0,65,30,21,0,69,10,71,0,17,34,45,0,77,18,41,0

%N a(n) = Sum_{d divides n} phi(n/d)*(-1)^bigomega(d).

%C a(n) = 0 if and only if n == 2 (mod 4). - _Robert Israel_, Jan 04 2017

%H Robert Israel, <a href="/A074722/b074722.txt">Table of n, a(n) for n = 1..10000</a>

%F Multiplicative with a(p^e) = 2*(-1)^(e+1)*((-p)^(e+1)-1)/(p+1)-p^e.

%F Dirichlet g.f.: zeta(2s)*zeta(s-1)/(zeta(s)^2). - _Benedict W. J. Irwin_, Jul 11 2018

%F Sum_{k=1..n} a(k) ~ n^2 / 5. - _Vaclav Kotesovec_, Feb 01 2019

%F a(n) = Sum_{k=1..n} (-1)^bigomega(gcd(n,k)). - _Ilya Gutkovskiy_, Feb 22 2020

%F Möbius transform of A206369: a(n) = Sum_{d|n} A008683(d) * A206369(n/d). - _Amiram Eldar_, Aug 28 2023

%p f:= proc(n) uses numtheory; local d;

%p add(phi(n/d)*(-1)^bigomega(d), d=divisors(n))

%p end proc:

%p map(f, [$1..100]); # _Robert Israel_, Jan 04 2017

%t f[d_] := EulerPhi[n/d] LiouvilleLambda[d]

%t Table[DivisorSum[n, f], {n, 1, 50}] (* _Benedict W. J. Irwin_, Jul 11 2018 *)

%t f[p_, e_] := 2*(-1)^(e + 1)*((-p)^(e + 1) - 1)/(p + 1) - p^e; a[1] = 1; a[n_] := Times @@ f @@@ FactorInteger[n]; Array[a, 100] (* _Amiram Eldar_, Oct 30 2022 *)

%o (PARI) a(n) = sumdiv(n, d, eulerphi(n/d)*(-1)^bigomega(d)); \\ _Michel Marcus_, Jul 11 2018

%Y Cf. A000010, A001222, A008683, A058026, A206369.

%K nonn,easy,mult

%O 1,4

%A _Vladeta Jovovic_, Sep 27 2002

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recents
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified March 28 10:55 EDT 2024. Contains 371241 sequences. (Running on oeis4.)