A000396 -id:A000396 - cp's OEIS frontend

A153006 Toothpick sequence starting at the outside corner of an infinite square from which protrudes a half toothpick.

0, 1, 3, 6, 9, 13, 20, 28, 33, 37, 44, 53, 63, 78, 100, 120, 129, 133, 140, 149, 159, 174, 196, 217, 231, 246, 269, 297, 332, 384, 448, 496, 513, 517, 524, 533, 543, 558, 580, 601, 615, 630, 653, 681, 716, 768, 832, 881, 903, 918, 941

Offset: 0

Omar E. Pol, Dec 17 2008, Dec 19 2008, Apr 28 2009

nonn

a(n) is the total number of integer toothpicks after n steps.
It appears that this sequence is related to triangular numbers, Mersenne primes and even perfect numbers. Conjecture: a(A000668(n))=A000217(A000668(n)). Conjecture: a(A000668(n))=A000396(n), assuming there are no odd perfect numbers.
The main entry for this sequence is A139250. See also A152980 (the first differences) and A147646.
The Mersenne prime conjectures are true, but aren't really about Mersenne primes. a(2^i-1) = 2^i (2^i-1)/2 for all i (whether or not i or 2^i-1 is prime). This follows from the formulas for A139250(2^i-1) and A139250(2^i). - David Applegate, May 11 2009
Then we can write a(A000225(k)) = A006516(k), for k > 0. - Omar E. Pol, May 23 2009
Equals A151550 convolved with [1, 2, 2, 2, ...]. (This is equivalent to the observation that the g.f. is x((1+x)/(1-x)) * Product_{n >= 1} (1 + x^(2^n-1) + 2*x^(2^n)).) Equivalently, equals A151555 convolved with A151575. - Gary W. Adamson, May 25 2009
It appears that a(n) is also 1/4 of the total path length of a toothpick structure as A139250 after n-th stage which is constructed following a special rule: toothpicks of the new generation have length 4 when are placed on the square grid (every toothpick has four components of length 1), but after every stage, one (or two) of the four components of every toothpick of the new generation is removed, if such component contains a endpoint of the toothpick and if such endpoint is touching the midpoint or the endpoint of another toothpick. The truncated endpoints of the toothpicks remain exposed forever. Note that there are three sizes of toothpicks in the structure: toothpicks of length 4, 3 and 2. a(n) is also 1/4 of the number of grid points that are covered after n-th stage, except the central point of the structure. A159795 gives the total path length and also the total number of components in the structure after n-th stage. - Omar E. Pol, Oct 24 2011

N. J. A. Sloane, Table of n, a(n) for n = 0..16387
David Applegate, The movie version
David Applegate, Omar E. Pol and N. J. A. Sloane, The Toothpick Sequence and Other Sequences from Cellular Automata, Congressus Numerantium, Vol. 206 (2010), 157-191. [There is a typo in Theorem 6: (13) should read u(n) = 4.3^(wt(n-1)-1) for n >= 2.]
Omar E. Pol, Illustration of initial terms
Omar E. Pol, Illustration of A153006(31) = 496
N. J. A. Sloane, Catalog of Toothpick and Cellular Automata Sequences in the OEIS
Index entries for sequences related to toothpick sequences
Index entries for sequences related to cellular automata

Cf. A000079, A000217, A000396, A000668, A139250, A139251, A139560, A152978, A152979, A152980, A153000, A153001, A153002.
Cf. A000225, A006516.
Cf. A153007, A078008, A151555.

G:=x*((1 + x)/(1 - x)) * mul( (1 + x^(2^n-1) + 2*x^(2^n)), n=1..20); # N. J. A. Sloane, May 20 2009

G.f.: x*((1 + x)/(1 - x)) * Product_{n >= 1} (1 + x^(2^n-1) + 2*x^(2^n)). - N. J. A. Sloane, May 20 2009

Edited by N. J. A. Sloane, Dec 19 2008

A071395 Primitive abundant numbers (abundant numbers all of whose proper divisors are deficient numbers).

Original entry on oeis.org

20, 70, 88, 104, 272, 304, 368, 464, 550, 572, 650, 748, 836, 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, 4030, 4070, 4095, 4216, 4288

Offset: 1

Joe McCauley (mccauley(AT)davesworld.net), Jun 12 2002

nonn

This is a subsequence of the primitive abundant number sequence A091191, since none of these numbers are a positive integer multiple of a perfect number (A000396). - Timothy L. Tiffin, Jul 15 2016
If the terms of this sequence are removed from A091191, then the resulting sequence will be A275082. - Timothy L. Tiffin, Jul 16 2016
Numbers n such that A294927(n) = 0 and A294937(n) = 1. - Antti Karttunen, Nov 14 2017

			20 is a term since 1, 2, 4, 5, and 10 (the proper divisors of 20) are all deficient numbers. - _Timothy L. Tiffin_, Jul 15 2016

Guy, R. K. Unsolved Problems in Number Theory, 2nd ed. New York: Springer-Verlag, p. 46, also section B2, 1994.

Donovan Johnson, Table of n, a(n) for n = 1..10000
Eric Weisstein's World of Mathematics, Primitive Abundant Number

Cf. A006038, A000396, A005100, A005101, subsequence of A091191, A275082.

Cf. A294927, A294937.

abundance:= proc(n) option remember;  numtheory:-sigma(n)-2*n end proc:
select(n -> abundance(n) > 0 and andmap(t -> abundance(t) < 0, numtheory:-divisors(n) minus {n}), [$1..10000]); # Robert Israel, Nov 15 2017

Select[Range@ 5000, DivisorSigma[1, #] > 2 # && Times @@ Boole@ Map[DivisorSigma[1, #] < 2 # &, Most@ Divisors@ #] == 1 &] (* Michael De Vlieger, Jul 16 2016 *)

isA071395(v) = {if (sigma(v) <= 2*v, return (0)); fordiv (v, d, if ((d != v) && (sigma(d) >= 2*d), return (0));); return (1);} \\ Michel Marcus, Mar 10 2013

Offset corrected by Donovan Johnson, Aug 28 2011

A111592 Admirable numbers. A number n is admirable if there exists a proper divisor d' of n such that sigma(n)-2d'=2n, where sigma(n) is the sum of all divisors of n.

Original entry on oeis.org

12, 20, 24, 30, 40, 42, 54, 56, 66, 70, 78, 84, 88, 102, 104, 114, 120, 138, 140, 174, 186, 222, 224, 234, 246, 258, 270, 282, 308, 318, 354, 364, 366, 368, 402, 426, 438, 464, 474, 476, 498, 532, 534, 582, 606, 618, 642, 644, 650, 654, 672, 678, 762, 786, 812

Offset: 1

Jason Earls, Aug 09 2005

All admirable numbers are abundant.
If 2^n-2^k-1 is an odd prime then m=2^(n-1)*(2^n-2^k-1) is in the sequence because 2^k is one of the proper divisors of m and sigma(m)-2m=(2^n-1)*(2^n-2^k)-2^n*(2^n-2^k-1)=2^k hence m=(sigma(m)-m)-2^k, namely m is an Admirable number. This is one of the results of the following theorem that I have found. Theorem: If 2^n-j-1 is an odd prime and m=2^(n-1)*(2^n-j-1) then sigma(m)-2m=j. The case j=0 is well known. - Farideh Firoozbakht, Jan 28 2006
In particular, these numbers have abundancy 2 to 3: 2 < sigma(n)/n <= 3. - Charles R Greathouse IV, Jan 30 2014
Subsequence of A083207. - Ivan N. Ianakiev, Mar 20 2017
The concept of admirable numbers was developed by educator Jerome Michael Sachs (1914-2012) for a television in-service training course in mathematics for elementary school teachers. - Amiram Eldar, Aug 22 2018
Odd terms are listed in A109729. For abundant nonsquares, it is equivalent to say sigma(n)/2 - n divides n. For squares, sigma(n)/2 - n is half-integer, but n could still be an integer multiple. This first occurs for n = m^2 with even m = 2^k*(2^(2*k+1)-1), k = 1, 2, 3, 6, ... (A146768), and odd m = 13167. - M. F. Hasler, Jan 26 2020

			12 = 1+3+4+6-2, 20 = 2+4+5+10-1, etc.

Charles R Greathouse IV, Table of n, a(n) for n = 1..10000
F. Firoozbakht, M. F. Hasler, Variations on Euclid's formula for Perfect Numbers, JIS 13 (2010) #10.3.1.
Giovanni Resta, admirable numbers
J. M. Sachs, Admirable Numbers and Compatible Pairs, The Arithmetic Teacher, Vol. 7, No. 6 (1960), pp. 293-295.
T. Trotter, Admirable Numbers

Subsequence of A005101 (abundant numbers).
Cf. A000396 (perfect numbers), A005100 (deficient numbers), A000203 (sigma), A061645.
Cf. A109729 (odd admirable numbers).

with(numtheory); isadmirable := proc(n) local b, d, S; b:=false; S:=divisors(n) minus {n}; for d in S do if sigma(n)-2*d=2*n then b:=true; break fi od; return b; end: select(proc(z) isadmirable(z) end, [$1..1000]); # Walter Kehowski, Aug 12 2005

fQ[n_] := Block[{d = Most[Divisors[n]], k = 1}, l = Length[d]; s = Plus @@ d; While[k < l && s - 2d[[k]] > n, k++ ]; If[k > l || s != n + 2d[[k]], False, True]]; Select[ Range[821], fQ[ # ] &] (* Robert G. Wilson v, Aug 13 2005 *)
Select[Range[812],MemberQ[Most[Divisors[#]],(DivisorSigma[1,#]-2*#)/2]&] (* Ivan N. Ianakiev, Mar 23 2017 *)

for(n=1,10^3,ap=sigma(n)-2*n;if(ap>0 && (ap%2)==0,d=ap/2;if(d!=n && (n%d)==0, print1(n",")))) \\ Herman Jamke (hermanjamke(AT)fastmail.fm), Mar 30 2008

is(n)=if(issquare(n)||issquare(n/2),0,my(d=sigma(n)/2-n); d>0 && d!=n && n%d==0) \\ Charles R Greathouse IV, Jun 21 2011

Better definition from Walter Kehowski, Aug 12 2005

A027687 4-perfect (quadruply-perfect or sous-triple) numbers: sum of divisors of n is 4n.

Original entry on oeis.org

30240, 32760, 2178540, 23569920, 45532800, 142990848, 1379454720, 43861478400, 66433720320, 153003540480, 403031236608, 704575228896, 181742883469056, 6088728021160320, 14942123276641920, 20158185857531904, 275502900594021408

Offset: 1

Jean-Yves Perrier (nperrj(AT)ascom.ch)

nonn

It is conjectured that there are only finitely many terms. - N. J. A. Sloane, Jul 22 2012
Odd perfect number (unlikely to exist) and infinitely many Mersenne primes will make the sequence infinite - take the product of the OPN and coprime EPNs.
Conjecture: A010888(a(n)) divides a(n). Tested for n up to 36 incl. - Ivan N. Ianakiev, Oct 31 2013
From Farideh Firoozbakht, Dec 26 2014: (Start)
Theorem: If k>1 and p=a(n)/2^(k-2)+1 is prime then for each positive integer m, 2^(k-1)*p^m is a solution to the equation sigma(phi(x))=2*x-2^k, which implies the equation has infinitely many solutions.
Proof: sigma(phi(2^(k-1)*p^m)) = sigma(2^(k-2)*(p-1)*p^(m-1)) = sigma(2^(k-2)*(p-1))*sigma(p^(m-1)) = sigma(a(n))*(p^m-1)/(p-1) = 4*a(n)*(p^m-1)/(p-1) = 2^k*(p^m-1) = 2*(2^(k-1)*p^m)-2^k.
It seems that for all such equations there exist such an infinite set of solutions. So I conjecture that the sequence is infinite! (End)
If 3 were prepended to this sequence, then it would be the sequence of integers k such that numerator(sigma(k)/k)=4. - Michel Marcus, Nov 22 2015

			From _Daniel Forgues_, May 09 2010: (Start)
30240 = 2^5*3^3*5*7
sigma(30240) = (2^6-1)/1*(3^4-1)/2*(5^2-1)/4*(7^2-1)/6
= (63)*(40)*(6)*(8)
= (7*3^2)*(2^3*5)*(2*3)*(2^3)
= 2^7*3^3*5*7
= (2^2) * (2^5*3^3*5*7)
= 4 * 30240 (End)

R. K. Guy, Unsolved Problems in Number Theory, B2.
James J. Tattersall, Elementary Number Theory in Nine Chapters, Cambridge University Press, 1999, page 143.

T. D. Noe, Table of n, a(n) for n=1..36 (complete sequence from Flammenkamp)
Abiodun E. Adeyemi, A Study of @-numbers, arXiv:1906.05798 [math.NT], 2019.
Kevin A. Broughan and Qizhi Zhou, Divisibility by 3 of even multiperfect numbers of abundancy 3 and 4, JIS 13 (2010) 10.1.5
Seth Colbert-Pollack, Judy Holdener, Emily Rachfal, and Yanqi Xu, A DIY Project: Construct Your Own Multiply Perfect Number!, Math Horizons, Vol. 28, pp. 20-23, February 2021.
Farideh Firoozbakht and Maxmilian F. Hasler, Variations on Euclid's formula for Perfect Numbers, JIS 13 (2010) #10.3.1.
Achim Flammenkamp, The Multiply Perfect Numbers Page
Shyam Sunder Gupta, Perfect, Multiply Perfect, and Sociable Numbers, Exploring the Beauty of Fascinating Numbers, Springer (2025) Ch. 6, 185-207.
Fred Helenius, Link to Glossary and Lists
Walter Nissen, Abundancy : Some Resources
Eric Weisstein's World of Mathematics, Multiperfect Number.
Eric Weisstein's World of Mathematics, Sous-Triple.
Wikipedia, Multiply perfect number

Cf. A000396, A005820, A007539, A046060, A046061.

AbundantQ[n_]:=DivisorSigma[1, n]==4*n;a={};Do[If[AbundantQ[n], AppendTo[a, n]], {n, 10^6}];a (* Vladimir Joseph Stephan Orlovsky, Aug 16 2008 *)

4 more terms from Labos Elemer

A246282 Numbers k for which A003961(k) > 2k; numbers n such that if n = Product_{k >= 1} (p_k)^(c_k), then Product_{k >= 1} (p_{k+1})^(c_k) > 2n, where p_k indicates the k-th prime, A000040(k).

Original entry on oeis.org

4, 6, 8, 9, 10, 12, 14, 15, 16, 18, 20, 21, 24, 27, 28, 30, 32, 35, 36, 39, 40, 42, 44, 45, 48, 49, 50, 52, 54, 56, 57, 60, 63, 64, 66, 68, 69, 70, 72, 75, 76, 78, 80, 81, 84, 88, 90, 91, 92, 96, 98, 99, 100, 102, 104, 105, 108, 110, 112, 114, 116, 117, 120, 124, 125, 126, 128, 130, 132, 135, 136, 138, 140, 144

Offset: 1

Antti Karttunen, Aug 24 2014

Numbers n such that A003961(n) > 2*n.
Numbers n such that A048673(n) > n.
The sequence grows as:
        a(10) = 18
       a(100) = 192
      a(1000) = 1830
     a(10000) = 18636
    a(100000) = 187350
   a(1000000) = 1865226
  a(10000000) = 18654333
and the powers of 10 occur at:
        a(5) = 10
       a(53) = 100
      a(536) = 1000
     a(5423) = 10000
    a(53290) = 100000
   a(535797) = 1000000
  a(5361886) = 10000000
suggesting that the ratio a(n)/n is converging to an constant and an arbitrary natural number is slightly more likely to be in this sequence than in the complement A246281. See also comments at A246351 and compare to quite a different ratio present in the "inverse" case A246362.
From Antti Karttunen, Aug 27 2020: (Start)
Any perfect number, including all odd perfect numbers (if such numbers exist), must occur in this sequence. See A286385 and A326042 for the reason why.
Like abundancy index (ratio A000203(n)/n), also ratio A003961(n)/n is multiplicative and always > 1 for all n > 1. Thus if the number has a proper divisor that is in this sequence, then the number itself also is. See A337372 for terms included here, but with no proper divisor in this sequence. (End)
For k >= 2, if m * A130789(k) is a term then m * A130789(k-1) is a term. - Peter Munn, Sep 01 2025
Could be called "primeshift-abundant numbers", in analogy with A005101. - Antti Karttunen, Sep 01 2025

			3 = p_2 (3 is the second prime, A000040(2)) is not a member, because p_3 = 5 (5 is the next prime after 3, A000040(3)) and 5/3 < 2.
4 = 2*2 = p_1 * p_1 is a member, as p_2 * p_2 = 3*3 = 9, and 9/4 > 2.
33 = 3*11 = p_2 * p_5 is not a member, as p_3 * p_6 = 5*13 = 65, and 65/33 < 2.
35 = 5*7 = p_3 * p_4 is a member, as p_4 * p_5 = 7*11 = 77, and 77/35 > 2.

Complement: A246281.
Setwise difference of A246352 and A048674.
Cf. A000040, A003961, A048673, A130789, A246362, A252742 (characteristic function), A286385, A326042, A337345.
Positions of positive terms in A252748 and in A337345.
Union of A337372 (primitive terms), A341610 (non-primitive terms).
Subsequences: A000396, (A005101 U A341614), (A023196 U A341615), A326134, A337373, A337374, A337378, A337381, A337384, A337386, A337543, A341611.
Cf. also A275717, A275718.

Select[Range[144], 2 # < Times @@ Map[#1^#2 & @@ # &, FactorInteger[#] /. {p_, e_} /; e > 0 :> {Prime[PrimePi@ p + 1], e}] - Boole[# == 1] &] (* Michael De Vlieger, Feb 22 2021 *)

A003961(n) = { my(f = factor(n)); for (i=1, #f~, f[i, 1] = nextprime(f[i, 1]+1)); factorback(f); }; \\ From A003961
isA246282(n) = (A003961(n) > (n+n));
n = 0; i = 0; while(i < 10000, n++; if(isA246282(n), i++; write("b246282.txt", i, " ", n)));

;; With Antti Karttunen's IntSeq-library.
(define A246282 (MATCHING-POS 1 1 (lambda (n) (> (A003961 n) (* 2 n)))))

;; With Antti Karttunen's IntSeq-library.
(define A246282 (MATCHING-POS 1 1 (lambda (n) (> (A048673 n) n))))

A new shorter version of name prepended by Antti Karttunen, Aug 27 2020

A336702 Numbers whose abundancy index is a power of 2.

Original entry on oeis.org

1, 6, 28, 496, 8128, 30240, 32760, 2178540, 23569920, 33550336, 45532800, 142990848, 1379454720, 8589869056, 43861478400, 66433720320, 137438691328, 153003540480, 403031236608, 704575228896, 181742883469056, 6088728021160320, 14942123276641920, 20158185857531904, 275502900594021408, 622286506811515392, 2305843008139952128

Offset: 1

Antti Karttunen, Aug 05 2020

nonn

Apart from missing 2, this sequence gives all numbers k such that the binary expansion of A156552(k) is a prefix of that of A156552(sigma(k)), that is, for k > 1, numbers k for which sigma(k) is a descendant of k in A005940-tree. This follows because of the two transitions x -> A005843(x) (doubling) and x -> A003961(x) (prime shift) used to generate descendants in A005940-tree, using A003961 at any step of the process will ruin the chances of encountering sigma(k) anywhere further down that subtree.
Proof: Any left child in A005940 (i.e., A003961(k) for k) is larger than sigma(k), for any k > 2 [see A286385 for a proof], and A003961(n) > n for all n > 1. Thus, apart from A003961(2) = 3 = sigma(2), A003961^t(k) > sigma(k), where A003961^t means t-fold application of prime shift, here with t >= 1. On the other hand, sigma(2n) > sigma(n) for all n, thus taking first some doubling steps before a run of one or more prime shift steps will not rescue us, as neither will taking further doubling steps after a bout of prime shifts.
The first terms of A325637 not included in this sequence are 154345556085770649600 and 9186050031556349952000, as they have abundancy index 6.
From Antti Karttunen, Nov 29 2021: (Start)
Odd terms of this sequence are given by the intersection of A349169 and A349174.
A064989 applied to the odd terms of this sequence gives the fixed points of A326042, i.e., the positions of zeros in A348736, and a subset of the positions of ones in A348941.
Odd terms of this sequence form a subsequence of A348943, but should occur neither in A348748 nor in A348749.
(End)

			For 30240, sigma(30240) = 120960 = 4*30240, therefore, as sigma(k)/k = 2^2, a power of two, 30240 is present.

Antti Karttunen, Table of n, a(n) for n = 1..769 (computed from the b-file of A007691 prepared by T. D. Noe, using Flammenkamp's data)
Index entries for sequences where odd perfect numbers must occur, if they exist at all
Index entries for sequences related to sigma(n)

Cf. A000203, A003961, A005843, A005940, A156552, A163511, A209229, A286385.
Cf. A000396, A027687 (subsequences).
Subsequence of A007691, and after 1, also subsequence of A325637.
Union with {2} gives the positions of zeros in A347381.
Cf. also A326042, A347391, A347392, A347393, A347394, A348736, A348748, A348749, A348941, A348943, A349169, A349174.

isA336702(n) = { my(r=sigma(n)/n); (1==denominator(r)&&!bitand(r, r-1)); }; \\ (Corrected) - Antti Karttunen, Aug 31 2021

A046060 5-multiperfect numbers.

Original entry on oeis.org

14182439040, 31998395520, 518666803200, 13661860101120, 30823866178560, 740344994887680, 796928461056000, 212517062615531520, 69357059049509038080, 87934476737668055040, 170206605192656148480, 1161492388333469337600, 1245087725796543283200, 1802582780370364661760

Offset: 1

Eric W. Weisstein

nonn

Conjectured finite and probably these are the only terms; cf. Flammenkamp's link. [Georgi Guninski, Jul 25 2012]

			From _Daniel Forgues_, May 09 2010: (Start)
14182439040 = 2^7*3^4*5*7*11^2*17*19
sigma(14182439040)
  = (2^8-1)/1*(3^5-1)/2*(5^2-1)/4*(7^2-1)/6*(11^3-1)/10*(17^2-1)/16*(19^2-1)/18
  = (255)*(121)*(6)*(8)*(133)*(18)*(20)
  = (3*5*17)*(11^2)*(2*3)*(2^3)*(7*19)*(2*3^2)*(2^2*5)
  = 2^7*3^4*5^2*7*11^2*17*19
  = (5) * (2^7*3^4*5*7*11^2*17*19)
  = 5 * 14182439040 (End)

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

T. D. Noe, Table of n, a(n) for n = 1..65 (complete sequence from Flammenkamp)
Farideh Firoozbakht and Maximilian F. Hasler, Variations on Euclid's formula for Perfect Numbers, JIS 13 (2010) #10.3.1.
Achim Flammenkamp, The Multiply Perfect Numbers Page.
Shyam Sunder Gupta, Perfect, Multiply Perfect, and Sociable Numbers, Exploring the Beauty of Fascinating Numbers, Springer (2025) Ch. 6, 185-207.
Fred Helenius, Link to Glossary and Lists.
Walter Nissen, Abundancy : Some Resources .
Eric Weisstein's World of Mathematics, Multiperfect Number.

Cf. A000396, A005820, A027687, A046061, A007539.

is(n)=sigma(n)==5*n \\ Charles R Greathouse IV, Apr 05 2013

A006038 Odd primitive abundant numbers.

Original entry on oeis.org

945, 1575, 2205, 3465, 4095, 5355, 5775, 5985, 6435, 6825, 7245, 7425, 8085, 8415, 8925, 9135, 9555, 9765, 11655, 12705, 12915, 13545, 14805, 15015, 16695, 18585, 19215, 19635, 21105, 21945, 22365, 22995, 23205, 24885, 25935, 26145, 26565, 28035, 28215

Offset: 1

N. J. A. Sloane

nonn

Dickson proves that there are only a finite number of odd primitive abundant numbers having n distinct prime factors. Sequence A188342 lists the smallest such numbers. - T. D. Noe, Mar 28 2011
Sequence A188439 sorts the numbers in this sequence by the number of distinct prime factors. Eight numbers have exactly three prime factors; 576 have exactly four prime factors. - T. D. Noe, Apr 04 2011
Any multiple of an abundant number (A005101) is again an abundant number. Primitive abundant numbers (A091191) are those not of this form, i.e., without an abundant proper divisor. We don't know any odd perfect number (A000396), so the (odd) terms here have only deficient proper divisors (A071395), and their prime factors p are less than sigma(n/p)/deficiency(n/p). See A005231 (odd abundant numbers) for an explanation why all terms have at least 3 distinct prime factors, and 5 prime factors when counted with multiplicity (A001222), whence a(1) = 3^3*5*7. All known terms are semiperfect (A005835, and thus in A006036): no odd weird number (A006037) is known, but if it exists, the smallest one is in this sequence. - M. F. Hasler, Jul 28 2016
So far, a(173) = 351351 is the only known term of A122036, i.e., which can't be written as sum of its proper divisors > 1. - M. F. Hasler, Jan 26 2020

N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

T. D. Noe, Table of n, a(n) for n = 1..10000
L. E. Dickson, Finiteness of the odd perfect and primitive abundant numbers with n distinct prime factors, American Journal of Mathematics 35 (1913), pp. 413-422.
R. K. Guy, Letter to N. J. A. Sloane with attachment, Jun. 1991
Jacob Liddy, An algorithm to determine all odd primitive abundant numbers with d prime divisors, Honors Research Projects (2018), 728.
Eric Weisstein's World of Mathematics, Primitive Abundant Number

Cf. A005101, A005231. Subsequence of A091191.

Cf. A000203, A027751, A379949 (subsequence of square terms).

a006038 n = a006038_list !! (n-1)
a006038_list = filter f [1, 3 ..] where
   f x = sum pdivs > x && all (<= 0) (map (\d -> a000203 d - 2 * d) pdivs)
         where pdivs = a027751_row x
-- Reinhard Zumkeller, Jan 31 2014

isA005101 := proc(n) is(numtheory[sigma](n) > 2*n ); end proc:
isA005100 := proc(n) is(numtheory[sigma](n) < 2*n ); end proc:
isA006038 := proc(n) local d; if type(n,'odd') and isA005101(n) then for d in numtheory[divisors](n) minus {1,n} do if not isA005100(d) then return false; end if; end do: return true;else false; end if; end proc:
n := 1 ; for a from 1 by 2 do if isA006038(a) then printf("%d %d\n",n,a) ; n := n+1 ; end if; end do: # R. J. Mathar, Mar 28 2011

t = {}; n = 1; While[Length[t] < 50, n = n + 2; If[DivisorSigma[1, n] > 2 n && Intersection[t, Divisors[n]] == {}, AppendTo[t, n]]]; t (* T. D. Noe, Mar 28 2011 *)

is(n)=n%2 && sumdiv(n,d,sigma(d,-1)>2)==1 \\ Charles R Greathouse IV, Jun 10 2013

is_A006038(n)=bittest(n,0) && sigma(n)>2*n && !for(i=1,#f=factor(n)[,1],sigma(n\f[i],-1)>2&&return) \\ More than 5 times faster. - M. F. Hasler, Jul 28 2016

A323244 a(1) = 0; and for n > 1, a(n) = A033879(A156552(n)).

Original entry on oeis.org

0, 1, 1, 2, 1, 4, 1, 6, 0, 5, 1, 10, 1, 16, 2, 6, 1, 12, 1, 18, -3, 18, 1, 22, -4, 46, 4, 22, 1, 10, 1, 30, 14, 82, -2, 14, 1, 256, -12, 22, 1, 36, 1, 66, 8, 226, 1, 46, -12, 19, 8, 130, 1, 28, -19, 70, -12, 748, 1, 42, 1, 1362, 16, 22, 10, 42, 1, 214, 254, 40, 1, 38, 1, 3838, 10, 406, -10, 106, 1, 78, -12, 5458, 1, 26, -72, 12250, -348, 30, 1, 12

Offset: 1

Antti Karttunen, Jan 10 2019

sign

After a(1) = 0, the other zeros occur for k >= 1, at A005940(1+A000396(k)), which, provided no odd perfect numbers exist, is equal to A324201(k) = A062457(A000043(k)): 9, 125, 161051, 410338673, ..., etc.

There are 2321 negative terms among the first 10000 terms.

Antti Karttunen, Table of n, a(n) for n = 1..10000 (based on Hans Havermann's factorization of A156552)
Index entries for sequences related to binary expansion of n
Index entries for sequences computed from indices in prime factorization
Index entries for sequences related to sigma(n)

Cf. A000043, A000396, A033879, A064989, A156552, A297112, A323240, A323243, A323245, A323248, A324115, A324051, A324103, A324396, A324398, A324543, A324713.
Cf. A324201 (positions of zeros, conjectured), A324551 (of negative terms), A324720 (of nonnegative terms), A324721 (of positive terms), A324731, A324732.
Cf. A329644 (Möbius transform).
Cf. A323174, A324055, A324185, A324546 for other permutations of deficiency, and also A324574, A324575, A324654.

Array[2 # - If[# == 0, 0, DivisorSigma[1, #]] &@ Floor@ Total@ Flatten@ MapIndexed[#1 2^(#2 - 1) &, Flatten[Table[2^(PrimePi@ #1 - 1), {#2}] & @@@ FactorInteger@ #]] &, 90] (* Michael De Vlieger, Apr 21 2019 *)

A064989(n) = {my(f); f = factor(n); if((n>1 && f[1,1]==2), f[1,2] = 0); for (i=1, #f~, f[i,1] = precprime(f[i,1]-1)); factorback(f)};
A156552(n) = if(1==n, 0, if(!(n%2), 1+(2*A156552(n/2)), 2*A156552(A064989(n))));
A323244(n) = if(1==n, 0, my(k=A156552(n)); (2*k)-sigma(k));

from sympy import divisor_sigma, primepi, factorint
def A323244(n): return (lambda n: (n<<1)-divisor_sigma(n))(sum((1< 1 else 0 # Chai Wah Wu, Mar 10 2023

a(n) = 2*A156552(n) - A323243(n).
a(1) = 0; and for n > 1, a(n) = A033879(A156552(n)).
a(n) = A323248(n) + A001222(n) = (A323247(n) - A323243(n)) + A001222(n).
From Antti Karttunen, Mar 12 2019 & Nov 23 2019: (Start)
a(n) = Sum_{d|n} (2*A297112(d) - A324543(d)) = Sum_{d|n} A329644(d).
A002487(a(n)) = A324115(n).
a(n) = A329638(n) - A329639(n).
a(n) = A329645(n) - A329646(n).
(End)

A263837 Non-abundant numbers (or nonabundant numbers): complement of A005101; numbers k for which sigma(k) <= 2*k.

Original entry on oeis.org

1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 13, 14, 15, 16, 17, 19, 21, 22, 23, 25, 26, 27, 28, 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, 71, 73, 74, 75, 76, 77, 79, 81, 82, 83, 85, 86, 87, 89, 91, 92, 93, 94, 95, 97, 98, 99

Offset: 1

N. J. A. Sloane, Oct 30 2015

nonn

For all n, A003961(a(n)) is in A005100. - Antti Karttunen, Aug 28 2020

Reinhard Zumkeller, Table of n, a(n) for n = 1..10000
Paul Pollack and Carl Pomerance, Some problems of Erdős on the sum-of-divisors function, For Richard Guy on his 99th birthday: May his sequence be unbounded, Trans. Amer. Math. Soc. Ser. B 3 (2016), 1-26.

Union of A000396 and A005100.
Cf. A005101 (complement), A023196, A294935 (characteristic function).
Cf. A000203, A001065, A003961.

a263837 n = a263837_list !! (n-1)
a263837_list = filter (\x -> a001065 x <= x) [1..]
-- Reinhard Zumkeller, Oct 31 2015

isA263837 := proc(n)
    if 2*n-numtheory[sigma](n) >=0 then
        true;
    else
        false;
    end if;
end proc:
A263837 := proc(n)
    option remember;
    local a;
    if n =1 then
        1;
    else
        for a from procname(n-1)+1 do
            if isA263837(a) then
                return a;
            end if;
        end do:
    end if;
end proc:
seq(A263837(n),n=1..100) ; # R. J. Mathar, Jun 06 2024

Select[Range[100], DivisorSigma[1, #] <= 2*# &] (* Amiram Eldar, Mar 14 2024 *)

isok(n) = sigma(n) <= 2*n; \\ Michel Marcus, Dec 27 2015

A001065(a(n)) <= a(n). - Reinhard Zumkeller, Oct 31 2015

Additional description to the definition added by Antti Karttunen, Aug 28 2020

cp's OEIS Frontend

A153006 Toothpick sequence starting at the outside corner of an infinite square from which protrudes a half toothpick.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Maple

Formula

Extensions

A071395 Primitive abundant numbers (abundant numbers all of whose proper divisors are deficient numbers).

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Maple

Mathematica

PARI

Extensions

A111592 Admirable numbers. A number n is admirable if there exists a proper divisor d' of n such that sigma(n)-2d'=2n, where sigma(n) is the sum of all divisors of n.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

Mathematica

PARI

PARI

Extensions

A027687 4-perfect (quadruply-perfect or sous-triple) numbers: sum of divisors of n is 4n.

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Mathematica

Extensions

A246282 Numbers k for which A003961(k) > 2*k; numbers n such that if n = Product_{k >= 1} (p_k)^(c_k), then Product_{k >= 1} (p_{k+1})^(c_k) > 2*n, where p_k indicates the k-th prime, A000040(k).

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

PARI

Scheme

Scheme

Extensions

A336702 Numbers whose abundancy index is a power of 2.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

PARI

A046060 5-multiperfect numbers.

Views

Author

Keywords

A246282 Numbers k for which A003961(k) > 2k; numbers n such that if n = Product_{k >= 1} (p_k)^(c_k), then Product_{k >= 1} (p_{k+1})^(c_k) > 2n, where p_k indicates the k-th prime, A000040(k).