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-10 of 31 results. Next

A000043 Mersenne exponents: primes p such that 2^p - 1 is prime. Then 2^p - 1 is called a Mersenne prime.

Original entry on oeis.org

2, 3, 5, 7, 13, 17, 19, 31, 61, 89, 107, 127, 521, 607, 1279, 2203, 2281, 3217, 4253, 4423, 9689, 9941, 11213, 19937, 21701, 23209, 44497, 86243, 110503, 132049, 216091, 756839, 859433, 1257787, 1398269, 2976221, 3021377, 6972593, 13466917, 20996011, 24036583, 25964951, 30402457, 32582657, 37156667, 42643801, 43112609, 57885161, 74207281
Offset: 1

Views

Author

Keywords

Comments

Equivalently, integers k such that 2^k - 1 is prime.
It is believed (but unproved) that this sequence is infinite. The data suggest that the number of terms up to exponent N is roughly K log N for some constant K.
Length of prime repunits in base 2.
The associated perfect number N=2^(p-1)*M(p) (=A019279*A000668=A000396), has 2p (=A061645) divisors with harmonic mean p (and geometric mean sqrt(N)). - Lekraj Beedassy, Aug 21 2004
In one of his first publications Euler found the numbers up to 31 but erroneously included 41 and 47.
Equals number of bits in binary expansion of n-th Mersenne prime (A117293). - Artur Jasinski, Feb 09 2007
Number of divisors of n-th even perfect number, divided by 2. Number of divisors of n-th even perfect number that are powers of 2. Number of divisors of n-th even perfect number that are multiples of n-th Mersenne prime A000668(n). - Omar E. Pol, Feb 24 2008
Number of divisors of n-th even superperfect number A061652(n). Numbers of divisors of n-th superperfect number A019279(n), assuming there are no odd superperfect numbers. - Omar E. Pol, Mar 01 2008
Differences between exponents when the even perfect numbers are represented as differences of powers of 2, for example: The 5th even perfect number is 33550336 = 2^25 - 2^12 then a(5)=25-12=13 (see A135655, A133033, A090748). - Omar E. Pol, Mar 01 2008
Number of 1's in binary expansion of n-th even perfect number (see A135650). Number of 1's in binary expansion of divisors of n-th even perfect number that are multiples of n-th Mersenne prime A000668(n) (see A135652, A135653, A135654, A135655). - Omar E. Pol, May 04 2008
Indices of the numbers A006516 that are also even perfect numbers. - Omar E. Pol, Aug 30 2008
Indices of Mersenne numbers A000225 that are also Mersenne primes A000668. - Omar E. Pol, Aug 31 2008
The (prime) number p appears in this sequence if and only if there is no prime q<2^p-1 such that the order of 2 modulo q equals p; a special case is that if p=4k+3 is prime and also q=2p+1 is prime then the order of 2 modulo q is p so p is not a term of this sequence. - Joerg Arndt, Jan 16 2011
Primes p such that sigma(2^p) - sigma(2^p-1) = 2^p-1. - Jaroslav Krizek, Aug 02 2013
Integers k such that every degree k irreducible polynomial over GF(2) is also primitive, i.e., has order 2^k-1. Equivalently, the integers k such that A001037(k) = A011260(k). - Geoffrey Critzer, Dec 08 2019
Conjecture: for k > 1, 2^k-1 is (a Mersenne) prime or k = 2^(2^m)+1 (is a Fermat number) if and only if (k-1)^(2^k-2) == 1 (mod (2^k-1)k^2). - Thomas Ordowski, Oct 05 2023
Conjecture: for p prime, 2^p-1 is (a Mersenne) prime or p = 2^(2^m)+1 (is a Fermat number) if and only if (p-1)^(2^p-2) == 1 (mod 2^p-1). - David Barina, Nov 25 2024
Already as of Dec. 2020, all exponents up to 10^8 had been verified, implying that 74207281, 77232917 and 82589933 are indeed the next three terms. As of today, all exponents up to 130439863 have been tested at least once, see the GIMPS Milestones Report. - M. F. Hasler, Apr 11 2025
On June 23. 2025 all exponents up to 74340751 have been verified, confirming that 74207281 is the exponent of the 49th Mersenne Prime. - Rodolfo Ruiz-Huidobro, Jun 23 2025

Examples

			Corresponding to the initial terms 2, 3, 5, 7, 13, 17, 19, 31 ... we get the Mersenne primes 2^2 - 1 = 3, 2^3 - 1 = 7, 2^5 - 1 = 31, 127, 8191, 131071, 524287, 2147483647, ... (see A000668).
		

References

  • T. M. Apostol, Introduction to Analytic Number Theory, Springer-Verlag, 1976, page 4.
  • J. Brillhart et al., Factorizations of b^n +- 1. Contemporary Mathematics, Vol. 22, Amer. Math. Soc., Providence, RI, 2nd edition, 1985; and later supplements.
  • Jan Gullberg, Mathematics from the Birth of Numbers, W. W. Norton & Co., NY & London, 1997, §3.2 Prime Numbers, p. 79.
  • R. K. Guy, Unsolved Problems in Number Theory, Section A3.
  • F. Lemmermeyer, Reciprocity Laws From Euler to Eisenstein, Springer-Verlag, 2000, p. 57.
  • Clifford A. Pickover, A Passion for Mathematics, Wiley, 2005; see p. 19.
  • Alfred S. Posamentier, Math Charmers, Tantalizing Tidbits for the Mind, Prometheus Books, NY, 2003, page 47.
  • N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence).
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).
  • James J. Tattersall, Elementary Number Theory in Nine Chapters, Cambridge University Press, 1999, pages 132-134.
  • B. Tuckerman, The 24th Mersenne prime, Notices Amer. Math. Soc., 18 (Jun, 1971), Abstract 684-A15, p. 608.

Crossrefs

Cf. A000668 (Mersenne primes).
Cf. A028335 (integer lengths of Mersenne primes).
Cf. A000225 (Mersenne numbers).
Cf. A001348 (Mersenne numbers with n prime).

Programs

  • Mathematica
    MersennePrimeExponent[Range[48]] (* Eric W. Weisstein, Jul 17 2017; updated Oct 21 2024 *)
  • PARI
    isA000043(n) = isprime(2^n-1) \\ Michael B. Porter, Oct 28 2009
    
  • PARI
    is(n)=my(h=Mod(2,2^n-1)); for(i=1, n-2, h=2*h^2-1); h==0||n==2 \\ Lucas-Lehmer test for exponent e. - Joerg Arndt, Jan 16 2011, and Charles R Greathouse IV, Jun 05 2013
    forprime(e=2,5000,if(is(e),print1(e,", "))); /* terms < 5000 */
    
  • Python
    from sympy import isprime, prime
    for n in range(1,100):
        if isprime(2**prime(n)-1):
            print(prime(n), end=', ') # Stefano Spezia, Dec 06 2018

Formula

a(n) = log((1/2)*(1+sqrt(1+8*A000396(n))))/log(2). - Artur Jasinski, Sep 23 2008 (under the assumption there are no odd perfect numbers, Joerg Arndt, Feb 23 2014)
a(n) = A000005(A061652(n)). - Omar E. Pol, Aug 26 2009
a(n) = A000120(A000396(n)), assuming there are no odd perfect numbers. - Omar E. Pol, Oct 30 2013

Extensions

Also in the sequence: p = 74207281. - Charles R Greathouse IV, Jan 19 2016
Also in the sequence: p = 77232917. - Eric W. Weisstein, Jan 03 2018
Also in the sequence: p = 82589933. - Gord Palameta, Dec 21 2018
a(46) = 42643801 and a(47) = 43112609, whose ordinal positions in the sequence are now confirmed, communicated by Eric W. Weisstein, Apr 12 2018
a(48) = 57885161, whose ordinal position in the sequence is now confirmed, communicated by Benjamin Przybocki, Jan 05 2022
Also in the sequence: p = 136279841. - Eric W. Weisstein, Oct 21 2024
As of Jan 31 2025, 48 terms are known, and are shown in the DATA section. Four additional numbers are known to be in the sequence, namely 74207281, 77232917, 82589933, and 136279841, but they may not be the next terms. See the GIMP website for the latest information. - N. J. A. Sloane, Jan 31 2025

A000396 Perfect numbers k: k is equal to the sum of the proper divisors of k.

Original entry on oeis.org

6, 28, 496, 8128, 33550336, 8589869056, 137438691328, 2305843008139952128, 2658455991569831744654692615953842176, 191561942608236107294793378084303638130997321548169216
Offset: 1

Views

Author

Keywords

Comments

A number k is abundant if sigma(k) > 2k (cf. A005101), perfect if sigma(k) = 2k (this sequence), or deficient if sigma(k) < 2k (cf. A005100), where sigma(k) is the sum of the divisors of k (A000203).
The numbers 2^(p-1)*(2^p - 1) are perfect, where p is a prime such that 2^p - 1 is also prime (for the list of p's see A000043). There are no other even perfect numbers and it is believed that there are no odd perfect numbers.
Numbers k such that Sum_{d|k} 1/d = 2. - Benoit Cloitre, Apr 07 2002
For number of divisors of a(n) see A061645(n). Number of digits in a(n) is A061193(n). - Lekraj Beedassy, Jun 04 2004
All terms other than the first have digital root 1 (since 4^2 == 4 (mod 6), we have, by induction, 4^k == 4 (mod 6), or 2*2^(2*k) = 8 == 2 (mod 6), implying that Mersenne primes M = 2^p - 1, for odd p, are of the form 6*t+1). Thus perfect numbers N, being M-th triangular, have the form (6*t+1)*(3*t+1), whence the property N mod 9 = 1 for all N after the first. - Lekraj Beedassy, Aug 21 2004
The earliest recorded mention of this sequence is in Euclid's Elements, IX 36, about 300 BC. - Artur Jasinski, Jan 25 2006
Theorem (Euclid, Euler). An even number m is a perfect number if and only if m = 2^(k-1)*(2^k-1), where 2^k-1 is prime. Euler's idea came from Euclid's Proposition 36 of Book IX (see Weil). It follows that every even perfect number is also a triangular number. - Mohammad K. Azarian, Apr 16 2008
Triangular numbers (also generalized hexagonal numbers) A000217 whose indices are Mersenne primes A000668, assuming there are no odd perfect numbers. - Omar E. Pol, May 09 2008, Sep 15 2013
If a(n) is even, then 2*a(n) is in A181595. - Vladimir Shevelev, Nov 07 2010
Except for a(1) = 6, all even terms are of the form 30*k - 2 or 45*k + 1. - Arkadiusz Wesolowski, Mar 11 2012
a(4) = A229381(1) = 8128 is the "Simpsons's perfect number". - Jonathan Sondow, Jan 02 2015
Theorem (Farideh Firoozbakht): If m is an integer and both p and p^k-m-1 are prime numbers then x = p^(k-1)*(p^k-m-1) is a solution to the equation sigma(x) = (p*x+m)/(p-1). For example, if we take m=0 and p=2 we get Euclid's result about perfect numbers. - Farideh Firoozbakht, Mar 01 2015
The cototient of the even perfect numbers is a square; in particular, if 2^p - 1 is a Mersenne prime, cototient(2^(p-1) * (2^p - 1)) = (2^(p-1))^2 (see A152921). So, this sequence is a subsequence of A063752. - Bernard Schott, Jan 11 2019
Euler's (1747) proof that all the even perfect number are of the form 2^(p-1)*(2^p-1) implies that their asymptotic density is 0. Kanold (1954) proved that the asymptotic density of odd perfect numbers is 0. - Amiram Eldar, Feb 13 2021
If k is perfect and semiprime, then k = 6. - Alexandra Hercilia Pereira Silva, Aug 30 2021
This sequence lists the fixed points of A001065. - Alois P. Heinz, Mar 10 2024

Examples

			6 is perfect because 6 = 1+2+3, the sum of all divisors of 6 less than 6; 28 is perfect because 28 = 1+2+4+7+14.
		

References

  • Tom M. Apostol, Introduction to Analytic Number Theory, Springer-Verlag, 1976, page 4.
  • Albert H. Beiler, Recreations in the Theory of Numbers, Dover, NY, 2d ed. 1966, pp. 11-23.
  • Stanley J. Bezuszka, Perfect Numbers (Booklet 3, Motivated Math. Project Activities), Boston College Press, Chestnut Hill MA, 1980.
  • John H. Conway and Richard K. Guy, The Book of Numbers, New York: Springer-Verlag, 1996. See pp. 136-137.
  • Euclid, Elements, Book IX, Section 36, about 300 BC.
  • Jan Gullberg, Mathematics from the Birth of Numbers, W. W. Norton & Co., NY & London, 1997, §3.3 Perfect and Amicable Numbers, pp. 82-83.
  • R. K. Guy, Unsolved Problems in Number Theory, Springer, 1st edition, 1981. See section B1.
  • G. H. Hardy and E. M. Wright, An Introduction to the Theory of Numbers. 3rd ed., Oxford Univ. Press, 1954, p. 239.
  • T. Koshy, "The Ends Of A Mersenne Prime And An Even Perfect Number", Journal of Recreational Mathematics, Baywood, NY, 1998, pp. 196-202.
  • Joseph S. Madachy, Madachy's Mathematical Recreations, New York: Dover Publications, Inc., 1979, p. 149 (First publ. by Charles Scribner's Sons, New York, 1966, under the title: Mathematics on Vacation).
  • Alfred S. Posamentier, Math Charmers, Tantalizing Tidbits for the Mind, Prometheus Books, NY, 2003, pages 46-48, 244-245.
  • Paulo Ribenboim, The Little Book of Bigger Primes, Springer-Verlag NY 2004. See pp. 83-87.
  • József Sándor and Borislav Crstici, Handbook of Number Theory, II, Springer Verlag, 2004.
  • N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence).
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).
  • Ian Stewart, L'univers des nombres, "Diviser Pour Régner", Chapter 14, pp. 74-81, Belin-Pour La Science, Paris, 2000.
  • James J. Tattersall, Elementary Number Theory in Nine Chapters, Cambridge University Press, 1999, chapter 4, pages 127-149.
  • Horace S. Uhler, On the 16th and 17th perfect numbers, Scripta Math., Vol. 19 (1953), pp. 128-131.
  • André Weil, Number Theory, An approach through history, From Hammurapi to Legendre, Birkhäuser, 1984, p. 6.
  • David Wells, The Penguin Dictionary of Curious and Interesting Numbers, pp. 107-110, Penguin Books, 1987.

Crossrefs

See A000043 for the current state of knowledge about Mersenne primes.
Cf. A228058 for Euler's criterion for odd terms.
Positions of 0's in A033879 and in A033880.
Cf. A001065.

Programs

  • Haskell
    a000396 n = a000396_list !! (n-1)
    a000396_list = [x | x <- [1..], a000203 x == 2 * x]
    -- Reinhard Zumkeller, Jan 20 2012
    
  • Mathematica
    Select[Range[9000], DivisorSigma[1,#]== 2*# &] (* G. C. Greubel, Oct 03 2017 *)
    PerfectNumber[Range[15]] (* Requires Mathematica version 10 or later *) (* Harvey P. Dale, Dec 10 2018 *)
  • PARI
    isA000396(n) = (sigma(n) == 2*n);
    
  • Python
    from sympy import divisor_sigma
    def ok(n): return n > 0 and divisor_sigma(n) == 2*n
    print([k for k in range(9999) if ok(k)]) # Michael S. Branicky, Mar 12 2022

Formula

The perfect number N = 2^(p-1)*(2^p - 1) is also multiplicatively p-perfect (i.e., A007955(N) = N^p), since tau(N) = 2*p. - Lekraj Beedassy, Sep 21 2004
a(n) = 2^A133033(n) - 2^A090748(n), assuming there are no odd perfect numbers. - Omar E. Pol, Feb 28 2008
a(n) = A000668(n)*(A000668(n)+1)/2, assuming there are no odd perfect numbers. - Omar E. Pol, Apr 23 2008
a(n) = A000217(A000668(n)), assuming there are no odd perfect numbers. - Omar E. Pol, May 09 2008
a(n) = Sum of the first A000668(n) positive integers, assuming there are no odd perfect numbers. - Omar E. Pol, May 09 2008
a(n) = A000384(A019279(n)), assuming there are no odd perfect numbers and no odd superperfect numbers. a(n) = A000384(A061652(n)), assuming there are no odd perfect numbers. - Omar E. Pol, Aug 17 2008
a(n) = A006516(A000043(n)), assuming there are no odd perfect numbers. - Omar E. Pol, Aug 30 2008
From Reikku Kulon, Oct 14 2008: (Start)
A144912(2, a(n)) = 1;
A144912(4, a(n)) = -1 for n > 1;
A144912(8, a(n)) = 5 or -5 for all n except 2;
A144912(16, a(n)) = -4 or -13 for n > 1. (End)
a(n) = A019279(n)*A000668(n), assuming there are no odd perfect numbers and odd superperfect numbers. a(n) = A061652(n)*A000668(n), assuming there are no odd perfect numbers. - Omar E. Pol, Jan 09 2009
a(n) = A007691(A153800(n)), assuming there are no odd perfect numbers. - Omar E. Pol, Jan 14 2009
Even perfect numbers N = K*A000203(K), where K = A019279(n) = 2^(p-1), A000203(A019279(n)) = A000668(n) = 2^p - 1 = M(p), p = A000043(n). - Lekraj Beedassy, May 02 2009
a(n) = A060286(A016027(n)), assuming there are no odd perfect numbers. - Omar E. Pol, Dec 13 2012
For n >= 2, a(n) = Sum_{k=1..A065549(n)} (2*k-1)^3, assuming there are no odd perfect numbers. - Derek Orr, Sep 28 2013
a(n) = A275496(2^((A000043(n) - 1)/2)) - 2^A000043(n), assuming there are no odd perfect numbers. - Daniel Poveda Parrilla, Aug 16 2016
a(n) = A156552(A324201(n)), assuming there are no odd perfect numbers. - Antti Karttunen, Mar 28 2019
a(n) = ((2^(A000043(n)))^3 - (2^(A000043(n)) - 1)^3 - 1)/6, assuming there are no odd perfect numbers. - Jules Beauchamp, Jun 06 2025

Extensions

I removed a large number of comments that assumed there are no odd perfect numbers. There were so many it was getting hard to tell which comments were true and which were conjectures. - N. J. A. Sloane, Apr 16 2023
Reference to Albert H. Beiler's book updated by Harvey P. Dale, Jan 13 2025

A051027 a(n) = sigma(sigma(n)) = sum of the divisors of the sum of the divisors of n.

Original entry on oeis.org

1, 4, 7, 8, 12, 28, 15, 24, 14, 39, 28, 56, 24, 60, 60, 32, 39, 56, 42, 96, 63, 91, 60, 168, 32, 96, 90, 120, 72, 195, 63, 104, 124, 120, 124, 112, 60, 168, 120, 234, 96, 252, 84, 224, 168, 195, 124, 224, 80, 128, 195, 171, 120, 360, 195, 360, 186, 234, 168, 480, 96
Offset: 1

Views

Author

Keywords

Examples

			a(2) = 4 because sigma(2)=1+2=3 and sigma(3)=1+3=4. - _Zak Seidov_, Aug 29 2012
		

References

  • József Sándor, On the composition of some arithmetic functions, Studia Univ. Babeș-Bolyai, Vol. 34, No. 1 (1989), pp. 7-14.
  • József Sándor and Borislav Crstici, Handbook of Number theory II, Kluwer Academic Publishers, 2004, Chapter 1, p. 39.

Crossrefs

Cf. A000203.

Programs

  • Maple
    with(numtheory): [seq(sigma(sigma(n)), n=1..100)];
  • Mathematica
    DivisorSigma[1,DivisorSigma[1,Range[100]]] (* Zak Seidov, Aug 29 2012 *)
  • PARI
    a(n)=sigma(sigma(n)); \\ Joerg Arndt, Feb 16 2014
    
  • Python
    from sympy import divisor_sigma as sigma
    def a(n): return sigma(sigma(n))
    print([a(n) for n in range(1, 62)]) # Michael S. Branicky, Dec 05 2021

Formula

a(n) = A000203(A000203(n)). - Zak Seidov, Aug 29 2012
a(p) = sigma(p+1) = A000203(p+1), for p prime. - Wesley Ivan Hurt, Feb 14 2014
a(n) = 2*n iff n = 2^q with M_(q+1) = 2^(q+1) - 1 is a Mersenne prime, hence iff n = 2^q with q in A090748. - Bernard Schott, Aug 08 2019
a(n) >= 2*n for even n, with equality only when n = 2^k and 2^(k+1) - 1 is prime (Sándor, 1989). - Amiram Eldar, Mar 09 2021

A138148 Cyclops numbers with binary digits only.

Original entry on oeis.org

0, 101, 11011, 1110111, 111101111, 11111011111, 1111110111111, 111111101111111, 11111111011111111, 1111111110111111111, 111111111101111111111, 11111111111011111111111, 1111111111110111111111111, 111111111111101111111111111, 11111111111111011111111111111
Offset: 0

Views

Author

Omar E. Pol, Mar 18 2008

Keywords

Comments

All members are palindromes A002113. The first five members are mentioned in A129868.
Also, binary representation of A129868.
a(A090748(n)) is equal to A138831(n), the n-th perfect number minus 1, written in base 2.
Except for the first term (replace 0 with 1) the binary representation of the n-th iteration of the elementary cellular automaton, Rule 219 starting with a single ON (black) cell. - Robert Price, Feb 21 2016
a(1) = 101 is only prime number in this sequence since a(n) = (10^(n+1)+1)*(10^n-1)/9. - Altug Alkan, May 11 2016

Examples

			n ........ a(n) .... A129868(n): value of a(n) read in base 2.
0 ......... 0 ......... 0
1 ........ 101 ........ 5
2 ....... 11011 ....... 27
3 ...... 1110111 ...... 119
4 ..... 111101111 ..... 495
5 .... 11111011111 .... 2015
6 ... 1111110111111 ... 8127
		

References

  • S. Wolfram, A New Kind of Science, Wolfram Media, 2002; p. 55.

Crossrefs

Cyclops numbers: A134808. Cf. A002113, A129868.
Cf. A002275 (repunits R_n = (10^n-1)/9), A011557 (10^n).

Programs

Formula

From Colin Barker, Feb 21 2013: (Start)
a(n) = (-1-9*10^n+10^(1+2*n))/9.
G.f.: x*(200*x-101) / ((x-1)*(10*x-1)*(100*x-1)). (End)
a(n) = 111*a(n-1) - 1110*a(n-2) + 1000*a(n-3) for n>2. - Wesley Ivan Hurt, Dec 08 2015
a(n) = A000533(n+1)*A002275(n). - Altug Alkan, May 12 2016
E.g.f.: (-1 - 9*exp(9*x) + 10*exp(99*x))*exp(x)/9. - Ilya Gutkovskiy, May 12 2016
a(n) = A002275(2n+1) - A011557(n). - M. F. Hasler, Feb 08 2020

Extensions

More terms from Omar E. Pol, Feb 09 2020

A019279 Superperfect numbers: numbers k such that sigma(sigma(k)) = 2*k where sigma is the sum-of-divisors function (A000203).

Original entry on oeis.org

2, 4, 16, 64, 4096, 65536, 262144, 1073741824, 1152921504606846976, 309485009821345068724781056, 81129638414606681695789005144064, 85070591730234615865843651857942052864
Offset: 1

Views

Author

Keywords

Comments

Let sigma_m(n) be result of applying sum-of-divisors function m times to n; call n (m,k)-perfect if sigma_m (n) = k*n; sequence gives (2,2)-perfect numbers.
Even values of these are 2^(p-1) where 2^p-1 is a Mersenne prime (A000043 and A000668). No odd superperfect numbers are known. Hunsucker and Pomerance checked that there are no odd ones below 7 * 10^24. - Jud McCranie, Jun 01 2000
The number of divisors of a(n) is equal to A000043(n), if there are no odd superperfect numbers. - Omar E. Pol, Feb 29 2008
The sum of divisors of a(n) is the n-th Mersenne prime A000668(n), provided that there are no odd superperfect numbers. - Omar E. Pol, Mar 11 2008
Largest proper divisor of A072868(n) if there are no odd superperfect numbers. - Omar E. Pol, Apr 25 2008
This sequence is a divisibility sequence if there are no odd superperfect numbers. - Charles R Greathouse IV, Mar 14 2012
For n>1, sigma(sigma(a(n))) + phi(phi(a(n))) = (9/4)*a(n). - Farideh Firoozbakht, Mar 02 2015
The term "super perfect number" was coined by Suryanarayana (1969). He and Kanold (1969) gave the general form of even superperfect numbers. - Amiram Eldar, Mar 08 2021

Examples

			sigma(sigma(4))=2*4, so 4 is in the sequence.
		

References

  • Dieter Bode, Über eine Verallgemeinerung der vollkommenen Zahlen, Dissertation, Braunschweig, 1971.
  • Richard K. Guy, Unsolved Problems in Number Theory, 3rd Edition, Springer, 2004, Section B9, pp. 99-100.
  • József Sándor, Dragoslav S. Mitrinovic and Borislav Crstici, Handbook of Number Theory I, Springer Science & Business Media, 2005, Chapter III, pp. 110-111.
  • József Sándor and Borislav Crstici, Handbook of Number theory II, Kluwer Academic Publishers, 2004, Chapter 1, pp. 38-42.
  • James J. Tattersall, Elementary Number Theory in Nine Chapters, Cambridge University Press, 1999, page 147.

Crossrefs

Programs

  • Mathematica
    sigma = DivisorSigma[1, #]&;
    For[n = 2, True, n++, If[sigma[sigma[n]] == 2 n, Print[n]]] (* Jean-François Alcover, Sep 11 2018 *)
  • PARI
    is(n)=sigma(sigma(n))==2*n \\ Charles R Greathouse IV, Nov 20 2012
    
  • Python
    from itertools import count, islice
    def A019279_gen(): # generator of terms
        return (n for n in count(1) if divisor_sigma(divisor_sigma(n)) == 2*n)
    A019279_list = list(islice(A019279_gen(),6)) # Chai Wah Wu, Feb 18 2022

Formula

a(n) = (1 + A000668(n))/2, if there are no odd superperfect numbers. - Omar E. Pol, Mar 11 2008
Also, if there are no odd superperfect numbers then a(n) = 2^A000043(n)/2 = A072868(n)/2 = A032742(A072868(n)). - Omar E. Pol, Apr 25 2008
a(n) = 2^A090748(n), if there are no odd superperfect numbers. - Ivan N. Ianakiev, Sep 04 2013

Extensions

a(8)-a(9) from Jud McCranie, Jun 01 2000
Corrected by Michel Marcus, Oct 28 2017

A061652 Even superperfect numbers: 2^(p-1) where 2^p-1 is a Mersenne prime (A000668).

Original entry on oeis.org

2, 4, 16, 64, 4096, 65536, 262144, 1073741824, 1152921504606846976, 309485009821345068724781056, 81129638414606681695789005144064, 85070591730234615865843651857942052864
Offset: 1

Views

Author

Jason Earls, Jun 16 2001

Keywords

Comments

It is conjectured that there are no odd superperfect numbers, in which case this coincides with A019279.
The number of divisors of a(n) is equal to A000043(n). - Omar E. Pol, Feb 29 2008
The sum of divisors of a(n) is equal to A000668(n), the n-th Mersenne prime. - Omar E. Pol, Mar 11 2008
Largest proper divisor of A072868(n). - Omar E. Pol, Apr 25 2008
Indices of hexagonal numbers (A000384) that are also even perfect numbers. [Omar E. Pol, Aug 26 2008]
Except for the first perfect number 6, this sequence is the greatest common divisor of a perfect number (A000396) and its arithmetic derivative (A003415). - Giorgio Balzarotti, Apr 21 2011
If n is in the sequence then n is a solution to the equation phi(sigma(x)) = 2x-2. It seems that there is no other solution to this equation. - Jahangeer Kholdi, Sep 09 2014
The sum of sums of elements of subsets of divisors of a(n), i.e. A229335(a(n)), is a perfect number (A000396). - Jaroslav Krizek, Nov 02 2017

References

  • James J. Tattersall, Elementary Number Theory in Nine Chapters, Cambridge University Press, 1999, page 147.

Crossrefs

Programs

  • Mathematica
    2^(Select[Range[512], PrimeQ[2^# - 1] &] - 1) (* Alonso del Arte, Apr 22 2011 *)
    2^(MersennePrimeExponent[Range[15]]-1) (* Requires Mathematica version 10 or later *) (* Harvey P. Dale, Jun 20 2021 *)
  • PARI
    forprime(p=2,1e3,if(ispseudoprime(2^p-1),print1(2^(p-1)", "))) \\ Charles R Greathouse IV, Mar 14 2012

Formula

a(n) = 2^(A090748(n)). - Lekraj Beedassy, Dec 07 2007
a(n) = (1 + A000668(n))/2. - Omar E. Pol, Mar 11 2008
a(n) = 2^A000043(n)/2 = A072868(n)/2 = A032742(A072868(n)). - Omar E. Pol, Apr 25 2008

A019278 Numbers j such that sigma(sigma(j)) = k*j for some k.

Original entry on oeis.org

1, 2, 4, 8, 15, 16, 21, 24, 42, 60, 64, 84, 160, 168, 240, 336, 480, 504, 512, 960, 1023, 1344, 1536, 4092, 4096, 10752, 13824, 16368, 29127, 32256, 32736, 47360, 57120, 58254, 61440, 65472, 65536, 86016, 116508, 217728, 262144, 331520, 343976, 466032, 550095
Offset: 1

Views

Author

Keywords

Comments

Let sigma^m (j) be the result of applying the sum-of-divisors function (A000203) m times to j; call j (m,k)-perfect if sigma^m (j) = k*j; then this is the sequence of (2,k)-perfect numbers.
From Michel Marcus, May 14 2016: (Start)
For these numbers, the quotient k = sigma(sigma(j))/j is an integer (see A098223). Then also k = (sigma(s)/s)*(sigma(j)/j) with s = sigma(j). That is, k = abundancy(s)*abundancy(j).
So looking at the abundancy of these terms may be interesting. Indeed we see that 459818240 and 51001180160 are actually 3-perfect numbers (A005820), and the reason they are here is that they are coprime to 3. So their sums of divisors are 4-perfect numbers (A027687), yielding q=12.
In a similar way, we can see that the 5-perfect numbers (A046060) that are coprime to 5 will be terms of this sequence with q=30. There are 20 such numbers, the smallest being 13188979363639752997731839211623940096. (End)
From Michel Marcus, May 15 2016: (Start)
It is also interesting to note that for a(2)=8, s=sigma(8)=15 is also a term. This happens to be the case for chains of several terms in a row:
8, 15, 24, 60, 168, 480 with k = 3,4,7,8,9,10;
512, 1023, 1536, 4092, 10752, 32736 with k = 3,4,7,8,9,10;
29127, 47360, 116508, 331520, 932064, 2983680 with k = 4,7,8,9,13,14;
1556480, 3932040, 14008320 with k = 9,13,14;
106151936, 251650560, 955367424 with k = 9,13,14;
312792480, 1505806848 with k = 19,20;
6604416000, 30834059256 with k = 19,20;
9623577600, 46566269568 with k = 19,20.
When j is a term, we can test if s=sigma(j) is also a term; this way we get 6 more terms: 572941926400, 845734196736, 1422976331052, 4010593484800, 11383810648416, 36095341363200.
And the corresponding chains are:
173238912000, 845734196736 with k = 19,20;
355744082763, 572941926400, 1422976331052, 4010593484800, 11383810648416, 36095341363200 with k = 4,7,8,9,13,14. (End)
From Altug Alkan, May 17 2016: (Start)
Here are additional chains for the above list:
57120, 217728 with k = 13,14;
343976, 710400 with k = 7,8;
1980342, 5621760 with k = 10,14;
4404480, 14913024 with k = 11,12;
238608384, 775898880 with k = 11,12. (End)
Currently, the coefficient pairs are [1, 1], [3, 4], [4, 7], [7, 8], [8, 9], [9, 10], [9, 13], [10, 14], [11, 12], [13, 14], [16, 17], [16, 21], [17, 18], [19, 20], [23, 24], [25, 26], [25, 31], [27, 28], [29, 30], [31, 32], [32, 33], [37, 38]. It is interesting to note that for some of them, the pair (s,t) also satisfies t=sigma(s). - Michel Marcus, Jul 03 2016; Sep 06 2016
Using these empirical pairs of coefficients in conjunction with the first comment allows us to determine whether some term is the sum of divisors of another yet unknown smaller term. - Michel Marcus, Jul 04 2016
For m in A090748 = A000043 - 1 and c in A205597 (= odd a(n)), c*2^m is in the sequence, unless 2^(m+1)-1 | sigma(c). Indeed, from sigma(x*y) = sigma(x)*sigma(y) for gcd(x,y) = 1, we get sigma(sigma(c*2^m)) = sigma(sigma(c))*2^(m+1), so c*2^m is in the sequence if sigma(sigma(c))/c = k/2 (where k can't be odd: A330598 has no odd c). - M. F. Hasler, Jan 06 2020

Crossrefs

For sigma see A000203 and A007691.
Cf. A205597 (odd terms), A323653 (those terms that are in A007691, i.e., for which sigma(n)/n is also an integer), A330598 (half-integer ratio).

Programs

  • Magma
    [m: m in [1..560000]| IsIntegral(DivisorSigma(1,DivisorSigma(1,m))/m)]; // Marius A. Burtea, Nov 16 2019
  • Mathematica
    Select[Range[100000], Mod[DivisorSigma[1, DivisorSigma[1, #]], #] == 0 &] (* Carl Najafi, Aug 22 2011 *)
  • PARI
    is_A019278(n)=sigma(sigma(n))%n==0 \\ M. F. Hasler, Jul 02 2016
    
  • Python
    from sympy.ntheory import divisor_sigma as D
    print([i for i in range(1, 10000) if D(D(i, 1), 1)%i==0]) # Indranil Ghosh, Mar 17 2017
    

Extensions

Simpler definition from M. F. Hasler, Jul 02 2016

A133033 Number of proper divisors of n-th even perfect number.

Original entry on oeis.org

3, 5, 9, 13, 25, 33, 37, 61, 121, 177, 213, 253, 1041, 1213, 2557, 4405, 4561, 6433, 8505, 8845, 19377, 19881, 22425, 39873, 43401, 46417, 88993, 172485, 221005, 264097, 432181, 1513677, 1718865, 2515573, 2796537, 5952441, 6042753, 13945185, 26933833, 41992021, 48073165, 51929901, 60804913
Offset: 1

Views

Author

Omar E. Pol, Oct 27 2007, Feb 23 2008, Apr 28 2009

Keywords

Comments

Perfect numbers: A000396(n) = 2^a(n) - 2^A090748(n), assuming there are no odd perfect numbers.
Also, a(n) is equal to the number of bits in A135650(n), the n-th even perfect number written in base 2.
These are odd numbers k such that 2^((k+1)/2) - 1 is (a Mersenne) prime. - Thomas Ordowski, Apr 20 2019

Crossrefs

Programs

  • Mathematica
    2 * MersennePrimeExponent[Range[48]] - 1 (* Amiram Eldar, Oct 18 2024 *)

Formula

a(n) = A061645(n) - 1.
a(n) = A000043(n) + A090748(n) = 2*A000043(n) - 1 = 2*A090748(n) + 1.

Extensions

a(39)-a(43) from Ivan Panchenko, Apr 12 2018

A135650 Even perfect numbers written in base 2.

Original entry on oeis.org

110, 11100, 111110000, 1111111000000, 1111111111111000000000000, 111111111111111110000000000000000, 1111111111111111111000000000000000000
Offset: 1

Views

Author

Omar E. Pol, Feb 21 2008, Feb 22 2008, Apr 28 2009

Keywords

Comments

The number of digits of a(n) is equal to 2*A000043(n)-1. The central digit is "1". The first digits are "1". The last digits are "0". The number of digits "1" is equal A000043(n). The number of digits "0" is equal A000043(n)-1.
The concatenation of digits "1" of a(n) gives the n-th Mersenne prime written in binary (see A117293(n)).
Also, the number of digits of a(n) is equal to A133033(n), the number of proper divisors of n-th even perfect number.

Examples

			a(3) = 111110000 because the 3rd even perfect number is 496 and 496 written in base 2 is 111110000. Note that 11111 is the 3rd Mersenne prime A000668(3) = 31 written in base 2.
		

Crossrefs

Programs

  • Mathematica
    Map[FromDigits[IntegerDigits[#, 2]] &, PerfectNumber[Range[8], "Even"]] (* Amiram Eldar, Oct 18 2024 *)

A120375 Integers k such that 2*5^k - 1 is prime.

Original entry on oeis.org

4, 6, 16, 24, 30, 54, 96, 178, 274, 1332, 2766, 3060, 4204, 17736, 190062, 223536, 260400, 683080
Offset: 1

Views

Author

Walter Kehowski, Jun 28 2006

Keywords

Comments

See comments for A057472. Examined in base 12, all n must be even and all primes must be 1-primes. For example, 1249 is 881 in base 12.
a(16) > 2*10^5. - Robert Price, Mar 14 2015

Examples

			a(1) = 4 since 2*5^4 - 1 = 1249 is the first prime.
		

Crossrefs

Integers k such that 2*b^k - 1 is prime: A090748 (b=2), A003307 (b=3), this sequence (b=5), A057472 (b=6), A002959 (b=7), A002957 (b=10), A120378 (b=11).
Primes of the form 2*b^k - 1: A000668 (b=2), A079363 (b=3), A120376 (b=5), A158795 (b=7), A055558 (b=10), A120377 (b=11).
Cf. also A000043, A002958.

Programs

  • Magma
    [n: n in [0..2800] |IsPrime(2*5^n - 1)]; // Vincenzo Librandi, Sep 23 2018
  • Maple
    for w to 1 do for k from 1 to 2000 do n:=2*5^k-1; if isprime(n) then printf("%d, %d ",k,n) fi od od;
  • Mathematica
    Select[Range[0, 100], PrimeQ[2*5^# - 1] &] (* Robert Price, Mar 14 2015 *)
  • PARI
    isok(k) = ispseudoprime(2*5^k-1); \\ Altug Alkan, Sep 22 2018
    

Formula

a(n) = 2*A002958(n).

Extensions

More terms from Ryan Propper, Mar 28 2007
a(14) from Herman Jamke (hermanjamke(AT)fastmail.fm), May 02 2007
a(15) from Robert Price, Mar 14 2015
a(16)-a(18) from Jorge Coveiro and Tyler NeSmith, Jun 14 2020
Showing 1-10 of 31 results. Next