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

A379029 Modified exponential abundant numbers: numbers k such that A241405(k) > 2*k.

Original entry on oeis.org

30, 42, 66, 70, 78, 102, 114, 120, 138, 150, 168, 174, 186, 210, 222, 246, 258, 270, 282, 294, 318, 330, 354, 366, 390, 402, 420, 426, 438, 462, 474, 498, 510, 534, 546, 570, 582, 606, 618, 630, 642, 654, 660, 678, 690, 714, 726, 750, 762, 770, 780, 786, 798, 822
Offset: 1

Views

Author

Amiram Eldar, Dec 14 2024

Keywords

Comments

All the squarefree abundant numbers (A087248) are terms since A241405(k) = A000203(k) for a squarefree number k.
If k is a term and m is coprime to k them k*m is also a term.
The numbers of terms that do no exceed 10^k, for k = 2, 3, ..., are 5, 67, 767, 7595, 76581, 764321, 7644328, 76468851, 764630276, ... . Apparently, the asymptotic density of this sequence exists and equals 0.07646... .

Crossrefs

Subsequence of A005101.
Subsequences: A034683, A087248, A379030, A379031.
Similar sequences: A064597, A129575, A129656, A292982, A348274, A348604.

Programs

  • Mathematica
    f[p_, e_] := DivisorSum[e + 1, p^(# - 1) &]; mesigma[1] = 1; mesigma[n_] := Times @@ f @@@ FactorInteger[n]; meAbQ[n_] := mesigma[n] > 2*n; Select[Range[1000], meAbQ]
  • PARI
    is(n) = {my(f=factor(n)); prod(i=1, #f~, sumdiv(f[i, 2]+1, d, f[i, 1]^(d-1))) > 2*n;}

A380932 Odd numbers k such that A380845(k) > 2*k.

Original entry on oeis.org

322245, 590205, 874665, 966735, 1934415, 2900205, 3224025, 3378375, 3869775, 4729725, 6081075, 6449625, 6818175, 7740495, 8783775, 8906625, 9029475, 9889425, 10135125, 10961685, 11609325, 11821425, 12900825, 13378365, 14189175, 15049125, 15481935, 15909075, 16253055
Offset: 1

Views

Author

Amiram Eldar, Feb 08 2025

Keywords

Comments

The odd terms in A380929.
Analogous to odd abundant numbers (A005231) with A380845 instead of A000203.

Examples

			322245 is a term since it is odd, and A380845(322245) = 679582 > 2 * 322245 = 644490.
		

Crossrefs

Intersection of A005408 and A380929.
Subsequence of A005231.

Programs

  • Mathematica
    q[k_] := Module[{h = DigitCount[k, 2, 1]}, DivisorSum[k, # &, DigitCount[#, 2, 1] == h &] > 2*k]; Select[Range[1,10^6,2], q]
  • PARI
    isok(k) = if(!(k % 2), 0, my(h = hammingweight(k)); sumdiv(k, d, d*(hammingweight(d) == h)) > 2*k);
Showing 1-2 of 2 results.