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.

A383693 Exponential unitary abundant numbers: numbers k such that A322857(k) > 2*k.

Original entry on oeis.org

900, 1764, 4356, 4500, 4900, 6084, 6300, 8820, 9900, 10404, 11700, 12348, 12996, 14700, 15300, 17100, 19044, 19404, 20700, 21780, 22500, 22932, 26100, 27900, 29988, 30276, 30420, 30492, 31500, 33300, 33516, 34596, 36900, 38700, 40572, 42300, 42588, 44100, 47700, 47916, 49284, 49500
Offset: 1

Views

Author

Amiram Eldar, May 05 2025

Keywords

Comments

First differs from its subsequence A383697 at n = 21.
All the terms are nonsquarefree numbers (A013929), since A322857(k) = k if k is a squarefree number (A005117).
If an exponential abundant number (A129575) is exponentially squarefree (A209061), then it is in this sequence. Terms of this sequence that are not exponentially squarefree are a(21) = 22500, a(77) = 86436, a(140) = 157500, etc..
The least odd term is a(202273) = 225450225, and the least term that is coprime to 6 is a(1.002..*10^18) = 1117347505588495206025.
The asymptotic density of this sequence is Sum_{n>=1} f(A383694(n)) = 0.00089722..., where f(n) = (6/(Pi^2*n))*Product_{prime p|n}(p/(p+1)).

Examples

			900 is a term since A322857(900) = 2160 > 2*900 = 1800.
		

Crossrefs

Subsequence of A013929 and A129575.
Subsequences: A383694, A383697, A383698.

Programs

  • Mathematica
    f[p_, e_] := DivisorSum[e, p^# &, GCD[#, e/#] == 1 &]; q[n_] := Times @@ f @@@ FactorInteger[n] > 2 n; Select[Range[50000], q]
  • PARI
    fun(p, e) = sumdiv(e, d, if(gcd(d, e/d) == 1, p^d));
    isok(k) = {my(f = factor(k)); prod(i = 1, #f~, fun(f[i, 1], f[i, 2])) > 2*k;}

A383698 Primitive exponential squarefree exponential abundant numbers: the powerful terms of A383697.

Original entry on oeis.org

900, 1764, 4356, 4500, 4900, 6084, 10404, 12348, 12996, 19044, 30276, 34596, 44100, 47916, 49284, 60516, 66564, 79092, 79524, 88200, 101124, 108900, 112500, 125316, 132300, 133956, 152100, 161604, 176868, 181476, 191844, 213444, 217800, 220500, 224676, 246924
Offset: 1

Views

Author

Amiram Eldar, May 06 2025

Keywords

Comments

Subsequence of A383694 and first differs from it at n = 11.
The least odd term is a(1345) = 225450225, and the least term that is coprime to 6 is 1117347505588495206025.
For squarefree numbers k, essigma(k) = k, where essigma is the sum of exponential squarefree exponential divisors function (A361174). Thus, if m is a term (essigma(m) > 2*m) and k is a squarefree number coprime to m, then essigma(k*m) = essigma(k) * essigma(m) = k * essigma(m) > 2*k*m, so k*m is an exponential squarefree exponential abundant number. Therefore, the sequence of exponential squarefree exponential abundant numbers (A383697) can be generated from this sequence by multiplying with coprime squarefree numbers.

Crossrefs

Intersection of A001694 and A383697.
Subsequence of A383694.

Programs

  • Mathematica
    fun[p_, e_] := DivisorSum[e, p^# &, SquareFreeQ[#] &]; q[n_] := Min[(f = FactorInteger[n])[[;; , 2]]] > 1 && Times @@ fun @@@ f > 2*n; Select[Range[250000], q]
  • PARI
    fun(p, e) = sumdiv(e, d, if(issquarefree(d), p^d, 0));
    isok(k) = {my(f = factor(k)); ispowerful(f) && prod(i = 1, #f~, fun(f[i, 1], f[i, 2])) > 2*k;}
Showing 1-2 of 2 results.