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 75 results. Next

A006037 Weird numbers: abundant (A005101) but not pseudoperfect (A005835).

Original entry on oeis.org

70, 836, 4030, 5830, 7192, 7912, 9272, 10430, 10570, 10792, 10990, 11410, 11690, 12110, 12530, 12670, 13370, 13510, 13790, 13930, 14770, 15610, 15890, 16030, 16310, 16730, 16870, 17272, 17570, 17990, 18410, 18830, 18970, 19390, 19670
Offset: 1

Views

Author

Keywords

Comments

OProject@Home in subproject Weird Engine calculates and stores the weird numbers.
There are no odd weird numbers < 10^17. - Robert A. Hearn (rah(AT)ai.mit.edu), May 25 2005
From Alois P. Heinz, Oct 30 2009: (Start)
The first weird number that has more than one decomposition of its divisors set into two subsets with equal sum (and thus is not a member of A083209) is 10430:
1+5+7+10+14+35+298+10430 = 2+70+149+745+1043+1490+2086+5215
2+70+298+10430 = 1+5+7+10+14+35+149+745+1043+1490+2086+5215. (End)
There are no odd weird numbers < 1.8*10^19. - Wenjie Fang, Sep 04 2013
S. Benkowski and P. Erdős (1974) proved that the asymptotic density W of weird numbers is positive. It can be shown that W < 0.0101 (see A005835). - Jaycob Coleman, Oct 26 2013
No odd weird number exists below 10^21. This search was done on the volunteer computing project yoyo@home. - Wenjie Fang, Feb 23 2014
No odd weird number with abundance less than 10^14 exists below 10^28. See Odd Weird Search link. - Wenjie Fang, Feb 25 2015
A weird number k multiplied by a prime p > sigma(k) is again weird. Primitive weird numbers (A002975) are those which are not a multiple of a smaller term, i.e., don't have a weird proper divisor. Sequence A065235 lists odd numbers that can be written in only one way as sum of their divisors, and A122036 lists those which are not in A136446, i.e., not sum of proper divisors > 1. - M. F. Hasler, Jul 30 2016

References

  • J.-M. De Koninck, Ces nombres qui nous fascinent, Entry 70, p. 24, Ellipses, Paris 2008.
  • R. K. Guy, Unsolved Problems in Number Theory, B2.
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).
  • David Wells, The Penguin Dictionary of Curious and Interesting Numbers. Penguin Books, NY, 1986, Revised edition 1987. See p. 129.

Crossrefs

Programs

  • Haskell
    a006037 n = a006037_list !! (n-1)
    a006037_list = filter ((== 0) . a210455) a005101_list
    -- Reinhard Zumkeller, Jan 21 2013
  • Maple
    isA006037 := proc(n)
        isA005101(n) and not isA005835(n) ;
    end proc:
    for n from 1 do
        if isA006037(n) then
            print(n);
        end if;
    end do: # R. J. Mathar, Jun 18 2015
  • Mathematica
    (* first do *) Needs["DiscreteMath`Combinatorica`"] (* then *) fQ[n_] := Block[{d, l, t, i}, If[ DivisorSigma[1, n] > 2n && Mod[n, 6] != 0, d = Take[Divisors[n], {1, -2}]; l = 2^Length[d]; t = Table[ NthSubset[j, d], {j, l - 1}]; i = 1; While[i < l && Plus @@ t[[i]] != n, i++ ]]; If[i == l, True, False]]; Select[ Range[ 20000], fQ[ # ] &] (* Robert G. Wilson v, May 20 2005 *)
  • PARI
    is_A006037(n,d=divisors(n),s=vecsum(d)-n,m=#d-1)={ m||return; while(d[m]>n, s-=d[m]; m--); d[m]n, is_A006037(n-d[m], d, s-d[m], m-1) && is_A006037(n, d, s-d[m], m-1), sM. F. Hasler, Mar 30 2008; improved and updated to current PARI syntax by M. F. Hasler, Jul 15 2016
    
  • PARI
    is_A006037(n, d=divisors(n)[^-1], s=vecsum(d))={s>n && !is_A005835(n,d,s)} \\ Equivalent but slightly faster than the self-contained version above.-- For efficiency, ensure that the argument is even or add "!bittest(n,0) && ..." to check this first. - M. F. Hasler, Jul 17 2016
    
  • PARI
    t=0; A006037=vector(100,i, until( is_A006037(t+=2),); t) \\ M. F. Hasler, Mar 30 2008
    

Extensions

More terms from Jud McCranie, Oct 21 2001

A136447 Numbers which are not pseudoperfect: Complement of A005835.

Original entry on oeis.org

1, 2, 3, 4, 5, 7, 8, 9, 10, 11, 13, 14, 15, 16, 17, 19, 21, 22, 23, 25, 26, 27, 29, 31, 32, 33, 34, 35, 37, 38, 39, 41, 43, 44, 45, 46, 47, 49, 50, 51, 52, 53, 55, 57, 58, 59, 61, 62, 63, 64, 65, 67, 68, 69, 70, 71, 73, 74, 75, 76, 77, 79, 81, 82, 83, 85, 86, 87, 89, 91, 92, 93, 94, 95
Offset: 1

Views

Author

N. J. A. Sloane, Apr 06 2008

Keywords

Comments

Contains the weird numbers 70, 836, ... (A006037), so is different from the set A005100 of deficient numbers.
A065205(a(n)) = A210455(a(n)) = 0. - Reinhard Zumkeller, Jan 21 2013

Crossrefs

Programs

  • Haskell
    a136447 n = a136447_list !! (n-1)
    a136447_list = filter ((== 0) . a210455) [1..]
    -- Reinhard Zumkeller, Jan 21 2013
  • Mathematica
    ppQ[n_] := AnyTrue[Subsets[Divisors[n] // Most], Total[#] == n&];
    Select[Range[100], Not[ppQ[#]]&] (* Jean-François Alcover, Aug 03 2018 *)

Extensions

Edited by M. F. Hasler, Jul 18 2016

A035480 Pseudoperfect numbers (A005835) that are not practical (second definition, A007620).

Original entry on oeis.org

102, 114, 138, 174, 186, 222, 246, 258, 282, 318, 350, 354, 366, 372, 402, 426, 438, 444, 474, 490, 492, 498, 516, 534, 550, 564, 572, 582, 606, 618, 636, 642, 650, 654, 678, 708, 732, 738, 748, 762, 770, 774, 786, 804, 822, 834, 846, 852, 876, 894, 906, 910, 940, 942, 945, 948, 954, 978, 996, 1002
Offset: 1

Views

Author

Keywords

Crossrefs

Numbers in A005835 but not in A007620.
Cf. A007621.

Programs

  • Mathematica
    Cases[ A005835, q_/;(Count[ CoefficientList[ Series[ Times@@( (1+z^#)& /@ Divisors[ q ] ), {z, 0, q} ], z ], 0 ]>0) ]

Extensions

Corrected by David W. Wilson.

A306987 Primitive abundant numbers (A071395) that are pseudoperfect (A005835).

Original entry on oeis.org

20, 88, 104, 272, 304, 368, 464, 550, 572, 650, 748, 945, 1184, 1312, 1376, 1430, 1504, 1575, 1696, 1870, 1888, 1952, 2002, 2090, 2205, 2210, 2470, 2530, 2584, 2990, 3128, 3190, 3230, 3410, 3465, 3496, 3770, 3944, 4070, 4095, 4216, 4288, 4408, 4510, 4544, 4672
Offset: 1

Views

Author

Amiram Eldar, Mar 18 2019

Keywords

Comments

By definition these numbers are also primitive pseudoperfect (A006036).
Benkoski and Erdős proved that this sequence is infinite, since it includes all the numbers of the form 2^k * p with p a prime such that 2^k < p < 2^(k+1).

Crossrefs

Programs

  • Mathematica
    paQ[n_]:=DivisorSigma[1, n] > 2n && Times @@ Boole@ Map[DivisorSigma[1, #] < 2 # &, Most@ Divisors@ n] == 1; psQ[n_]:=Module[{d= Most[Divisors[n] ]}, SeriesCoefficient[Series[Product[1 + x^d[[i]], {i, Length[d]}], {x, 0, n}], n] > 0]; Select[Range[5000], paQ[#]&&psQ[#]&] (* after Michael De Vlieger at A071395 and T. D. Noe at A005835 *)

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

A005101 Abundant numbers (sum of divisors of m exceeds 2m).

Original entry on oeis.org

12, 18, 20, 24, 30, 36, 40, 42, 48, 54, 56, 60, 66, 70, 72, 78, 80, 84, 88, 90, 96, 100, 102, 104, 108, 112, 114, 120, 126, 132, 138, 140, 144, 150, 156, 160, 162, 168, 174, 176, 180, 186, 192, 196, 198, 200, 204, 208, 210, 216, 220, 222, 224, 228, 234, 240, 246, 252, 258, 260, 264, 270
Offset: 1

Views

Author

Keywords

Comments

A number m is abundant if sigma(m) > 2m (this sequence), perfect if sigma(m) = 2m (cf. A000396), or deficient if sigma(m) < 2m (cf. A005100), where sigma(m) is the sum of the divisors of m (A000203).
While the first even abundant number is 12 = 2^2*3, the first odd abundant is 945 = 3^3*5*7, the 232nd abundant number!
It appears that for m abundant and > 23, 2*A001055(m) - A101113(m) is NOT 0. - Eric Desbiaux, Jun 01 2009
If m is a term so is every positive multiple of m. "Primitive" terms are in A091191.
If m=6k (k>=2), then sigma(m) >= 1 + k + 2*k + 3*k + 6*k > 12*k = 2*m. Thus all such m are in the sequence.
According to Deléglise (1998), the abundant numbers have natural density 0.2474 < A(2) < 0.2480. Thus the n-th abundant number is asymptotic to 4.0322*n < n/A(2) < 4.0421*n. - Daniel Forgues, Oct 11 2015
From Bob Selcoe, Mar 28 2017 (prompted by correspondence with Peter Seymour): (Start)
Applying similar logic as the proof that all multiples of 6 >= 12 appear in the sequence, for all odd primes p:
i) all numbers of the form j*p*2^k (j >= 1) appear in the sequence when p < 2^(k+1) - 1;
ii) no numbers appear when p > 2^(k+1) - 1 (i.e., are deficient and are in A005100);
iii) when p = 2^(k+1) - 1 (i.e., perfect numbers, A000396), j*p*2^k (j >= 2) appear.
Note that redundancies are eliminated when evaluating p only in the interval [2^k, 2^(k+1)].
The first few even terms not of the forms i or iii are {70, 350, 490, 550, 572, 650, 770, ...}. (End)

References

  • L. E. Dickson, Theorems and tables on the sum of the divisors of a number, Quart. J. Pure Appl. Math., Vol. 44 (1913), pp. 264-296.
  • Richard K. Guy, Unsolved Problems in Number Theory, 3rd Edition, Springer, 2004, Section B2, pp. 74-84.
  • Clifford A. Pickover, A Passion for Mathematics, Wiley, 2005; see p. 59.
  • 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, page 128.

Crossrefs

Cf. A005231 and A006038 (odd abundant numbers).
Cf. A094268 (n consecutive abundant numbers).
Cf. A173490 (even abundant numbers).
Cf. A001065.
Cf. A000396 (perfect numbers).
Cf. A302991.

Programs

  • Haskell
    a005101 n = a005101_list !! (n-1)
    a005101_list = filter (\x -> a001065 x > x) [1..]
    -- Reinhard Zumkeller, Nov 01 2015, Jan 21 2013
    
  • Maple
    with(numtheory): for n from 1 to 270 do if sigma(n)>2*n then printf(`%d,`,n) fi: od:
    isA005101 := proc(n)
        simplify(numtheory[sigma](n) > 2*n) ;
    end proc: # R. J. Mathar, Jun 18 2015
    A005101 := proc(n)
        option remember ;
        local a ;
        if n =1 then
            12 ;
        else
            a := procname(n-1)+1 ;
            while numtheory[sigma](a) <= 2*a do
                a := a+1 ;
            end do ;
            a ;
        end if ;
    end proc: # R. J. Mathar, Oct 11 2017
  • Mathematica
    abQ[n_] := DivisorSigma[1, n] > 2n; A005101 = Select[ Range[270], abQ[ # ] &] (* Robert G. Wilson v, Sep 15 2005 *)
    Select[Range[300], DivisorSigma[1, #] > 2 # &] (* Vincenzo Librandi, Oct 12 2015 *)
  • PARI
    isA005101(n) = (sigma(n) > 2*n) \\ Michael B. Porter, Nov 07 2009
    
  • Python
    from sympy import divisors
    def ok(n): return sum(divisors(n)) > 2*n
    print(list(filter(ok, range(1, 271)))) # Michael S. Branicky, Aug 29 2021
    
  • Python
    from sympy import divisor_sigma
    from itertools import count, islice
    def A005101_gen(startvalue=1): return filter(lambda n:divisor_sigma(n) > 2*n, count(max(startvalue, 1))) # generator of terms >= startvalue
    A005101_list = list(islice(A005101_gen(), 20)) # Chai Wah Wu, Jan 14 2022

Formula

a(n) is asymptotic to C*n with C=4.038... (Deléglise, 1998). - Benoit Cloitre, Sep 04 2002
A005101 = { n | A033880(n) > 0 }. - M. F. Hasler, Apr 19 2012
A001065(a(n)) > a(n). - Reinhard Zumkeller, Nov 01 2015

A083207 Zumkeller or integer-perfect numbers: numbers n whose divisors can be partitioned into two disjoint sets with equal sum.

Original entry on oeis.org

6, 12, 20, 24, 28, 30, 40, 42, 48, 54, 56, 60, 66, 70, 78, 80, 84, 88, 90, 96, 102, 104, 108, 112, 114, 120, 126, 132, 138, 140, 150, 156, 160, 168, 174, 176, 180, 186, 192, 198, 204, 208, 210, 216, 220, 222, 224, 228, 234, 240, 246, 252, 258, 260, 264, 270, 272
Offset: 1

Views

Author

Reinhard Zumkeller, Apr 22 2003

Keywords

Comments

The 229026 Zumkeller numbers less than 10^6 have a maximum difference of 12. This leads to the conjecture that any 12 consecutive numbers include at least one Zumkeller number. There are 1989 odd Zumkeller numbers less than 10^6; they are exactly the odd abundant numbers that have even abundance, A174865. - T. D. Noe, Mar 31 2010
For k >= 0, numbers of the form 18k + 6 and 18k + 12 are terms (see Remark 2.3. in Somu et al., 2023). Corollary: The maximum difference between any two consecutive terms is at most 12. - Ivan N. Ianakiev, Jan 02 2024
All 205283 odd abundant numbers less than 10^8 that have even abundance (see A174865) are Zumkeller numbers. - T. D. Noe, Nov 14 2010
Except for 1 and 2, all primorials (A002110) are Zumkeller numbers (follows from Fact 6 in the Rao/Peng paper). - Ivan N. Ianakiev, Mar 23 2016
Supersequence of A111592 (follows from Fact 3 in the Rao/Peng paper). - Ivan N. Ianakiev, Mar 20 2017
Conjecture: Any 4 consecutive terms include at least one number k such that sigma(k)/2 is also a Zumkeller number (verified for the first 10^5 Zumkeller numbers). - Ivan N. Ianakiev, Apr 03 2017
LeVan studied these numbers using the equivalent definition of numbers n such that n = Sum_{d|n, dA180332) "minimal integer-perfect numbers". - Amiram Eldar, Dec 20 2018
The numbers 3 * 2^k for k > 0 are all Zumkeller numbers: half of one such partition is {3*2^k, 3*2^(k-2), ...}, replacing 3 with 2 if it appears. With this and the lemma that the product of a Zumkeller number and a number coprime to it is again a Zumkeller number (see A179527), we have that all numbers divisible by 6 but not 9 (or numbers congruent to 6 or 12 modulo 18) are Zumkeller numbers, proving that the difference between consecutive Zumkeller numbers is at most 12. - Charlie Neder, Jan 15 2019
Improvements on the previous comment: 1) For every integer q > 0, every odd integer r > 0 and every integer s > 0 relatively prime to 6, the integer 2^q*3^r*s is a Zumkeller number, and therefore 2) there exist Zumkeller numbers divisible by 9 (such as 54, 90, 108, 126, etc.). - Ivan N. Ianakiev, Jan 16 2020
Conjecture: If d > 1, d|k and tau(d)*sigma(d) = k, then k is a Zumkeller number (cf. A331668). - Ivan N. Ianakiev, Apr 24 2020
This sequence contains A378541, the intersection of the practical numbers (A005153) with numbers with even sum of divisors (A028983). - David A. Corneth, Nov 03 2024
Sequence gives the positions of even terms in A119347, and correspondingly, of odd terms in A308605. - Antti Karttunen, Nov 29 2024
If s = sigma(m) is odd and p > s then m*p is not in the sequence. - David A. Corneth, Dec 07 2024

Examples

			Given n = 48, we can partition the divisors thus: 1 + 3 + 4 + 6 + 8 + 16 + 24 = 2 + 12 + 48, therefore 48 is a term (A083206(48) = 5).
From _David A. Corneth_, Dec 04 2024: (Start)
30 is in the sequence. sigma(30) = 72. So we look for distinct divisors of 30 that sum to 72/2 = 36. That set or its complement contains 30. The other divisors in that set containing 30 sum to 36 - 30 = 6. So we look for some distinct proper divisors of 30 that sum to 6. That is from the divisors of {1, 2, 3, 5, 6, 10, 15}. It turns out that both 1+2+3 and 6 satisfy this condition. So 36 is in the sequence.
25 is not in the sequence as sigma(25) = 31 which is odd so the sum of two equal integers cannot be the sum of divisors of 25.
33 is not in the sequence as sigma(33) = 48 < 2*33. So is impossible to have a partition of the set of divisors into two disjoint set the sum of each of them sums to 48/2 = 24 as one of them contains 33 > 24 and any other divisors are nonnegative. (End)
		

References

  • Marijo O. LeVan, Integer-perfect numbers, Journal of Natural Sciences and Mathematics, Vol. 27, No. 2 (1987), pp. 33-50.
  • Marijo O. LeVan, On the order of nu(n), Journal of Natural Sciences and Mathematics, Vol. 28, No. 1 (1988), pp. 165-173.
  • J. Sandor and B. Crstici, Handbook of Number Theory, II, Springer Verlag, 2004, chapter 1.10, pp. 53-54.

Crossrefs

Positions of nonzero terms in A083206, positions of 0's in A103977 and in A378600.
Positions of even terms in A119347, of odd terms in A308605.
Complement of A083210.
Subsequence of A023196 and of A028983.
Union of A353061 and A378541.
Conjectured subsequences: A007691, A331668 (after their initial 1's), A351548 (apart from 0-terms).
Cf. A174865 (Odd abundant numbers whose abundance is even).
Cf. A204830, A204831 (equal sums of 3 or 4 disjoint subsets).
Cf. A000203, A005101, A005153 (practical numbers), A005835, A027750, A048055, A083206, A083208, A083211, A171641, A175592, A179527 (characteristic function), A221054.

Programs

  • Haskell
    a083207 n = a083207_list !! (n-1)
    a083207_list = filter (z 0 0 . a027750_row) $ [1..] where
       z u v []     = u == v
       z u v (p:ps) = z (u + p) v ps || z u (v + p) ps
    -- Reinhard Zumkeller, Apr 18 2013
    
  • Maple
    with(numtheory): with(combstruct):
    is_A083207 := proc(n) local S, R, Found, Comb, a, s; s := sigma(n);
    if not(modp(s, 2) = 0 and n * 2 <= s) then return false fi;
    S := s / 2 - n; R := select(m -> m <= S, divisors(n)); Found := false;
    Comb := iterstructs(Combination(R)):
    while not finished(Comb) and not Found do
       Found := add(a, a = nextstruct(Comb)) = S
    od; Found end:
    A083207_list := upto -> select(is_A083207, [$1..upto]):
    A083207_list(272); # Peter Luschny, Dec 14 2009, updated Aug 15 2014
  • Mathematica
    ZumkellerQ[n_] := Module[{d=Divisors[n], t, ds, x}, ds = Plus@@d; If[Mod[ds, 2] > 0, False, t = CoefficientList[Product[1 + x^i, {i, d}], x]; t[[1 + ds/2]] > 0]]; Select[Range[1000], ZumkellerQ] (* T. D. Noe, Mar 31 2010 *)
    znQ[n_]:=Length[Select[{#,Complement[Divisors[n],#]}&/@Most[Rest[ Subsets[ Divisors[ n]]]],Total[#[[1]]]==Total[#[[2]]]&]]>0; Select[Range[300],znQ] (* Harvey P. Dale, Dec 26 2022 *)
  • PARI
    part(n,v)=if(n<1, return(n==0)); forstep(i=#v,2,-1,if(part(n-v[i],v[1..i-1]), return(1))); n==v[1]
    is(n)=my(d=divisors(n),s=sum(i=1,#d,d[i])); s%2==0 && part(s/2-n,d[1..#d-1]) \\ Charles R Greathouse IV, Mar 09 2014
    
  • PARI
    \\ See Corneth link
    
  • Python
    from sympy import divisors
    from sympy.combinatorics.subsets import Subset
    for n in range(1,10**3):
        d = divisors(n)
        s = sum(d)
        if not s % 2 and max(d) <= s/2:
            for x in range(1,2**len(d)):
                if sum(Subset.unrank_binary(x,d).subset) == s/2:
                    print(n,end=', ')
                    break
    # Chai Wah Wu, Aug 13 2014
    
  • Python
    from sympy import divisors
    import numpy as np
    A083207 = []
    for n in range(2,10**3):
        d = divisors(n)
        s = sum(d)
        if not s % 2 and 2*n <= s:
            d.remove(n)
            s2, ld = int(s/2-n), len(d)
            z = np.zeros((ld+1,s2+1),dtype=int)
            for i in range(1,ld+1):
                y = min(d[i-1],s2+1)
                z[i,range(y)] = z[i-1,range(y)]
                z[i,range(y,s2+1)] = np.maximum(z[i-1,range(y,s2+1)],z[i-1,range(0,s2+1-y)]+y)
                if z[i,s2] == s2:
                    A083207.append(n)
                    break
    # Chai Wah Wu, Aug 19 2014
    
  • Sage
    def is_Zumkeller(n):
        s = sigma(n)
        if not (2.divides(s) and n*2 <= s): return False
        S = s // 2 - n
        R = (m for m in divisors(n) if m <= S)
        return any(sum(c) == S for c in Combinations(R))
    A083207_list = lambda lim: [n for n in (1..lim) if is_Zumkeller(n)]
    print(A083207_list(272)) # Peter Luschny, Sep 03 2018

Formula

A083206(a(n)) > 0.
A083208(n) = A083206(a(n)).
A179529(a(n)) = 1. - Reinhard Zumkeller, Jul 19 2010

Extensions

Name improved by T. D. Noe, Mar 31 2010
Name "Zumkeller numbers" added by N. J. A. Sloane, Jul 08 2010

A005231 Odd abundant numbers (odd numbers m whose sum of divisors exceeds 2m).

Original entry on oeis.org

945, 1575, 2205, 2835, 3465, 4095, 4725, 5355, 5775, 5985, 6435, 6615, 6825, 7245, 7425, 7875, 8085, 8415, 8505, 8925, 9135, 9555, 9765, 10395, 11025, 11655, 12285, 12705, 12915, 13545, 14175, 14805, 15015, 15435, 16065, 16695, 17325, 17955
Offset: 1

Views

Author

Keywords

Comments

While the first even abundant number is 12 = 2^2*3, the first odd abundant is 945 = 3^3*5*7, the 232nd abundant number.
Schiffman notes that 945+630k is in this sequence for all k < 52. Most of the first initial terms are of the form. Among the 1996 terms below 10^6, 1164 terms are of that form, and only 26 terms are not divisible by 5, cf. A064001. - M. F. Hasler, Jul 16 2016
From M. F. Hasler, Jul 28 2016: (Start)
Any multiple of an abundant number is again abundant, see A006038 for primitive terms, i.e., those which are not a multiple of an earlier term.
An odd abundant number must have at least 3 distinct prime factors, and 5 prime factors when counted with multiplicity (A001222), whence a(1) = 3^3*5*7. To see this, write the relative abundancy A(N) = sigma(N)/N = sigma[-1](N) as A(Product p_i^e_i) = Product (p_i-1/p_i^e_i)/(p_i-1) < Product p_i/(p_i-1).
See A115414 for terms not divisible by 3, A064001 for terms not divisible by 5, A112640 for terms coprime to 5*7, and A047802 for other generalizations.
As of today, we don't know of a difference between this set S of odd abundant numbers and the set S' of odd semiperfect numbers: Elements of S' \ S would be perfect (A000396), and elements of S \ S' would be weird (A006037), but no odd weird or perfect number is known. (End)
For any term m in this sequence, A064989(m) is also an abundant number (in A005101), and for any term x in A115414, A064989(x) is in this sequence. If there are no odd perfect numbers, then applying A064989 to these terms and sorting into ascending order gives A337386. - Antti Karttunen, Aug 28 2020
There exist infinitely many terms m such that 2*m+1 is also a term. An example of such a term is given by m = 985571808130707987847768908867571007187. - Max Alekseyev, Nov 16 2023

References

  • W. Dunham, Euler: The Master of Us All, The Mathematical Association of America Inc., Washington, D.C., 1999, p. 13.
  • R. K. Guy, Unsolved Problems in Number Theory, B2.
  • James J. Tattersall, Elementary Number Theory in Nine Chapters, Cambridge University Press, 1999, page 128.

Crossrefs

Programs

  • Maple
    A005231 := proc(n) option remember ; local a ; if n = 1 then 945 ; else for a from procname(n-1)+2 by 2 do if numtheory[sigma](a) > 2*a then return a; end if; end do: end if; end proc: # R. J. Mathar, Mar 20 2011
  • Mathematica
    fQ[n_] := DivisorSigma[1, n] > 2n; Select[1 + 2Range@ 9000, fQ] (* Robert G. Wilson v, Mar 20 2011 *)
  • PARI
    je=[]; forstep(n=1,15000,2, if(sigma(n)>2*n, je=concat(je,n))); je
    
  • PARI
    is_A005231(n)={bittest(n,0)&&sigma(n)>2*n} \\ M. F. Hasler, Jul 28 2016
    
  • PARI
    list(lim)=my(v=List()); forfactored(n=945,lim\1, if(n[2][1,1]>2 && sigma(n,-1)>2, listput(v,n[1]))); Vec(v) \\ Charles R Greathouse IV, Apr 21 2022

Formula

a(n) ~ k*n for some constant k (perhaps around 500). - Charles R Greathouse IV, Apr 21 2022
482.8 < k < 489.8 (based on density bounds by Kobayashi et al., 2009). - Amiram Eldar, Jul 17 2022

Extensions

More terms from James Sellers

A033630 Number of partitions of n into distinct divisors of n.

Original entry on oeis.org

1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 3, 1, 1, 1, 1, 1, 3, 1, 2, 1, 1, 1, 6, 1, 1, 1, 2, 1, 4, 1, 1, 1, 1, 1, 8, 1, 1, 1, 4, 1, 3, 1, 1, 1, 1, 1, 11, 1, 1, 1, 1, 1, 4, 1, 3, 1, 1, 1, 35, 1, 1, 1, 1, 1, 3, 1, 1, 1, 1, 1, 32, 1, 1, 1, 1, 1, 2, 1, 7, 1, 1, 1, 26, 1, 1, 1, 2, 1, 24, 1, 1, 1, 1, 1, 22, 1, 1, 1, 3
Offset: 0

Views

Author

Keywords

Examples

			a(12) = 3 because we have the partitions [12], [6, 4, 2], and [6, 3, 2, 1].
		

Crossrefs

Programs

  • Haskell
    a033630 0 = 1
    a033630 n = p (a027750_row n) n where
       p _  0 = 1
       p [] _ = 0
       p (d:ds) m = if d > m then 0 else p ds (m - d) + p ds m
    -- Reinhard Zumkeller, Feb 23 2014, Apr 04 2012, Oct 27 2011
  • Maple
    with(numtheory): a:=proc(n) local div, g, gser: div:=divisors(n): g:=product(1+x^div[j],j=1..tau(n)): gser:=series(g,x=0,105): coeff(gser,x^n): end: seq(a(n),n=1..100); # Emeric Deutsch, Mar 30 2006
    # second Maple program:
    with(numtheory):
    a:= proc(n) local b, l; l:= sort([(divisors(n))[]]):
          b:= proc(m, i) option remember; `if`(m=0, 1, `if`(i<1, 0,
                 b(m, i-1)+`if`(l[i]>m, 0, b(m-l[i], i-1))))
              end; forget(b):
          b(n, nops(l))
        end:
    seq(a(n), n=0..100); # Alois P. Heinz, Feb 05 2014
  • Mathematica
    A033630 = Table[SeriesCoefficient[Series[Times@@((1 + z^#) & /@ Divisors[n]), {z, 0, n}], n ], {n, 512}] (* Wouter Meeussen *)
    A033630[n_] := f[n, n, 1]; f[n_, m_, k_] := f[n, m, k] = If[k <= m, f[n, m, k + 1] + f[n, m - k, k + 1] * Boole[Mod[n, k] == 0], Boole[m == 0]]; Array[A033630, 101, 0] (* Jean-François Alcover, Jul 29 2015, after Reinhard Zumkeller *)

Formula

a(n) = A065205(n) + 1.
a(A005100(n)) = 1; a(A005835(n)) > 1. - Reinhard Zumkeller, Mar 02 2007
a(n) = f(n, n, 1) with f(n, m, k) = if k <= m then f(n, m, k + 1) + f(n, m - k, k + 1)*0^(n mod k) else 0^m. - Reinhard Zumkeller, Dec 11 2009
a(n) = [x^n] Product_{d|n} (1 + x^d). - Ilya Gutkovskiy, Jul 26 2017
a(n) = 1 if n is deficient (A005100) or weird (A006037). a(n) = 2 if n is perfect (A000396). - Alonso del Arte, Sep 24 2017

Extensions

More terms from Reinhard Zumkeller, Apr 21 2003

A023196 Nondeficient numbers: numbers k such that sigma(k) >= 2k; union of A000396 and A005101.

Original entry on oeis.org

6, 12, 18, 20, 24, 28, 30, 36, 40, 42, 48, 54, 56, 60, 66, 70, 72, 78, 80, 84, 88, 90, 96, 100, 102, 104, 108, 112, 114, 120, 126, 132, 138, 140, 144, 150, 156, 160, 162, 168, 174, 176, 180, 186, 192, 196, 198, 200, 204, 208, 210, 216, 220, 222, 224, 228, 234, 240, 246, 252
Offset: 1

Views

Author

Keywords

Comments

Also called the non-deficient numbers.
If k is a term, so is every positive multiple of k. The "primitive" terms form A006039.
The sequence of numbers k that give local minima for A004125, i.e., such that A004125(k-1) > A004125(k) and A004125(k) < A004125(k+1) coincides with this sequence for the first 1014 terms. Then there appears 4095 which is a term of A023196 but is not a local minimum. - Max Alekseyev, Jan 26 2005
Also, union of pseudoperfect and weird numbers. Cf. A005835, A006037. - Franklin T. Adams-Watters, Mar 29 2006
k is in the sequence if and only if A004125(k-1) > A004125(k). - Jaycob Coleman, Mar 31 2014
Like the abundant numbers, this sequence has density between 0.2474 and 0.2480, see A005101. - Charles R Greathouse IV, Nov 30 2022

Crossrefs

Programs

Showing 1-10 of 75 results. Next