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

A096399 Numbers k such that both k and k+1 are abundant.

Original entry on oeis.org

5775, 5984, 7424, 11024, 21735, 21944, 26144, 27404, 39375, 43064, 49664, 56924, 58695, 61424, 69615, 70784, 76544, 77175, 79695, 81080, 81675, 82004, 84524, 84644, 89775, 91664, 98175, 103455, 104895, 106784, 109395, 111824, 116655, 116864, 120015, 121904, 122264
Offset: 1

Views

Author

John L. Drost, Aug 06 2004

Keywords

Comments

Numbers k such that both sigma(k) > 2k and sigma(k+1) > 2*(k+1).
Numbers k such that both k and k+1 are in A005101.
Set difference of sequences A103289 and {2^m-1} for m in A103291.
The numbers of terms not exceeding 10^k, for k = 4, 5, ..., are 3, 27, 357, 3723, 36640, 365421, 3665799, 36646071, ... . Apparently, the asymptotic density of this sequence exists and equals 0.000366... . - Amiram Eldar, Sep 02 2022

Examples

			sigma(5775) = sigma(3*5*5*7*11) = 11904 > 2*5775.
sigma(5776) = sigma(2*2*2*2*19*19) = 11811 > 2*5776.
		

Crossrefs

Programs

  • Mathematica
    fQ[n_] := DivisorSigma[1, n] > 2 n; Select[ Range@ 117000, fQ[ # ] && fQ[ # + 1] &] (* Robert G. Wilson v, Jun 11 2010 *)
    Select[Partition[Select[Range[120000], DivisorSigma[1, #] > 2 # &], 2, 1], Differences@ # == {1} &][[All, 1]] (* Michael De Vlieger, May 20 2017 *)
  • PARI
    for(i=1,1000000,if(sigma(i)>2*i && sigma(i+1)>2*(i+1),print(i))); \\ Max Alekseyev, Jan 28 2005

Formula

a(n) = A005101(A169822(n)). - Amiram Eldar, Mar 01 2025

Extensions

Two further terms from Max Alekseyev, Jan 28 2005
Entry revised by N. J. A. Sloane, Dec 03 2006
Edited by T. D. Noe, Nov 15 2010

A103291 Numbers k such that sigma(2^k-1) >= 2*(2^k-1)-1, i.e., the number 2^k-1 is perfect, abundant, or least deficient.

Original entry on oeis.org

1, 12, 24, 36, 40, 48, 60, 72, 80, 84, 90, 96, 108, 120, 132, 140, 144, 156, 160, 168, 180, 192, 200, 204, 210, 216, 220, 228, 240, 252, 264, 270, 276, 280, 288, 300, 312, 320, 324, 330, 336, 348, 360, 372, 384, 396, 400, 408, 420, 432, 440, 444, 450, 456, 468
Offset: 1

Views

Author

Max Alekseyev, Jan 28 2005

Keywords

Comments

Is there an odd term besides 1? Numbers 2^a(i)-1 form set difference of sequences A103289 and A096399.
Odd terms > 1 exist, but there are none < 10^7. If k > 1 is an odd term, then 2^k-1 must have more than 900000 distinct prime factors and all of them must be members of A014663. - David Wasserman, Apr 15 2008

Crossrefs

Programs

  • PARI
    for(i=1,1000,n=2^i-1;if(sigma(n)>=2*n-1,print1(i, ", ")));

Formula

Numbers k such that 2^k-1 is in A103288.

Extensions

More terms from David Wasserman, Apr 15 2008

A103292 Numbers k such that sigma(2^k - 1) >= 2(2^k - 1), i.e., the number 2^k - 1 is perfect or abundant.

Original entry on oeis.org

12, 24, 36, 40, 48, 60, 72, 80, 84, 90, 96, 108, 120, 132, 140, 144, 156, 160, 168, 180, 192, 200, 204, 210, 216, 220, 228, 240, 252, 264, 270, 276, 280, 288, 300, 312, 320, 324, 330, 336, 348, 360, 372, 384, 396, 400, 408, 420, 432, 440, 444, 450, 456, 468, 480, 492, 504
Offset: 1

Views

Author

Max Alekseyev, Jan 28 2005

Keywords

Comments

Numbers k that 2^k - 1 is in A023196.
Are there any odd terms? This is a subsequence of A103291. Is the number 1 the only term where they differ? This is so if there is no least deficient number of the form 2^n-1 besides 1.
For each n in the sequence, 2n is also in the sequence: sigma[2^(2n)-1] = sigma[(2^n+1)(2^n-1)] >= (2^n+1)*sigma(2^n-1) because for each divisor d|2^n-1 there is (at least) the divisor (2^n+1)d |[(2^n+1)(2^n-1)]. Inserting sigma(2^n-1) >=2(2^n-1) yields (2^n+1)*sigma(2^n-1)>=(2^n+1)*2*(2^n-1)=2*[2^(2n)-1] qed. - R. J. Mathar, Aug 07 2007
From David Wasserman, May 16 2008: (Start)
Odd members exist. One such n is the lcm of the first 4416726 members of A139686, which has 6864499 digits. To show that n is a member, it's not necessary to exactly compute sigma(2^n-1).
The function f(x) = sigma(x)/x is multiplicative and has the property that for any a, b > 1, f(ab) > f(a). So it suffices to find some y such that f(y) >= 2 and y divides 2^n-1. In this case, y is the product of the first 4416726 members of A014663 and has 35260810 digits. (A014663(4416726) = 278379727.)
To see that this works, note that if a divides b, then 2^a-1 divides 2^b-1. For 1 <= i <= 4416726, A014663(i) divides 2^A139686(i)-1 by definition and A139686(i) divides n, so 2^A139686(i)-1 divides 2^n-1 and therefore A014663(i) divides 2^n-1. Then we can compute that f(y) = Product_{i = 1..4416726} (1 + 1/A014663(i)) is > 2.
The members of A014663 are the only primes that can divide 2^n-1 with n odd. Any powers of these primes are also possible divisors.
By including powers, we can construct a much smaller y. I found a y with 7057382 digits, omega(y) = 969004 and bigomega(y) = 969440. This y is close to the minimum possible. The least n such that y divides 2^n-1 is an odd number with 1472897 digits.
However, minimizing y is not the way to minimize n. We can get a smaller n by skipping primes p such that the order of 2 mod p is divisible by a large prime. This increases the number and size of the prime factors needed to make f(y) >= 2 and the time needed to find them.
The least odd n that I've found has 28375 digits. The corresponding y has 305621222 digits, omega(y) = 31903142 and bigomega(y) = 32796897. To find these prime factors, I searched up to A014663(96433108) = 7154804519.
I believe that the smallest odd member has between 10000 and 20000 digits, but the largest lower bound I can prove has 8 digits: f(p^i) is bounded above by 1 + 1/(p-1) and Product_{i=1..c} (1 + 1/(A014663(i)-1)) < 2 if c < 968858, so y must be at least Product_{i=1..968858} A014663(i), which has 7054790 digits.
Then n must be large enough that 2^n-1 >= y, yielding a lower bound of 23435503. I don't see any way to increase this significantly. (End)

Crossrefs

Programs

  • PARI
    for(i=1,1000,n=2^i-1;if(sigma(n)>=2*n,print(i)));

Extensions

Extended to a(32) by R. J. Mathar, Aug 07 2007
Terms from a(33) onwards from David Wasserman, May 16 2008

A387654 Numbers k such that each of k and k+1 is either a practical number (A005153) or an almost practical number (A174533).

Original entry on oeis.org

1, 4095, 5775, 5984, 11024, 21735, 21944, 26144, 39375, 49664, 58695, 61424, 69615, 76544, 79695, 89775, 91664, 98175, 104895, 106784, 111824, 116655, 116864, 121904, 134504, 135135, 144584, 152775, 155295, 160544, 165375, 170624, 174824, 180495, 185535, 192464
Offset: 1

Views

Author

Amiram Eldar, Sep 05 2025

Keywords

Comments

Differs from A103289 by not having the terms 7424, 27404, 43064, 56924, 70784, ... . The first 344 terms of this sequence are in A103289. Is this sequence a subsequence of A103289?
Differs from A096399 by not having the terms 7424, 27404, 43064, 56924, 70784, ... . The first 342 terms after 1 and 4095 are in A096399. Is this sequence \ {1, 4095} a subsequence of A096399?
Terms k such that both k and k+1 are almost practical numbers are in A387653.
The only pair of consecutive integers that are both practical is 1 and 2, since 1 is the only odd practical number.
All the rest are pairs in which one member (the odd member) is almost practical and the second member (the even member) is practical.
Are there 3 consecutive numbers that are all either practical or almost practical? There are none below 2.8*10^6.

Crossrefs

A387653 is subsequence.

Programs

  • Mathematica
    q[n_] := q[n] = Module[{d = Divisors[n], c, x}, c = CoefficientList[Product[1 + x^i, {i, d}], x]; MemberQ[{0, 2}, Total[d] - Count[Rest[c], _?(# > 0 &)]]];
    Select[Range[40000], q[#] && q[#+1] &] (* warning: a slow program *)
  • PARI
    isp(k) = {my(d = divisors(k), nd = #d, s = vecsum(d), p = prod(i = 1, nd, 1 + 'x^d[i])); abs(#select(x -> x > 0, Col(p)) - s) == 1;}
    list(kmax) = {my(is1 = 1, is2); for(k = 2, kmax, is2 = isp(k); if(is1 && is2, print1(k-1, ", ")); is1 = is2);} \\ warning: a slow program
Showing 1-4 of 4 results.