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-3 of 3 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

A344995 Numbers k such that A051709(k)/A173557(k) is a positive natural number and a divisor of k.

Original entry on oeis.org

4, 6, 12, 28, 42, 312, 456, 496, 588, 828, 1080, 1216, 1242, 1377, 1560, 4560, 5964, 7320, 7480, 7584, 8128, 11400, 13728, 14784, 23760, 33462, 59400, 59520, 124020, 147840, 188600, 277648, 321000, 543552, 1288224, 1510272, 1596048, 1964544, 2038140, 3323736, 3611520, 3780672, 3909816, 6137440, 9034032, 10783890
Offset: 1

Views

Author

Antti Karttunen, Jun 05 2021

Keywords

Crossrefs

Subsequence of A344994.
Cf. also A344755, A345002.

Programs

  • PARI
    A051709(n) = ((sigma(n) + eulerphi(n)) - (2*n));
    A173557(n) = factorback(apply(p -> p-1, factor(n)[, 1]));
    isA344995(n) = { my(u=A051709(n),t=A173557(n),r=u/t); ((u>0)&&(1==denominator(r)&&!(n%r))); };

A345054 Odd numbers k such that A173557(k) divides nonzero A051709(k).

Original entry on oeis.org

27, 243, 1377, 2187, 3125, 19683, 28125, 55233, 68445, 177147, 195625, 203125, 239805, 253125, 453125, 823543, 907137, 1323297, 1378125, 1464561, 1594323, 1953125, 2278125, 3341637, 3572829, 5255361, 5877117, 9034497, 9819837, 11701053, 14348907, 17453125, 19460393, 20503125, 22209633, 26010621, 30074733, 44910045
Offset: 1

Views

Author

Antti Karttunen, Jun 06 2021

Keywords

Comments

Question: Are there any common terms with A345051?

Crossrefs

Odd terms in A344994.

Programs

  • PARI
    A051709(n) = ((sigma(n) + eulerphi(n)) - (2*n));
    A173557(n) = factorback(apply(p -> p-1, factor(n)[, 1]));
    isA345054(n) = if(!(n%2),0,my(u=A051709(n)); ((u>0)&&(0==(u%A173557(n)))));
Showing 1-3 of 3 results.