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

A308053 Coreful abundant numbers: numbers k such that csigma(k) > 2*k, where csigma(k) is the sum of the coreful divisors of k (A057723).

Original entry on oeis.org

72, 108, 144, 200, 216, 288, 324, 360, 400, 432, 504, 540, 576, 600, 648, 720, 756, 784, 792, 800, 864, 900, 936, 972, 1000, 1008, 1080, 1152, 1188, 1200, 1224, 1296, 1368, 1400, 1404, 1440, 1512, 1568, 1584, 1600, 1620, 1656, 1728, 1764, 1800, 1836, 1872, 1936
Offset: 1

Views

Author

Amiram Eldar, May 10 2019

Keywords

Comments

Analogous to A005101 as A307958 is analogous to A000396.
The asymptotic density of this sequence is Sum_{n>=1} f(A356871(n)) = 0.0262215..., where f(n) = (6/(Pi^2*n))*Product_{prime p|n}(p/(p+1)). - Amiram Eldar, Sep 02 2022

Examples

			72 is in the sequence since its coreful divisors are 6, 12, 18, 24, 36, 72, whose sum is 168 > 2 * 72.

Links

Crossrefs

Cf. A005101, A057723, A307958.
A339940 and A356871 are subsequences.
Subsequence of A129575.

Programs

  • Mathematica
    f[p_, e_] := (p^(e+1)-1)/(p-1)-1; a[1]=1; a[n_] := Times @@ (f @@@ FactorInteger[n]); s={}; Do[If[a[n] > 2n, AppendTo[s, n]], {n, 1, 2000}]; s
  • PARI
    rad(n) = factorback(factorint(n)[, 1]); \\ A007947
s(n) = rad(n)*sigma(n/rad(n)); \\ A057723
isok(k) = s(k) > 2*k; \\ Michel Marcus, May 11 2019
  • PARI
    isok(k) = {my(f=factor(k)); prod(i = 1, #f~, (f[i,1]^(f[i,2]+1)-1) / (f[i,1]-1)-1) > 2*k}; \\ Amiram Eldar, Sep 02 2022

A356871 Primitive coreful abundant numbers (second definition): coreful abundant numbers (A308053) that are powerful numbers (A001694).

Original entry on oeis.org

72, 108, 144, 200, 216, 288, 324, 400, 432, 576, 648, 784, 800, 864, 900, 972, 1000, 1152, 1296, 1568, 1600, 1728, 1764, 1800, 1936, 1944, 2000, 2304, 2592, 2700, 2704, 2744, 2916, 3136, 3200, 3456, 3528, 3600, 3872, 3888, 4000, 4356, 4500, 4608, 4900, 5000, 5184
Offset: 1

Views

Author

Amiram Eldar, Sep 02 2022

Keywords

Comments

For squarefree numbers k, csigma(k) = k, where csigma(k) is the sum of the coreful divisors of k (A057723). Thus, if m is a term (csigma(m) > 2*m) and k is a squarefree number coprime to k, then csigma(k*m) = csigma(k) * csigma(m) = k * csigma(m) > 2*k*m, so k*m is a coreful abundant number. Therefore, the sequence of coreful abundant numbers (A308053) can be generated from this sequence by multiplying with coprime squarefree numbers. The asymptotic density of the coreful abundant numbers can be calculated from this sequence (see comment in A308053).

Examples

			72 is a term since csigma(72) = 168 > 2 * 72, and 72 = 2^3 * 3^2 is powerful.

Links

Crossrefs

Intersection of A001694 and A308053.
A339940 is a subsequence.
Cf. A057723.
Similar sequences: A307959, A328136.

Programs

  • Mathematica
    f[p_, e_] := (p^(e+1)-1)/(p-1)-1; s[1] = 1; s[n_] := If[AllTrue[(fct = FactorInteger[n])[[;;, 2]], #>1 &], Times @@ f @@@ fct, 0]; seq={}; Do[If[s[n] > 2*n, AppendTo[seq, n]], {n, 1, 5000}]; seq

A372299 Primitive infinitary abundant numbers (definition 2): infinitary abundant numbers (A129656) having no proper infinitary divisors that are infinitary abundant numbers.

Original entry on oeis.org

24, 30, 40, 42, 54, 56, 66, 70, 72, 78, 88, 96, 102, 104, 114, 138, 150, 174, 186, 222, 246, 258, 282, 294, 318, 354, 366, 402, 420, 426, 438, 474, 486, 498, 534, 540, 582, 606, 618, 642, 654, 660, 678, 726, 756, 762, 780, 786, 822, 834, 894, 906, 924, 942, 945, 960, 978, 990
Offset: 1

Views

Author

Amiram Eldar, Apr 25 2024

Keywords

Examples

			24 is a term since it is an infinitary abundant number and none of its proper infinitary divisors, {1, 2, 3, 4, 6, 8, 12}, are infinitary abundant numbers.
The least infinitary abundant number that is not primitive is 120. It has 3 infinitary divisors, 24, 30, and 40, that are also infinitary abundant numbers.

Links

Crossrefs

Subsequence of A129656.
A372298 is a subsequence.
Cf. A077609, A129657.
Similar sequences: A091191, A302574, A339940.

Programs

  • Mathematica
    f[p_, e_] := Module[{b = IntegerDigits[e, 2]}, m = Length[b]; Product[If[b[[j]] > 0, 1 + p^(2^(m - j)), 1], {j, 1, m}]];
isigma[1] = 1; isigma[n_] := Times @@ f @@@ FactorInteger[n]; iabQ[n_] := isigma[n] > 2*n; idivs[1] = {1};
idivs[n_] := Sort@ Flatten@ Outer[Times, Sequence @@ (FactorInteger[n] /. {p_, e_Integer} :> p^Select[Range[0, e], BitOr[e, #] == e &])];
q[n_] := Module[{d = idivs[n]}, Total[d] > 2*n && AllTrue[Most[d], !iabQ[#] &]]; Select[Range[1000], q]
  • PARI
    isidiv(d, f) = {if (d==1, return (1)); for (k=1, #f~, bne = binary(f[k, 2]); bde = binary(valuation(d, f[k, 1])); if (#bde < #bne, bde = concat(vector(#bne-#bde), bde)); for (j=1, #bne, if (! bne[j] && bde[j], return (0)); ); ); return (1); }
idivs(n) = {my(f = factor(n), d = divisors(f), idiv = []); for (k=1, #d, if (isidiv(d[k], f), idiv = concat(idiv, d[k])); ); idiv; } \\ Michel Marcus at A077609
isigma(n) = {my(f = factor(n), b); prod(i=1, #f~, b = binary(f[i, 2]); prod(k=1, #b, if(b[k], 1+f[i, 1]^(2^(#b-k)), 1)))} ;
is(n) = isigma(n) > 2*n && select(x -> x < n && isigma(x) > 2*x, idivs(n)) == [];
Showing 1-3 of 3 results.

Started in 1964 by Neil J. A. Sloane | Maintained by The OEIS Foundation Inc.

Content is available under The OEIS End-User License Agreement.