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
-
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
5400 is a term since csigma(5400) = 16380 > 3 * 5400.
-
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
A372298
Primitive infinitary abundant numbers (definition 1): infinitary abundant numbers (A129656) whose all proper infinitary divisors are infinitary deficient numbers.
Original entry on oeis.org
40, 56, 70, 72, 88, 104, 756, 924, 945, 1092, 1188, 1344, 1386, 1428, 1430, 1596, 1638, 1760, 1870, 2002, 2016, 2080, 2090, 2142, 2176, 2210, 2394, 2432, 2470, 2530, 2584, 2720, 2750, 2944, 2990, 3040, 3128, 3190, 3200, 3230, 3250, 3400, 3410, 3496, 3712, 3770
Offset: 1
40 is a term since it is an infinitary abundant number and all its proper infinitary divisors, {1, 2, 4, 5, 8, 10, 20}, are infinitary deficient numbers.
24 and 30, which are infinitary abundant numbers, are not primitive, because they are divisible by 6 which is an infinitary perfect number.
-
f[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 @@ f @@@ FactorInteger[n]; idefQ[n_] := isigma[n] < 2*n; idivs[1] = {1};
idivs[n_] := Sort@ Flatten@ Outer[Times, Sequence @@ (FactorInteger[n] /. {p_, e_Integer} :> p^Select[Range[0, e], BitOr[e, #] == e &])];
q[n_] := Module[{d = idivs[n]}, Total[d] > 2*n && AllTrue[Most[d], idefQ]]; Select[Range[4000], q]
-
isidiv(d, f) = {if (d==1, return (1)); for (k=1, #f~, bne = binary(f[k, 2]); bde = binary(valuation(d, f[k, 1])); if (#bde < #bne, bde = concat(vector(#bne-#bde), bde)); for (j=1, #bne, if (! bne[j] && bde[j], return (0)); ); ); return (1); }
idivs(n) = {my(f = factor(n), d = divisors(f), idiv = []); for (k=1, #d, if (isidiv(d[k], f), idiv = concat(idiv, d[k])); ); idiv; } \\ Michel Marcus at A077609
isigma(n) = {my(f = factor(n), b); prod(i=1, #f~, b = binary(f[i, 2]); prod(k=1, #b, if(b[k], 1+f[i, 1]^(2^(#b-k)), 1)))} ;
is(n) = isigma(n) > 2*n && select(x -> x < n && isigma(x) >= 2*x, idivs(n)) == [];
A383699
Primitive exponential 3-abundant numbers: the powerful terms of A328135.
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
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.
-
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]
A307111
a(n) is the least primitive n-abundant number k with the largest possible abundancy index sigma(k)/k.
Original entry on oeis.org
3465, 6930, 19399380, 8172244080
Offset: 2
3465 is in the sequence since it is the primitive abundant (A071395) number with the largest possible abundancy index among the primitive abundant numbers: sigma(3465)/3465 = 832/385 = 2.161003... The abundancy indices of the next terms are 1248/385 = 3.241558..., 193536/46189 = 4.190088..., 642816/124729 = 5.153701...
Showing 1-5 of 5 results.
Comments