A334974 Infinitary admirable numbers: numbers k such that there is a proper infinitary divisor d of k such that isigma(k) - 2*d = 2*k, where isigma is the sum of infinitary divisors function (A049417).
24, 30, 40, 42, 54, 56, 66, 70, 78, 88, 96, 102, 104, 114, 120, 138, 150, 174, 186, 222, 246, 258, 270, 282, 294, 318, 354, 360, 366, 402, 420, 426, 438, 474, 486, 498, 534, 540, 582, 606, 618, 630, 642, 654, 660, 678, 726, 762, 780, 786, 822, 834, 894, 906, 942
Offset: 1
Keywords
Examples
150 is in the sequence since 150 = 1 + 2 + 3 - 6 + 25 + 50 + 75 is the sum of its proper infinitary divisors with one of them, 6, taken with a minus sign.
Links
- Amiram Eldar, Table of n, a(n) for n = 1..10000
Crossrefs
Programs
-
Mathematica
fun[p_, e_] := Module[{b = IntegerDigits[e, 2], m}, 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]; infDivQ[n_, 1] = True; infDivQ[n_, d_] := BitAnd[IntegerExponent[n, First /@ (f = FactorInteger[d])], (e = Last /@ f)] == e; infAdmQ[n_] := (ab = isigma[n] - 2 n) > 0 && EvenQ[ab] && ab/2 < n && Divisible[n, ab/2] && infDivQ[n, ab/2]; Select[Range[1000], infAdmQ]
Comments