OFFSET
1,1
EXAMPLE
MATHEMATICA
f[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 @@ f @@@ FactorInteger[n]; q[n_] := DivisorSigma[1, n] - isigma[n] > n; seq = {}; q1 = q[1]; Do[q2 = q[n]; If[q1 && q2, AppendTo[seq, n-1]]; q1=q2 , {n, 2, 10^8}]; seq
CROSSREFS
KEYWORD
nonn
AUTHOR
Amiram Eldar, Oct 09 2021
STATUS
approved