A336253
Exponential barely deficient numbers: exponential deficient numbers whose exponential abundancy is closer to 2 than that of any smaller exponential deficient number.
Original entry on oeis.org
1, 4, 72, 100, 144, 3528, 12100, 15876, 24336, 441000, 1334025, 2205000, 5664400, 24206400, 71267364, 151880976, 3252372552, 9346201200, 13319078472, 26828235000, 347372082000, 1851803856100, 2260121356900, 3198696480100, 5202286387272, 10330374528100, 16316106062400
Offset: 1
4 is a term since it is exponential deficient, and esigma(4)/4 = 3/2 is higher than esigma(k)/k for all the exponential deficient numbers k < 4.
-
fun[p_, e_] := DivisorSum[e, p^# &]; esigma[1] = 1; esigma[n_] := Times @@ fun @@@ FactorInteger[n]; rm = 0; s={}; Do[r = esigma[n]/n; If[r >= 2, Continue[]]; If[r > rm, rm = r; AppendTo[s, n]], {n, 1, 10^6}]; s
A336252
Infinitary barely deficient numbers: infinitary deficient numbers whose infinitary abundancy is closer to 2 than that of any smaller infinitary deficient number.
Original entry on oeis.org
1, 2, 8, 84, 110, 128, 1155, 3680, 6490, 8200, 8648, 12008, 18632, 32768, 724000, 1495688, 2095208, 3214090, 3477608, 3660008, 5076008, 12026888, 16102808, 26347688, 29322008, 33653888, 73995392, 615206030, 815634435, 2147483648, 42783299288, 80999455688
Offset: 1
8 is a term since it is infinitary deficient (A129657), and isigma(8)/8 = 15/8 is higher than isigma(k)/k for all the infinitary deficient numbers k < 8.
-
fun[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 @@ fun @@@ FactorInteger[n]; seq = {}; r = 0; Do[s = isigma[n]/n; If[s < 2 && s > r, AppendTo[seq, n]; r = s], {n, 1, 10^6}]; seq
A336672
Unitary barely 3-deficient numbers: numbers m such that usigma(k)/k < usigma(m)/m < 3 for all numbers k < m, where usigma is the sum of unitary divisors function (A034448).
Original entry on oeis.org
1, 2, 6, 30, 210, 2310, 110670, 182910, 898590, 22851570, 26266170, 45255210, 64124970, 265402410, 1374105810, 1631268870
Offset: 1
-
usigma[1] = 1; usigma[n_] := Times @@ (1 + Power @@@ FactorInteger[n]); s = {}; rm = 0; Do[r = usigma[n]/n; If[r < 3 && r > rm, rm = r; AppendTo[s, n]], {n, 1, 10^5}]; s
Showing 1-3 of 3 results.
Comments