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.

A328135 Exponential 3-abundant numbers: numbers m such that esigma(m) >= 3m, where esigma(m) is the sum of exponential divisors of m (A051377).

Original entry on oeis.org

901800900, 1542132900, 1926332100, 2153888100, 2690496900, 2822796900, 3942584100, 4487660100, 4600908900, 5127992100, 6267888900, 6742052100, 7162236900, 7305120900, 8421732900, 8969984100, 9866448900, 10203020100, 10718460900, 11723411700, 11787444900, 12528324900
Offset: 1

Views

Author

Amiram Eldar, Oct 04 2019

Keywords

Comments

Aiello et al. found bounds on e-multiperfect numbers, i.e., numbers m such that esigma(m) = k * m for k > 2: 2 * 10^7 for k = 3, and 10^85, 10^320, and 10^1210 for k = 4, 5, and 6. The data of this sequence raise the bound for exponential 3-perfect numbers to 3 * 10^10.
The least odd term is (59#/2)^2 = 924251841031287598942273821762233522616225. The least term which is coprime to 6 is (239#/6)^2 = 3.135... * 10^190.
The least exponential 4-abundant number (esigma(m) >= 4m) is (31#)^2 = 40224510201185827416900. In general, the least exponential k-abundant number (esigma(m) >= k*m), for k > 2, is (A002110(A072986(k)))^2.
The asymptotic density of this sequence is Sum_{n>=1} f(A383699(n)) = 1.325...*10^(-9), where f(n) = (6/(Pi^2*n))*Product_{prime p|n}(p/(p+1)). - Amiram Eldar, May 06 2025

Crossrefs

Subsequence of A129575.
A383699 is a subsequence.
Cf. A023197, A307112, A285615 (unitary), A293187 (bi-unitary), A300664 (infinitary).

Programs

  • Mathematica
    f[p_, e_] := DivisorSum[e, p^# &]; esigma[1] = 1; esigma[n_] := Times @@ f @@@ FactorInteger[n]; Select[Range[10^10], esigma[#] >= 3 # &]

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

Original entry on oeis.org

5400, 7200, 10800, 14400, 16200, 18000, 21168, 21600, 27000, 28800, 32400, 36000, 37800, 42336, 43200, 48600, 50400, 54000, 56448, 57600, 59400, 63504, 64800, 70200, 72000, 75600, 79200, 81000, 84672, 86400, 88200, 90000, 91800, 93600, 97200, 98784, 100800, 102600
Offset: 1

Views

Author

Amiram Eldar, Dec 28 2020

Keywords

Comments

A coreful divisor d of a number k is a divisor with the same set of distinct prime factors as k, or rad(d) = rad(k), where rad(k) is the largest squarefree divisor of k (A007947).
Analogous to A068403 as A308053 is analogous to A005101.
Apparently, the least odd term in this sequence is 3^4 * 5^3 * 7^3 * 11^2 * 13^2 * 17^2 * 19^2 * 23^2 * 29^2 = 3296233276111741840875.
The asymptotic density of this sequence is Sum_{n>=1} f(A364991(n)) = 0.0004006..., where f(n) = (6/(Pi^2*n)) * Product_{prime p|n} (p/(p+1)). - Amiram Eldar, Aug 15 2023

Examples

			5400 is a term since csigma(5400) = 16380 > 3 * 5400.
		

Crossrefs

Subsequence of A308053.
Cf. A007947, A057723, A364991 (primitive terms).
Similar sequences: A068403, A285615, A293187, A300664, A307112, A328135.

Programs

  • Mathematica
    f[p_, e_] := (p^(e + 1) - 1)/(p - 1) - 1; s[1] = 1; s[n_] := Times @@ (f @@@ FactorInteger[n]); Select[Range[10^5], s[#] > 3*# &]
  • PARI
    s(n) = {my(f = factor(n)); prod(i = 1, #f~, sigma(f[i, 1]^f[i, 2]) - 1);}
    is(n) = s(n) > 3*n; \\ Amiram Eldar, Aug 15 2023

A380930 Numbers k such that A380845(k) > 3*k.

Original entry on oeis.org

1080, 2160, 3600, 4320, 7200, 7440, 8640, 11340, 13608, 14400, 14880, 15120, 17280, 18600, 22680, 22860, 27216, 28800, 29760, 30240, 30480, 31752, 33264, 34020, 34560, 37200, 41664, 45360, 45720, 45900, 51408, 53340, 54432, 57600, 59520, 60480, 60960, 61200, 63504
Offset: 1

Views

Author

Amiram Eldar, Feb 08 2025

Keywords

Comments

Analogous to 3-abundant numbers (A068403) with A380845 instead of A000203.

Crossrefs

Subsequence of A068403 and A380929.
Subsequences: A380848, A380931.
Similar sequences: A285615, A293187, A300664, A328135, A340109.

Programs

  • Mathematica
    q[k_] := Module[{h = DigitCount[k, 2, 1]}, DivisorSum[k, # &, DigitCount[#, 2, 1] == h &] > 3*k]; Select[Range[64000], q]
  • PARI
    isok(k) = {my(h = hammingweight(k)); sumdiv(k, d, d*(hammingweight(d) == h)) > 3*k;}

Formula

1080 is a term since A380845(1080) = 3330 > 3 * 1080 = 3240.
Showing 1-3 of 3 results.