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-8 of 8 results.

A054024 Sum of the divisors of n reduced modulo n.

Original entry on oeis.org

0, 1, 1, 3, 1, 0, 1, 7, 4, 8, 1, 4, 1, 10, 9, 15, 1, 3, 1, 2, 11, 14, 1, 12, 6, 16, 13, 0, 1, 12, 1, 31, 15, 20, 13, 19, 1, 22, 17, 10, 1, 12, 1, 40, 33, 26, 1, 28, 8, 43, 21, 46, 1, 12, 17, 8, 23, 32, 1, 48, 1, 34, 41, 63, 19, 12, 1, 58, 27, 4, 1, 51, 1, 40, 49, 64, 19, 12, 1, 26, 40
Offset: 1

Views

Author

Asher Auel, Jan 19 2000

Keywords

Comments

If a(n) = 0, then n is a multiply-perfect number (A007691). - Alonso del Arte, Mar 30 2014

Examples

			a(12) = 4 because sigma(12) = 28 and 28 == 4 (mod 12).
a(13) = 1 because 13 is prime.
a(14) = 10 because sigma(14) = 24 and 24 == 10 (mod 14).
		

Crossrefs

Cf. A000203 (sigma), A005114 (untouchable numbers), A007691 (positions of 0's), A045768, A045769, A088834, A045770, A076496, A159907.

Programs

Formula

a(n) = sigma(n) mod n.
a(p) = 1 for p prime.

A045770 Numbers k such that sigma(k) == 8 (mod k).

Original entry on oeis.org

1, 7, 10, 49, 56, 368, 836, 11096, 17816, 45356, 77744, 91388, 128768, 254012, 388076, 2087936, 2291936, 13174976, 29465852, 35021696, 45335936, 120888092, 260378492, 381236216, 775397948, 3381872252, 4856970752, 6800228816, 8589344768, 44257207676, 114141404156, 1461083549696, 1471763808896, 2199013818368
Offset: 1

Views

Author

Keywords

Comments

Every number of the form 2^(j-1)*(2^j - 9), where 2^j - 9 is prime, is a term. - Jon E. Schoenfield, Jun 02 2019
If m is a term of A045768 with gcd(m,3) = 1 and sigma(m) = 3*q*m + 2 for some integer q, then 3*m is a term of this sequence since sigma(3*m) = 4*q*(3*m) + 8. Some other large terms: 36893488108764397568, 877615520070055755776, 1700388548189538291286016, 85954979333046510417991676, 2081228720695521934665574252544. - Max Alekseyev, May 25 2025

Crossrefs

Programs

  • Maple
    q:= k-> nops(map(x-> x mod k, {8, numtheory[sigma](k)}))=1:
    select(q, [$1..100000])[];  # Alois P. Heinz, Apr 07 2025
  • Mathematica
    Select[Range[1000000], Mod[DivisorSigma[1, #] - 8, #] == 0 &] (* Pontus von Brömssen, Apr 07 2025 *)
  • PARI
    isok(k) = Mod(sigma(k),k) == 8; \\ Pontus von Brömssen, Apr 07 2025

Extensions

a(18)-a(26) from T. D. Noe, Apr 06 2011
Initial term 1 added and a(27)-a(31) from Donovan Johnson, Mar 01 2012
a(32)-a(34) from Giovanni Resta, Apr 02 2014
Term a(2)=7 inserted by Pontus von Brömssen, Apr 07 2025

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

A067702 Numbers k such that sigma(k) == 0 (mod k+2).

Original entry on oeis.org

12, 70, 88, 180, 1888, 4030, 5830, 32128, 521728, 1848964, 8378368, 34359083008, 66072609790, 549753192448, 259708613909470, 2251799645913088, 9223372026117357568
Offset: 1

Views

Author

Benoit Cloitre, Feb 05 2002

Keywords

Comments

If 2^i-5 is prime for i > 2 then let x = (2^i-5)*2^(i-1). Then sigma(x)=2*(x+2), so x is in the sequence. There are other terms that are not of this form. - Jud McCranie, Jan 12 2019
Contains terms of A088832, terms m of A045769 with (sigma(m)-4)/m = 2, terms m of A088834 with (sigma(m)-6)/m = 3, terms m of A045770 with (sigma(m)-8)/m = 4, terms m of A076496 with (sigma(m)-12)/m = 6. - Max Alekseyev, Sep 04 2025

Examples

			sigma(180) = 546 = 3(180+2), so 180 is in the sequence.
		

Crossrefs

Contains subsequence A088832.

Programs

  • Mathematica
    Select[Range[84*10^5],Divisible[DivisorSigma[1,#],#+2]&] (* Harvey P. Dale, May 11 2018 *)
  • PARI
    isok(n) = sigma(n) % (n+2) == 0; \\ Michel Marcus, Nov 22 2013

Extensions

a(9)-a(11) from Michel Marcus, Nov 22 2013
a(12)-a(13) from Jud McCranie, Jan 12 2019
a(14) from Jud McCranie, Jan 13 2019
a(15) from Jud McCranie, Dec 02 2019
a(16)-a(17) from Max Alekseyev, Sep 04 2025

A274551 Numbers k such that sigma(k) == 0 (mod k+3).

Original entry on oeis.org

4, 8925, 32445, 442365
Offset: 1

Views

Author

Paolo P. Lava, Jun 28 2016

Keywords

Comments

a(5) > 10^8 if it exists. - Felix Fröhlich, Jul 01 2016
No more terms < 6.5*10^14. - Jud McCranie, Dec 02 2019
No more terms < 2.7*10^15. - Jud McCranie, Jul 27 2025

Examples

			sigma(4) mod (4+3) = 7 mod 7 = 0.
		

Crossrefs

Programs

  • Magma
    [n: n in [1..2*10^6] | SumOfDivisors(n) mod (n+3) eq 0 ]; // Vincenzo Librandi, Jul 02 2016
  • Mathematica
    Select[Range[10^6], Mod[DivisorSigma[1, #], # + 3] == 0 &] (* Michael De Vlieger, Jul 01 2016 *)
  • PARI
    is(n) = Mod(sigma(n), n+3)==0 \\ Felix Fröhlich, Jul 01 2016
    

A076496 Numbers k such that sigma(k) == 12 (mod k).

Original entry on oeis.org

1, 6, 11, 24, 30, 42, 54, 66, 78, 102, 114, 121, 138, 174, 186, 222, 246, 258, 282, 304, 318, 354, 366, 402, 426, 438, 474, 498, 534, 582, 606, 618, 642, 654, 678, 762, 780, 786, 822, 834, 894, 906, 942, 978, 1002, 1038, 1074, 1086, 1146, 1158, 1182, 1194, 1266, 1338
Offset: 1

Views

Author

Labos Elemer, Oct 21 2002

Keywords

Examples

			6*p is a solution if p > 3 is prime, since sigma(6*p) = 1 + 2 + 3 + 6 + p + 2*p + 3*p + 6*p = 12*(p+1) = 2*6*p + 12 = 2*k + 12. These are "regular" solutions. Also k = 121, 304 are "singular" solutions. See other remainders in cross-references.
		

Crossrefs

Cf. A141545 (a subsequence).

Programs

  • Mathematica
    Select[Range[2000], Mod[DivisorSigma[1, #] - 12, #] == 0 &] (* Vincenzo Librandi, Mar 11 2014, corrected by Amiram Eldar, Jan 04 2023 *)
  • PARI
    isok(k) = Mod(sigma(k), k) == 12; \\ Michel Marcus, Jan 04 2023

Extensions

Initial term 1 added by Vincenzo Librandi, Mar 11 2014
Terms 6 and 11 inserted by Michel Marcus, Jan 04 2023

A045769 Numbers k such that sigma(k) == 4 (mod k).

Original entry on oeis.org

1, 3, 9, 12, 70, 88, 1888, 4030, 5830, 32128, 521728, 1848964, 8378368, 34359083008, 66072609790, 549753192448, 259708613909470, 2251799645913088, 9223372026117357568
Offset: 1

Views

Author

Keywords

Comments

Every number of the form 2^(j-1)*(2^j - 5), where 2^j - 5 is prime, is a term. See A059608. - Jon E. Schoenfield, Jun 02 2019

Crossrefs

Contains subsequence A088832.

Programs

  • PARI
    isok(k) = Mod(sigma(k), k) == 4; \\ Michel Marcus, Jan 04 2023

Extensions

a(13) from Harvey P. Dale, Mar 20 2011
Initial term 1 inserted and a(14)-a(16) from Donovan Johnson, Mar 01 2012
Term 3 inserted by Michel Marcus, Jan 04 2023
a(18) from Jon E. Schoenfield confirmed, and a(17), a(19) added by Max Alekseyev, Jun 08 2025

A301859 Abundant numbers whose abundance is a perfect number.

Original entry on oeis.org

48, 2002, 2632, 4540, 5170, 6952, 8925, 29056, 32445, 32980, 88330, 133042, 174856, 189472, 280228, 442365, 518368, 566752, 892552, 1266952, 2030368, 2052256, 2218450, 3959752, 4120672, 4558936, 5568448, 9071752, 15921112, 38551936, 65969536, 70114936, 88149352, 97364848
Offset: 1

Views

Author

Waldemar Puszkarz, Mar 27 2018

Keywords

Comments

There are 34 terms up to 10^8. The abundance of odd terms (only 3 terms) is 6 (see also A087167). The abundance of even terms is 28, 496, 8128, and 33550336 (for 97364848). There exist deficient numbers whose abundance is a perfect number in absolute terms, e.g., 7, 29, 62.

Examples

			48 is a term as it is abundant and its abundance, sigma(48)-2*48 = 28, is the second perfect number.
		

Crossrefs

Cf. A005101 (abundant numbers), A033880 (abundance), A000396 (perfect numbers), A087167, A088834, A088012, A077374 (sequences related to the odd terms of this sequence).

Programs

  • Mathematica
    Select[Range[10^8], PerfectNumberQ[DivisorSigma[1,# ]-2#]&]
  • PARI
    for(n=1,10^8, a=sigma(n)-2*n; a>0&&sigma(a)==2*a&&print1(n ","))
Showing 1-8 of 8 results.