cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

Showing 1-2 of 2 results.

A051709 a(n) = sigma(n) + phi(n) - 2n.

Original entry on oeis.org

0, 0, 0, 1, 0, 2, 0, 3, 1, 2, 0, 8, 0, 2, 2, 7, 0, 9, 0, 10, 2, 2, 0, 20, 1, 2, 4, 12, 0, 20, 0, 15, 2, 2, 2, 31, 0, 2, 2, 26, 0, 24, 0, 16, 12, 2, 0, 44, 1, 13, 2, 18, 0, 30, 2, 32, 2, 2, 0, 64, 0, 2, 14, 31, 2, 32, 0, 22, 2, 28, 0, 75, 0, 2, 14, 24, 2, 36, 0
Offset: 1

Views

Author

Keywords

Comments

Sigma is the sum of divisors (A000203), and phi is the Euler totient function (A000010). - Michael B. Porter, Jul 05 2013
Because sigma and phi are multiplicative functions, it is easy to show that (1) if a(n)=0, then n is prime or 1 and (2) if a(n)=2, then n is the product of two distinct prime numbers. Note that a(n) is the n-th term of the Dirichlet series whose generating function is given below. Using the generating function, it is theoretically possible to compute a(n). Hence a(n)=0 could be used as a primality test and a(n)=2 could be used as a test for membership in P2 (A006881). - T. D. Noe, Aug 01 2002
It appears that a(n) - A002033(n) = zeta(s-1) * (zeta(s) - 2 + 1/zeta(s)) + 1/(zeta(s)-2). - Eric Desbiaux, Jul 04 2013
a(n) = 1 if and only if n = prime(k)^2 (n is in A001248). It seems that a(n) = k has only finitely many solutions for k >= 3. - Jianing Song, Jun 27 2021

Examples

			a(5) = sigma(5) + phi(5) - 2*5 = 6 + 4 - 10 = 0.
		

Crossrefs

Cf. A278373 (range of this sequence), A056996 (numbers not present).
Cf. also A344753, A345001 (analogous sequences).

Programs

  • Mathematica
    Table[DivisorSigma[1,n]+EulerPhi[n]-2n,{n,80}] (* Harvey P. Dale, Apr 08 2015 *)
  • PARI
    a(n)=sigma(n)+eulerphi(n)-2*n \\ Charles R Greathouse IV, Jul 05 2013
    
  • PARI
    A051709(n) = -sumdiv(n,d,(dAntti Karttunen, Mar 02 2018

Formula

Dirichlet g.f.: zeta(s-1) * (zeta(s) - 2 + 1/zeta(s)). - T. D. Noe, Aug 01 2002
From Antti Karttunen, Mar 02 2018: (Start)
a(n) = A001065(n) - A051953(n). [Difference between the sum of proper divisors of n and their Moebius-transform.]
a(n) = -Sum_{d|n, dA008683(n/d)*A001065(d). (End)
Sum_{k=1..n} a(k) = (3/(Pi^2) + Pi^2/12 - 1) * n^2 + O(n*log(n)). - Amiram Eldar, Dec 03 2023

A072779 a(n) = sigma_2(n) + phi(n) * sigma(n).

Original entry on oeis.org

2, 8, 18, 35, 50, 74, 98, 145, 169, 202, 242, 322, 338, 394, 452, 589, 578, 689, 722, 882, 884, 970, 1058, 1330, 1271, 1354, 1540, 1722, 1682, 1876, 1922, 2373, 2180, 2314, 2452, 3003, 2738, 2890, 3044, 3650, 3362, 3652, 3698, 4242, 4238, 4234, 4418
Offset: 1

Views

Author

T. D. Noe, Jul 15 2002

Keywords

Comments

This sequence is interesting because (1) a(n) >= 2 n^2, with equality only when n is prime (or 1) and (2) a(n) = 2 + 2*n^2 if and only if n is the product of two distinct primes. Note for twin primes: let n = m^2 - 1, then m-1 and m+1 are twin primes if and only if a(n) = 2 + 2*n^2. Note for the Goldbach conjecture: let n = m^2 - r^2, then m-r and m+r are primes that add to 2m if and only if a(n) = 2 + 2*n^2. See A072780 for a(n) - 2*n^2.

Crossrefs

Programs

  • Haskell
    a072779 n = a001157 n + (a000203 n) * (a000010 n)
    -- Reinhard Zumkeller, Jan 15 2013
    
  • Mathematica
    Table[DivisorSigma[2, n]+EulerPhi[n]DivisorSigma[1, n], {n, 100}]
  • PARI
    a(n)=sigma(n,2)+eulerphi(n)*sigma(n) \\ Charles R Greathouse IV, May 15 2013

Formula

a(n) = A001157(n) + A000203(n)*A000010(n). - Reinhard Zumkeller, Jan 15 2013
Sum_{k=1..n} a(k) ~ c * n^3 / 3, where c = zeta(3) + Product_{p prime} (1 - 1/(p^2*(p+1))) = A002117 + A065465 = 2.083570742884... . - Amiram Eldar, Dec 03 2023
Showing 1-2 of 2 results.