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

A063734 Square abundant numbers.

Original entry on oeis.org

36, 100, 144, 196, 324, 400, 576, 784, 900, 1296, 1600, 1764, 1936, 2304, 2500, 2704, 2916, 3136, 3600, 4356, 4624, 4900, 5184, 5776, 6084, 6400, 7056, 7744, 8100, 8464, 9216, 9604, 10000, 10404, 10816, 11025, 11664, 12100, 12544, 12996, 13456, 14400
Offset: 1

Views

Author

Jason Earls, Aug 13 2001

Keywords

Crossrefs

Intersection of A000290 and A005101.
Subsequence of A363169.
Cf. A381738.

Programs

  • Mathematica
    Select[Range[150]^2,DivisorSigma[1,#]>2#&] (* Harvey P. Dale, Aug 14 2014 *)
  • PARI
    j=[]; for(n=1,25000, if(sigma(n)>(2*n),a=n; if(issquare(a),j=concat(j,a)))); j
    
  • PARI
    { n=0; for (m=1, 10^9, if (issquare(m) && sigma(m)>(2*m), write("b063734.txt", n++, " ", m); if (n==1000, break)) ) } \\ Harry J. Smith, Aug 29 2009
    
  • PARI
    { n=0; for (m=1, 10^9, s=m^2; if (sigma(s)>(2*s), write("b063734.txt", n++, " ", s); if (n==1000, break)) ) } \\ Harry J. Smith, Aug 29 2009

Formula

a(n) = A381738(n)^2. - Amiram Eldar, Mar 07 2025

A381741 Squarefree numbers k such that k^2 is abundant, and d^2 is nonabundant for any proper divisor d of k.

Original entry on oeis.org

6, 10, 14, 105, 286, 374, 418, 442, 506, 2145, 2805, 3135, 3315, 3705, 3795, 4485, 4785, 4845, 5115, 5655, 6045, 6105, 6765, 7095, 7755, 8745, 9735, 10065, 11362, 14326, 14858, 15314, 17342, 18278, 18538, 18734, 19778, 20026, 20254, 21242, 22126, 22678, 23218
Offset: 1

Views

Author

Amiram Eldar, Mar 06 2025

Keywords

Comments

The primitive terms of A381740. Each term of A381740 is a multiple of a term in this sequence.

Crossrefs

Intersection of A005117 and A381739.
Subsequence of A381738 and A381740.

Programs

  • Mathematica
    q[k_] := SquareFreeQ[k] && DivisorSigma[-1, k^2] > 2 && AllTrue[Divisors[k], DivisorSigma[-1, #^2] <= 2 || # == k &]; Select[Range[24000], q]
  • PARI
    is1(k) = {my(f = factor(k)); if(!issquarefree(f), 0, prod(i = 1, #f~, f[i,2] *= 2); sigma(f, -1) > 2);}
    isok(k) = if(!is1(k), 0, fordiv(k, d, if(d < k && is1(d), return(0))); 1);

A381739 Number k such that k^2 is abundant, and d^2 is nonabundant for any proper divisor d of k.

Original entry on oeis.org

6, 10, 14, 44, 52, 68, 76, 92, 105, 116, 124, 286, 296, 328, 344, 374, 376, 418, 424, 442, 472, 488, 495, 506, 536, 568, 584, 585, 632, 664, 712, 776, 808, 824, 856, 872, 904, 1016, 2096, 2145, 2192, 2224, 2384, 2416, 2512, 2608, 2672, 2768, 2805, 2864, 2896, 3056
Offset: 1

Views

Author

Amiram Eldar, Mar 05 2025

Keywords

Comments

The primitive terms of A381738. Each term of A381738 is a multiple of a term in this sequence.

Crossrefs

Subsequence of A381738.
A381741 is a subsequence.

Programs

  • Mathematica
    q[k_] := DivisorSigma[-1, k^2] > 2 && AllTrue[Divisors[k], DivisorSigma[-1, #^2] <= 2 || # == k &]; Select[Range[3200], q]
  • PARI
    is1(k) = {my(f = factor(k)); prod(i = 1, #f~, f[i,2] *= 2); sigma(f, -1) > 2;}
    isok(k) = if(!is1(k), 0, fordiv(k, d, if(d < k && is1(d), return(0))); 1);

A381740 Squarefree numbers k such that k^2 is abundant.

Original entry on oeis.org

6, 10, 14, 30, 42, 66, 70, 78, 102, 105, 110, 114, 130, 138, 154, 170, 174, 182, 186, 190, 210, 222, 230, 238, 246, 258, 266, 282, 286, 290, 310, 318, 322, 330, 354, 366, 370, 374, 390, 402, 406, 410, 418, 426, 430, 434, 438, 442, 462, 470, 474, 498, 506, 510
Offset: 1

Views

Author

Amiram Eldar, Mar 05 2025

Keywords

Comments

The least odd term is a(105) = 105.
The least term that is coprime to 6 is a(3637276) = 37182145.
If k is a squarefree number that is divisible by 6, 10 or 14, then it is a term. Therefore a lower bound for the asymptotic density of this sequence is 29/(192*zeta(2)) = 0.0918... .
The numbers of terms that do not exceed 10^k, for k = 1, 2, ..., are 2, 8, 99, 972, 9826, 97610, 979190, 9770801, 97650638, 976893969, ... . Apparently, the asymptotic density of this sequence exists and equals 0.0976... .
If k is a term then any multiple of k that is squarefree is a term. The primitive terms are in A381741.

Crossrefs

Intersection of A005117 and A381738.
Subsequences: A087248, A381741.
Cf. A013661.

Programs

  • Mathematica
    Select[Range[600], SquareFreeQ[#] && DivisorSigma[-1, #^2] > 2 &]
  • PARI
    isok(k) = {my(f = factor(k)); if(!issquarefree(f), 0, prod(i = 1, #f~, f[i,2] *= 2); sigma(f, -1) > 2);}

A381742 Numbers k such that k^2 is abundant but d*k is nonabundant for any proper divisor d of k.

Original entry on oeis.org

14, 124, 585, 1016, 16748, 32085, 33892, 37882, 39962, 41925, 46665, 121605, 134589, 181305, 212175, 388455, 495465, 522488, 524224, 544065, 839865, 1061565, 1152921, 1165515, 1243275, 1247103, 1335411, 1676829, 1943638, 2151075, 2290869, 2478075, 2625514, 2673998
Offset: 1

Views

Author

Amiram Eldar, Mar 06 2025

Keywords

Comments

Numbers k such that k^2 is primitive abundant number (A091191).
If p is an odd Mersenne exponent (A174265), then 2^((p-1)/2) * (2^p-1) is a term.

Crossrefs

Subsequence of A381738.
A379950 is a subsequence.

Programs

  • Mathematica
    q[k_] := DivisorSigma[-1, k^2] > 2 &&  AllTrue[Divisors[k], DivisorSigma[-1, #*k] <= 2 || # == k &]; Select[Range[10^6], q]
  • PARI
    isok(k) = fordiv(k, d, if(d < k && sigma(d*k, -1) > 2, return(0))); sigma(k^2, -1) > 2;

A383390 Numbers k such that k^2 and (k+1)^2 are both abundant numbers.

Original entry on oeis.org

104, 495, 584, 735, 944, 1155, 1364, 1484, 2144, 2204, 2415, 2624, 2924, 2925, 3135, 3255, 3794, 3795, 4304, 4484, 4784, 4844, 5264, 5355, 5445, 5564, 5565, 5655, 5775, 5984, 6104, 6764, 7424, 7455, 7664, 7755, 7875, 8084, 8294, 8295, 8414, 8415, 8924, 9009, 9344, 9944, 9975
Offset: 1

Views

Author

Amiram Eldar, Apr 25 2025

Keywords

Comments

The numbers of terms that do not exceed 10^k, for k = 3, 4, ..., are 5, 47, 459, 4655, 46733, 460693, 4612685, 46177602, ... . Apparently, the asymptotic density of this sequence exists and equals 0.00461... .

Crossrefs

Subsequence of A381738.
A383391 and A096399 are subsequences.

Programs

  • Mathematica
    Select[Range[10000], DivisorSigma[-1, #^2] > 2 && DivisorSigma[-1, (#+1)^2] > 2 &]
  • PARI
    is1(k) = {my(f = factor(k)); prod(i = 1, #f~, f[i,2] *= 2); sigma(f, -1) > 2;}
    list(lim) = {my(q1 = is1(1), q2); for(k = 2, lim, q2 = is1(k); if(q1 && q2, print1(k-1, ", ")); q1 = q2);}

A383391 Numbers k such that k^2, (k+1)^2 and (k+2)^2 are all abundant numbers.

Original entry on oeis.org

2924, 3794, 5564, 8294, 8414, 10064, 13454, 19304, 22154, 22814, 35684, 39974, 40544, 40754, 41768, 46214, 49994, 52064, 56264, 60884, 63854, 65624, 68354, 68474, 69068, 70244, 78974, 84824, 88604, 92168, 93224, 95354, 100694, 102464, 106028, 107084, 111110, 111824
Offset: 1

Views

Author

Amiram Eldar, Apr 25 2025

Keywords

Comments

Are there numbers k such that k^2, (k+1)^2, (k+2)^2 and (k+3)^2 are all abundant numbers? There are none below 2.5*10^10.
Are there odd terms in this sequence? There are none below 2.5*10^10.
From David A. Corneth, Apr 26 2025: (Start)
If it exists then it is at least sqrt(A002110(24)/2 * 155925 + 1) - 1 ~= 4.3*10^19.
Proof: Exactly two of k, k+1, k+2 and k+3 are odd. Those two are coprime and differ by 2. Let them be m and m+2. Then sigma(m) > 2*m and sigma(m+2) > 2*(m+2). As they are coprime we have sigma(m*(m+2)) > 2*m*2*(m+2) so for a lower bound we look for the smallest odd t that sigma(t) > 4*t. The partial product of p / (p-1) for odd primes p first exceeds 4 when odd primes <= 79 are multiplied so t is divisible by 3 * 5 * 7 * ... * 79. A small search of multiples of this number gives A002110(24)/2 * 155925.
k * (k + 2) >= A002110(24)/2 * 155925 so k * (k + 2) + 1 = (k + 1)^2 >= A002110(24)/2 * 155925 + 1. Taking square roots on both sides and keeping the positive root gives the desired lower bound. (End)
From Yifan Xie, Apr 30 2025: (Start)
Both types of numbers exist, but the constructed ones are too large to be displayed here. For numbers k such that k^2, (k+1)^2, (k+2)^2 and (k+3)^2 are all abundant numbers, choose 4 disjoint subsets of the primes P_1, P_2, P_3 and P_4, and let the product of elements in P_i divide k+i-1. This is achievable because of the Chinese remainder theorem. If P_i contains p_1, ..., p_k, then sigma((k+i-1)^2)/(k+i-1)^2 >= Product_{i=1..k} (p_i+1)/p_i.
We are able to make the right hand side larger than 2 for each i because the infinite product Product_{p is prime} (p+1)/p = Product_{p is prime} (1+1/p) = Sum_{k is squarefree} 1/k diverges, since the squarefree numbers have asymptotic density 6/Pi^2.
For odd terms in this sequence, we can use only the odd primes to construct 3 prime subsets instead, and add a constraint that k == 1 (mod 2) after which the Chinese remainder theorem still applies. (End)

Crossrefs

Subsequence of A381738 and A383390.

Programs

  • Mathematica
    q[n_] := q[n] = DivisorSigma[-1, n^2] > 2; Select[Range[120000], AllTrue[# + {0, 1, 2}, q] &]
  • PARI
    is1(k) = {my(f = factor(k)); prod(i = 1, #f~, f[i,2] *= 2); sigma(f, -1) > 2;}
    list(lim) = {my(q1 = is1(1), q2 = is1(2), q3); for(k = 3, lim, q3 = is1(k); if(q1 && q2 && q3, print1(k-2, ", ")); q1 = q2; q2 = q3);}

A381743 The number of divisors d of n such that d*n is abundant.

Original entry on oeis.org

0, 0, 0, 0, 0, 3, 0, 0, 0, 2, 0, 6, 0, 1, 0, 0, 0, 6, 0, 6, 0, 0, 0, 8, 0, 0, 0, 5, 0, 8, 0, 0, 0, 0, 0, 9, 0, 0, 0, 8, 0, 8, 0, 4, 0, 0, 0, 10, 0, 3, 0, 4, 0, 8, 0, 8, 0, 0, 0, 12, 0, 0, 0, 0, 0, 8, 0, 2, 0, 8, 0, 12, 0, 0, 0, 2, 0, 8, 0, 10, 0, 0, 0, 12, 0, 0
Offset: 1

Views

Author

Amiram Eldar, Mar 06 2025

Keywords

Comments

The indices of records of this sequence are the highly composite numbers (A002182) without 2 and 4.

Crossrefs

Programs

  • Mathematica
    a[n_] := DivisorSum[n, 1 &, DivisorSigma[-1, #*n] > 2 &]; Array[a, 100]
  • PARI
    a(n) = sumdiv(n, d, sigma(d*n, -1) > 2);

Formula

a(n) = A000005(n) if and only if n is abundant (A005101).
a(n) = A000005(n) - 1 if and only if n is either 1, a perfect number (A000396), or in A341358.
a(n) >= 1 if and only if n is in A381738.
a(n) = 1 if and only if n is in A381742.
Showing 1-8 of 8 results.