A053650 - OEIS

A053650

Cototient function of n^2.

0, 2, 3, 8, 5, 24, 7, 32, 27, 60, 11, 96, 13, 112, 105, 128, 17, 216, 19, 240, 189, 264, 23, 384, 125, 364, 243, 448, 29, 660, 31, 512, 429, 612, 385, 864, 37, 760, 585, 960, 41, 1260, 43, 1056, 945, 1104, 47, 1536, 343, 1500, 969, 1456, 53, 1944, 825, 1792, 1197

(list; graph; refs; listen; history; text; internal format)

OFFSET

1,2

COMMENTS

Seems to be invertible like n*Phi(n). Compare with A002618, A038040.

LINKS

Reinhard Zumkeller, Table of n, a(n) for n = 1..10000

FORMULA

a(n) = n*(n - phi(n)) = n^2 - n*phi(n) = Cototient(n^2) = A051953(A000290(n)).

a(n) = n^2 - A002618(n).

For p prime, Cototient(p)=1 and a(p)=p.

a(n) = n*cototient(n) = n*A051953(n). - Omar E. Pol, Nov 22 2012

Dirichlet g.f.: zeta(s-2)*(1 - 1/zeta(s-1)). - Ilya Gutkovskiy, Jul 26 2016

Sum_{k=1..n} a(k) ~ c * n^3 / 3, where c = 1 - 6/Pi^2 (A229099). - Amiram Eldar, Dec 15 2023

MATHEMATICA

Table[n(n-EulerPhi[n]), {n, 60}] (* Michael De Vlieger, Jul 26 2016 *)

PROG

(PARI) a(n) = n^2 - eulerphi(n^2) \\ Michel Marcus, Jul 27 2013

(Haskell)

a053650 = a051953 . a000290 -- Reinhard Zumkeller, Jan 21 2014

(Magma) [n*(n-EulerPhi(n)): n in [1..60]]; // Vincenzo Librandi, Jul 27 2016

(Sage) [n*(n - euler_phi(n)) for n in (1..60)] # G. C. Greubel, May 18 2019

(GAP) List([1..60], n-> n*(n- Phi(n)) ); # G. C. Greubel, May 18 2019

CROSSREFS

Cf. A000005, A038040.

Cf. A001248, A002618, A053650, A053192, A053193, A036689, A229099.

Sequence in context: A062956 A333375 A289667 * A119794 A117987 A376016

Adjacent sequences: A053647 A053648 A053649 * A053651 A053652 A053653

KEYWORD

nonn,look,easy

AUTHOR

Labos Elemer, Feb 18 2000

STATUS

approved