OFFSET
1,1
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..10000
MATHEMATICA
esigma[n_] := Times @@ (Sum[First[#]^d, {d, Divisors[Last[#]]}] &) /@ FactorInteger[n]; q[n_] := DivisorSigma[1, n] - esigma[n] > n; Select[Range[1, 3*10^6], q[#] && q[#+1] &]
CROSSREFS
KEYWORD
nonn
AUTHOR
Amiram Eldar, Oct 25 2021
STATUS
approved