A328135 -id:A328135 - cp's OEIS frontend

A383699 Primitive exponential 3-abundant numbers: the powerful terms of A328135.

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

Offset: 1

Amiram Eldar, May 06 2025

nonn

Subsequence of A328135 and first differ from it at n = 25: A328135(25) = 15330615300 is not a term of this sequence.

For squarefree numbers k, esigma(k) = k, where esigma is the sum of exponential divisors function (A051377). Thus, if m is a term (esigma(m) >= 3*m) and k is a squarefree number coprime to m, then esigma(k*m) = esigma(k) * esigma(m) = k * esigma(m) >= 3*k*m, so k*m is an exponential 3-abundant number. Therefore, the sequence of exponential 3-abundant numbers (A328135) can be generated from this sequence by multiplying with coprime squarefree numbers.

			901800900 is a term since esigma(901800900) = 2905943040 > 3 * 901800900 = 2705402700, and 901800900 = 2^2 * 3^2 * 5^2 * 7^2 * 11^2 * 13^2 is a powerful number.

Amiram Eldar, Table of n, a(n) for n = 1..10136 (terms below 10^16)

Intersection of A001694 and A328135.
Subsequence of A328136.
Cf. A051377, A307112.

pows[max_] := Union[Flatten[Table[i^2*j^3, {j, 1, Surd[max, 3]}, {i, 1, Sqrt[max/j^3]}]]];
f[p_, e_] := DivisorSum[e, p^# &]; esigma[1] = 1; esigma[n_] := Times @@ f @@@ FactorInteger[n];
seq[max_] := Select[pows[max], esigma[#] >= 3 # &]; seq[10^10]

A323332 The Dedekind psi function values of the powerful numbers, A001615(A001694(n)).

Original entry on oeis.org

1, 6, 12, 12, 24, 30, 36, 48, 72, 56, 96, 144, 108, 180, 216, 132, 150, 192, 288, 182, 336, 360, 432, 360, 324, 384, 576, 306, 648, 392, 380, 672, 720, 864, 672, 792, 900, 768, 552, 1152, 750, 1296, 1080, 1092, 972, 1344, 1440, 870, 1728, 2160, 992, 1584

Offset: 1

Amiram Eldar, Jan 11 2019

nonn

The sum of the reciprocals of all the terms of this sequence is Pi^2/6 (A013661).

The asymptotic density of a sequence S that possesses the property that an integer k is a term if and only if its powerful part, A057521(k) is a term, is (1/zeta(2)) * Sum_{n>=1, A001694(n) is a term of S} 1/a(n). Examples for such sequences are the e-perfect numbers (A054979), the exponential abundant numbers (A129575), and other sequences listed in the Crossrefs section. - Amiram Eldar, May 06 2025

Amiram Eldar, Table of n, a(n) for n = 1..10000
Eckford Cohen, A property of Dedekind's psi-function, Proceedings of the American Mathematical Society, Vol. 12, No. 6 (1961), p. 996.

Cf. A001615, A001694, A013661, A082695, A112526.

Sequences whose density can be calculated using this sequence: A054979, A129575, A307958, A308053, A321147, A322858, A323310, A328135, A339936, A340109, A364990, A382061, A383693, A383695, A383697.

psi[1]=1; psi[n_] := n * Times@@(1+1/Transpose[FactorInteger[n]][[1]]); psi /@ Join[{1}, Select[Range@ 1200, Min@ FactorInteger[#][[All, 2]] > 1 &]] (* after T. D. Noe at A001615 and Harvey P. Dale at A001694 *)

from math import isqrt, prod
from sympy import mobius, integer_nthroot, primefactors
def A323332(n):
    def squarefreepi(n): return int(sum(mobius(k)*(n//k**2) for k in range(1, isqrt(n)+1)))
    def bisection(f,kmin=0,kmax=1):
        while f(kmax) > kmax: kmax <<= 1
        while kmax-kmin > 1:
            kmid = kmax+kmin>>1
            if f(kmid) <= kmid:
                kmax = kmid
            else:
                kmin = kmid
        return kmax
    def f(x):
        c, l = n+x-squarefreepi(integer_nthroot(x,3)[0]), 0
        j = isqrt(x)
        while j>1:
            k2 = integer_nthroot(x//j**2,3)[0]+1
            w = squarefreepi(k2-1)
            c -= j*(w-l)
            l, j = w, isqrt(x//k2**3)
        return c+l
    a = primefactors(m:=bisection(f,n,n))
    return m*prod(p+1 for p in a)//prod(a) # Chai Wah Wu, Sep 14 2024

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

Amiram Eldar, Dec 28 2020

nonn

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

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

Amiram Eldar, Table of n, a(n) for n = 1..10000

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

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*# &]

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

Amiram Eldar, Feb 08 2025

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

Amiram Eldar, Table of n, a(n) for n = 1..10000

Cf. A000203, A380845, A380846, A380847.
Subsequence of A068403 and A380929.
Subsequences: A380848, A380931.
Similar sequences: A285615, A293187, A300664, A328135, A340109.

q[k_] := Module[{h = DigitCount[k, 2, 1]}, DivisorSum[k, # &, DigitCount[#, 2, 1] == h &] > 3*k]; Select[Range[64000], q]

isok(k) = {my(h = hammingweight(k)); sumdiv(k, d, d*(hammingweight(d) == h)) > 3*k;}

1080 is a term since A380845(1080) = 3330 > 3 * 1080 = 3240.

cp's OEIS Frontend

A383699 Primitive exponential 3-abundant numbers: the powerful terms of A328135.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

A323332 The Dedekind psi function values of the powerful numbers, A001615(A001694(n)).

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Mathematica

Python

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).

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

PARI

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

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Mathematica

PARI

Formula