OFFSET
1,2
COMMENTS
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..10000
MATHEMATICA
lps = Cases[Import["https://oeis.org/A025487/b025487.txt", "Table"], {_, _}][[;; , 2]]; f[p_, e_] := e/(e + 1); d[1] = 1; d[n_] := Denominator[Plus @@ f @@@ FactorInteger[n]]; Select[lps, d[#] == 1 &]
CROSSREFS
KEYWORD
nonn
AUTHOR
Amiram Eldar, Jul 01 2021
STATUS
approved