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.

A045768 Numbers k such that sigma(k) == 2 (mod k).

Original entry on oeis.org

1, 20, 104, 464, 650, 1952, 130304, 522752, 8382464, 134193152, 549754241024, 8796086730752, 140737463189504, 144115187270549504
Offset: 1

Views

Author

Keywords

Comments

Equivalently, Chowla function of k is congruent to 1 (mod k).
If p=2^i-3 is prime, then 2^(i-1)*p is a term of the sequence. 650 is in the sequence, but is not of this form.
Terms k from a(2) to a(14) satisfy sigma(k) = 2*k + 2, implying that sigma(k) == 0 (mod k+1). It is not known if this holds in general, for there might be solutions of sigma(k)=3k+2 or 4k+2 or ... (Comments from Jud McCranie and Dean Hickerson, updated by Jon E. Schoenfield, Sep 25 2021 and by Max Alekseyev, May 23 2025).
k | sigma(k) produces the multiperfect numbers (A007691). It is an open question whether k | sigma(k) - 1 iff k is a prime or 1. It is not known if there exist solutions to sigma(k) = 2k+1.
Sequence also gives the nonprime solutions to sigma(k) == 0 (mod k+1), k > 1. - Benoit Cloitre, Feb 05 2002
Sequence seems to give nonprime k such that the numerator of the sum of the reciprocals of the divisors of k equals k+1 (nonprime k such that A017665(k)=k+1). - Benoit Cloitre, Apr 04 2002
For k > 1, composite numbers k such that A108775(k) = floor(sigma(k)/k) = sigma(k) mod k = A054024(k). Complement of primes (A000040) with respect to A230606. There are no numbers k > 2 such that sigma(x) = k*(x+1) has a solution. - Jaroslav Krizek, Dec 05 2013

Examples

			sigma(650) = 1302 == 2 (mod 650).
		

References

  • R. K. Guy, Unsolved Problems in Number Theory, B2.

Crossrefs

Numbers k such that A054013(k)=1.

Programs

  • Mathematica
    Do[If[Mod[DivisorSigma[1, n]-2, n]==0, Print[n]], {n, 1, 10^8}]
    Join[{1}, Select[Range[8000000], Mod[DivisorSigma[1, #], #]==2 &]] (* Vincenzo Librandi, Mar 11 2014 *)
  • PARI
    is(n)=sigma(n)%n==2 || n==1 \\ Charles R Greathouse IV, Mar 09 2014

Extensions

More terms from Jud McCranie, Dec 22 1999.
a(11) from Donovan Johnson, Mar 01 2012
a(12) from Giovanni Resta, Apr 02 2014
a(13) from Jud McCranie, Jun 02 2019
Edited and a(14) from Jon E. Schoenfield confirmed by Max Alekseyev, May 23 2025