A033880 -id:A033880 - cp's OEIS frontend

A357976 Numbers with a divisor having the same sum of prime indices as their quotient.

1, 4, 9, 12, 16, 25, 30, 36, 40, 48, 49, 63, 64, 70, 81, 84, 90, 100, 108, 112, 120, 121, 144, 154, 160, 165, 169, 192, 196, 198, 210, 220, 225, 252, 256, 264, 270, 273, 280, 286, 289, 300, 324, 325, 336, 351, 352, 360, 361, 364, 390, 400, 432, 441, 442, 448

Offset: 1

Gus Wiseman, Oct 26 2022

nonn

A prime index of n is a number m such that prime(m) divides n. The multiset of prime indices of n is row n of A112798.

			The terms together with their prime indices begin:
   1: {}
   4: {1,1}
   9: {2,2}
  12: {1,1,2}
  16: {1,1,1,1}
  25: {3,3}
  30: {1,2,3}
  36: {1,1,2,2}
  40: {1,1,1,3}
  48: {1,1,1,1,2}
  49: {4,4}
For example, 40 has factorization 8*5, and both factors have the same sum of prime indices 3, so 40 is in the sequence.

Robert Israel, Table of n, a(n) for n = 1..10000

The partitions with these Heinz numbers are counted by A002219.
A subset of A300061.
The squarefree case is A357854, counted by A237258.
Positions of nonzero terms in A357879.
A001222 counts prime factors, distinct A001221.
A056239 adds up prime indices, row sums of A112798.
Cf. A033879, A033880, A064914, A181819, A213086, A235130, A237194, A276107, A300273, A321144, A357975.

filter:= proc(n) local F,s,t,i,R;
  F:= ifactors(n)[2];
  F:= map(t -> [numtheory:-pi(t[1]),t[2]], F);
  s:= add(t[1]*t[2],t=F)/2;
  if not s::integer then return false fi;
  try
  R:= Optimization:-Maximize(0, [add(F[i][1]*x[i],i=1..nops(F)) = s, seq(x[i]<= F[i][2],i=1..nops(F))], assume=nonnegint, depthlimit=20);
  catch "no feasible integer point found; use feasibilitytolerance option to adjust tolerance": return false;
  end try;
  true
end proc:
filter(1):= true:
select(filter, [$1..1000]); # Robert Israel, Oct 26 2023

sumprix[n_]:=Total[Cases[FactorInteger[n],{p_,k_}:>k*PrimePi[p]]];
Select[Range[100],MemberQ[sumprix/@Divisors[#],sumprix[#]/2]&]

A181595 Abundant numbers n for which the abundance d = sigma(n) - 2*n is a proper divisor, that is, 0 < d < n and d | n.

Original entry on oeis.org

12, 18, 20, 24, 40, 56, 88, 104, 196, 224, 234, 368, 464, 650, 992, 1504, 1888, 1952, 3724, 5624, 9112, 11096, 13736, 15376, 15872, 16256, 17816, 24448, 28544, 30592, 32128, 77744, 98048, 122624, 128768, 130304, 174592, 396896, 507392, 521728, 522752, 537248

Offset: 1

Vladimir Shevelev, Nov 01 2010

nonn

Named near-perfect numbers by sequence author.
Union of this sequence and A005820 is A153501.
Every even perfect number n = 2^(p-1)*(2^p-1), p and 2^p-1 prime, of A000396 generates three entries: 2*n, 2^p*n and (2^p-1)*n.
Every number M=2^(t-1)*P, where P is a prime of the form 2^t-2^k-1, is an entry for which (2^k)|M and sigma(M)-2^k=2*M (see A181701).
Conjecture 1: For every k>=1, there exist infinitely many entries m for which (2^k)|m and sigma(m)-2^k = 2*m.
Conjecture 2. All entries are even. [Proved to be false, see below. (Ed.)]
Conjecture 3. If the suitable (according to the definition) divisor d of an entry is not a power of 2, then it is not suitable divisor for any other entry.
Conjecture 4. If a suitable divisor for an even entry is odd, then it is a Mersenne prime (A000043).
If Conjectures 3 and 4 are true, then an entry with odd suitable divisor has the form 2^(p-1)*(2^p-1)^2, where p and 2^p-1 are primes. - Vladimir Shevelev, Nov 08 2010 to Dec 16 2010
The only odd term in this sequence < 2*10^12 is 173369889. - Donovan Johnson, Feb 15 2012
173369889 remains only odd term up to 1.4*10^19. - Peter J. C. Moses, Mar 05 2012
These numbers are obviously pseudoperfect (A005835) since they are equal to the sum of all the proper divisors except the one that is the same as the abundance. - Alonso del Arte, Jul 16 2012

			The abundance of 12 is A033880(12) = 4, which is a proper divisor of 12, so 12 is in the sequence.

Donovan Johnson, Table of n, a(n) for n = 1..200
Hùng Việt Chu, Divisibility of Divisor Functions of Even Perfect Numbers, J. Int. Seq., Vol. 24 (2021), Article 21.3.4.
Yanbin Li and Qunying Liao, A class of new near-perfect numbers, J. Korean Math. Soc. 52 (2015), No. 4, pp. 751-763.
Paul Pollack and Vladimir Shevelev, On perfect and near-perfect numbers, J. Number Theory 132 (2012), pp. 3037-3046. arXiv preprint, arXiv:1011.6160 [math.NT], 2010-2012.
X.-Z. Ren, Y.-G. Chen, On near-perfect numbers with two distinct prime factors, Bulletin of the Australian Mathematical Society, No 3 (2013), available on CJO2013. doi:10.1017/S0004972713000178.
M. Tang, X. Z. Ren and M. Li, On near-perfect and deficient-perfect numbers, Colloq. Math. 133 (2013), 221-226.

Cf. A000396, A005101, A153501, A005820.

q:= n-> (t-> t>0 and tAlois P. Heinz, May 11 2023

Select[Range[550000], 0 < (d = DivisorSigma[1, #] - 2*#) < # && Divisible[#, d] &] (* Amiram Eldar, May 12 2023 *)

is_A181595(n)=my(d=sigma(n)-2*n); (d>0) && (dA181595(n)&&print1(n","))  \\ M. F. Hasler, Apr 14 2012; corrected by Michel Marcus, May 12 2023

Definition shortened, entries checked by R. J. Mathar, Nov 17 2010

A088833 Numbers n whose abundance is 8: sigma(n) - 2n = 8.

Original entry on oeis.org

56, 368, 836, 11096, 17816, 45356, 77744, 91388, 128768, 254012, 388076, 2087936, 2291936, 13174976, 29465852, 35021696, 45335936, 120888092, 260378492, 381236216, 775397948, 3381872252, 4856970752, 6800228816, 8589344768, 44257207676, 114141404156

Offset: 1

Labos Elemer, Oct 28 2003

nonn

A subset of A045770.
If p=2^m-9 is prime (m is in the sequence A059610) then n=2^(m-1)*p is in the sequence. See comment lines of the sequence A088831. 56, 368, 128768, 2087936 & 8589344768 are of the mentioned form. - Farideh Firoozbakht, Feb 15 2008
a(28) > 10^12. - Donovan Johnson, Dec 08 2011
a(31) > 10^13. - Giovanni Resta, Mar 29 2013
a(38) > 10^18. - Hiroaki Yamanouchi, Aug 23 2018
Any term x of this sequence can be combined with any term y of A125247 to satisfy the property (sigma(x)+sigma(y))/(x+y) = 2, which is a necessary (but not sufficient) condition for two numbers to be amicable. - Timothy L. Tiffin, Sep 13 2016

			Except first 4 terms of A045770 (1, 7, 10, and 49) are here: abundances = {-1,-6,-2,-41,8,8,8,8,8,8,8,8,8,8,8,8,8}.

Giovanni Resta and Hiroaki Yamanouchi, Table of n, a(n) for n = 1..37 (terms a(1)-a(30) from _Giovanni Resta_)

Cf. A033880, A045668, A045669, A045770, A088831, A088832, A059610, A125247 (deficiency 8).

Do[If[DivisorSigma[1,n]==2n+8,Print[n]],{n,100000000}] (* Farideh Firoozbakht, Feb 15 2008 *)

is(n)=sigma(n)==2*n+8 \\ Charles R Greathouse IV, Feb 21 2017

a(14)-a(17) from Farideh Firoozbakht, Feb 15 2008
a(18)-a(25) from Donovan Johnson, Dec 23 2008
a(26)-a(27) from Donovan Johnson, Dec 08 2011

A357854 Squarefree numbers with a divisor having the same sum of prime indices as their quotient.

Original entry on oeis.org

1, 30, 70, 154, 165, 210, 273, 286, 390, 442, 462, 561, 595, 646, 714, 741, 858, 874, 910, 1045, 1155, 1173, 1254, 1326, 1330, 1334, 1495, 1653, 1771, 1794, 1798, 1870, 1938, 2139, 2145, 2294, 2415, 2465, 2470, 2530, 2622, 2639, 2730, 2926, 2945, 2958, 3034

Offset: 1

Gus Wiseman, Oct 27 2022

nonn

A prime index of n is a number m such that prime(m) divides n. The multiset of prime indices of n is row n of A112798.

			The terms together with their prime indices begin:
     1: {}
    30: {1,2,3}
    70: {1,3,4}
   154: {1,4,5}
   165: {2,3,5}
   210: {1,2,3,4}
   273: {2,4,6}
   286: {1,5,6}
   390: {1,2,3,6}
For example, 210 has factorization 14*15, and both factors have the same sum of prime indices 5, so 210 is in the sequence.

The partitions with these Heinz numbers are counted by A237258.
A subset of A319241, squarefree case of A300061.
Squarefree positions of nonzero terms in A357879.
This is the squarefree case of A357976, counted by A002219.
A001222 counts prime factors, distinct A001221.
A056239 adds up prime indices, row sums of A112798.
Cf. A033879, A033880, A064914, A181819, A235130, A237194, A276107, A300273, A321144, A357975, A357976.

sumprix[n_]:=Total[Cases[FactorInteger[n],{p_,k_}:>k*PrimePi[p]]];
Select[Range[1000],SquareFreeQ[#]&&MemberQ[sumprix/@Divisors[#],sumprix[#]/2]&]

A077374 Odd numbers m whose abundance by absolute value is at most 10, that is, -10 <= sigma(m) - 2m <= 10.

Original entry on oeis.org

1, 3, 5, 7, 9, 11, 15, 21, 315, 1155, 8925, 32445, 442365, 815634435

Offset: 1

Jason Earls, Nov 30 2002

Apart from {1, 3, 5, 7, 9, 11, 15, 21, 315}, subset of A088012. Probably finite. - Charles R Greathouse IV, Mar 28 2011
a(15) > 10^13. - Giovanni Resta, Mar 29 2013
The abundance of the given terms a(1..14) is: (-1, -2, -4, -6, -5, -10, -6, -10, -6, -6, 6, 6, 6, -6). See also A171929, A188263 and A188597 for numbers with abundancy sigma(n)/n close to 2. - M. F. Hasler, Feb 21 2017
a(15) > 10^22. - Wenjie Fang, Jul 13 2017

			sigma(32445) = 64896 and 32445*2 = 64890, which makes the odd number 32445 six away from perfection: A(32445) = 6 and hence in this sequence.

Eric Weisstein's World of Mathematics, Abundance.

Contains odd terms of: A191363, A088831, A125246, A088832, A141548, A087167, A125247, A088833, A101223, A223609.

Cf. A033879, A033880, A000203, A088012.

Select[Range[1, 10^6, 2], -10 <= DivisorSigma[1, #] - 2 # <= 10 &] (* Michael De Vlieger, Feb 22 2017 *)

forstep(n=1,442365,2,if(abs(sigma(n)-2*n)<=10,print1(n,",")))

a(14) from Farideh Firoozbakht, Jan 12 2004

A141548 Numbers n whose deficiency is 6.

Original entry on oeis.org

7, 15, 52, 315, 592, 1155, 2102272, 815634435

Offset: 1

Vladimir Joseph Stephan Orlovsky, Aug 16 2008

a(9) > 10^12. - Donovan Johnson, Dec 08 2011
a(9) > 10^13. - Giovanni Resta, Mar 29 2013
a(9) > 10^18. - Hiroaki Yamanouchi, Aug 21 2018
For all k in A059242, the number m = 2^(k-1)*(2^k+5) is in this sequence. This yields further terms 2^46*(2^47+5), 2^52*(2^53+5), 2^140*(2^141+5), ... All even terms known so far and the initial 7 = 2^0*(2^1+5) are of this form. All odd terms beyond a(2) are of the form a(n) = a(k)*p*q, k < n. We have proved that there is no further term of this form with the a(k) given so far. - M. F. Hasler, Apr 23 2015
A term n of this sequence multiplied by a prime p not dividing it is abundant if and only if p < sigma(n)/6 = n/3-1. For the even terms 592 and 2102272, there is such a prime near this limit (191 resp. 693571) such that n*p is a primitive weird number, cf. A002975. For a(3)=52, the largest such prime, 11, is already too small. Odd weird numbers do not exist within these limits. - M. F. Hasler, Jul 19 2016
Any term x of this sequence can be combined with any term y of A087167 to satisfy the property (sigma(x)+sigma(y))/(x+y) = 2, which is a necessary (but not sufficient) condition for two numbers to be amicable. - Timothy L. Tiffin, Sep 13 2016

			a(1) = 7, since 2*7 - sigma(7) = 14 - 8 = 6. - _Timothy L. Tiffin_, Sep 13 2016

Gianluca Amato, Maximilian Hasler, Giuseppe Melfi, Maurizio Parton, Primitive weird numbers having more than three distinct prime factors, Riv. Mat. Univ. Parma, 7(1), (2016), 153-163, arXiv:1803.00324 [math.NT], 2018.

Cf. A087485 (odd terms).
Cf. A000203, A033880, A005100; A191363 (deficiency 2), A125246 (deficiency 4), A141548 (deficiency 6), A125247 (deficiency 8), A101223 (deficiency 10), A141549 (deficiency 12), A141550 (deficiency 14), A125248 (deficiency 16), A223608 (deficiency 18), A223607 (deficiency 20).
Cf. A087167 (abundance 6).

[n: n in [1..9*10^6] | (SumOfDivisors(n)-2*n) eq -6]; // Vincenzo Librandi, Sep 14 2016

lst={};Do[If[n==Plus@@Divisors[n]-n+6,AppendTo[lst,n]],{n,10^4}];Print[lst];
Select[Range[1, 10^8], DivisorSigma[1, #] - 2 # == - 6 &] (* Vincenzo Librandi, Sep 14 2016 *)

is(n)=sigma(n)==2*n-6 \\ Charles R Greathouse IV, Apr 23 2015, corrected by M. F. Hasler, Jul 18 2016

a(8) from Donovan Johnson, Dec 08 2011

A326057 a(n) = gcd(A003961(n)-2n, A003961(n)-sigma(n)), where A003961(n) is fully multiplicative function with a(prime(k)) = prime(k+1).

Original entry on oeis.org

1, 1, 1, 1, 1, 3, 3, 1, 1, 1, 1, 1, 3, 1, 1, 1, 1, 3, 3, 1, 1, 1, 1, 3, 1, 1, 1, 43, 1, 3, 5, 1, 1, 1, 1, 1, 3, 1, 1, 1, 1, 3, 3, 1, 1, 5, 1, 1, 1, 1, 1, 1, 1, 3, 19, 1, 1, 1, 1, 3, 5, 1, 1, 1, 1, 3, 3, 5, 7, 1, 1, 3, 1, 1, 1, 1, 1, 3, 3, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 3, 5, 1, 1, 1, 1, 3, 3, 1, 1, 1, 1, 3, 3, 1, 1

Offset: 1

Antti Karttunen, Jun 06 2019

nonn

Terms a(n) larger than 1 and equal to A252748(n) occur at n = 6, 28, 69, 91, 496, ..., see A326134. See also A349753.

Records 1, 3, 43, 45, 2005, 79243, ... occur at n = 1, 6, 28, 360, 496, 8128, ...

Cf. A000203, A000360, A003961, A033879, A033880, A252748, A286385, A326134.

Cf. also A009194, A325385, A325813, A325975, A326046, A326047, A326048, A326056, A326060, A326062, A349753.

Array[GCD[#3 - #1, #3 - #2] & @@ {2 #, DivisorSigma[1, #], Times @@ Map[#1^#2 & @@ # &, FactorInteger[#] /. {p_, e_} /; e > 0 :> {Prime[PrimePi@ p + 1], e}] - Boole[# == 1]} &, 78] (* 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
A252748(n) = (A003961(n) - (2*n));
A286385(n) = (A003961(n) - sigma(n));
A326057(n) = gcd(A252748(n), A286385(n));

a(n) = gcd(A252748(n), A286385(n)) = gcd(A003961(n) - 2n, A003961(n) - A000203(n)).

a(n) = gcd(A252748(n), A033879(n)) = gcd(A286385(n), A033879(n)). [Also A033880 can be used] - Antti Karttunen, May 06 2024

A191363 Numbers m such that sigma(m) = 2*m - 2.

Original entry on oeis.org

3, 10, 136, 32896, 2147516416

Offset: 1

Luis H. Gallardo, May 31 2011

Let k be a nonnegative integer such that F(k) = 2^(2^k) + 1 is prime (a Fermat prime A019434), then m = (F(k)-1)*F(k)/2 appears in the sequence.
Conjecture: a(1)=3 is the only odd term of the sequence.
Conjecture: All terms of the sequence are of the above form derived from Fermat primes.
The sequence has 5 (known) terms in common with sequences A055708 (k-1 | sigma(k)) and A056006 (k | sigma(k)+2) since {a(n)} is a subsequence of both.
The first five terms of the sequence are respectively congruent to 3, 4, 4, 4, 4 modulo 6.
After a(5) there are no further terms < 8*10^9.
Up to m = 1312*10^8 there are no further terms in the class congruent to 4 modulo 6.
a(6) > 10^12. - Donovan Johnson, Dec 08 2011
a(6) > 10^13. - Giovanni Resta, Mar 29 2013
a(6) > 10^18. - Hiroaki Yamanouchi, Aug 21 2018
See A125246 for numbers with deficiency 4, i.e., sigma(m) = 2*m - 4, and A141548 for numbers with deficiency 6. - M. F. Hasler, Jun 29 2016 and Jul 17 2016
A term m of this sequence multiplied by a prime p not dividing it is abundant if and only if p < m-1. For each of a(2..5) there is such a prime near this limit (here: 7, 127, 30197, 2147483647) such that a(k)*p is a primitive weird number, cf. A002975. - M. F. Hasler, Jul 19 2016
Any term m of this sequence can be combined with any term j of A088831 to satisfy the property (sigma(m) + sigma(j))/(m+j) = 2, which is a necessary (but not sufficient) condition for two numbers to be amicable. [Proof: If m = a(n) and j = A088831(k), then sigma(m) = 2m-2 and sigma(j) = 2j+2. Thus, sigma(m) + sigma(j) = (2m-2) + (2j+2) = 2m + 2j = 2(m+j), which implies that (sigma(m) + sigma(j))/(m+j) = 2(m+j)/(m+j) = 2.] - Timothy L. Tiffin, Sep 13 2016
At least the first five terms are a subsequence of A295296 and of A295298. - David A. Corneth, Antti Karttunen, Nov 26 2017
Conjectures: all terms are second hexagonal numbers (A014105). There are no terms with middle divisors. - Omar E. Pol, Oct 31 2018
The symmetric representation of sigma(m) of each of the 5 numbers in the sequence consists of 2 parts of width 1 that meet at the diagonal (subsequence of A246955). - Hartmut F. W. Hoft, Mar 04 2022
The first five terms coincide with the sum of two successive terms of A058891. The same is not true for a(6), if such exists. - Omar E. Pol, Mar 03 2023

			For n=1, a(1) = 3 since sigma(3) = 4 = 2*3 - 2.

Gianluca Amato, Maximilian Hasler, Giuseppe Melfi, and Maurizio Parton, Primitive weird numbers having more than three distinct prime factors, Riv. Mat. Univ. Parma, 7(1), (2016), 153-163, arXiv:1803.00324 [math.NT], 2018.

Cf. A000203, A002975, A056006, A055708, A088831 (abundance 2).
Cf. A033880, A125246 (deficiency 4), A141548 (deficiency 6), A125247 (deficiency 8), A125248 (deficiency 16).
Cf. A295296, A295298.
Cf. A014105, A237593, A246955.
Cf. A058891.

[n: n in [1..9*10^6] | (SumOfDivisors(n)-2*n) eq -2]; // Vincenzo Librandi, Sep 15 2016

ok[n_] := DivisorSigma[1,n] == 2*n-2; Select[ Table[ 2^(2^k-1) * (2^(2^k)+1), {k, 0, 5}], ok] (* Jean-François Alcover, Sep 14 2011, after conjecture *)
Select[Range[10^6], DivisorSigma[1, #] == 2 # - 2 &] (* Michael De Vlieger, Sep 14 2016 *)

zp(a,b) = {my(c,c1,s); c = a; c1 = 2*c-2;
while(c

a(k)=(2^2^k+1)<<(2^k-1) \\ For k<6. - M. F. Hasler, Jul 27 2016

a(n) = (A019434(n)-1)*A019434(n)/2 for all terms known so far. - M. F. Hasler, Jun 29 2016

A125246 Numbers m whose abundance sigma(m) - 2m = -4. Numbers whose deficiency is 4.

Original entry on oeis.org

5, 14, 44, 110, 152, 884, 2144, 8384, 18632, 116624, 8394752, 15370304, 73995392, 536920064, 2147581952, 34360131584, 27034175140420610, 36028797421617152, 576460753914036224

Offset: 1

Jason G. Wurtzel, Nov 25 2006

a(17) > 10^12. - Donovan Johnson, Dec 08 2011
a(17) > 10^13. - Giovanni Resta, Mar 29 2013
a(17) <= b(28) = 36028797421617152 ~ 3.6*10^16, since b(k) := 2^(k-1)*(2^k+3) is in this sequence for all k in A057732, i.e., whenever 2^k+3 is prime, and 28 = A057732(11). Further terms of this form are b(30), b(55), b(67), b(84), ... The only terms not of the form b(k), below 10^13, are {110, 884, 18632, 116624, 15370304, 73995392}. - M. F. Hasler, Apr 27 2015, edited on Jul 17 2016
See A191363 for numbers with deficiency 2, and A141548 for numbers with deficiency 6. - M. F. Hasler, Jun 29 2016 and Jul 17 2016
A term of this sequence multiplied with a prime p not dividing it is abundant if and only if p < sigma(a(n))/4. For each of a(2..16) there is such a prime, near this limit, such that a(n)*p is a primitive weird number, cf. A002975. - M. F. Hasler, Jul 17 2016
Any term x of this sequence can be combined with any term y of A088832 to satisfy the property (sigma(x)+sigma(y))/(x+y) = 2, which is a necessary (but not sufficient) condition for two numbers to be amicable. - Timothy L. Tiffin, Sep 13 2016
Is 5 the only odd number in this sequence? Is it possible to prove this? - M. F. Hasler, Feb 22 2017
a(20) > 10^18. - Hiroaki Yamanouchi, Aug 21 2018
If m is an even term, then (m-2)/2 is a term of A067680. - Jinyuan Wang, Apr 08 2020

			The abundance of 5 = (1+5)-10 = -4.
More generally, whenever p = 2^k + 3 is prime (as p = 5 for k = 1), then A(2^(k-1)*p) = (2^k-1)*(p+1) - 2^k*p = 2^k - p - 1 = -4.

Cf. A033880, A057732, A067680, A191363, A141548, A125247, A125248, A088832 (abundance 4).

[n: n in [1..9*10^6] | (SumOfDivisors(n)-2*n) eq -4]; // Vincenzo Librandi, Sep 15 2016

Select[Range[10^7], DivisorSigma[1, #] - 2 # == -4 &] (* Michael De Vlieger, Jul 18 2016 *)

for(n=1,1000000,if(((sigma(n)-2*n)==-4),print1(n,",")))

a(11) to a(14) from Klaus Brockhaus, Nov 29 2006
a(15)-a(16) from Donovan Johnson, Dec 23 2008
a(17)-a(19) from Hiroaki Yamanouchi, Aug 21 2018

A125247 Numbers n whose abundance sigma(n) - 2n = -8. Numbers n whose deficiency is 8.

Original entry on oeis.org

22, 130, 184, 1012, 2272, 18904, 33664, 70564, 85936, 100804, 391612, 527872, 1090912, 17619844, 2147713024, 6800695312, 34360655872, 549759483904, 1661355408388, 28502765343364, 82994670582016, 99249696661504, 120646991405056, 431202442356004, 952413274955776

Offset: 1

Jason G. Wurtzel, Nov 25 2006

a(19) > 10^12. - Donovan Johnson, Dec 08 2011
a(20) > 10^13. - Giovanni Resta, Mar 29 2013
a(30) > 10^18. - Hiroaki Yamanouchi, Aug 21 2018
a(20) <= 36028797958488064 ~ 3.6*10^16. Indeed, if k is in A057195 then 2^(k-1)*A168415(k) is in this sequence, and k=28 yields this upper bound for a(20) which is in any case a term of this sequence. - M. F. Hasler, Apr 27 2015
If n is in this sequence and p a prime not dividing n, then np is abundant if and only if p < sigma(n)/8 = n/4-1. For all n=a(k) except {22, 70564, 100804, 17619844}, there is such a p near this limit, such that n*p is a primitive weird number (A002975; in A258882 for the terms mentioned in the preceding comment). - M. F. Hasler, Jul 20 2016
Any term x of this sequence can be combined with any term y of A088833 to satisfy the property (sigma(x)+sigma(y))/(x+y) = 2, which is a necessary (but not sufficient) condition for two numbers to be amicable. - Timothy L. Tiffin, Sep 13 2016
Is there any odd number in this sequence? Is it possible to prove the contrary? - M. F. Hasler, Feb 22 2017

			The abundance of 22 = (1+2+11+22)-44 = -8

Hiroaki Yamanouchi, Table of n, a(n) for n = 1..29

Cf. A033880, A088833 (abundance 8).

[n: n in [1..2*10^7] | (DivisorSigma(1,n)-2*n) eq - 8]; // Vincenzo Librandi, Jul 22 2016

Select[Range[10^6], DivisorSigma[1, #] - 2 # == -8 &] (* Michael De Vlieger, Jul 21 2016 *)

for(n=1,1000000,if(((sigma(n)-2*n)==-8),print1(n,",")))

a(13)-a(15) from Klaus Brockhaus, Nov 29 2006
a(16)-a(17) from Donovan Johnson, Dec 23 2008
a(18) from Donovan Johnson, Dec 08 2011
a(19) from Giovanni Resta, Mar 29 2013
a(20)-a(25) from Hiroaki Yamanouchi, Aug 21 2018

cp's OEIS Frontend

A357976 Numbers with a divisor having the same sum of prime indices as their quotient.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

Mathematica

A181595 Abundant numbers n for which the abundance d = sigma(n) - 2*n is a proper divisor, that is, 0 < d < n and d | n.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

Mathematica

PARI

Extensions

A088833 Numbers n whose abundance is 8: sigma(n) - 2n = 8.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

PARI

Extensions

A357854 Squarefree numbers with a divisor having the same sum of prime indices as their quotient.

Views

Author

Keywords

Comments

Examples

Crossrefs

Programs

Mathematica

A077374 Odd numbers m whose abundance by absolute value is at most 10, that is, -10 <= sigma(m) - 2m <= 10.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

PARI

Extensions

A141548 Numbers n whose deficiency is 6.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Magma

Mathematica

PARI

Extensions

A326057 a(n) = gcd(A003961(n)-2n, A003961(n)-sigma(n)), where A003961(n) is fully multiplicative function with a(prime(k)) = prime(k+1).

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Mathematica