A069249 a(n) = n^2 - phi(n)*sigma(n).

0, 1, 1, 2, 1, 12, 1, 4, 3, 28, 1, 32, 1, 52, 33, 8, 1, 90, 1, 64, 57, 124, 1, 96, 5, 172, 9, 112, 1, 324, 1, 16, 129, 292, 73, 204, 1, 364, 177, 160, 1, 612, 1, 256, 153, 532, 1, 320, 7, 640, 297, 352, 1, 756, 145, 256, 369, 844, 1, 912, 1, 964, 225, 32, 193, 1476, 1, 592

Offset

1,4

Comments

Always >0 for n>0. a(n)=1 if n is prime.

If p is a prime and k is a natural number then a(p^k)=p^(k-1) because a(p^k)=(p^k)^2-sigma(p^k)*phi(p^k) =p^(2k)-(p-1)*p^(k-1)*(p^(k+1)-1)/(p-1)=p^(k-1). If n is a composite number then a(n)>1 and a(1)=0, so n is prime iff a(n)=1. - Farideh Firoozbakht, Nov 15 2005

Links

T. D. Noe, Table of n, a(n) for n = 1..10000

Formula

a(n) = n^2-A062354(n). - R. J. Mathar, Oct 01 2011

Sum_{k=1..n} a(k) ~ c * n^3 / 3, where c = 1 - A065465 = 0.118486... . - Amiram Eldar, Dec 04 2023

Example

sigma(10) = 18; phi(10) = 4; 10^2 - sigma(10)*phi(10) = 28. sigma(p) = p+1; phi(p) = p-1; p^2 - (p+1)(p-1) = 1. [From Walter Nissen, Aug 29 2009]

Mathematica

Table[n^2-EulerPhi[n]DivisorSigma[1, n], {n, 70}] (* Harvey P. Dale, Oct 22 2016 *)

Prog

(PARI) a(n)=n^2-eulerphi(n)*sigma(n) \\ Charles R Greathouse IV, Nov 27 2013

Crossrefs

Cf. A000010, A000203, A065465, A164875, A164876.

Sequence in context: A005730 A112284 A167401 * A128247 A161150 A163088

Adjacent sequences: A069246 A069247 A069248 * A069250 A069251 A069252

Keyword

easy,nonn

Author

Benoit Cloitre, Apr 13 2002

Status

approved