OFFSET
1,2
LINKS
MATHEMATICA
Join[{1}, Select[Range[2, 15000], IntegerQ[Plus @@ (1/PrimePi[#[[1]]] & /@ FactorInteger[#])] &]]
PROG
(PARI) isok(k) = my(f=factor(k)); denominator(sum(i=1, #f~, 1/primepi(f[i, 1]))) == 1; \\ Michel Marcus, Mar 22 2025
CROSSREFS
KEYWORD
nonn
AUTHOR
Ilya Gutkovskiy, Mar 22 2025
STATUS
approved
