OFFSET
1,1
COMMENTS
LINKS
MATHEMATICA
Select[Range[2500], If[# == 1, 1, DivisorSigma[1, # ]] >= 2# &@ Apply[Times, FactorInteger[#] /. {p_, e_} /; e > 0 :> Prime[PrimePi@ p + 1]^e] &] (* Michael De Vlieger, Aug 27 2020 *)
PROG
CROSSREFS
KEYWORD
nonn
AUTHOR
Antti Karttunen, Aug 27 2020
STATUS
approved